helios

Each planet's position comes from solving Kepler's equation M = E − e sin E for eccentric anomaly E (Newton's method, ~5 iterations), then rotating into ellipse coordinates. No mean orientations, no inclinations, no perturbations — the goal is the gestalt of orbital motion, not an ephemeris.

This is the JS prototype. The real version (Rust → WASM) will add outer planets, eccentric anomaly via WASM SIMD, and time-of-night sky views.