Iterative and Direct Determinations of the Gravitational Potential of Axisymmetric Mass Distributions in Spherical Coordinates