Mathematical Algorithm for Solving Two–Body Problem

The classical two–body problem is the dynamical system which describes the motion of two objects. Of course, it is the simplest and only integrable system within frame of the classical Newtonian potential between two–bodies is applied. There are many applications in mathematical astronomy or in engineering and physical sciences can be analyzed within the frame work of two–body systems, it can be used in both quantum mechanics (particles motion) and celestial mechanics (Stars motion)

The multivariate of the perturbed forces in inner and outer space change the two–body from a simple and integrable system to one is complex and is not integrable. Thereby the analytical solutions will be invalid in most real applications. The analysis of two–body problem under the effect of many perturbed forces have received a comprehensive an extended study in the literature space dynamics. For example, the perturbed two–body problem by radiation pressure force and many types of drag forces or both together has been investigated [6,9,10,11,12]. Furthermore, the analytical solutions of the satellite motion within frame of non–sphericity and zonal harmonics perturbations effect have been studied by [7, 8]. [4] have studied also the dynamics of anisotropic Kepler problem with small anisotropy. They proved that at every energy level the anisotropic dynamical system has two periodic orbits.

The perturbed two–body system is not limited to the perturbations of either radiation pressure and drag force or the non–sphericity of the bodies. But the other bodies can be taken as a perturbed forces for the motion of two–body system [5]. Another consequence of two–body problem is that the three-body problem can be reduced to the system of two–body, if we consider either the third body has mass equals zero or if it moves to infinity while the other primary does not. There are a considerable work have been constructed to find the periodic solution for many perturbed two–body systems, see for details [1, 2, 13]. In addition the two–body problem within frame of the corrections law of Newtonian inverse square law of gravitation can be studied as a perturbed model [3].

The motion of planets and asteroids is a very important field of interest in astronomy, and space dynamics. It can be formulated as a dynamical system of differential equations based on the Newton’s laws of gravitation [14]. The law states that everybody attracts every other body along a line intersecting the two bodies with a force equals to F = Gm₁m₂/r², where F is the force between the two bodies, G is the gravitational constant, m₁ is the mass o f the first body, m₂ is the mass of the second body, and r the distance between the centres of masses of the bodies [15].

Additionally, the Newtonian gravitation can be extended to N–bodies by simply summing the forces [16,17,18]. Moreover, various analytical and numerical methods have been used to examine such problems of celestial mechanics comprising of the two–body, three–body and the generalized N–body problem, see [19,20,21,22,23,24,25,26]. But we aim in this paper to computationally tackle the gravitational two–body problem with the application of the power series method [14].

Vector two–body equation

We start off this section by first introducing the equation of the gravitational N–body problem as [16] (1) ${\ddot{r}}_{i} = - \frac{G}{r_{i}^{3}} (m_{0} + m_{i}) {\ddot{r}}_{i} - G \sum_{j = 1, j \neq i}^{N} m_{j} {\frac{r_{i} - r_{j}}{r_{i j}^{3}} + \frac{r_{j}}{r_{j}^{3}}}, i = 1, 2, ..., N$ {{\bf{\ddot r}}_i} = - {G \over {r_i^3}}({m_0} + {m_i}){{\bf{\ddot r}}_i} - G\sum\limits_{j = 1,j \ne i}^N {m_j}\{{{{{\bf{r}}_i} - {{\bf{r}}_j}} \over {r_{ij}^3}} + {{{{\bf{r}}_j}} \over {r_j^3}}\},i = 1,2,...,N where ${\ddot{r}}_{i}$ {{\bf{\ddot r}}_i} is the position vector of m_i relative to m₀, m_i and m₀ are masses of i^th body and central body, respectively and G is the universal constant of gravitation. Note that if all masses are equal zero except m₀, m_i, then Eq. (1) becomes (2) $\ddot{r} + \frac{μ}{r^{3}} r = o$ {\bf{\ddot r}} + {\mu \over {{r^3}}}{\bf{r}} = {\bf{o}} where μ = G(m₀ + m_i). The above equation is called the classical vector two–body equation.

Power series solution

3.1

Lagrange’s fundamental invariants

In favour of Eq. (2), the Lagrange’s fundamental invariants ε, λ and ψ are defined by [14, 27, 28] (3) $\begin{array}{l} ε = & \frac{μ}{r^{3}} \\ λ = & \frac{1}{r^{2}} (r \cdot v) \\ ψ = & \frac{1}{r^{2}} (v \cdot v) \end{array}$ \matrix{{\varepsilon =} \hfill & {{\kern 1pt} {\mu \over {{r^3}}}} \hfill \cr {\lambda =} \hfill & {{\kern 1pt} {1 \over {{r^2}}}({\bf{r}} \cdot {\bf{v}})} \hfill \cr {\psi =} \hfill & {{\kern 1pt} {1 \over {{r^2}}}({\bf{v}} \cdot {\bf{v}})} \hfill \cr} where $v = \dot{r}$ {\bf{v}} = {\bf{\dot r}} .

