Numerical solution for solving procedure for 3D motions near libration points in the Circular Restricted Three Body Problem (CR3BP)

In a recent paper in Astrophysics and Space Science Vol. 364 no. 11 (2019), S. Ershkov & D. Leschenko presented a new solving procedure for Euler-Poisson equations for solving momentum equations of the CR3BP near libration points for uniformly rotating planets having inclined orbits in the solar system with respect to the orbit of the Earth. The system of equations of the CR3BP has been explored with regard to the existence of an analytic way of presentation of the approximated solution in the vicinity of libration points. A new and elegant ansatz has been suggested in their publication whereby, in solving, the momentum equation is reduced to a system two coupled Riccati ODEs. In this paper, we presented a numerical solution of such coupled Riccati ODEs using Mathematica software package.