double f(double x, double y, double z, double t) {
double r500044 = x;
double r500045 = 1.0;
double r500046 = r500044 + r500045;
double r500047 = sqrt(r500046);
double r500048 = sqrt(r500044);
double r500049 = r500047 - r500048;
double r500050 = y;
double r500051 = r500050 + r500045;
double r500052 = sqrt(r500051);
double r500053 = sqrt(r500050);
double r500054 = r500052 - r500053;
double r500055 = r500049 + r500054;
double r500056 = z;
double r500057 = r500056 + r500045;
double r500058 = sqrt(r500057);
double r500059 = sqrt(r500056);
double r500060 = r500058 - r500059;
double r500061 = r500055 + r500060;
double r500062 = t;
double r500063 = r500062 + r500045;
double r500064 = sqrt(r500063);
double r500065 = sqrt(r500062);
double r500066 = r500064 - r500065;
double r500067 = r500061 + r500066;
return r500067;
}