Tidal torque locks a planet's rotation to its mean orbital motion, so its spin rate is constant: the body-fixed reference angle advances as the mean anomaly M = n·t, where n = 2π/T is the mean motion. But Kepler's second law says the planet does not sweep out true orbital angle at a constant rate on an eccentric orbit — it moves faster near periapsis and slower near apoapsis. The actual direction to the star, the true anomaly ν, is found from the eccentric anomaly E via Kepler's equation:
M = E − e·sin(E) (Kepler's equation, solved for E by Newton iteration)
tan(ν/2) = √((1+e)/(1−e)) · tan(E/2)
Because the body's spin angle tracks M exactly while the true sub-stellar direction tracks ν, the two drift apart and back together once per orbit — the forced libration in longitude, ψ(t) = ν(t) − M(t). This 2D view drops the spatial orbit picture entirely and instead plots the two angles themselves: a dual clock dial where the red hand (M) turns at perfectly constant speed and the green hand (ν) visibly speeds up and slows down, a scrolling strip chart of ψ(t) so the oscillation period and amplitude are readable directly, and a closed polar loop of ψ against M that traces out a bigger figure-eight as e grows.
ψ(t) ≈ 2e·sin(M) + (5/4)e²·sin(2M) + … (leading-order libration series)
- Eccentricity slider — sets e directly; the loop and strip-chart amplitude update live.
- Time speed — scales simulated time; the green hand visibly races ahead near periapsis.
- Clear trace — wipes the accumulated polar loop and strip chart so a new orbit builds up from scratch.