double f(double c_p, double c_n, double t, double s) {
double r252535 = 1.0;
double r252536 = s;
double r252537 = -r252536;
double r252538 = exp(r252537);
double r252539 = r252535 + r252538;
double r252540 = r252535 / r252539;
double r252541 = c_p;
double r252542 = pow(r252540, r252541);
double r252543 = r252535 - r252540;
double r252544 = c_n;
double r252545 = pow(r252543, r252544);
double r252546 = r252542 * r252545;
double r252547 = t;
double r252548 = -r252547;
double r252549 = exp(r252548);
double r252550 = r252535 + r252549;
double r252551 = r252535 / r252550;
double r252552 = pow(r252551, r252541);
double r252553 = r252535 - r252551;
double r252554 = pow(r252553, r252544);
double r252555 = r252552 * r252554;
double r252556 = r252546 / r252555;
return r252556;
}