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 r174067 = 1.0;
        double r174068 = s;
        double r174069 = -r174068;
        double r174070 = exp(r174069);
        double r174071 = r174067 + r174070;
        double r174072 = r174067 / r174071;
        double r174073 = c_p;
        double r174074 = pow(r174072, r174073);
        double r174075 = r174067 - r174072;
        double r174076 = c_n;
        double r174077 = pow(r174075, r174076);
        double r174078 = r174074 * r174077;
        double r174079 = t;
        double r174080 = -r174079;
        double r174081 = exp(r174080);
        double r174082 = r174067 + r174081;
        double r174083 = r174067 / r174082;
        double r174084 = pow(r174083, r174073);
        double r174085 = r174067 - r174083;
        double r174086 = pow(r174085, r174076);
        double r174087 = r174084 * r174086;
        double r174088 = r174078 / r174087;
        return r174088;
}