double f(double x, double y, double z, double t) {
double r521179 = x;
double r521180 = 1.0;
double r521181 = r521179 + r521180;
double r521182 = sqrt(r521181);
double r521183 = sqrt(r521179);
double r521184 = r521182 - r521183;
double r521185 = y;
double r521186 = r521185 + r521180;
double r521187 = sqrt(r521186);
double r521188 = sqrt(r521185);
double r521189 = r521187 - r521188;
double r521190 = r521184 + r521189;
double r521191 = z;
double r521192 = r521191 + r521180;
double r521193 = sqrt(r521192);
double r521194 = sqrt(r521191);
double r521195 = r521193 - r521194;
double r521196 = r521190 + r521195;
double r521197 = t;
double r521198 = r521197 + r521180;
double r521199 = sqrt(r521198);
double r521200 = sqrt(r521197);
double r521201 = r521199 - r521200;
double r521202 = r521196 + r521201;
return r521202;
}