pub fn reduced_mass(mass_a: f64, mass_b: f64) -> Option<f64>
Computes reduced mass using μ = (m1 * m2) / (m1 + m2).
μ = (m1 * m2) / (m1 + m2)