double f(double c_p, double c_n, double t, double s) {
double r280409 = 1.0;
double r280410 = s;
double r280411 = -r280410;
double r280412 = exp(r280411);
double r280413 = r280409 + r280412;
double r280414 = r280409 / r280413;
double r280415 = c_p;
double r280416 = pow(r280414, r280415);
double r280417 = r280409 - r280414;
double r280418 = c_n;
double r280419 = pow(r280417, r280418);
double r280420 = r280416 * r280419;
double r280421 = t;
double r280422 = -r280421;
double r280423 = exp(r280422);
double r280424 = r280409 + r280423;
double r280425 = r280409 / r280424;
double r280426 = pow(r280425, r280415);
double r280427 = r280409 - r280425;
double r280428 = pow(r280427, r280418);
double r280429 = r280426 * r280428;
double r280430 = r280420 / r280429;
return r280430;
}