1 2 3 4
#[inline] pub(crate) fn msfn(sinphi: f64, cosphi: f64, es: f64) -> f64 { cosphi / (1. - es * sinphi * sinphi).sqrt() }