double f(double x, double y, double z, double t) {
double r587330 = x;
double r587331 = 1.0;
double r587332 = r587330 + r587331;
double r587333 = sqrt(r587332);
double r587334 = sqrt(r587330);
double r587335 = r587333 - r587334;
double r587336 = y;
double r587337 = r587336 + r587331;
double r587338 = sqrt(r587337);
double r587339 = sqrt(r587336);
double r587340 = r587338 - r587339;
double r587341 = r587335 + r587340;
double r587342 = z;
double r587343 = r587342 + r587331;
double r587344 = sqrt(r587343);
double r587345 = sqrt(r587342);
double r587346 = r587344 - r587345;
double r587347 = r587341 + r587346;
double r587348 = t;
double r587349 = r587348 + r587331;
double r587350 = sqrt(r587349);
double r587351 = sqrt(r587348);
double r587352 = r587350 - r587351;
double r587353 = r587347 + r587352;
return r587353;
}