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 r8624049 = 1.0;
        double r8624050 = s;
        double r8624051 = -r8624050;
        double r8624052 = exp(r8624051);
        double r8624053 = r8624049 + r8624052;
        double r8624054 = r8624049 / r8624053;
        double r8624055 = c_p;
        double r8624056 = pow(r8624054, r8624055);
        double r8624057 = r8624049 - r8624054;
        double r8624058 = c_n;
        double r8624059 = pow(r8624057, r8624058);
        double r8624060 = r8624056 * r8624059;
        double r8624061 = t;
        double r8624062 = -r8624061;
        double r8624063 = exp(r8624062);
        double r8624064 = r8624049 + r8624063;
        double r8624065 = r8624049 / r8624064;
        double r8624066 = pow(r8624065, r8624055);
        double r8624067 = r8624049 - r8624065;
        double r8624068 = pow(r8624067, r8624058);
        double r8624069 = r8624066 * r8624068;
        double r8624070 = r8624060 / r8624069;
        return r8624070;
}