double f(double c_p, double c_n, double t, double s) {
double r2257827 = 1.0;
double r2257828 = s;
double r2257829 = -r2257828;
double r2257830 = exp(r2257829);
double r2257831 = r2257827 + r2257830;
double r2257832 = r2257827 / r2257831;
double r2257833 = c_p;
double r2257834 = pow(r2257832, r2257833);
double r2257835 = r2257827 - r2257832;
double r2257836 = c_n;
double r2257837 = pow(r2257835, r2257836);
double r2257838 = r2257834 * r2257837;
double r2257839 = t;
double r2257840 = -r2257839;
double r2257841 = exp(r2257840);
double r2257842 = r2257827 + r2257841;
double r2257843 = r2257827 / r2257842;
double r2257844 = pow(r2257843, r2257833);
double r2257845 = r2257827 - r2257843;
double r2257846 = pow(r2257845, r2257836);
double r2257847 = r2257844 * r2257846;
double r2257848 = r2257838 / r2257847;
return r2257848;
}