double f(double lambda1, double lambda2, double phi1, double phi2) {
double r11227349 = lambda1;
double r11227350 = lambda2;
double r11227351 = r11227349 - r11227350;
double r11227352 = sin(r11227351);
double r11227353 = phi2;
double r11227354 = cos(r11227353);
double r11227355 = r11227352 * r11227354;
double r11227356 = phi1;
double r11227357 = cos(r11227356);
double r11227358 = sin(r11227353);
double r11227359 = r11227357 * r11227358;
double r11227360 = sin(r11227356);
double r11227361 = r11227360 * r11227354;
double r11227362 = cos(r11227351);
double r11227363 = r11227361 * r11227362;
double r11227364 = r11227359 - r11227363;
double r11227365 = atan2(r11227355, r11227364);
return r11227365;
}