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 r141451 = 1.0;
        double r141452 = s;
        double r141453 = -r141452;
        double r141454 = exp(r141453);
        double r141455 = r141451 + r141454;
        double r141456 = r141451 / r141455;
        double r141457 = c_p;
        double r141458 = pow(r141456, r141457);
        double r141459 = r141451 - r141456;
        double r141460 = c_n;
        double r141461 = pow(r141459, r141460);
        double r141462 = r141458 * r141461;
        double r141463 = t;
        double r141464 = -r141463;
        double r141465 = exp(r141464);
        double r141466 = r141451 + r141465;
        double r141467 = r141451 / r141466;
        double r141468 = pow(r141467, r141457);
        double r141469 = r141451 - r141467;
        double r141470 = pow(r141469, r141460);
        double r141471 = r141468 * r141470;
        double r141472 = r141462 / r141471;
        return r141472;
}