double f(double c_p, double c_n, double t, double s) {
double r7354086 = 1.0;
double r7354087 = s;
double r7354088 = -r7354087;
double r7354089 = exp(r7354088);
double r7354090 = r7354086 + r7354089;
double r7354091 = r7354086 / r7354090;
double r7354092 = c_p;
double r7354093 = pow(r7354091, r7354092);
double r7354094 = r7354086 - r7354091;
double r7354095 = c_n;
double r7354096 = pow(r7354094, r7354095);
double r7354097 = r7354093 * r7354096;
double r7354098 = t;
double r7354099 = -r7354098;
double r7354100 = exp(r7354099);
double r7354101 = r7354086 + r7354100;
double r7354102 = r7354086 / r7354101;
double r7354103 = pow(r7354102, r7354092);
double r7354104 = r7354086 - r7354102;
double r7354105 = pow(r7354104, r7354095);
double r7354106 = r7354103 * r7354105;
double r7354107 = r7354097 / r7354106;
return r7354107;
}