Numerical solution for 3D motions in the Circular Restricted Three Body Problem

Abstract

Recently Ershkov & D. Leschenko [1] 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.