Lagrange’s invariants in Eqs.(3) satisfy the following differential equations [14, 27, 28] $\begin{array}{r} \dot{ε} + 3 ε λ = 0 \\ \dot{λ} + 2 λ^{2} + (ε - ψ) = 0 \\ \dot{ψ} + 2 λ ψ + 2 λ ε = 0 \end{array}$ \matrix{\hfill {\dot \varepsilon + 3\varepsilon \lambda = 0} \cr \hfill {\dot \lambda + 2{\lambda^2} + (\varepsilon - \psi) = 0} \cr \hfill {\dot \psi + 2\lambda \psi + 2\lambda \varepsilon = 0} \cr}

3.2

Basic differential equations

Here, to analyze the gravitational two-body problem by formulate and examine the following nonlinear differential equations [14, 27, 28] (4) $\begin{array}{r} \ddot{q} + ε q = 0 \\ \dot{ε} + 3 ε λ = 0 \\ \dot{λ} + 2 λ^{2} + (ε - ψ) = 0 \\ \dot{ψ} + 2 λ ψ + 2 λ ε = 0 \end{array}$ \matrix{\hfill {\ddot q + \varepsilon q = 0} \cr \hfill {\dot \varepsilon + 3\varepsilon \lambda = 0} \cr \hfill {\dot \lambda + 2{\lambda^2} + (\varepsilon - \psi) = 0} \cr \hfill {\dot \psi + 2\lambda \psi + 2\lambda \varepsilon = 0} \cr}

In analyzing the previous system of nonlinear differential equations, with employ the power series method [16]. Thus, leads to the power series expansion of the four dependent variables in Eqs. (4) as follows: (5) $\begin{array}{l} q (t) & = \sum_{n = 0}^{\infty} q_{n} {(t - t_{0})}^{n} \\ ε (t) & = \sum_{n = 0}^{\infty} ε_{n} {(t - t_{0})}^{n} \\ λ (t) & = \sum_{n = 0}^{\infty} λ_{n} {(t - t_{0})}^{n} \\ ψ (t) & = \sum_{n = 0}^{\infty} ψ_{n} {(t - t_{0})}^{n} \end{array}$ \matrix{{q(t)} \hfill & {= \sum\limits_{n = 0}^\infty {q_n}{{(t - {t_0})}^n}} \hfill \cr {\varepsilon (t)} \hfill & {= \sum\limits_{n = 0}^\infty {\varepsilon_n}{{(t - {t_0})}^n}} \hfill \cr {\lambda (t)} \hfill & {= \sum\limits_{n = 0}^\infty {\lambda_n}{{(t - {t_0})}^n}} \hfill \cr {\psi (t)} \hfill & {= \sum\limits_{n = 0}^\infty {\psi_n}{{(t - {t_0})}^n}} \hfill \cr}

Upon substitution of Eqs. (5) to Eqs. (4) and then solve for the coefficients of q_n, ε_n, λ_n and ψ_n, we get the following recurrence relations (6) $\begin{array}{l} q_{n + 2} & = - \frac{1}{(n + 1) (n + 2)} \sum_{p = 0}^{n} ε_{p} q_{n - p} \\ ε_{n + 1} & = - \frac{3}{(n + 1)} \sum_{p = 0}^{n} ε_{p} λ_{n - p} \\ λ_{n + 1} & = \frac{1}{(n + 1)} (ψ_{n} - ε_{n} - 2 \sum_{p = 0}^{n} λ_{p} λ_{n - p}) \\ ψ_{n + 1} & = - \frac{2}{(n + 1)} \sum_{p = 0}^{n} λ_{p} (ε_{n - p} + ψ_{n - p}) \end{array}$ \matrix{{{q_{n + 2}}} \hfill & {= - {1 \over {(n + 1)(n + 2)}}\sum\limits_{p = 0}^n {\varepsilon_p}{q_{n - p}}} \hfill \cr {{\varepsilon_{n + 1}}} \hfill & {= - {3 \over {(n + 1)}}\sum\limits_{p = 0}^n {\varepsilon_p}{\lambda_{n - p}}} \hfill \cr {{\lambda_{n + 1}}} \hfill & {= {1 \over {(n + 1)}}\left({{\psi_n} - {\varepsilon_n} - 2\sum\limits_{p = 0}^n {\lambda_p}{\lambda_{n - p}}} \right)} \hfill \cr {{\psi_{n + 1}}} \hfill & {= - {2 \over {(n + 1)}}\sum\limits_{p = 0}^n {\lambda_p}\left({{\varepsilon_{n - p}} + {\psi_{n - p}}} \right)} \hfill \cr}

