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 r12955225 = R;
        double r12955226 = 2.0;
        double r12955227 = phi1;
        double r12955228 = phi2;
        double r12955229 = r12955227 - r12955228;
        double r12955230 = r12955229 / r12955226;
        double r12955231 = sin(r12955230);
        double r12955232 = pow(r12955231, r12955226);
        double r12955233 = cos(r12955227);
        double r12955234 = cos(r12955228);
        double r12955235 = r12955233 * r12955234;
        double r12955236 = lambda1;
        double r12955237 = lambda2;
        double r12955238 = r12955236 - r12955237;
        double r12955239 = r12955238 / r12955226;
        double r12955240 = sin(r12955239);
        double r12955241 = r12955235 * r12955240;
        double r12955242 = r12955241 * r12955240;
        double r12955243 = r12955232 + r12955242;
        double r12955244 = sqrt(r12955243);
        double r12955245 = 1.0;
        double r12955246 = r12955245 - r12955243;
        double r12955247 = sqrt(r12955246);
        double r12955248 = atan2(r12955244, r12955247);
        double r12955249 = r12955226 * r12955248;
        double r12955250 = r12955225 * r12955249;
        return r12955250;
}