pub fn planet_position(planet: Planet, jd: f64) -> PlanetPositionExpand description
Compute heliocentric position of a planet at a given Julian Date.
Uses mean Keplerian elements with secular perturbations. Returns full 3D ecliptic coordinates (x, y, z). Accuracy: ~1° for outer planets, ~2° for inner planets over centuries.