geocentricLatitudeDifference function

double geocentricLatitudeDifference(
  1. double phi
)

Geographic latitude minus geocentric latitude (radians).

Implementation

double geocentricLatitudeDifference(double phi) {
  return secToRad(692.73 * math.sin(2 * phi) - 1.16 * math.sin(4 * phi));
}