Function apoapsis_speed 

pub fn apoapsis_speed(
    mu: f64,
    periapsis_radius: f64,
    apoapsis_radius: f64,
) -> Option<f64>

Computes speed at apoapsis using the vis-viva equation.

Returns None when the apsides are invalid or when the derived speed is invalid.