pub fn circular_orbit_state(mu: Mu, radius: Km) -> PropState
Create initial state for a circular orbit in the XY plane.
Position: (r, 0, 0) Velocity: (0, v_circ, 0) where v_circ = sqrt(μ/r)