double f(double x, double y, double z, double t) {
double r589995 = x;
double r589996 = 1.0;
double r589997 = r589995 + r589996;
double r589998 = sqrt(r589997);
double r589999 = sqrt(r589995);
double r590000 = r589998 - r589999;
double r590001 = y;
double r590002 = r590001 + r589996;
double r590003 = sqrt(r590002);
double r590004 = sqrt(r590001);
double r590005 = r590003 - r590004;
double r590006 = r590000 + r590005;
double r590007 = z;
double r590008 = r590007 + r589996;
double r590009 = sqrt(r590008);
double r590010 = sqrt(r590007);
double r590011 = r590009 - r590010;
double r590012 = r590006 + r590011;
double r590013 = t;
double r590014 = r590013 + r589996;
double r590015 = sqrt(r590014);
double r590016 = sqrt(r590013);
double r590017 = r590015 - r590016;
double r590018 = r590012 + r590017;
return r590018;
}