double f(double c_p, double c_n, double t, double s) {
double r218859 = 1.0;
double r218860 = s;
double r218861 = -r218860;
double r218862 = exp(r218861);
double r218863 = r218859 + r218862;
double r218864 = r218859 / r218863;
double r218865 = c_p;
double r218866 = pow(r218864, r218865);
double r218867 = r218859 - r218864;
double r218868 = c_n;
double r218869 = pow(r218867, r218868);
double r218870 = r218866 * r218869;
double r218871 = t;
double r218872 = -r218871;
double r218873 = exp(r218872);
double r218874 = r218859 + r218873;
double r218875 = r218859 / r218874;
double r218876 = pow(r218875, r218865);
double r218877 = r218859 - r218875;
double r218878 = pow(r218877, r218868);
double r218879 = r218876 * r218878;
double r218880 = r218870 / r218879;
return r218880;
}