pub fn vis_viva(mu: Mu, r: Km, a: Km) -> KmPerSec
Vis-viva equation: v = sqrt(μ * (2/r - 1/a))
Calculates orbital velocity at distance r from the central body, for an orbit with semi-major axis a.