double f(double c_p, double c_n, double t, double s) {
double r7112258 = 1.0;
double r7112259 = s;
double r7112260 = -r7112259;
double r7112261 = exp(r7112260);
double r7112262 = r7112258 + r7112261;
double r7112263 = r7112258 / r7112262;
double r7112264 = c_p;
double r7112265 = pow(r7112263, r7112264);
double r7112266 = r7112258 - r7112263;
double r7112267 = c_n;
double r7112268 = pow(r7112266, r7112267);
double r7112269 = r7112265 * r7112268;
double r7112270 = t;
double r7112271 = -r7112270;
double r7112272 = exp(r7112271);
double r7112273 = r7112258 + r7112272;
double r7112274 = r7112258 / r7112273;
double r7112275 = pow(r7112274, r7112264);
double r7112276 = r7112258 - r7112274;
double r7112277 = pow(r7112276, r7112267);
double r7112278 = r7112275 * r7112277;
double r7112279 = r7112269 / r7112278;
return r7112279;
}