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 r4253224 = R;
        double r4253225 = 2.0;
        double r4253226 = phi1;
        double r4253227 = phi2;
        double r4253228 = r4253226 - r4253227;
        double r4253229 = r4253228 / r4253225;
        double r4253230 = sin(r4253229);
        double r4253231 = pow(r4253230, r4253225);
        double r4253232 = cos(r4253226);
        double r4253233 = cos(r4253227);
        double r4253234 = r4253232 * r4253233;
        double r4253235 = lambda1;
        double r4253236 = lambda2;
        double r4253237 = r4253235 - r4253236;
        double r4253238 = r4253237 / r4253225;
        double r4253239 = sin(r4253238);
        double r4253240 = r4253234 * r4253239;
        double r4253241 = r4253240 * r4253239;
        double r4253242 = r4253231 + r4253241;
        double r4253243 = sqrt(r4253242);
        double r4253244 = 1.0;
        double r4253245 = r4253244 - r4253242;
        double r4253246 = sqrt(r4253245);
        double r4253247 = atan2(r4253243, r4253246);
        double r4253248 = r4253225 * r4253247;
        double r4253249 = r4253224 * r4253248;
        return r4253249;
}