double f(double c_p, double c_n, double t, double s) {
double r9630464 = 1.0;
double r9630465 = s;
double r9630466 = -r9630465;
double r9630467 = exp(r9630466);
double r9630468 = r9630464 + r9630467;
double r9630469 = r9630464 / r9630468;
double r9630470 = c_p;
double r9630471 = pow(r9630469, r9630470);
double r9630472 = r9630464 - r9630469;
double r9630473 = c_n;
double r9630474 = pow(r9630472, r9630473);
double r9630475 = r9630471 * r9630474;
double r9630476 = t;
double r9630477 = -r9630476;
double r9630478 = exp(r9630477);
double r9630479 = r9630464 + r9630478;
double r9630480 = r9630464 / r9630479;
double r9630481 = pow(r9630480, r9630470);
double r9630482 = r9630464 - r9630480;
double r9630483 = pow(r9630482, r9630473);
double r9630484 = r9630481 * r9630483;
double r9630485 = r9630475 / r9630484;
return r9630485;
}