“…Indeed, around the poles of the star the cylindrical coordinate ρ is such that ρ = r sin θ r. The results obtained in Section 3.1 can then easily be adapted to such a case (as done in Ref. [23] for the inside of a sphere) by expanding here, instead, the non-relativistic Newtonian potential around the poles of the star into that of a simple harmonic oscillator. Denoting now the radius of the star by R, we have at the poles, r = R 2 + ρ 2 , with ρ R. Therefore, we have, up to a constant, −GM/r ∼ 1 2 GMρ 2 /R 3 .…”