double f(double x, double y, double z, double t) {
double r23095158 = x;
double r23095159 = 1.0;
double r23095160 = r23095158 + r23095159;
double r23095161 = sqrt(r23095160);
double r23095162 = sqrt(r23095158);
double r23095163 = r23095161 - r23095162;
double r23095164 = y;
double r23095165 = r23095164 + r23095159;
double r23095166 = sqrt(r23095165);
double r23095167 = sqrt(r23095164);
double r23095168 = r23095166 - r23095167;
double r23095169 = r23095163 + r23095168;
double r23095170 = z;
double r23095171 = r23095170 + r23095159;
double r23095172 = sqrt(r23095171);
double r23095173 = sqrt(r23095170);
double r23095174 = r23095172 - r23095173;
double r23095175 = r23095169 + r23095174;
double r23095176 = t;
double r23095177 = r23095176 + r23095159;
double r23095178 = sqrt(r23095177);
double r23095179 = sqrt(r23095176);
double r23095180 = r23095178 - r23095179;
double r23095181 = r23095175 + r23095180;
return r23095181;
}