“…To avoid singularity issues, a set of nonsingular orbital elements is used in this work [27], where Ω is the right ascension of the ascending node, ω is the argument of perigee, e is the eccentricity, i is the inclination, M is the mean anomaly, and θ G is the Greenwich sidereal time, measured from the instantaneous equinox. When low-thrust propulsion is considered, the equations of motions are given as [16] ( ) ( ) where v represents the velocity of near-circular orbits, v ¼na, n is the mean motion, Ω E is the Earth's rotation rate, l is the right ascension of the satellite, Ω ω = + + l M .…”