double f(double c_p, double c_n, double t, double s) {
double r6952338 = 1.0;
double r6952339 = s;
double r6952340 = -r6952339;
double r6952341 = exp(r6952340);
double r6952342 = r6952338 + r6952341;
double r6952343 = r6952338 / r6952342;
double r6952344 = c_p;
double r6952345 = pow(r6952343, r6952344);
double r6952346 = r6952338 - r6952343;
double r6952347 = c_n;
double r6952348 = pow(r6952346, r6952347);
double r6952349 = r6952345 * r6952348;
double r6952350 = t;
double r6952351 = -r6952350;
double r6952352 = exp(r6952351);
double r6952353 = r6952338 + r6952352;
double r6952354 = r6952338 / r6952353;
double r6952355 = pow(r6952354, r6952344);
double r6952356 = r6952338 - r6952354;
double r6952357 = pow(r6952356, r6952347);
double r6952358 = r6952355 * r6952357;
double r6952359 = r6952349 / r6952358;
return r6952359;
}