\nabla \phi = g \nabla r^{-1} =
\frac{gmM}{r^2} \frac{(x, y, z)}{\|\overrightarrow{r}\|}