double f(double c_p, double c_n, double t, double s) {
double r206517 = 1.0;
double r206518 = s;
double r206519 = -r206518;
double r206520 = exp(r206519);
double r206521 = r206517 + r206520;
double r206522 = r206517 / r206521;
double r206523 = c_p;
double r206524 = pow(r206522, r206523);
double r206525 = r206517 - r206522;
double r206526 = c_n;
double r206527 = pow(r206525, r206526);
double r206528 = r206524 * r206527;
double r206529 = t;
double r206530 = -r206529;
double r206531 = exp(r206530);
double r206532 = r206517 + r206531;
double r206533 = r206517 / r206532;
double r206534 = pow(r206533, r206523);
double r206535 = r206517 - r206533;
double r206536 = pow(r206535, r206526);
double r206537 = r206534 * r206536;
double r206538 = r206528 / r206537;
return r206538;
}