double f(double c_p, double c_n, double t, double s) {
double r180283 = 1.0;
double r180284 = s;
double r180285 = -r180284;
double r180286 = exp(r180285);
double r180287 = r180283 + r180286;
double r180288 = r180283 / r180287;
double r180289 = c_p;
double r180290 = pow(r180288, r180289);
double r180291 = r180283 - r180288;
double r180292 = c_n;
double r180293 = pow(r180291, r180292);
double r180294 = r180290 * r180293;
double r180295 = t;
double r180296 = -r180295;
double r180297 = exp(r180296);
double r180298 = r180283 + r180297;
double r180299 = r180283 / r180298;
double r180300 = pow(r180299, r180289);
double r180301 = r180283 - r180299;
double r180302 = pow(r180301, r180292);
double r180303 = r180300 * r180302;
double r180304 = r180294 / r180303;
return r180304;
}