double f(double x, double y, double z, double t) {
double r574057 = x;
double r574058 = 1.0;
double r574059 = r574057 + r574058;
double r574060 = sqrt(r574059);
double r574061 = sqrt(r574057);
double r574062 = r574060 - r574061;
double r574063 = y;
double r574064 = r574063 + r574058;
double r574065 = sqrt(r574064);
double r574066 = sqrt(r574063);
double r574067 = r574065 - r574066;
double r574068 = r574062 + r574067;
double r574069 = z;
double r574070 = r574069 + r574058;
double r574071 = sqrt(r574070);
double r574072 = sqrt(r574069);
double r574073 = r574071 - r574072;
double r574074 = r574068 + r574073;
double r574075 = t;
double r574076 = r574075 + r574058;
double r574077 = sqrt(r574076);
double r574078 = sqrt(r574075);
double r574079 = r574077 - r574078;
double r574080 = r574074 + r574079;
return r574080;
}