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 r8853297 = 1.0;
        double r8853298 = s;
        double r8853299 = -r8853298;
        double r8853300 = exp(r8853299);
        double r8853301 = r8853297 + r8853300;
        double r8853302 = r8853297 / r8853301;
        double r8853303 = c_p;
        double r8853304 = pow(r8853302, r8853303);
        double r8853305 = r8853297 - r8853302;
        double r8853306 = c_n;
        double r8853307 = pow(r8853305, r8853306);
        double r8853308 = r8853304 * r8853307;
        double r8853309 = t;
        double r8853310 = -r8853309;
        double r8853311 = exp(r8853310);
        double r8853312 = r8853297 + r8853311;
        double r8853313 = r8853297 / r8853312;
        double r8853314 = pow(r8853313, r8853303);
        double r8853315 = r8853297 - r8853313;
        double r8853316 = pow(r8853315, r8853306);
        double r8853317 = r8853314 * r8853316;
        double r8853318 = r8853308 / r8853317;
        return r8853318;
}