double f(double c_p, double c_n, double t, double s) {
double r105275 = 1.0;
double r105276 = s;
double r105277 = -r105276;
double r105278 = exp(r105277);
double r105279 = r105275 + r105278;
double r105280 = r105275 / r105279;
double r105281 = c_p;
double r105282 = pow(r105280, r105281);
double r105283 = r105275 - r105280;
double r105284 = c_n;
double r105285 = pow(r105283, r105284);
double r105286 = r105282 * r105285;
double r105287 = t;
double r105288 = -r105287;
double r105289 = exp(r105288);
double r105290 = r105275 + r105289;
double r105291 = r105275 / r105290;
double r105292 = pow(r105291, r105281);
double r105293 = r105275 - r105291;
double r105294 = pow(r105293, r105284);
double r105295 = r105292 * r105294;
double r105296 = r105286 / r105295;
return r105296;
}