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 r8073333 = 1.0;
        double r8073334 = s;
        double r8073335 = -r8073334;
        double r8073336 = exp(r8073335);
        double r8073337 = r8073333 + r8073336;
        double r8073338 = r8073333 / r8073337;
        double r8073339 = c_p;
        double r8073340 = pow(r8073338, r8073339);
        double r8073341 = r8073333 - r8073338;
        double r8073342 = c_n;
        double r8073343 = pow(r8073341, r8073342);
        double r8073344 = r8073340 * r8073343;
        double r8073345 = t;
        double r8073346 = -r8073345;
        double r8073347 = exp(r8073346);
        double r8073348 = r8073333 + r8073347;
        double r8073349 = r8073333 / r8073348;
        double r8073350 = pow(r8073349, r8073339);
        double r8073351 = r8073333 - r8073349;
        double r8073352 = pow(r8073351, r8073342);
        double r8073353 = r8073350 * r8073352;
        double r8073354 = r8073344 / r8073353;
        return r8073354;
}