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 r120833 = 1.0;
        double r120834 = s;
        double r120835 = -r120834;
        double r120836 = exp(r120835);
        double r120837 = r120833 + r120836;
        double r120838 = r120833 / r120837;
        double r120839 = c_p;
        double r120840 = pow(r120838, r120839);
        double r120841 = r120833 - r120838;
        double r120842 = c_n;
        double r120843 = pow(r120841, r120842);
        double r120844 = r120840 * r120843;
        double r120845 = t;
        double r120846 = -r120845;
        double r120847 = exp(r120846);
        double r120848 = r120833 + r120847;
        double r120849 = r120833 / r120848;
        double r120850 = pow(r120849, r120839);
        double r120851 = r120833 - r120849;
        double r120852 = pow(r120851, r120842);
        double r120853 = r120850 * r120852;
        double r120854 = r120844 / r120853;
        return r120854;
}