pub fn rectifying_radius(n: f64) -> f64
Compute the rectifying radius = quarter meridian / (pi/2 * a). The accuracy of this series is the same as those used for the computation of the auxiliary latitudes.