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 r3202053 = 1.0;
        double r3202054 = s;
        double r3202055 = -r3202054;
        double r3202056 = exp(r3202055);
        double r3202057 = r3202053 + r3202056;
        double r3202058 = r3202053 / r3202057;
        double r3202059 = c_p;
        double r3202060 = pow(r3202058, r3202059);
        double r3202061 = r3202053 - r3202058;
        double r3202062 = c_n;
        double r3202063 = pow(r3202061, r3202062);
        double r3202064 = r3202060 * r3202063;
        double r3202065 = t;
        double r3202066 = -r3202065;
        double r3202067 = exp(r3202066);
        double r3202068 = r3202053 + r3202067;
        double r3202069 = r3202053 / r3202068;
        double r3202070 = pow(r3202069, r3202059);
        double r3202071 = r3202053 - r3202069;
        double r3202072 = pow(r3202071, r3202062);
        double r3202073 = r3202070 * r3202072;
        double r3202074 = r3202064 / r3202073;
        return r3202074;
}