Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r5818342 = 1.0;
        double r5818343 = s;
        double r5818344 = -r5818343;
        double r5818345 = exp(r5818344);
        double r5818346 = r5818342 + r5818345;
        double r5818347 = r5818342 / r5818346;
        double r5818348 = c_p;
        double r5818349 = pow(r5818347, r5818348);
        double r5818350 = r5818342 - r5818347;
        double r5818351 = c_n;
        double r5818352 = pow(r5818350, r5818351);
        double r5818353 = r5818349 * r5818352;
        double r5818354 = t;
        double r5818355 = -r5818354;
        double r5818356 = exp(r5818355);
        double r5818357 = r5818342 + r5818356;
        double r5818358 = r5818342 / r5818357;
        double r5818359 = pow(r5818358, r5818348);
        double r5818360 = r5818342 - r5818358;
        double r5818361 = pow(r5818360, r5818351);
        double r5818362 = r5818359 * r5818361;
        double r5818363 = r5818353 / r5818362;
        return r5818363;
}