From the potential theory, the gravitational potential (any object has axial
symmetry \(m_{0}\) ) experienced by the satellite \(m\)  will be controlled by (see Murray and Dermott [1]})
\(V=-\dfrac{Gm_{0}m}{r_{0}}\left[1-\sum _{n=2}^{+\infty }J_{n}\left(\dfrac{R}{r_{0}}\right)p_{n}(\sin (\delta ))\right]\) 
