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 r4208895 = 1.0;
        double r4208896 = s;
        double r4208897 = -r4208896;
        double r4208898 = exp(r4208897);
        double r4208899 = r4208895 + r4208898;
        double r4208900 = r4208895 / r4208899;
        double r4208901 = c_p;
        double r4208902 = pow(r4208900, r4208901);
        double r4208903 = r4208895 - r4208900;
        double r4208904 = c_n;
        double r4208905 = pow(r4208903, r4208904);
        double r4208906 = r4208902 * r4208905;
        double r4208907 = t;
        double r4208908 = -r4208907;
        double r4208909 = exp(r4208908);
        double r4208910 = r4208895 + r4208909;
        double r4208911 = r4208895 / r4208910;
        double r4208912 = pow(r4208911, r4208901);
        double r4208913 = r4208895 - r4208911;
        double r4208914 = pow(r4208913, r4208904);
        double r4208915 = r4208912 * r4208914;
        double r4208916 = r4208906 / r4208915;
        return r4208916;
}