double f(double c_p, double c_n, double t, double s) {
double r192585 = 1.0;
double r192586 = s;
double r192587 = -r192586;
double r192588 = exp(r192587);
double r192589 = r192585 + r192588;
double r192590 = r192585 / r192589;
double r192591 = c_p;
double r192592 = pow(r192590, r192591);
double r192593 = r192585 - r192590;
double r192594 = c_n;
double r192595 = pow(r192593, r192594);
double r192596 = r192592 * r192595;
double r192597 = t;
double r192598 = -r192597;
double r192599 = exp(r192598);
double r192600 = r192585 + r192599;
double r192601 = r192585 / r192600;
double r192602 = pow(r192601, r192591);
double r192603 = r192585 - r192601;
double r192604 = pow(r192603, r192594);
double r192605 = r192602 * r192604;
double r192606 = r192596 / r192605;
return r192606;
}