double f(double c_p, double c_n, double t, double s) {
double r3576354 = 1.0;
double r3576355 = s;
double r3576356 = -r3576355;
double r3576357 = exp(r3576356);
double r3576358 = r3576354 + r3576357;
double r3576359 = r3576354 / r3576358;
double r3576360 = c_p;
double r3576361 = pow(r3576359, r3576360);
double r3576362 = r3576354 - r3576359;
double r3576363 = c_n;
double r3576364 = pow(r3576362, r3576363);
double r3576365 = r3576361 * r3576364;
double r3576366 = t;
double r3576367 = -r3576366;
double r3576368 = exp(r3576367);
double r3576369 = r3576354 + r3576368;
double r3576370 = r3576354 / r3576369;
double r3576371 = pow(r3576370, r3576360);
double r3576372 = r3576354 - r3576370;
double r3576373 = pow(r3576372, r3576363);
double r3576374 = r3576371 * r3576373;
double r3576375 = r3576365 / r3576374;
return r3576375;
}