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 r150162 = 1.0;
        double r150163 = s;
        double r150164 = -r150163;
        double r150165 = exp(r150164);
        double r150166 = r150162 + r150165;
        double r150167 = r150162 / r150166;
        double r150168 = c_p;
        double r150169 = pow(r150167, r150168);
        double r150170 = r150162 - r150167;
        double r150171 = c_n;
        double r150172 = pow(r150170, r150171);
        double r150173 = r150169 * r150172;
        double r150174 = t;
        double r150175 = -r150174;
        double r150176 = exp(r150175);
        double r150177 = r150162 + r150176;
        double r150178 = r150162 / r150177;
        double r150179 = pow(r150178, r150168);
        double r150180 = r150162 - r150178;
        double r150181 = pow(r150180, r150171);
        double r150182 = r150179 * r150181;
        double r150183 = r150173 / r150182;
        return r150183;
}