3.2.15 · D1Orbital Mechanics & Astrodynamics

Foundations — Solving Kepler's equation — Newton-Raphson iteration

3,009 words14 min readBack to topic

Before you can read that central equation, every letter in it must mean something visual to you. This page introduces the letters one at a time, in build-up order, and only after all of them are defined do we finally write the equation down. Read top to bottom — each symbol leans on the one before.


1. The ellipse and its two secret numbers: and

An ellipse is a circle that has been stretched more in one direction than the other. Picture squashing a round balloon between two hands.

The picture: stand at the centre of the ellipse. Walk to the far end — that distance is . Turn 90° and walk to the near edge — that distance is . Because an ellipse is squashed, (for a circle, ).

Figure 1 (below) draws exactly this: a violet ellipse with the magenta arrow marking along the long axis and the orange arrow marking along the short axis, plus the centre and the off-centre focus we meet next.

Why the topic needs it: the whole orbit lives inside the box set by and . The clever trick that derives Kepler's equation is to compare the ellipse to a perfect circle of radius drawn around it — so and are the first tools we reach for.

Figure — Solving Kepler's equation — Newton-Raphson iteration

2. The focus and the offset — where the Sun sits

An ellipse has a special point called the focus. The Sun (or Earth, for a satellite) sits at one focus — not at the centre. This off-centre-ness is the source of all the difficulty.

The picture: in Figure 1 the centre is the navy dot at the geometric middle; slide toward the "fat" end by a distance (the dashed navy segment) and you reach the orange star — that is where the Sun sits. If the Sun is dead centre and the orbit is a perfect circle.

Why the topic needs it: because the Sun is offset by , the body is sometimes close to it and sometimes far. Close = fast, far = slow. That uneven speed is exactly why we cannot compute position from time with simple algebra.


3. Eccentricity — the "how squashed" dial

We now package the squash into a single dimensionless number.

Rearranging gives — the focal offset written using . Combined with the triangle law from the last section, this also yields the handy relation (substitute and solve for ).

Read the dial:

  • → focus at centre → a perfect circle (no squash).
  • near (e.g. for Earth) → almost a circle.
  • near (e.g. for some comets) → a long, thin, cigar-shaped orbit.
  • → the ellipse breaks open into a parabola (not an orbit anymore).

4. Angles measured in radians — the unit that makes addable

The topic's central equation will add an angle to a sine of that angle. For that to make sense, the angle must be measured in radians, not degrees.

The picture: take a circle of radius . Lay a piece of string of length along its rim. The angle from the centre spanning that string is exactly one radian ().

Why the topic needs it: all three anomalies we meet next are radian angles, and every / in the topic assumes radians.


5. The three anomalies: , , — three ways to say "where"

"Anomaly" is an old astronomy word that just means angle telling you where the body is. There are three of them because there are three natural reference points.

Figure 2 (below) shows all three at once: the magenta ellipse (real orbit) inside its violet dashed auxiliary circle. The body is the magenta dot; its "shadow" pushed straight up onto the circle is the violet dot. The violet angle is measured at the navy centre, and the orange angle is measured at the orange Sun — you can literally see them pointing in slightly different directions.

Figure — Solving Kepler's equation — Newton-Raphson iteration

Why the topic needs it: the mission is time position, i.e. . is easy from time, is easy from , but is the transcendental knot the whole topic exists to untie.


6. Mass, gravity, and the parameter — the strength of the pull

Before we can turn time into an angle, we need one number that captures how strongly the central body pulls.

The picture: think of as a single "gravity dial" for the Sun. You never need and separately in orbit work — they always travel together as this one product, so astronomers bundle them into once and for all.

Why the topic needs it: the steady rate at which the mean angle grows (defined next) depends on how strong gravity is, so must be on the table before we can write that rate.


7. Mean motion , period , and periapsis time — turning time into

The picture: if a full lap ( radians) takes time , then each second earns radians. That rate is . From it, the elapsed radians since the body last passed closest approach at time . More in Mean Motion and Orbital Period.

The link between , and hidden in is exactly Kepler's third law.


8. The moving tools: , , the derivative , and the root

The final ingredients are the mathematical machines the topic runs.

Why enters Kepler? In the derivation the triangle from centre-to-focus-to-circle-point has height . That height is what pulls away from the steady — so the correction term ends up being .


9. Assembling the central equation and its root function

Now that , and all have pictures, we may finally write the topic's central relation:

We cannot rearrange this to get alone, so we hunt for the that makes it balance. Move everything to one side:

Why the derivative at all? Newton-Raphson uses this slope to predict where the curve will cross zero: follow the tangent line down to the axis and land near the root. Figure 3 (below) shows exactly this — the violet curve , its magenta root, an orange guess, and the dashed orange tangent whose slope is sliding the guess toward the crossing. Compare this with the slower slope-free cousin in Fixed-Point Iteration for Kepler.

Figure — Solving Kepler's equation — Newton-Raphson iteration

How it all feeds the topic

ellipse axes a and b

c squared = a squared minus b squared

focus offset c

eccentricity e

radians

mean anomaly M

eccentric anomaly E

true anomaly nu

periapsis time t_p

gravitational parameter mu

period T and mean motion n

Kepler equation M = E minus e sin E

sine and cosine on circle

root function f of E

derivative f prime

Newton-Raphson update


Equipment checklist

Test yourself — cover the right side and answer out loud.

What do and measure on an ellipse?
Half the longest () and half the shortest () diameters — the semi-major and semi-minor axes.
How are , , and related?
By the right-triangle law ; the focus offset, minor axis, and major axis form a right triangle with as hypotenuse.
Where does the Sun sit on the orbit?
At one focus, offset from the centre by — never at the centre.
What does the eccentricity tell you?
How squashed the orbit is: is a circle, near is a long thin ellipse; it scales every correction term.
Why must every angle be in radians?
Because Kepler's equation adds a raw angle to a sine, and only in radians does an angle behave like a plain addable number.
What is the difference between , , and ?
= steady fake-body angle (a proxy for time); = angle to the auxiliary-circle shadow at the centre; = the real angle at the Sun.
Why can plain give the wrong , and how do you fix it?
Because only returns and cannot tell top half from bottom; use atan2 or carry the sign of so and stay in the same half-plane.
What is and why is it there?
The clock reading at closest approach; measuring makes exactly at perihelion.
What is the gravitational parameter ?
The product of Newton's constant and the central mass — a single "gravity dial" that sets how fast the body must move.
How do you get from a time ?
with mean motion .
On a unit circle, what are and ?
The height and the rightward reach of the point at angle ; each lies in .
What is a "root" of , and what is ?
A root is an input making ; is the local slope of , used by Newton to aim at the root.

Ready? Head back to the main topic and every symbol should now feel like an old friend.