pub fn enfn(n: f64) -> Vec<f64>
meridional distance for ellipsoid and inverse using 6th-order expansion in the third flattening n. This gives full double precision accuracy for |f| <= 1/150.