Polynomial Long Division and Megan Moroney

A brief clip from Megan Moroney’s video “I’m Not Pretty” correctly uses polynomial long division to establish that $2x+3$ is a factor of $2x^4+5x^3+7x^2+16x+15$ . Even more amazingly, the fact that the remainder is $0$ actually fits artistically with the video.

And while I have her music on my mind, I can’t resist sharing her masterpiece “Tennessee Orange” and its playful commentary on the passion of college football fans.

Bracket Symbols

Source: https://xkcd.com/2954/

Recipe Relativity

Source: https://xkcd.com/2768/

Lagrange Points and Polynomial Equations: Index

I’m doing something that I should have done a long time ago: collecting a series of posts into one single post. The links below show my series on Lagrange points.

Part 1: Introduction

Part 2: Finding $L_1$

Part 3: Finding $L_2$

Part 4: A more analytical way of finding $L_1$ and $L_2$

Part 5: Historical anecdote: Lagrange made a careless mistake!

Confirming Einstein’s General Theory of Relativity with Calculus: Index

I’m doing something that I should have done a long time ago: collecting a series of posts into one single post. The links below show my series on general relativity and the precession of Mercury’s orbit.

Part 1: Introduction

Part 1a: Introduction
Part 1b: Observed precession of Mercury
Part 1c: Outline of argument

Part 2: Precession, polar coordinates, and conic sections

Part 2a: Graphically exploring precession
Part 2b: Polar coordinates and ellipses
Part 2c: Polar coordinates, circles, and parabolas
Part 2d: Polar coordinates and hyperbolas

Part 3: Method of successive approximations

Part 4: Principles from physics

Part 4a: Angular momentum
Part 4b: Acceleration in polar coordinates
Part 4c: Newton’s Second Law and Newton’s Law of Gravitation

Part 5: Orbits under Newtonian mechanics

Part 5a: Confirmation of solution
Part 5b: Derivation with calculus
Part 5c: Derivation with differential equations and the method of undetermined coefficients
Part 5d: Derivation with differential equations and variation of parameters

Part 6: Orbits under general relativity

Part 6a: New differential equation under general relativity
Part 6b: Confirmation of solution
Part 6c: Derivation with variation of parameters
Parts 6d, 6e, 6f, 6g, 6h, 6i, 6j: Rationale for the method of undetermined coefficients
Part 6k: Derivation with undetermined coefficients

Part 7: Computing precession

Parts 7a, 7b, 7c, 7d: Predicting precession
Part 7e: Computing precession

Part 8: Second- and third-order solutions with the method of successive approximations

Part 9: Pedagogical thoughts

Earlier this year, I presented these ideas for the UNT Math Department’s Undergraduate Mathematics Colloquium Series. The video of my lecture is below.

Lagrange Points and Polynomial Equations: Part 5

This series was motivated by a terrific article that I read in the American Mathematical Monthly about Lagrange points, which are (from Wikipedia) “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 article points out a delicious historical factoid: Lagrange had a slight careless mistake in his derivation!

From the article:

Equation (d) would be just the tool to use to determine where to locate the JWST [James Webb Space Telescope, which is now in orbit about $L_2$ ], except for one thing: Lagrange got it wrong!… Do you see it? His algebra in converting $1 - \displaystyle \frac{1}{(m-1)^3}$ to common denominator form is incorrect… Fortunately, at some point in the two-and-a-half centuries between Lagrange’s work and the launch of JWST, this error has been recognized and corrected.

This little historical anecdote illustrates that, despite our best efforts, even the best of us are susceptible to careless mistakes. The simplification should have been

$q' = \displaystyle \left[ 1 - \frac{1}{(m-1)^3} \right] \cdot \frac{1}{r^3}$

$= \displaystyle \frac{(m-1)^3 - 1}{(m-1)^3} \cdot \frac{1}{r^3}$

$= \displaystyle \frac{m^3 - 3m^2 + 3m - 1 - 1}{(m-1)^3} \cdot \frac{1}{r^3}$

$= \displaystyle \frac{m^3 - 3m^2 + 3m - 2}{(m-1)^3} \cdot \frac{1}{r^3}$ .

(Parenthetically, The article also notes a clear but unintended typesetting error, as the correct but smudged exponent of 3 in the first equation became an incorrect exponent of 2 in the second.)

Lagrange Points and Polynomial Equations: Part 4

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.

As we’ve seen, the positions of $L_1$ and $L_2$ can be found by numerically solving the fifth-order polynomial equations

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

and

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

respectively. In these equations, $\mu = \displaystyle \frac{m_2}{m_1+m_2}$ where $m_1$ is the mass of the Sun and $m_2$ is the mass of Earth. Also, $t$ is the distance from the Earth to $L_1$ or $L_2$ measured as a proportion of the distance from the Sun to Earth.

We’ve also seen that, for the Sun and Earth, $mu \approx 3.00346 \times 10^{-6}$ , and numerically solving the above quintics yields $t \approx 0.000997$ for $L_1$ and $t \approx 0.01004$ for $L_2$ . In other words, $L_1$ and $L_2$ are approximately the same distance from Earth but in opposite directions.

There’s a good reason why the positive real roots of these two similar quintics are almost equal. We know that $t$ will be a lot closer to 0 than 1 because, for gravity to balance, the Lagrange points have to be a lot closer to Earth than the Sun. For this reason, the terms $\mu t^2$ and $2\mu t$ will be a lot smaller than $\mu$ , and so those two terms can be safely ignored in a first-order approximation. Also, the terms $t^5$ and $(3-\mu)t^4$ will be a lot smaller than $(3-2\mu)t^3$ , and so those two terms can also be safely ignored in a first-order approximation. Furthermore, since $\mu$ is also close to 0, the coefficient $(3-2\mu)$ can be safely replaced by just $3$ .

