Lagrange Points and Polynomial Equations: Part 3

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.

The position of L_2 can be found by numerically solving the fifth-order polynomial equation

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

-(m_2+3m_3)t^2-(2m_2+3m_3)t-(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 t is the distance from the Earth to L_2 measured as a proportion of the distance from the Sun to Earth. In other words, if the distance from the Sun to Earth is 1 unit, then the distance from the Earth to L_2 is t units. The above equation is derived using principles from physics which are not elaborated upon here.

We notice that the coefficients of t^5, t^4, and t^3 are all positive, while the coefficients of t^2, t, 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+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.

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 + (3-\mu) 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.01004, so that L_2 is located 1.004\% of the distance from the Earth to the Sun in the direction away from the Sun.

One thought on “Lagrange Points and Polynomial Equations: Part 3

Leave a comment

This site uses Akismet to reduce spam. Learn how your comment data is processed.