double f(double c_p, double c_n, double t, double s) {
double r92107 = 1.0;
double r92108 = s;
double r92109 = -r92108;
double r92110 = exp(r92109);
double r92111 = r92107 + r92110;
double r92112 = r92107 / r92111;
double r92113 = c_p;
double r92114 = pow(r92112, r92113);
double r92115 = r92107 - r92112;
double r92116 = c_n;
double r92117 = pow(r92115, r92116);
double r92118 = r92114 * r92117;
double r92119 = t;
double r92120 = -r92119;
double r92121 = exp(r92120);
double r92122 = r92107 + r92121;
double r92123 = r92107 / r92122;
double r92124 = pow(r92123, r92113);
double r92125 = r92107 - r92123;
double r92126 = pow(r92125, r92116);
double r92127 = r92124 * r92126;
double r92128 = r92118 / r92127;
return r92128;
}