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 r204663 = 1.0;
        double r204664 = s;
        double r204665 = -r204664;
        double r204666 = exp(r204665);
        double r204667 = r204663 + r204666;
        double r204668 = r204663 / r204667;
        double r204669 = c_p;
        double r204670 = pow(r204668, r204669);
        double r204671 = r204663 - r204668;
        double r204672 = c_n;
        double r204673 = pow(r204671, r204672);
        double r204674 = r204670 * r204673;
        double r204675 = t;
        double r204676 = -r204675;
        double r204677 = exp(r204676);
        double r204678 = r204663 + r204677;
        double r204679 = r204663 / r204678;
        double r204680 = pow(r204679, r204669);
        double r204681 = r204663 - r204679;
        double r204682 = pow(r204681, r204672);
        double r204683 = r204680 * r204682;
        double r204684 = r204674 / r204683;
        return r204684;
}