double f(double c_p, double c_n, double t, double s) {
double r6262961 = 1.0;
double r6262962 = s;
double r6262963 = -r6262962;
double r6262964 = exp(r6262963);
double r6262965 = r6262961 + r6262964;
double r6262966 = r6262961 / r6262965;
double r6262967 = c_p;
double r6262968 = pow(r6262966, r6262967);
double r6262969 = r6262961 - r6262966;
double r6262970 = c_n;
double r6262971 = pow(r6262969, r6262970);
double r6262972 = r6262968 * r6262971;
double r6262973 = t;
double r6262974 = -r6262973;
double r6262975 = exp(r6262974);
double r6262976 = r6262961 + r6262975;
double r6262977 = r6262961 / r6262976;
double r6262978 = pow(r6262977, r6262967);
double r6262979 = r6262961 - r6262977;
double r6262980 = pow(r6262979, r6262970);
double r6262981 = r6262978 * r6262980;
double r6262982 = r6262972 / r6262981;
return r6262982;
}