pub fn energy_drift(states: &[PropState], mu: Mu) -> (f64, f64)
Compute maximum relative energy drift during a propagation.
Returns (max_relative_error, final_relative_error). Useful for validating integrator accuracy.