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 r5859976 = 1.0;
        double r5859977 = s;
        double r5859978 = -r5859977;
        double r5859979 = exp(r5859978);
        double r5859980 = r5859976 + r5859979;
        double r5859981 = r5859976 / r5859980;
        double r5859982 = c_p;
        double r5859983 = pow(r5859981, r5859982);
        double r5859984 = r5859976 - r5859981;
        double r5859985 = c_n;
        double r5859986 = pow(r5859984, r5859985);
        double r5859987 = r5859983 * r5859986;
        double r5859988 = t;
        double r5859989 = -r5859988;
        double r5859990 = exp(r5859989);
        double r5859991 = r5859976 + r5859990;
        double r5859992 = r5859976 / r5859991;
        double r5859993 = pow(r5859992, r5859982);
        double r5859994 = r5859976 - r5859992;
        double r5859995 = pow(r5859994, r5859985);
        double r5859996 = r5859993 * r5859995;
        double r5859997 = r5859987 / r5859996;
        return r5859997;
}