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 r142262 = 1.0;
        double r142263 = s;
        double r142264 = -r142263;
        double r142265 = exp(r142264);
        double r142266 = r142262 + r142265;
        double r142267 = r142262 / r142266;
        double r142268 = c_p;
        double r142269 = pow(r142267, r142268);
        double r142270 = r142262 - r142267;
        double r142271 = c_n;
        double r142272 = pow(r142270, r142271);
        double r142273 = r142269 * r142272;
        double r142274 = t;
        double r142275 = -r142274;
        double r142276 = exp(r142275);
        double r142277 = r142262 + r142276;
        double r142278 = r142262 / r142277;
        double r142279 = pow(r142278, r142268);
        double r142280 = r142262 - r142278;
        double r142281 = pow(r142280, r142271);
        double r142282 = r142279 * r142281;
        double r142283 = r142273 / r142282;
        return r142283;
}