Cannot sample enough valid points. (more)

\[0 \lt c_p \land 0 \lt c_n\]
\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}
double f(double c_p, double c_n, double t, double s) {
        double r5822445 = 1.0;
        double r5822446 = s;
        double r5822447 = -r5822446;
        double r5822448 = exp(r5822447);
        double r5822449 = r5822445 + r5822448;
        double r5822450 = r5822445 / r5822449;
        double r5822451 = c_p;
        double r5822452 = pow(r5822450, r5822451);
        double r5822453 = r5822445 - r5822450;
        double r5822454 = c_n;
        double r5822455 = pow(r5822453, r5822454);
        double r5822456 = r5822452 * r5822455;
        double r5822457 = t;
        double r5822458 = -r5822457;
        double r5822459 = exp(r5822458);
        double r5822460 = r5822445 + r5822459;
        double r5822461 = r5822445 / r5822460;
        double r5822462 = pow(r5822461, r5822451);
        double r5822463 = r5822445 - r5822461;
        double r5822464 = pow(r5822463, r5822454);
        double r5822465 = r5822462 * r5822464;
        double r5822466 = r5822456 / r5822465;
        return r5822466;
}