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 r5676597 = R;
        double r5676598 = 2.0;
        double r5676599 = phi1;
        double r5676600 = phi2;
        double r5676601 = r5676599 - r5676600;
        double r5676602 = r5676601 / r5676598;
        double r5676603 = sin(r5676602);
        double r5676604 = pow(r5676603, r5676598);
        double r5676605 = cos(r5676599);
        double r5676606 = cos(r5676600);
        double r5676607 = r5676605 * r5676606;
        double r5676608 = lambda1;
        double r5676609 = lambda2;
        double r5676610 = r5676608 - r5676609;
        double r5676611 = r5676610 / r5676598;
        double r5676612 = sin(r5676611);
        double r5676613 = r5676607 * r5676612;
        double r5676614 = r5676613 * r5676612;
        double r5676615 = r5676604 + r5676614;
        double r5676616 = sqrt(r5676615);
        double r5676617 = 1.0;
        double r5676618 = r5676617 - r5676615;
        double r5676619 = sqrt(r5676618);
        double r5676620 = atan2(r5676616, r5676619);
        double r5676621 = r5676598 * r5676620;
        double r5676622 = r5676597 * r5676621;
        return r5676622;
}