// Reduced or Parametric Latitude
   Latitude Ellipsoid::
Eta(const Latitude & lat) const {
   // NOTE: sqrt(1-e^2) == b/a
   return Latitude(atan(b/a*tan(lat*deg2rad))*rad2deg);
}  // Eta


// Radius of Ellipsoid at Latitude.
   double Ellipsoid::
Radius(const Latitude & lat) const {
   Latitude eta = Eta(lat);
   double x = a * cos(eta*deg2rad);
   double y = b * sin(eta*deg2rad);
   return sqrt(x*x + y*y);
}  // Radius