Back

Lagrange points of the circular restricted three-body problem

What you are seeing: two heavy bodies on a circular orbit (yellow: the primary, like the Sun; blue: the secondary, like Jupiter or the Moon), in a frame that rotates with them so they sit still. A small test particle (red) is dropped at a chosen starting location with a chosen velocity and the playground integrates its motion. Trails fade with age. The five labeled black dots are the Lagrange points: locations where the test particle can in principle sit still in the rotating frame.

Equations of motion in the synodic frame : x¨2y˙x=(1μ)(x+μ)r13μ(x1+μ)r23\ddot x - 2\dot y - x = -\frac{(1 - \mu)(x + \mu)}{r_1^3} - \frac{\mu (x - 1 + \mu)}{r_2^3}, y¨+2x˙y=(1μ)yr13μyr23\ddot y + 2\dot x - y = -\frac{(1 - \mu) y}{r_1^3} - \frac{\mu y}{r_2^3}. Mass parameter μ=m2/(m1+m2)\mu = m_2 / (m_1 + m_2). L1, L2, L3 are on the xx-axis and always unstable. L4 and L5 form equilateral triangles with the primaries; they are linearly stable if μ<μR0.0385\mu \lt \mu_R \approx 0.0385 (Routh's stability criterion). The famous Trojan asteroid clouds live around Jupiter's L4 and L5 (μ103\mu \approx 10^{-3}).

Click anywhere on the plot to drop a test particle there at zero velocity. Use the buttons to drop one exactly at L4 or L5 with a tiny perturbation.

Figure 1. Lagrange points and test-particle orbits in the synodic frame of the circular restricted three-body problem. Method: velocity-Verlet with predictor-corrector for the Coriolis term, dt = 0.002.
mu0.01215
speed1.0

WHAT TO TRY

  • Vary each control and watch the rail readouts respond.
  • Compare the diagnostic plot against the live scene.