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 r10217363 = 1.0;
        double r10217364 = s;
        double r10217365 = -r10217364;
        double r10217366 = exp(r10217365);
        double r10217367 = r10217363 + r10217366;
        double r10217368 = r10217363 / r10217367;
        double r10217369 = c_p;
        double r10217370 = pow(r10217368, r10217369);
        double r10217371 = r10217363 - r10217368;
        double r10217372 = c_n;
        double r10217373 = pow(r10217371, r10217372);
        double r10217374 = r10217370 * r10217373;
        double r10217375 = t;
        double r10217376 = -r10217375;
        double r10217377 = exp(r10217376);
        double r10217378 = r10217363 + r10217377;
        double r10217379 = r10217363 / r10217378;
        double r10217380 = pow(r10217379, r10217369);
        double r10217381 = r10217363 - r10217379;
        double r10217382 = pow(r10217381, r10217372);
        double r10217383 = r10217380 * r10217382;
        double r10217384 = r10217374 / r10217383;
        return r10217384;
}