double f(double c_p, double c_n, double t, double s) {
double r187973 = 1.0;
double r187974 = s;
double r187975 = -r187974;
double r187976 = exp(r187975);
double r187977 = r187973 + r187976;
double r187978 = r187973 / r187977;
double r187979 = c_p;
double r187980 = pow(r187978, r187979);
double r187981 = r187973 - r187978;
double r187982 = c_n;
double r187983 = pow(r187981, r187982);
double r187984 = r187980 * r187983;
double r187985 = t;
double r187986 = -r187985;
double r187987 = exp(r187986);
double r187988 = r187973 + r187987;
double r187989 = r187973 / r187988;
double r187990 = pow(r187989, r187979);
double r187991 = r187973 - r187989;
double r187992 = pow(r187991, r187982);
double r187993 = r187990 * r187992;
double r187994 = r187984 / r187993;
return r187994;
}