Cannot sample enough valid points. (more)

\[R \cdot \left(2 \cdot \tan^{-1}_* \frac{\sqrt{{\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)}}{\sqrt{1 - \left({\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right)}}\right)\]
R \cdot \left(2 \cdot \tan^{-1}_* \frac{\sqrt{{\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)}}{\sqrt{1 - \left({\left(\sin \left(\frac{\phi_1 - \phi_2}{2}\right)\right)}^{2} + \left(\left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right) \cdot \sin \left(\frac{\lambda_1 - \lambda_2}{2}\right)\right)}}\right)
double f(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r1733530 = R;
        double r1733531 = 2.0;
        double r1733532 = phi1;
        double r1733533 = phi2;
        double r1733534 = r1733532 - r1733533;
        double r1733535 = r1733534 / r1733531;
        double r1733536 = sin(r1733535);
        double r1733537 = pow(r1733536, r1733531);
        double r1733538 = cos(r1733532);
        double r1733539 = cos(r1733533);
        double r1733540 = r1733538 * r1733539;
        double r1733541 = lambda1;
        double r1733542 = lambda2;
        double r1733543 = r1733541 - r1733542;
        double r1733544 = r1733543 / r1733531;
        double r1733545 = sin(r1733544);
        double r1733546 = r1733540 * r1733545;
        double r1733547 = r1733546 * r1733545;
        double r1733548 = r1733537 + r1733547;
        double r1733549 = sqrt(r1733548);
        double r1733550 = 1.0;
        double r1733551 = r1733550 - r1733548;
        double r1733552 = sqrt(r1733551);
        double r1733553 = atan2(r1733549, r1733552);
        double r1733554 = r1733531 * r1733553;
        double r1733555 = r1733530 * r1733554;
        return r1733555;
}