double f(double x, double y, double z, double t) {
double r525247 = x;
double r525248 = 1.0;
double r525249 = r525247 + r525248;
double r525250 = sqrt(r525249);
double r525251 = sqrt(r525247);
double r525252 = r525250 - r525251;
double r525253 = y;
double r525254 = r525253 + r525248;
double r525255 = sqrt(r525254);
double r525256 = sqrt(r525253);
double r525257 = r525255 - r525256;
double r525258 = r525252 + r525257;
double r525259 = z;
double r525260 = r525259 + r525248;
double r525261 = sqrt(r525260);
double r525262 = sqrt(r525259);
double r525263 = r525261 - r525262;
double r525264 = r525258 + r525263;
double r525265 = t;
double r525266 = r525265 + r525248;
double r525267 = sqrt(r525266);
double r525268 = sqrt(r525265);
double r525269 = r525267 - r525268;
double r525270 = r525264 + r525269;
return r525270;
}