double f(double lambda1, double lambda2, double phi1, double phi2) {
double r14104100 = lambda1;
double r14104101 = lambda2;
double r14104102 = r14104100 - r14104101;
double r14104103 = sin(r14104102);
double r14104104 = phi2;
double r14104105 = cos(r14104104);
double r14104106 = r14104103 * r14104105;
double r14104107 = phi1;
double r14104108 = cos(r14104107);
double r14104109 = sin(r14104104);
double r14104110 = r14104108 * r14104109;
double r14104111 = sin(r14104107);
double r14104112 = r14104111 * r14104105;
double r14104113 = cos(r14104102);
double r14104114 = r14104112 * r14104113;
double r14104115 = r14104110 - r14104114;
double r14104116 = atan2(r14104106, r14104115);
return r14104116;
}