double f(double x, double y, double z, double t) {
double r341224 = x;
double r341225 = 1.0;
double r341226 = r341224 + r341225;
double r341227 = sqrt(r341226);
double r341228 = sqrt(r341224);
double r341229 = r341227 - r341228;
double r341230 = y;
double r341231 = r341230 + r341225;
double r341232 = sqrt(r341231);
double r341233 = sqrt(r341230);
double r341234 = r341232 - r341233;
double r341235 = r341229 + r341234;
double r341236 = z;
double r341237 = r341236 + r341225;
double r341238 = sqrt(r341237);
double r341239 = sqrt(r341236);
double r341240 = r341238 - r341239;
double r341241 = r341235 + r341240;
double r341242 = t;
double r341243 = r341242 + r341225;
double r341244 = sqrt(r341243);
double r341245 = sqrt(r341242);
double r341246 = r341244 - r341245;
double r341247 = r341241 + r341246;
return r341247;
}