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 r181704 = 1.0;
        double r181705 = s;
        double r181706 = -r181705;
        double r181707 = exp(r181706);
        double r181708 = r181704 + r181707;
        double r181709 = r181704 / r181708;
        double r181710 = c_p;
        double r181711 = pow(r181709, r181710);
        double r181712 = r181704 - r181709;
        double r181713 = c_n;
        double r181714 = pow(r181712, r181713);
        double r181715 = r181711 * r181714;
        double r181716 = t;
        double r181717 = -r181716;
        double r181718 = exp(r181717);
        double r181719 = r181704 + r181718;
        double r181720 = r181704 / r181719;
        double r181721 = pow(r181720, r181710);
        double r181722 = r181704 - r181720;
        double r181723 = pow(r181722, r181713);
        double r181724 = r181721 * r181723;
        double r181725 = r181715 / r181724;
        return r181725;
}