Let $\mathbf{r} = \langle x, y, z \rangle$ be the position vector in $\mathbb{R}^3$ and $r = \|\mathbf{r}\|$. Compute the general divergence $\nabla \cdot (r^n \mathbf{r})$ for any real number $n$.