Also with the recurrence relations above taking the following starting values $q_{0} \equiv q (t_{0}), q_{1} \equiv \dot{q} (t_{0}), ε_{0} \equiv ε (t_{0}), λ_{0} \equiv λ (t_{0}), ψ_{0} \equiv ψ (t_{0})$ {q_0} \equiv q({t_0}),\quad\quad{q_1} \equiv \dot q({t_0}),\quad\quad{\varepsilon_0} \equiv \varepsilon ({t_0}),\quad\quad{\lambda_0} \equiv \lambda ({t_0}),\quad\quad{\psi_0} \equiv \psi ({t_0})

Furthermore, the values of q₀ and q₁ are known from the formulated problem, while ε₀, λ₀ and ψ₀ to be determined from r₀ ≡ r(t₀) and v₀ ≡ v(t₀) as follows (7) $\begin{array}{l} ε_{0} = & \frac{1}{r_{0}^{3}} μ \\ λ_{0} = & \frac{1}{r_{0}^{2}} (x_{0} {\dot{x}}_{0} + y_{0} {\dot{y}}_{0} + z_{0} {\dot{z}}_{0}) \\ ψ_{0} = & \frac{1}{r_{0}^{2}} ({\dot{x}}_{0}^{2} + {\dot{y}}_{0}^{2} + {\dot{z}}_{0}^{2}) \end{array}$ \matrix{{{\varepsilon_0} =} \hfill & {{1 \over {r_0^3}}\mu} \hfill \cr {{\lambda_0} =} \hfill & {{1 \over {r_0^2}}({x_0}{{\dot x}_0} + {y_0}{{\dot y}_0} + {z_0}{{\dot z}_0})} \hfill \cr {{\psi_0} =} \hfill & {{1 \over {r_0^2}}(\dot x_0^2 + \dot y_0^2 + \dot z_0^2)} \hfill \cr} with $r_{0} = \sqrt{x_{0}^{2} + y_{0}^{2} + z_{0}^{2}}$ {r_0} = \sqrt {x_0^2 + y_0^2 + z_0^2} .

In the next section with a help of Eqs. (6, 7), the computational values of q_n are determined for n = 9 using a designed algorithm which is called QPS, see Appendix A for algorithm.

3.3

f and g functions

Here, two Lagrangian coefficients functions f and g [14,27] introduced instead of one that led to the following position vector (8) $r = f r_{0} + g v_{0}$ {\bf{r}} = f{{\bf{r}}_0} + g{{\bf{v}}_0} where r₀ and v₀ are the initial values of the position and velocity vectors. Note also that the velocity vector is given by (9) $v = \dot{f} r_{0} + \dot{g} v_{0}$ {\bf{v}} = \dot f{{\bf{r}}_0} + \dot g{{\bf{v}}_0}

