double f(double c_p, double c_n, double t, double s) {
double r178429 = 1.0;
double r178430 = s;
double r178431 = -r178430;
double r178432 = exp(r178431);
double r178433 = r178429 + r178432;
double r178434 = r178429 / r178433;
double r178435 = c_p;
double r178436 = pow(r178434, r178435);
double r178437 = r178429 - r178434;
double r178438 = c_n;
double r178439 = pow(r178437, r178438);
double r178440 = r178436 * r178439;
double r178441 = t;
double r178442 = -r178441;
double r178443 = exp(r178442);
double r178444 = r178429 + r178443;
double r178445 = r178429 / r178444;
double r178446 = pow(r178445, r178435);
double r178447 = r178429 - r178445;
double r178448 = pow(r178447, r178438);
double r178449 = r178446 * r178448;
double r178450 = r178440 / r178449;
return r178450;
}