double f(double x, double y, double z, double t) {
double r511323 = x;
double r511324 = 1.0;
double r511325 = r511323 + r511324;
double r511326 = sqrt(r511325);
double r511327 = sqrt(r511323);
double r511328 = r511326 - r511327;
double r511329 = y;
double r511330 = r511329 + r511324;
double r511331 = sqrt(r511330);
double r511332 = sqrt(r511329);
double r511333 = r511331 - r511332;
double r511334 = r511328 + r511333;
double r511335 = z;
double r511336 = r511335 + r511324;
double r511337 = sqrt(r511336);
double r511338 = sqrt(r511335);
double r511339 = r511337 - r511338;
double r511340 = r511334 + r511339;
double r511341 = t;
double r511342 = r511341 + r511324;
double r511343 = sqrt(r511342);
double r511344 = sqrt(r511341);
double r511345 = r511343 - r511344;
double r511346 = r511340 + r511345;
return r511346;
}