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 r8543641 = 1.0;
        double r8543642 = s;
        double r8543643 = -r8543642;
        double r8543644 = exp(r8543643);
        double r8543645 = r8543641 + r8543644;
        double r8543646 = r8543641 / r8543645;
        double r8543647 = c_p;
        double r8543648 = pow(r8543646, r8543647);
        double r8543649 = r8543641 - r8543646;
        double r8543650 = c_n;
        double r8543651 = pow(r8543649, r8543650);
        double r8543652 = r8543648 * r8543651;
        double r8543653 = t;
        double r8543654 = -r8543653;
        double r8543655 = exp(r8543654);
        double r8543656 = r8543641 + r8543655;
        double r8543657 = r8543641 / r8543656;
        double r8543658 = pow(r8543657, r8543647);
        double r8543659 = r8543641 - r8543657;
        double r8543660 = pow(r8543659, r8543650);
        double r8543661 = r8543658 * r8543660;
        double r8543662 = r8543652 / r8543661;
        return r8543662;
}