Lagrange Points and Polynomial Equations: Part 2

From Wikipedia, Lagrange points are points of equilibrium for small-mass objects under the gravitational influence of two massive orbiting bodies. There are five such points in the Sun-Earth system, called $L_1$ , $L_2$ , $L_3$ , $L_4$ , and $L_5$ .

The stable equilibrium points $L_4$ and $L_5$ are easiest to explain: they are the corners of equilateral triangles in the plane of Earth’s orbit. The points $L_1$ and $L_2$ are also equilibrium points, but they are unstable. Nevertheless, they have practical applications for spaceflight.

We begin with $L_1$ , whose position can be found by numerically solving the fifth-order polynomial equation

$(m_1+m_3)x^5+(3m_1+2m_3)x^4+(3m_1+m_3)x^3$

$-(3m_2+m_3)x^2-(3m_2+m_3)x-(m_2+m_3)=0$ .

In this equation, $m_1$ is the mass of the Sun, $m_2$ is the mass of Earth, $m_3$ is the mass of the spacecraft, and $x$ is the distance from the Earth to $L_1$ measured as a proportion of the distance from the Sun to $L_1$ . In other words, if the distance from the Sun to $L_1$ is 1 unit, then the distance from the Earth to $L_1$ is $x$ units. The above equation is derived using principles from physics which are not elaborated upon here.

We notice that the coefficients of $x^5$ , $x^4$ , and $x^3$ are all positive, while the coefficients of $x^2$ , $x$ , and the constant term are all negative. Therefore, since there is only one change in sign, this equation has only one positive real root by Descartes’ Rule of Signs.

Since $m_3$ is orders of magnitude smaller than both $m_1$ and $m_2$ , this may safely approximated by

$m_1 x^5 + 3m_1 x^4 + 3m_1 x^3 - 3m_2 x^2 - 3m_2x - m_2=0$ .

Unfortunately, the unit $x$ is not as natural for Earth-bound observers as $t$ , the proportion of the distance of $L_1$ to Earth as a proportion of the distance from the Sun to Earth. Since $L_1$ is between the Sun and Earth, the distance from the Sun to Earth is $x+1$ units, so that $t = x/(x+1)$ . We then solve for $x$ in terms of $t$ (just like finding an inverse function):

$t = \displaystyle \frac{x}{x+1}$

$t(x+1) = x$

$tx + t = x$

$t = x - tx$

$t= x(1-t)$

$\displaystyle \frac{t}{1-t} = x$ .

Substituting into the above equation, we find an equation for $t$ :

$\displaystyle \frac{m_1t^5}{(1-t)^5} + \frac{3m_1t^4}{(1-t)^4} + \frac{3m_1t^3}{(1-t)^3} - \frac{3m_2t^2}{(1-t)^2} - \frac{3m_2t}{1-t} - m_2=0$

$m_1t^5 + 3m_1t^4(1-t) + 3m_1t^3(1-t)^2 - 3m_2t^2(1-t)^3 - 3m_2t(1-t)^4 - m_2(1-t)^5=0$

Expanding, we find

$m_1 t^5 + 3m_1 (t^4 - t^5) + 3m_1 (t^3-2t^4+t^5) - 3m_2 (t^2-3t^3+3t^4-t^5)$

$-3m_2(t - 4t^2 + 6t^3 - 4t^4 + t^5) - m_2(1 - 5t + 10t^2 - 10 t^3 + 5t^4 + t^5) = 0$

Collecting like terms, we find

$(m_1 - 3m_1 + 3m_1 + 3m_2 - 3m_2 + m_2)t^5 + (3m_1-6m_1-9m_2+12m_2-5m_2)t^4$

$+ (3m_1+9m_2-18m_2+10m_2)t^3 + (-3m_2+12m_2-10m_2) t^2$

$+ (-3m_2+5m_2)t - m_2 = 0$ ,

$(m_1+m_2) t^5 - (3m_1 +2m_2) t^4 + (3m_1 + m_2) t^3 - m_2 t^2 + 2m_2 t- m_2 = 0$ .

Again, this equation has only one positive real root since the original quintic in $x$ only had one positive real root. This new equation can be rewritten as

$t^5 - \displaystyle \frac{3m_1 + 2m_2}{m_1 + m_2} t^4 + \frac{3m_1+m_2}{m_1+m_2} t^3 - \frac{m_2}{m_1+m_2} t^2 + \frac{2m_2}{m_1+m_2} t- \frac{m_2}{m_1+m_2} = 0$ .

If we define

$\mu = \displaystyle \frac{m_2}{m_1+m_2}$ ,

we see that

$\displaystyle \frac{3m_1 + 2m_2}{m_1 + m_2} = \frac{3m_1 + 3m_2}{m_1 + m_2} - \frac{m_2}{m_1 + m_2} = 3 - \mu$

and

$\displaystyle \frac{3m_1 + m_2}{m_1 + m_2} = \frac{3m_1 + 3m_2}{m_1 + m_2} - \frac{2m_2}{m_1 + m_2} = 3 - 2\mu$ ,

so that the equation may be written as

$t^5 + (\mu-3) t^4 + (3-2\mu) - \mu t^2 + 2\mu t - \mu = 0$ ,

matching the equation found at Wikipedia.

For the Sun and Earth, $m_1 \approx 1.9885 \times 10^{30} ~ \hbox{kg}$ and $m_2 \approx 5.9724 \times 10^{24} ~ \hbox{kg}$ , so that

$\mu = \displaystyle \frac{5.9724 \times 10^{24}}{1.9885 \times 10^{30} + 5.9724 \times 10^{24}} \approx 3.00346 \times 10^{-6}$ .

This yields a quintic equation that is hopeless to solve using standard techniques from Precalculus, but the root can be found graphically by seeing where the function crosses the $x-$ axis (or, in this case, the $t-$ axis):

As it turns out, the root is $t \approx 0.00997$ , so that $L_1$ is located $0.997\%$ of the distance from the Earth to the Sun in the direction of the Sun.

Lagrange Points and Polynomial Equations: Part 2

Published by John Quintanilla

One thought on “Lagrange Points and Polynomial Equations: Part 2”

Leave a comment Cancel reply

Share this:

Related

Published by John Quintanilla

One thought on “Lagrange Points and Polynomial Equations: Part 2”

Leave a comment Cancel reply