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 r98648 = 1.0;
        double r98649 = s;
        double r98650 = -r98649;
        double r98651 = exp(r98650);
        double r98652 = r98648 + r98651;
        double r98653 = r98648 / r98652;
        double r98654 = c_p;
        double r98655 = pow(r98653, r98654);
        double r98656 = r98648 - r98653;
        double r98657 = c_n;
        double r98658 = pow(r98656, r98657);
        double r98659 = r98655 * r98658;
        double r98660 = t;
        double r98661 = -r98660;
        double r98662 = exp(r98661);
        double r98663 = r98648 + r98662;
        double r98664 = r98648 / r98663;
        double r98665 = pow(r98664, r98654);
        double r98666 = r98648 - r98664;
        double r98667 = pow(r98666, r98657);
        double r98668 = r98665 * r98667;
        double r98669 = r98659 / r98668;
        return r98669;
}