A Series Solution Approach to the Circular Restricted Gravitational Three-Body Dynamical Problem

The present manuscript examines the circular restricted gravitational three-body problem (CRGTBP) by the introduction of a new approach through the power series method. In addition, certain computational algorithms with the aid of Mathematica software are specifically designed for the problem. The algorithms or rather mathematical modules are established to determine the velocity and position of the third body’s motion. In fact, the modules led to ac-curate results and thus proved the new approach to be efficient.


Introduction
The circular restricted gravitational three-body problem (CRGTBP) is a special case of the gravitational three-body problem which is one of the most important n-body problems. The CRGTBP consists of three bodies; the two bodies of them (which are referred to as the primaries or the primary and the secondary) are moving in a circular orbit about their common mass's center under the influence of their mutual gravitation and the third body is infinitesimal mass moving under the gravitational significance of the two masses where its mass cannot influence the two masses. Furthermore, the third body has a common plane of movement as defined by both the primary and secondary [1]- [6]. The problem is tackled by studying the motion of the third body assuming full knowledge of the Many physicists and mathematicians have used the method of power series to solve a variety of unsolvable differential equations. The method has been used significantly to solve the celestial mechanic problems as better exactness of inspectional data is required for relatively optimal solutions to the bodies' dynamical equations of motion. Numerous scientists have used this method including Saad et al. [9] that used it to get recurrent algorithm for comets under non-gravitational motion. Rabe [10] used a new computational and iteration method to determine a series of periodic Trojan orbits in the restricted problem of three bodies. Deprit and Price [11] used the numerical methods to compute characteristic exponents in the planar restricted problem of three bodies. Sharaf et al. [12] found symbolic solution of the three-dimensional restricted three-body problem and applied it for any given set of initial values.
However, we aim in this paper to employ an approach based on the power series to establish an algorithm or mathematical module using Mathematica to tackle this important problem of circular restricted gravitational dynamical problem. More specifically, we determine the components of the velocity and position vectors of the third body with regard to the CRGTBP.

Dynamical Equations of Motion of CRGTBP
The dynamical equations of motion for the CRGTBP are given by the following coupled first-order system as follows [12] ( ) ( ) ( ) 1 , , , x y z , are the coordinates of the primary, the secondary and the third body, respectively. More, µ is the primary's 1 G = is the unit of the force of a gravitational constant; with 1 r denoting the distance of the primary and 2 r that of the secondary both to the third body [ .
Using Broucke's method [14], the following system of first-order differential Journal of Applied Mathematics and Physics equations [12] is thus obtained

Solution by Power Series
The equations given in (3) In these power series, the first coefficients are given by the known initial values of ( ) The following recurrence relations for the remaining coefficients are found using power series [12] Journal of Applied Mathematics and Physics

* Purpose
To generate the components of position and velocity for the third body at any time.

* Output
The components of position and velocity for the third body at any time.    It is remarkable to mention here that the accuracy is increased with an increase in the number of terms of the power series, but after NN = 50, the accuracy is fixed. So, we aren't needed to the more of terms of the power series.

Conclusion
In conclusion, the analytical solutions to the circular restricted gravitational three-body problem (CRGTBP) are determined via the application of the power series method. Also, a module or algorithm is specifically designed and implemented via the help of Mathematica software to find the components of the velocity and position vectors for the third body. Finally, the proposed methodology via the devised module has worked accurately and resulted in reliable results as shown.