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 r4580230 = 1.0;
        double r4580231 = s;
        double r4580232 = -r4580231;
        double r4580233 = exp(r4580232);
        double r4580234 = r4580230 + r4580233;
        double r4580235 = r4580230 / r4580234;
        double r4580236 = c_p;
        double r4580237 = pow(r4580235, r4580236);
        double r4580238 = r4580230 - r4580235;
        double r4580239 = c_n;
        double r4580240 = pow(r4580238, r4580239);
        double r4580241 = r4580237 * r4580240;
        double r4580242 = t;
        double r4580243 = -r4580242;
        double r4580244 = exp(r4580243);
        double r4580245 = r4580230 + r4580244;
        double r4580246 = r4580230 / r4580245;
        double r4580247 = pow(r4580246, r4580236);
        double r4580248 = r4580230 - r4580246;
        double r4580249 = pow(r4580248, r4580239);
        double r4580250 = r4580247 * r4580249;
        double r4580251 = r4580241 / r4580250;
        return r4580251;
}