3.2.29 · D1Orbital Mechanics & Astrodynamics

Foundations — Gauss's method for Lambert's problem

2,337 words11 min readBack to topic

This page is the toolbox. Before you can read Lambert's problem or the parent Gauss's method note, you need every symbol below. We build each from nothing: plain words → a picture → why the topic needs it. Do not skip; every later symbol leans on an earlier one.


1. A position vector — "an arrow to where something is"

The central body sits at a special point called the focus. Every orbit we discuss has the central body at one focus of its ellipse — that is Kepler's first law, and we simply take it as given here.

Figure — Gauss's method for Lambert's problem
  • The picture: look at figure s01. The orange dot is the focus. Two magenta arrows, and , reach out to two positions of the spacecraft — an earlier one and a later one.
  • Why the topic needs it: Lambert's problem is stated in terms of these two arrows. They are the known data. Everything else is squeezed out of them.

2. The angle between the two arrows — the transfer angle

The symbol (Greek "nu") is the true anomaly — the angle of the spacecraft measured around its orbit. So (the change in true anomaly) and are the same thing: the swept angle.

  • The picture: in figure s01 the violet wedge between the two magenta arrows is .
  • Why the topic needs it: tells you how "spread out" the two points are around the central body. A tiny means they are almost in line; a large means the spacecraft went most of the way round.

3. The dot product — how we compute from the arrows

We know the two arrows, but not yet the angle between them. The tool that extracts an angle from two arrows is the dot product.

The (cosine) here is the cosine function: for an angle, it is a number between and that says how much two arrows "agree" in direction. (same direction), (perpendicular), (opposite).


4. The chord and the triangle

Draw a straight line directly connecting the two spacecraft positions. That straight segment is the chord .

Figure — Gauss's method for Lambert's problem
  • The picture: figure s02. The focus, tip of , and tip of form a triangle. The two arrows are two sides; the chord is the third side (the orange straight line).
  • Why we need : it is a pure-geometry summary of "how far apart the two endpoints are." The chord appears in every serious Lambert solver.
  • Why the topic needs the triangle: it is the flat stand-in for the curved region the planet actually sweeps. Gauss compares the curved region to this flat one.

5. Sector vs triangle — the picture behind the ratio

The planet does not travel along the straight chord. It travels along the curved orbit. So the region it truly sweeps out is bounded by the two arrows and the curve — a "pie slice." That is the elliptic sector.

Figure — Gauss's method for Lambert's problem
  • The picture: figure s03. The violet region is the true elliptic sector (bounded by the curved orbit). The magenta dashed region is the flat triangle (bounded by the chord). The sector is always a little fatter than the triangle because the curve bulges outward.
  • Why we need both: Gauss's whole trick is one number, When the orbit barely bends (small ), the curve is almost the chord, sector triangle, so . That is the sensible first guess the algorithm starts from.

6. Kepler's second law — why area equals time

Why does area matter at all? Because of Kepler's second law:

Here is the specific angular momentum — a fixed number for a given orbit measuring how fast area is swept. It relates to orbit shape through (both symbols defined next).


7. — the strength of gravity

  • Why we need it: the shape of an orbit and the speed along it both depend on how hard the central body pulls. In canonical units the parent note sets to keep the arithmetic clean.
  • Where it enters: inside the time constant — bigger gravity pulls the planet round faster, so it changes the time budget.

8. — the semi-latus rectum (the orbit's "width" number)

  • Why we need it: once Gauss's iteration lands on the right geometry, is the number we extract, and from come the Lagrange f and g functions that hand back the velocities. It is the "answer knob" of the whole method.

9. Eccentric anomaly and the auxiliary variables ,

The true anomaly from section 2 is the real, physical angle. But solving the timing cleanly needs a different bookkeeping angle, the eccentric anomaly — the angle you'd measure if you stretched the ellipse into a circle. It appears in Kepler's equation, the equation that links angle to time on an orbit.

  • The picture: back in figure s03, a bigger bend of the curve = a bigger = a fatter sector = a bigger . They rise together.

10. The known geometry constants and

These two are just pre-computed numbers built from the data, so the iteration doesn't recompute them each loop.

Notice in both denominators. When (the two points nearly opposite the focus), and — so and blow up to infinity. This is exactly where Gauss's method fails, and why near-antipodal transfers need the Universal variable formulation or the Izzo Lambert solver instead.


Prerequisite map

Position vectors r1 and r2

Dot product gives cos theta

Transfer angle theta

Law of cosines gives chord c

Triangle area

Sector vs triangle ratio y

Keplers second law

Time equals area over h

Gravity parameter mu

Constants ell and m

Two Gauss equations

Eccentric anomaly and x

Semi latus rectum p

Lagrange f and g give velocities


Equipment checklist

Test yourself — cover the right side and see if each item is second nature before you read the parent note.

What does a bold mean versus a plain ?
Bold = the arrow from focus to object (has direction); plain = just its length (a number).
Which tool extracts the angle between two arrows, and what does it return?
The dot product; , from which and then follow.
Why can't alone give the transfer angle?
It only returns , so it can't distinguish short-way from long-way; you need a prograde/retrograde flag.
What is the chord and which law computes it?
The straight-line distance between the two positions; the law of cosines, .
What does the ratio compare, and what makes it near 1?
Curved elliptic sector area ÷ flat triangle area; near 1 when the orbit barely bends (small ).
Why is area the same as time here?
Kepler's second law: the radius sweeps equal areas in equal times, so swept area .
What is and where does it enter?
The gravitational parameter ; it sits inside the time constant .
Does enter linearly or squared?
Squared — matching Kepler's third-law scaling.
What happens to and as , and why?
They blow up because in their denominators; Gauss's method fails there.
What is used for once found?
It fixes the orbit's shape and feeds the Lagrange functions that recover .