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 r98500 = 1.0;
        double r98501 = s;
        double r98502 = -r98501;
        double r98503 = exp(r98502);
        double r98504 = r98500 + r98503;
        double r98505 = r98500 / r98504;
        double r98506 = c_p;
        double r98507 = pow(r98505, r98506);
        double r98508 = r98500 - r98505;
        double r98509 = c_n;
        double r98510 = pow(r98508, r98509);
        double r98511 = r98507 * r98510;
        double r98512 = t;
        double r98513 = -r98512;
        double r98514 = exp(r98513);
        double r98515 = r98500 + r98514;
        double r98516 = r98500 / r98515;
        double r98517 = pow(r98516, r98506);
        double r98518 = r98500 - r98516;
        double r98519 = pow(r98518, r98509);
        double r98520 = r98517 * r98519;
        double r98521 = r98511 / r98520;
        return r98521;
}