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 r8863671 = 1.0;
        double r8863672 = s;
        double r8863673 = -r8863672;
        double r8863674 = exp(r8863673);
        double r8863675 = r8863671 + r8863674;
        double r8863676 = r8863671 / r8863675;
        double r8863677 = c_p;
        double r8863678 = pow(r8863676, r8863677);
        double r8863679 = r8863671 - r8863676;
        double r8863680 = c_n;
        double r8863681 = pow(r8863679, r8863680);
        double r8863682 = r8863678 * r8863681;
        double r8863683 = t;
        double r8863684 = -r8863683;
        double r8863685 = exp(r8863684);
        double r8863686 = r8863671 + r8863685;
        double r8863687 = r8863671 / r8863686;
        double r8863688 = pow(r8863687, r8863677);
        double r8863689 = r8863671 - r8863687;
        double r8863690 = pow(r8863689, r8863680);
        double r8863691 = r8863688 * r8863690;
        double r8863692 = r8863682 / r8863691;
        return r8863692;
}