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 r173708 = 1.0;
        double r173709 = s;
        double r173710 = -r173709;
        double r173711 = exp(r173710);
        double r173712 = r173708 + r173711;
        double r173713 = r173708 / r173712;
        double r173714 = c_p;
        double r173715 = pow(r173713, r173714);
        double r173716 = r173708 - r173713;
        double r173717 = c_n;
        double r173718 = pow(r173716, r173717);
        double r173719 = r173715 * r173718;
        double r173720 = t;
        double r173721 = -r173720;
        double r173722 = exp(r173721);
        double r173723 = r173708 + r173722;
        double r173724 = r173708 / r173723;
        double r173725 = pow(r173724, r173714);
        double r173726 = r173708 - r173724;
        double r173727 = pow(r173726, r173717);
        double r173728 = r173725 * r173727;
        double r173729 = r173719 / r173728;
        return r173729;
}