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 r113752 = 1.0;
        double r113753 = s;
        double r113754 = -r113753;
        double r113755 = exp(r113754);
        double r113756 = r113752 + r113755;
        double r113757 = r113752 / r113756;
        double r113758 = c_p;
        double r113759 = pow(r113757, r113758);
        double r113760 = r113752 - r113757;
        double r113761 = c_n;
        double r113762 = pow(r113760, r113761);
        double r113763 = r113759 * r113762;
        double r113764 = t;
        double r113765 = -r113764;
        double r113766 = exp(r113765);
        double r113767 = r113752 + r113766;
        double r113768 = r113752 / r113767;
        double r113769 = pow(r113768, r113758);
        double r113770 = r113752 - r113768;
        double r113771 = pow(r113770, r113761);
        double r113772 = r113769 * r113771;
        double r113773 = r113763 / r113772;
        return r113773;
}