docs: private static double polarTriangleArea(double tan1, double lng1, double tan2, double lng2) {
          double deltaLng = lng1 - lng2;
          double t = tan1 * tan2;
          return 2 * atan2(t * sin(deltaLng), 1 + t * cos(deltaLng));
      }