double f(double c_p, double c_n, double t, double s) {
double r217259 = 1.0;
double r217260 = s;
double r217261 = -r217260;
double r217262 = exp(r217261);
double r217263 = r217259 + r217262;
double r217264 = r217259 / r217263;
double r217265 = c_p;
double r217266 = pow(r217264, r217265);
double r217267 = r217259 - r217264;
double r217268 = c_n;
double r217269 = pow(r217267, r217268);
double r217270 = r217266 * r217269;
double r217271 = t;
double r217272 = -r217271;
double r217273 = exp(r217272);
double r217274 = r217259 + r217273;
double r217275 = r217259 / r217274;
double r217276 = pow(r217275, r217265);
double r217277 = r217259 - r217275;
double r217278 = pow(r217277, r217268);
double r217279 = r217276 * r217278;
double r217280 = r217270 / r217279;
return r217280;
}