double f(double c_p, double c_n, double t, double s) {
double r7602253 = 1.0;
double r7602254 = s;
double r7602255 = -r7602254;
double r7602256 = exp(r7602255);
double r7602257 = r7602253 + r7602256;
double r7602258 = r7602253 / r7602257;
double r7602259 = c_p;
double r7602260 = pow(r7602258, r7602259);
double r7602261 = r7602253 - r7602258;
double r7602262 = c_n;
double r7602263 = pow(r7602261, r7602262);
double r7602264 = r7602260 * r7602263;
double r7602265 = t;
double r7602266 = -r7602265;
double r7602267 = exp(r7602266);
double r7602268 = r7602253 + r7602267;
double r7602269 = r7602253 / r7602268;
double r7602270 = pow(r7602269, r7602259);
double r7602271 = r7602253 - r7602269;
double r7602272 = pow(r7602271, r7602262);
double r7602273 = r7602270 * r7602272;
double r7602274 = r7602264 / r7602273;
return r7602274;
}