double f(double c_p, double c_n, double t, double s) {
double r216346 = 1.0;
double r216347 = s;
double r216348 = -r216347;
double r216349 = exp(r216348);
double r216350 = r216346 + r216349;
double r216351 = r216346 / r216350;
double r216352 = c_p;
double r216353 = pow(r216351, r216352);
double r216354 = r216346 - r216351;
double r216355 = c_n;
double r216356 = pow(r216354, r216355);
double r216357 = r216353 * r216356;
double r216358 = t;
double r216359 = -r216358;
double r216360 = exp(r216359);
double r216361 = r216346 + r216360;
double r216362 = r216346 / r216361;
double r216363 = pow(r216362, r216352);
double r216364 = r216346 - r216362;
double r216365 = pow(r216364, r216355);
double r216366 = r216363 * r216365;
double r216367 = r216357 / r216366;
return r216367;
}