Cannot sample enough valid points. (more)

\[0.0 \lt c_p \land 0.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 r141500 = 1.0;
        double r141501 = s;
        double r141502 = -r141501;
        double r141503 = exp(r141502);
        double r141504 = r141500 + r141503;
        double r141505 = r141500 / r141504;
        double r141506 = c_p;
        double r141507 = pow(r141505, r141506);
        double r141508 = r141500 - r141505;
        double r141509 = c_n;
        double r141510 = pow(r141508, r141509);
        double r141511 = r141507 * r141510;
        double r141512 = t;
        double r141513 = -r141512;
        double r141514 = exp(r141513);
        double r141515 = r141500 + r141514;
        double r141516 = r141500 / r141515;
        double r141517 = pow(r141516, r141506);
        double r141518 = r141500 - r141516;
        double r141519 = pow(r141518, r141509);
        double r141520 = r141517 * r141519;
        double r141521 = r141511 / r141520;
        return r141521;
}