pub fn mlfn(phi: f64, sphi: f64, cphi: f64, en: &[f64]) -> f64
meridional distance for ellipsoid and inverse using 6th-order expansion in the third flattening n.