double f(double c_p, double c_n, double t, double s) {
double r6552256 = 1.0;
double r6552257 = s;
double r6552258 = -r6552257;
double r6552259 = exp(r6552258);
double r6552260 = r6552256 + r6552259;
double r6552261 = r6552256 / r6552260;
double r6552262 = c_p;
double r6552263 = pow(r6552261, r6552262);
double r6552264 = r6552256 - r6552261;
double r6552265 = c_n;
double r6552266 = pow(r6552264, r6552265);
double r6552267 = r6552263 * r6552266;
double r6552268 = t;
double r6552269 = -r6552268;
double r6552270 = exp(r6552269);
double r6552271 = r6552256 + r6552270;
double r6552272 = r6552256 / r6552271;
double r6552273 = pow(r6552272, r6552262);
double r6552274 = r6552256 - r6552272;
double r6552275 = pow(r6552274, r6552265);
double r6552276 = r6552273 * r6552275;
double r6552277 = r6552267 / r6552276;
return r6552277;
}