^{1}

^{*}

^{2}

^{*}

In this paper of the series, the equations of motion for the spatial circular restricted three-body problem in sidereal spherical coordinates system were established. Initial value procedure that can be used to compute both the spherical and Cartesian sidereal coordinates and velocities was also developed. The application of the procedure was illustrated by numerical example and graphical representations of the variations of the two sidereal coordinate systems.

In a previous communication to this journal [

If two of the bodies, say m_{1} and m_{2} in the three-body problem move in circular, coplanar orbits about their common center of mass and the mass say m_{3} of the third body is too small to affect the motion of the other bodies, the problem of the motion of the third body is called the circular, restricted, three body problem. The two revolving bodies are called the primaries; their masses are arbitrary but have such internal mass distributions that they may be considered point masses.

The equations of motion of the third body in a dimensionless sidereal (inertial) coordinate system with the mean motion, are [

where is given as

denotes the mass of the smaller primary when the total mass of the primaries has been normalized to unity.

and are the distances of the third body from the primaries which are located at;, these coordinates are functions of the time t and are given as

Corresponding to the Cartesian sidereal coordinate system, the coordinate system related to the system by certain transformation, is also called sidereal coordinate system. In this respect the system of Equation (7) is called sidereal spherical coordinate system.

In what follows we shall establish, the differential equations for the spatial circular restricted three bodyproblem in sidereal spherical coordinate system.

where

From Equation (7) we have

Differentiating the first and the third of Equation (10) and the third of Equation (7) with respect to the time t we get:

where and are given in terms of and from the previous equations.

The kinetic energy of a particle of unit mass in the spherical coordinate system is

By using the transformation equations (Equations (7)), the gravitational potential V could be expressed in term of.

Consequently, we deduce for the equations of motion in sidereal spherical coordinate system, the forms

where are given as

, and; can be computed from Equation (7), while, and can be computed from Equations (1)-(3), so we get

where

,

,

,

.

In what follows, we shall establish a procedure that can be used to compute (say) both:

1) The spherical sidereal coordinates and velocities, and 2) The Cartesian sidereal coordinates and velocities.

So, such procedure is a double usefulness computational algorithm, for which a differential solver can be used for the spherical sidereal six order system to obtain. While the Cartesian sidereal coordinates and velocities are obtained by the substitutions in the direct transformation formulae (Equations (7) and (8)), rather than solving the six order system of Equations (1), (2) and (3). By this way, great time can be saved.

This initial value procedure using sidereal spherical coordinate system will be described through its basic points, input, output and computational steps.

Input: 1) at2) the final time

3)

Output: 1)

2)

1) Using the given values at and the inverse transformations to compute the initial values.

2) Using the partial derivatives (functions of) to construct the analytical forms of equations of motion as first order system.

3) Using the initial conditions from step 1 to solve numerically the above differential system of step 2 for , (note that).

4) Using from step 3 and the direct transformations of Equations (7) and (8) to compute numerically and .

5) End.

Consider the initial values

applying the above procedure we get the results as displayed in Tables 1 and 2.

In this paper of the series, the equations of motion for the spatial circular restricted three-body problem in sidereal spherical coordinates system were established. Initial value procedure that can be used to compute both the spherical and Cartesian sidereal coordinates and velocities was also developed. The application of the procedure

was illustrated by numerical example and graphical representations of the variations of the two sidereal coordinate systems.