double f(double x, double y, double z, double t) {
double r11227328 = x;
double r11227329 = 1.0;
double r11227330 = r11227328 + r11227329;
double r11227331 = sqrt(r11227330);
double r11227332 = sqrt(r11227328);
double r11227333 = r11227331 - r11227332;
double r11227334 = y;
double r11227335 = r11227334 + r11227329;
double r11227336 = sqrt(r11227335);
double r11227337 = sqrt(r11227334);
double r11227338 = r11227336 - r11227337;
double r11227339 = r11227333 + r11227338;
double r11227340 = z;
double r11227341 = r11227340 + r11227329;
double r11227342 = sqrt(r11227341);
double r11227343 = sqrt(r11227340);
double r11227344 = r11227342 - r11227343;
double r11227345 = r11227339 + r11227344;
double r11227346 = t;
double r11227347 = r11227346 + r11227329;
double r11227348 = sqrt(r11227347);
double r11227349 = sqrt(r11227346);
double r11227350 = r11227348 - r11227349;
double r11227351 = r11227345 + r11227350;
return r11227351;
}