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 r107209 = 1.0;
        double r107210 = s;
        double r107211 = -r107210;
        double r107212 = exp(r107211);
        double r107213 = r107209 + r107212;
        double r107214 = r107209 / r107213;
        double r107215 = c_p;
        double r107216 = pow(r107214, r107215);
        double r107217 = r107209 - r107214;
        double r107218 = c_n;
        double r107219 = pow(r107217, r107218);
        double r107220 = r107216 * r107219;
        double r107221 = t;
        double r107222 = -r107221;
        double r107223 = exp(r107222);
        double r107224 = r107209 + r107223;
        double r107225 = r107209 / r107224;
        double r107226 = pow(r107225, r107215);
        double r107227 = r107209 - r107225;
        double r107228 = pow(r107227, r107218);
        double r107229 = r107226 * r107228;
        double r107230 = r107220 / r107229;
        return r107230;
}