double f(double c_p, double c_n, double t, double s) {
double r7836846 = 1.0;
double r7836847 = s;
double r7836848 = -r7836847;
double r7836849 = exp(r7836848);
double r7836850 = r7836846 + r7836849;
double r7836851 = r7836846 / r7836850;
double r7836852 = c_p;
double r7836853 = pow(r7836851, r7836852);
double r7836854 = r7836846 - r7836851;
double r7836855 = c_n;
double r7836856 = pow(r7836854, r7836855);
double r7836857 = r7836853 * r7836856;
double r7836858 = t;
double r7836859 = -r7836858;
double r7836860 = exp(r7836859);
double r7836861 = r7836846 + r7836860;
double r7836862 = r7836846 / r7836861;
double r7836863 = pow(r7836862, r7836852);
double r7836864 = r7836846 - r7836862;
double r7836865 = pow(r7836864, r7836855);
double r7836866 = r7836863 * r7836865;
double r7836867 = r7836857 / r7836866;
return r7836867;
}