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 r130706 = 1.0;
        double r130707 = s;
        double r130708 = -r130707;
        double r130709 = exp(r130708);
        double r130710 = r130706 + r130709;
        double r130711 = r130706 / r130710;
        double r130712 = c_p;
        double r130713 = pow(r130711, r130712);
        double r130714 = r130706 - r130711;
        double r130715 = c_n;
        double r130716 = pow(r130714, r130715);
        double r130717 = r130713 * r130716;
        double r130718 = t;
        double r130719 = -r130718;
        double r130720 = exp(r130719);
        double r130721 = r130706 + r130720;
        double r130722 = r130706 / r130721;
        double r130723 = pow(r130722, r130712);
        double r130724 = r130706 - r130722;
        double r130725 = pow(r130724, r130715);
        double r130726 = r130723 * r130725;
        double r130727 = r130717 / r130726;
        return r130727;
}