pub fn solve_kepler_equation(m_anom: f64, e: f64) -> f64
Solve Kepler’s equation M = E - e*sin(E) for eccentric anomaly E.