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 r100661 = 1.0;
        double r100662 = s;
        double r100663 = -r100662;
        double r100664 = exp(r100663);
        double r100665 = r100661 + r100664;
        double r100666 = r100661 / r100665;
        double r100667 = c_p;
        double r100668 = pow(r100666, r100667);
        double r100669 = r100661 - r100666;
        double r100670 = c_n;
        double r100671 = pow(r100669, r100670);
        double r100672 = r100668 * r100671;
        double r100673 = t;
        double r100674 = -r100673;
        double r100675 = exp(r100674);
        double r100676 = r100661 + r100675;
        double r100677 = r100661 / r100676;
        double r100678 = pow(r100677, r100667);
        double r100679 = r100661 - r100677;
        double r100680 = pow(r100679, r100670);
        double r100681 = r100678 * r100680;
        double r100682 = r100672 / r100681;
        return r100682;
}