double f(double c_p, double c_n, double t, double s) {
double r182740 = 1.0;
double r182741 = s;
double r182742 = -r182741;
double r182743 = exp(r182742);
double r182744 = r182740 + r182743;
double r182745 = r182740 / r182744;
double r182746 = c_p;
double r182747 = pow(r182745, r182746);
double r182748 = r182740 - r182745;
double r182749 = c_n;
double r182750 = pow(r182748, r182749);
double r182751 = r182747 * r182750;
double r182752 = t;
double r182753 = -r182752;
double r182754 = exp(r182753);
double r182755 = r182740 + r182754;
double r182756 = r182740 / r182755;
double r182757 = pow(r182756, r182746);
double r182758 = r182740 - r182756;
double r182759 = pow(r182758, r182749);
double r182760 = r182757 * r182759;
double r182761 = r182751 / r182760;
return r182761;
}