double f(double c_p, double c_n, double t, double s) {
double r7342720 = 1.0;
double r7342721 = s;
double r7342722 = -r7342721;
double r7342723 = exp(r7342722);
double r7342724 = r7342720 + r7342723;
double r7342725 = r7342720 / r7342724;
double r7342726 = c_p;
double r7342727 = pow(r7342725, r7342726);
double r7342728 = r7342720 - r7342725;
double r7342729 = c_n;
double r7342730 = pow(r7342728, r7342729);
double r7342731 = r7342727 * r7342730;
double r7342732 = t;
double r7342733 = -r7342732;
double r7342734 = exp(r7342733);
double r7342735 = r7342720 + r7342734;
double r7342736 = r7342720 / r7342735;
double r7342737 = pow(r7342736, r7342726);
double r7342738 = r7342720 - r7342736;
double r7342739 = pow(r7342738, r7342729);
double r7342740 = r7342737 * r7342739;
double r7342741 = r7342731 / r7342740;
return r7342741;
}