Moreover, it is remarkable that Eqs. (8, 9) satisfies the first equation in Eq. (4) with the following initial conditions [14, 27] $\begin{array}{l} q (t_{0}) = {\begin{matrix} 1 & q = f \\ 0 & q = g \end{matrix} \\ \frac{d q (t_{0})}{d t} = {\begin{matrix} 0 & q = f \\ 1 & q = g \end{matrix} \end{array}$ \matrix{{q({t_0}) = \left\{{\matrix{1 & {q = f} \cr 0 & {q = g} \cr}} \right.} \hfill \cr\cr {{{dq({t_0})} \over {dt}} = \left\{{\matrix{0 & {q = f} \cr 1 & {q = g} \cr}} \right.} \hfill \cr}

As in the above, using the power series method to get the relation recursively. The power series expansion of these functions take the form [27] $\begin{array}{l} f = \sum_{n = 0}^{\infty} f_{n} {(t - t_{0})}^{n} \\ g = \sum_{n = 0}^{\infty} g_{n} {(t - t_{0})}^{n} \end{array}$ \matrix{{f = \sum\limits_{n = 0}^\infty {f_n}{{(t - {t_0})}^n}} \hfill \cr {g = \sum\limits_{n = 0}^\infty {g_n}{{(t - {t_0})}^n}} \hfill \cr} together with the starting values [14, 27] $\begin{array}{l} f_{0} = 1 & f_{1} = 0 \\ g_{0} = 0 & g_{1} = 1 \end{array}$ \matrix{{{f_0} = 1} \hfill & {} \quad\quad \hfill & {{f_1} = 0} \hfill \cr {{g_0} = 0} \hfill & {} \hfill & {{g_1} = 1} \hfill \cr}

Also in the next section, analytical computations determine for the values of f_n and g_n for n = 10 using a designed algorithm called FGPS, see Appendix B for the algorithm.

Results and Discussion

In this chapter, computational values determine for the values q_{n +2} for n=9, f_n and g_n for n = 10 using Mathematica software with designed computational algorithms QPS and FGPS are given in Appendix A and Appendix B, respectively.

Furthermore, a mathematical modules establish to find the scalar values of the position and velocity vectors of the two-body motion via Mathematica software with the designed computational algorithm RVPS is given in Appendix C then we apply this module for different kinds of orbits (elliptic, parabolic and hyperbolic).

4.1

Values of q_n

With the application of QPS algorithm for n = 9, the following symbolic expression of the q’s coefficients are $\begin{array}{l} q_{0} = q_{0} \\ q_{1} = q_{1} \\ q_{2} = \frac{1}{2} (- q_{0}) ε_{0} \\ q_{3} = \frac{1}{6} (- ε_{0}) (q_{1} - (3 q_{0}) λ_{0}) \end{array}$ \matrix{{{{\rm{q}}_0} = {q_0}} \hfill \cr {{{\rm{q}}_1} = {q_1}} \hfill \cr {{{\rm{q}}_2} = {1 \over 2}\left({- {q_0}} \right){\varepsilon_0}} \hfill \cr {{{\rm{q}}_3} = {1 \over 6}\left({- {\varepsilon_0}} \right)\left({{q_1} - \left({3{q_0}} \right){\lambda_0}} \right)} \hfill \cr} $\begin{array}{l} q_{4} = \frac{1}{24} ε_{0} ((6 q_{1}) λ_{0} + q_{0} (3 (ψ_{0} - 5 λ_{0}^{2}) - 2 ε_{0})) \\ q_{5} = \frac{1}{120} ε_{0} (((15 q_{0}) λ_{0}) (7 λ_{0}^{2} + 2 ε_{0} - 3 ψ_{0}) + q_{1} (9 (ψ_{0} - 5 λ_{0}^{2}) - 8 ε_{0})) \\ \begin{array}{l} q_{6} = \frac{1}{720} ε_{0} (((30 q_{1}) λ_{0}) (14 λ_{0}^{2} + 5 ε_{0} - 6 ψ_{0}) - q_{0} (22 ε_{0}^{2} + (6 ε_{0}) (70 λ_{0}^{2} - 11 ψ_{0}) \\ + 45 (21 λ_{0}^{4} + ψ_{0}^{2} - (14 λ_{0}^{2}) ψ_{0}))) \end{array} \end{array}$ \matrix{{{q_4} = {1 \over {24}}{\varepsilon_0}\left({\left({6{q_1}} \right){\lambda_0} + {q_0}\left({3\left({{\psi_0} - 5\lambda_0^2} \right) - 2{\varepsilon_0}} \right)} \right)} \hfill \cr {{q_5} = {1 \over {120}}{\varepsilon_0}\left({\left({\left({15{q_0}} \right){\lambda_0}} \right)\left({7\lambda_0^2 + 2{\varepsilon_0} - 3{\psi_0}} \right) + {q_1}\left({9\left({{\psi_0} - 5\lambda_0^2} \right) - 8{\varepsilon_0}} \right)} \right)} \hfill \cr \matrix{{q_6} = {1 \over {720}}{\varepsilon_0}(((30{q_1}){\lambda_0})(14\lambda_0^2 + 5{\varepsilon_0} - 6{\psi_0}) - {q_0}(22\varepsilon_0^2 + (6{\varepsilon_0})(70\lambda_0^2 - 11{\psi_0}) \hfill \cr \quad \quad + 45(21\lambda_0^4 + \psi_0^2 - (14\lambda_0^2){\psi_0}))) \hfill \cr} \hfill \cr} $\begin{array}{l} q_{7} = \frac{1}{5040} (ε_{0} (((63 q_{0}) λ_{0}) (165 λ_{0}^{4} + (100 ε_{0}) λ_{0}^{2} + 12 ε_{0}^{2} + 25 ψ_{0}^{2} - (6 (25 λ_{0}^{2} \\ + 6 ε_{0})) ψ_{0}) - q_{1} (172 ε_{0}^{2} + (36 ε_{0}) (70 λ_{0}^{2} - 11 ψ_{0}) + 225 (21 λ_{0}^{4} + ψ_{0}^{2} \\ - (14 λ_{0}^{2}) ψ_{0})))) \\ q_{8} = \frac{1}{40320} (ε_{0} (((126 q_{1}) λ_{0}) (495 λ_{0}^{4} + (350 ε_{0}) λ_{0}^{2} + 52 ε_{0}^{2} + 75 ψ_{0}^{2} - (18 (25 λ_{0}^{2} \\ + 7 ε_{0})) ψ_{0}) + q_{0} (- 584 ε_{0}^{3} + (36 ε_{0}^{2}) (73 ψ_{0} - 560 λ_{0}^{2}) - (54 ε_{0}) (1925 λ_{0}^{4} \\ + 67 ψ_{0}^{2} - (1120 λ_{0}^{2}) ψ_{0}) - 315 (429 λ_{0}^{6} - 5 ψ_{0}^{3} + (135 λ_{0}^{2}) ψ_{0}^{2} \\ - (495 λ_{0}^{4}) ψ_{0})))) \\ q_{9} = \frac{1}{362880} (ε_{0} (((15 q_{0}) λ_{0}) (2368 ε_{0}^{3} + (444 ε_{0}^{2}) (77 λ_{0}^{2} - 24 ψ_{0}) \\ + (18 ε_{0}) (7007 λ_{0}^{4} + 827 ψ_{0}^{2} - (5698 λ_{0}^{2}) ψ_{0}) + 189 (715 λ_{0}^{6} - 35 ψ_{0}^{3} \\ + (385 λ_{0}^{2}) ψ_{0}^{2} - (1001 λ_{0}^{4}) ψ_{0})) - q_{1} (7136 ε_{0}^{3} + (108 ε_{0}^{2}) (1785 λ_{0}^{2} \\ - 232 ψ_{0}) + (432 ε_{0}) (1925 λ_{0}^{4} + 67 ψ_{0}^{2} - (1120 λ_{0}^{2}) ψ_{0}) + 2205 (429 λ_{0}^{6} \\ - 5 ψ_{0}^{3} + (135 λ_{0}^{2}) ψ_{0}^{2} - (495 λ_{0}^{4}) ψ_{0})))) \end{array}$ \matrix{{{q_7} = {1 \over {5040}}({\varepsilon_0}(((63{q_0}){\lambda_0})(165\lambda_0^4 + (100{\varepsilon_0})\lambda_0^2 + 12\varepsilon_0^2 + 25\psi_0^2 - (6(25\lambda_0^2} \hfill \cr {\quad + 6{\varepsilon_0})){\psi_0}) - {q_1}(172\varepsilon_0^2 + (36{\varepsilon_0})(70\lambda_0^2 - 11{\psi_0}) + 225(21\lambda_0^4 + \psi_0^2} \hfill \cr {\quad - (14\lambda_0^2){\psi_0}))))} \hfill \cr {{q_8} = {1 \over {40320}}({\varepsilon_0}(((126{q_1}){\lambda_0})(495\lambda_0^4 + (350{\varepsilon_0})\lambda_0^2 + 52\varepsilon_0^2 + 75\psi_0^2 - (18(25\lambda_0^2} \hfill \cr {\quad + 7{\varepsilon_0})){\psi_0}) + {q_0}(- 584\varepsilon_0^3 + (36\varepsilon_0^2)(73{\psi_0} - 560\lambda_0^2) - (54{\varepsilon_0})(1925\lambda_0^4} \hfill \cr {\quad + 67\psi_0^2 - (1120\lambda_0^2){\psi_0}) - 315(429\lambda_0^6 - 5\psi_0^3 + (135\lambda_0^2)\psi_0^2} \hfill \cr {\quad - (495\lambda_0^4){\psi_0}))))} \hfill \cr {{q_9} = {1 \over {362880}}({\varepsilon_0}(((15{q_0}){\lambda_0})(2368\varepsilon_0^3 + (444\varepsilon_0^2)(77\lambda_0^2 - 24{\psi_0})} \hfill \cr {\quad + (18{\varepsilon_0})(7007\lambda_0^4 + 827\psi_0^2 - (5698\lambda_0^2){\psi_0}) + 189(715\lambda_0^6 - 35\psi_0^3} \hfill \cr {\quad + (385\lambda_0^2)\psi_0^2 - (1001\lambda_0^4){\psi_0})) - {q_1}(7136\varepsilon_0^3 + (108\varepsilon_0^2)(1785\lambda_0^2} \hfill \cr {\quad - 232{\psi_0}) + (432{\varepsilon_0})(1925\lambda_0^4 + 67\psi_0^2 - (1120\lambda_0^2){\psi_0}) + 2205(429\lambda_0^6} \hfill \cr {\quad - 5\psi_0^3 + (135\lambda_0^2)\psi_0^2 - (495\lambda_0^4){\psi_0}))))} \hfill \cr}

4.2

Values of f_n

With the application of FGPS for n = 10; T = 1, 2, the following symbolic expressions of the f’s are $\begin{array}{l} f_{0} = 1 \\ f_{1} = 0 \\ f_{2} = - \frac{ε_{0}}{2} \\ f_{3} = \frac{ε_{0} λ_{0}}{2} \end{array}$ \matrix{{{f_0} = 1} \hfill \cr {{f_1} = 0} \hfill \cr {{f_2} = - {{{\varepsilon_0}} \over 2}} \hfill \cr {{f_3} = {{{\varepsilon_0}{\lambda_0}} \over 2}} \hfill \cr} $\begin{array}{l} f_{4} = \frac{1}{24} (- ε_{0}) (15 λ_{0}^{2} + 2 ε_{0} - 3 ψ_{0}) \\ f_{5} = \frac{1}{8} (ε_{0} λ_{0}) (7 λ_{0}^{2} + 2 ε_{0} - 3 ψ_{0}) \\ f_{6} = \frac{1}{720} (- ε_{0}) (22 ε_{0}^{2} + (6 ε_{0}) (70 λ_{0}^{2} - 11 ψ_{0}) + 45 (21 λ_{0}^{4} + ψ_{0}^{2} - (14 λ_{0}^{2}) ψ_{0})) \end{array}$ \matrix{{{f_4} = {1 \over {24}}\left({- {\varepsilon_0}} \right)\left({15\lambda_0^2 + 2{\varepsilon_0} - 3{\psi_0}} \right)} \hfill \cr {{f_5} = {1 \over 8}\left({{\varepsilon_0}{\lambda_0}} \right)\left({7\lambda_0^2 + 2{\varepsilon_0} - 3{\psi_0}} \right)} \hfill \cr {{f_6} = {1 \over {720}}\left({- {\varepsilon_0}} \right)\left({22\varepsilon_0^2 + \left({6{\varepsilon_0}} \right)\left({70\lambda_0^2 - 11{\psi_0}} \right) + 45\left({21\lambda_0^4 + \psi_0^2 - \left({14\lambda_0^2} \right){\psi_0}} \right)} \right)} \hfill \cr} $\begin{array}{l} f_{7} = \frac{1}{80} (ε_{0} λ_{0}) (165 λ_{0}^{4} + (100 ε_{0}) λ_{0}^{2} + 12 ε_{0}^{2} + 25 ψ_{0}^{2} - (6 (25 λ_{0}^{2} + 6 ε_{0})) ψ_{0}) \\ f_{8} = - \frac{1}{40320} (ε_{0} (584 ε_{0}^{3} + (36 ε_{0}^{2}) (560 λ_{0}^{2} - 73 ψ_{0}) + (54 ε_{0}) (1925 λ_{0}^{4} + 67 ψ_{0}^{2} \\ - (1120 λ_{0}^{2}) ψ_{0}) + 315 (429 λ_{0}^{6} - 5 ψ_{0}^{3} + (135 λ_{0}^{2}) ψ_{0}^{2} - (495 λ_{0}^{4}) ψ_{0}))) \\ f_{9} = \frac{1}{24192} ((ε_{0} λ_{0}) (2368 ε_{0}^{3} + (444 ε_{0}^{2}) (77 λ_{0}^{2} - 24 ψ_{0}) + (18 ε_{0}) (7007 λ_{0}^{4} \\ + 827 ψ_{0}^{2} - (5698 λ_{0}^{2}) ψ_{0}) + 189 (715 λ_{0}^{6} - 35 ψ_{0}^{3} + (385 λ_{0}^{2}) ψ_{0}^{2} \\ - (1001 λ_{0}^{4}) ψ_{0}))) \\ f_{10} = - \frac{1}{3628800} (ε_{0} (28384 ε_{0}^{4} + (48 ε_{0}^{3}) (31735 λ_{0}^{2} - 3548 ψ_{0}) \\ + (54 ε_{0}^{2}) (245245 λ_{0}^{4} + 6559 ψ_{0}^{2} - (126940 λ_{0}^{2}) ψ_{0}) + (90 ε_{0}) (420420 λ_{0}^{6} \\ - 3461 ψ_{0}^{3} + (107514 λ_{0}^{2}) ψ_{0}^{2} - (441441 λ_{0}^{4}) ψ_{0}) + 14175 (2431 λ_{0}^{8} + 7 ψ_{0}^{4} \\ - (308 λ_{0}^{2}) ψ_{0}^{3} + (2002 λ_{0}^{4}) ψ_{0}^{2} - (4004 λ_{0}^{6}) ψ_{0}))) \end{array}$ \matrix{{{f_7} = {1 \over {80}}\left({{\varepsilon_0}{\lambda_0}} \right)\left({165\lambda_0^4 + \left({100{\varepsilon_0}} \right)\lambda_0^2 + 12\varepsilon_0^2 + 25\psi_0^2 - \left({6\left({25\lambda_0^2 + 6{\varepsilon_0}} \right)} \right){\psi_0}} \right)} \hfill \cr {{f_8} = - {1 \over {40320}}({\varepsilon_0}(584\varepsilon_0^3 + (36\varepsilon_0^2)(560\lambda_0^2 - 73{\psi_0}) + (54{\varepsilon_0})(1925\lambda_0^4 + 67\psi_0^2} \hfill \cr {\quad - (1120\lambda_0^2){\psi_0}) + 315(429\lambda_0^6 - 5\psi_0^3 + (135\lambda_0^2)\psi_0^2 - (495\lambda_0^4){\psi_0})))} \hfill \cr {{f_9} = {1 \over {24192}}(({\varepsilon_0}{\lambda_0})(2368\varepsilon_0^3 + (444\varepsilon_0^2)(77\lambda_0^2 - 24{\psi_0}) + (18{\varepsilon_0})(7007\lambda_0^4} \hfill \cr {\quad + 827\psi_0^2 - (5698\lambda_0^2){\psi_0}) + 189(715\lambda_0^6 - 35\psi_0^3 + (385\lambda_0^2)\psi_0^2} \hfill \cr {\quad - (1001\lambda_0^4){\psi_0})))} \hfill \cr {{f_{10}} = - {1 \over {3628800}}({\varepsilon_0}(28384\varepsilon_0^4 + (48\varepsilon_0^3)(31735\lambda_0^2 - 3548{\psi_0})} \hfill \cr {\quad + (54\varepsilon_0^2)(245245\lambda_0^4 + 6559\psi_0^2 - (126940\lambda_0^2){\psi_0}) + (90{\varepsilon_0})(420420\lambda_0^6} \hfill \cr {\quad - 3461\psi_0^3 + (107514\lambda_0^2)\psi_0^2 - (441441\lambda_0^4){\psi_0}) + 14175(2431\lambda_0^8 + 7\psi_0^4} \hfill \cr {\quad - (308\lambda_0^2)\psi_0^3 + (2002\lambda_0^4)\psi_0^2 - (4004\lambda_0^6){\psi_0})))} \hfill \cr}

4.3

Values of g_n

Again with the application of FGPS for n = 10; T = 1, 2, the following symbolic expressions of g’s coefficients are $\begin{array}{l} g_{0} = 0 \\ g_{1} = 1 \\ g_{2} = 0 \\ g_{3} = - \frac{ε_{0}}{6} \\ g_{4} = \frac{ε_{0} λ_{0}}{4} \end{array}$ \matrix{{{g_0} = 0} \hfill \cr {{g_1} = 1} \hfill \cr {{g_2} = 0} \hfill \cr {{g_3} = - {{{\varepsilon_0}} \over 6}} \hfill \cr {{g_4} = {{{\varepsilon_0}{\lambda_0}} \over 4}} \hfill \cr} $\begin{array}{l} g_{5} = \frac{1}{120} (- ε_{0}) (45 λ_{0}^{2} + 8 ε_{0} - 9 ψ_{0}) \\ g_{6} = \frac{1}{24} (ε_{0} λ_{0}) (14 λ_{0}^{2} + 5 ε_{0} - 6 ψ_{0}) \\ g_{7} = - \frac{ε_{0} (172 ε_{0}^{2} + (36 ε_{0}) (70 λ_{0}^{2} - 11 ψ_{0}) + 225 (21 λ_{0}^{4} + ψ_{0}^{2} - (14 λ_{0}^{2}) ψ_{0}))}{5040} \end{array}$ \matrix{{{g_5} = {1 \over {120}}\left({- {\varepsilon_0}} \right)\left({45\lambda_0^2 + 8{\varepsilon_0} - 9{\psi_0}} \right)} \hfill \cr {{g_6} = {1 \over {24}}\left({{\varepsilon_0}{\lambda_0}} \right)\left({14\lambda_0^2 + 5{\varepsilon_0} - 6{\psi_0}} \right)} \hfill \cr {{g_7} = - {{{\varepsilon_0}\left({172\varepsilon_0^2 + \left({36{\varepsilon_0}} \right)\left({70\lambda_0^2 - 11{\psi_0}} \right) + 225\left({21\lambda_0^4 + \psi_0^2 - \left({14\lambda_0^2} \right){\psi_0}} \right)} \right)} \over {5040}}} \hfill \cr} $\begin{array}{l} g_{8} = \frac{1}{320} (ε_{0} λ_{0}) (495 λ_{0}^{4} + (350 ε_{0}) λ_{0}^{2} + 52 ε_{0}^{2} + 75 ψ_{0}^{2} - (18 (25 λ_{0}^{2} + 7 ε_{0})) ψ_{0}) \\ g_{9} = - \frac{1}{362880} (ε_{0} (7136 ε_{0}^{3} + (108 ε_{0}^{2}) (1785 λ_{0}^{2} - 232 ψ_{0}) + (432 ε_{0}) (1925 λ_{0}^{4} \\ + 67 ψ_{0}^{2} - (1120 λ_{0}^{2}) ψ_{0}) + 2205 (429 λ_{0}^{6} - 5 ψ_{0}^{3} + (135 λ_{0}^{2}) ψ_{0}^{2} \\ - (495 λ_{0}^{4}) ψ_{0}))) \\ g_{10} = \frac{1}{120960} ((ε_{0} λ_{0}) (15220 ε_{0}^{3} + (12 ε_{0}^{2}) (14938 λ_{0}^{2} - 4647 ψ_{0}) + (81 ε_{0}) (7007 λ_{0}^{4} \\ + 827 ψ_{0}^{2} - (5698 λ_{0}^{2}) ψ_{0}) + 756 (715 λ_{0}^{6} - 35 ψ_{0}^{3} + (385 λ_{0}^{2}) ψ_{0}^{2} \\ - (1001 λ_{0}^{4}) ψ_{0}))) \end{array}$ \matrix{{{g_8} = {1 \over {320}}\left({{\varepsilon_0}{\lambda_0}} \right)\left({495\lambda_0^4 + \left({350{\varepsilon_0}} \right)\lambda_0^2 + 52\varepsilon_0^2 + 75\psi_0^2 - \left({18\left({25\lambda_0^2 + 7{\varepsilon_0}} \right)} \right){\psi_0}} \right)} \hfill \cr {{g_9} = - {1 \over {362880}}({\varepsilon_0}(7136\varepsilon_0^3 + (108\varepsilon_0^2)(1785\lambda_0^2 - 232{\psi_0}) + (432{\varepsilon_0})(1925\lambda_0^4} \hfill \cr {\quad\quad + 67\psi_0^2 - (1120\lambda_0^2){\psi_0}) + 2205(429\lambda_0^6 - 5\psi_0^3 + (135\lambda_0^2)\psi_0^2} \hfill \cr {\quad\quad - (495\lambda_0^4){\psi_0})))} \hfill \cr {{g_{10}} = {1 \over {120960}}(({\varepsilon_0}{\lambda_0})(15220\varepsilon_0^3 + (12\varepsilon_0^2)(14938\lambda_0^2 - 4647{\psi_0}) + (81{\varepsilon_0})(7007\lambda_0^4} \hfill \cr {\quad\quad + 827\psi_0^2 - (5698\lambda_0^2){\psi_0}) + 756(715\lambda_0^6 - 35\psi_0^3 + (385\lambda_0^2)\psi_0^2} \hfill \cr {\quad\quad - (1001\lambda_0^4){\psi_0})))} \hfill \cr}

4.4

Applications on orbits

With the application of RVPS algorithm on different kinds of orbits with the initial values components for the position and velocity vectors is given by then, we get the final scalar values of the position and velocity vectors

Conclusion

The dynamical system of the two–body problem is addressed, by applying the mathematical modules, which are established to find the scalar values of the position and velocity vectors of the two–body motion. The Application of this module for different kinds of orbits (elliptic, parabolic and hyperbolic) leads to accurate results, that prove module efficiency.

eISSN:: 2444-8656
Idioma:: Inglés

Calendario de la edición:: Volume Open
Temas de la revista:: Life Sciences, other, Mathematics, Applied Mathematics, General Mathematics, Physics

RSS Feed de revista

Mathematical Algorithm for Solving Two–Body Problem

Publicado en línea: 28 oct 2020

Páginas: 217 - 228

Recibido: 06 mar 2020

Aceptado: 05 abr 2020

DOI: https://doi.org/10.2478/amns.2020.2.00039

Palabras claveTwo–body problem, Power Series Methods, Mathematica software

© 2020 M. H. Alghamdi et al., published by Sciendo

This work is licensed under the Creative Commons Attribution 4.0 International License.

Palabras clave
Two–body problem, Power Series Methods, Mathematica software