double f(double c_p, double c_n, double t, double s) {
double r84803 = 1.0;
double r84804 = s;
double r84805 = -r84804;
double r84806 = exp(r84805);
double r84807 = r84803 + r84806;
double r84808 = r84803 / r84807;
double r84809 = c_p;
double r84810 = pow(r84808, r84809);
double r84811 = r84803 - r84808;
double r84812 = c_n;
double r84813 = pow(r84811, r84812);
double r84814 = r84810 * r84813;
double r84815 = t;
double r84816 = -r84815;
double r84817 = exp(r84816);
double r84818 = r84803 + r84817;
double r84819 = r84803 / r84818;
double r84820 = pow(r84819, r84809);
double r84821 = r84803 - r84819;
double r84822 = pow(r84821, r84812);
double r84823 = r84820 * r84822;
double r84824 = r84814 / r84823;
return r84824;
}