Consequently, the solution of both quintic equations should be close to the solution of the cubic equation

$3t^3 - \mu = 0$ ,

which is straightforward to solve:

$3t^3 = \mu$

$t^3 = \displaystyle \frac{\mu}{3}$

$t = \displaystyle \sqrt[3]{ \frac{\mu}{3} }$ .

If $\mu = 3.00346 \times 10^{-6}$ , we obtain $t \approx 0.010004$ , which is indeed reasonably close to the actual solutions for $L_1$ and $L_2$ . Indeed, this may be used as the first approximation in Newton’s method to quickly numerically evaluate the actual solutions of the two quintic polynomials.

Lagrange Points and Polynomial Equations: Part 3

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.

Lagrange Points and Polynomial Equations: Part 2

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}$ .

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 1

I recently read a terrific article in the American Mathematical Monthly about Lagrange points, which are (from Wikipedia) “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$ .

To describe these Lagrange points, I can do no better than the estimable Isaac Asimov. I quote from his essay “Colonizing the Heavens” from his book The Beginning and the End, which was published in 1977. I read the book over and over again as a boy in the mid-1980s. (Asimov’s essay originally concerned the Earth-Moon system; in the quote below, I changed the words to apply to the Sun-Earth system.)

Imagine the Sun at zenith, exactly overhead. Trace a line due eastward from the Sun down to the horizon. Two-thirds of the way along that line, one-third of the way up from the horizon, is one of those places. Trace another line westward away from the Sun down to the horizon. Two-thirds of the way along that line, one-third of the way up from the horizon, is another of those places.

Put an object in either place and it will form an equilateral triangle with the Sun and Earth…

What is so special about those places? Back in 1772, the astronomer Joseph Louis Lagrange showed that in those places any object remained stationary with respect to the Sun. As the Earth moved about the Sun, any object in either of those places would also move about the Sun in such a way as to keep perfect step with the Earth. The competing gravities of the Sun and Earth would keep it where it was. If anything happened to push it out of place it would promptly move back, wobbling back and forth a bit (“librating”) as it did so. The two places are called “Lagrangian points” or “libration points.”

Lagrange discovered five such places altogether, but three of them are of no importance since they don’t represent stable conditions. An object in those three places, once pushed out of place, would continue to drift outward and would never return.

The last paragraph of the above quote represents a rare failure of imagination by Asimov, who wrote prolifically about the future of spaceflight. Points $L_4$ and $L_5$ are indeed stable equilibria, and untold science fiction stories have placed spacecraft or colonies at these locations. (The rest of Asimov’s essay speculates about using these points in the Earth-Moon system for space colonization.) However, while the points $L_1$ and $L_2$ are unstable equilibria, they do have practical applications for spacecraft that can perform minor course corrections to stay in position. (The point $L_3$ is especially unstable to outside gravitational influences and thus seems unsuitable for spacecraft.) Again from Wikipedia,

Sun–Earth L₁ is suited for making observations of the Sun–Earth system. Objects here are never shadowed by Earth or the Moon and, if observing Earth, always view the sunlit hemisphere… Solar and heliospheric missions currently located around L₁ include the Solar and Heliospheric Observatory, Wind, Aditya-L1 Mission and the Advanced Composition Explorer. Planned missions include the Interstellar Mapping and Acceleration Probe(IMAP) and the NEO Surveyor.

Sun–Earth L₂ is a good spot for space-based observatories. Because an object around L₂ will maintain the same relative position with respect to the Sun and Earth, shielding and calibration are much simpler… The James Webb Space Telescope was positioned in a halo orbit about L₂ on January 24, 2022.

Earth–Moon L₁ allows comparatively easy access to Lunar and Earth orbits with minimal change in velocity and this has as an advantage to position a habitable space station intended to help transport cargo and personnel to the Moon and back. The SMART-1 Mission passed through the L₁ Lagrangian Point on 11 November 2004 and passed into the area dominated by the Moon’s gravitational influence.

Earth–Moon L₂ has been used for a communications satellite covering the Moon’s far side, for example, Queqiao, launched in 2018, and would be “an ideal location” for a propellant depot as part of the proposed depot-based space transportation architecture.

While the locations $L_4$ and $L_5$ are easy to describe, the precise locations of $L_1$ and $L_2$ are found by numerically solving a fifth-order polynomial equation. This was news to me when I read that article from the American Mathematical Monthly. While I had read years ago that finding the positions of the other three Lagrange points wasn’t simple, I did not realize that it was no more complicated that numerically finding the roots of a polynomial.

The above article from the American Mathematical Monthly concludes…

[t]he mathematical tools that Lagrange uses to arrive at a solution to this three-body problem lie entirely within the scope of modern courses in algebra, trigonometry, and first-semester calculus. But surely no ordinary person could have pursued the many extraordinarily complicated threads in his work to their ends, let alone woven them together into a magnificent solution to the problem as he has done. Lagrange noted in the introduction to his essay, “This research is really no more than for pure curiosity …” If only he could have watched on Christmas Day as the James Webb Space Telescope began its journey to the Lagrange point $L_2$ !

In this short series, we discuss the polynomial equations for finding $L_1$ and $L_2$ .