double f(double c_p, double c_n, double t, double s) {
double r178532 = 1.0;
double r178533 = s;
double r178534 = -r178533;
double r178535 = exp(r178534);
double r178536 = r178532 + r178535;
double r178537 = r178532 / r178536;
double r178538 = c_p;
double r178539 = pow(r178537, r178538);
double r178540 = r178532 - r178537;
double r178541 = c_n;
double r178542 = pow(r178540, r178541);
double r178543 = r178539 * r178542;
double r178544 = t;
double r178545 = -r178544;
double r178546 = exp(r178545);
double r178547 = r178532 + r178546;
double r178548 = r178532 / r178547;
double r178549 = pow(r178548, r178538);
double r178550 = r178532 - r178548;
double r178551 = pow(r178550, r178541);
double r178552 = r178549 * r178551;
double r178553 = r178543 / r178552;
return r178553;
}