pub fn ultimate_floor(alpha: f64, beta: f64, delta: f64, mu: f64) -> Result<f64>
ISS energy floor V_inf = delta^2 / (2 (alpha - beta)^2 mu).
V_inf = delta^2 / (2 (alpha - beta)^2 mu)
Requires alpha > beta and mu > 0 (PSP-8 System 2).
alpha > beta
mu > 0