double f(double x, double y, double z, double t) {
double r1017473 = x;
double r1017474 = 1.0;
double r1017475 = r1017473 + r1017474;
double r1017476 = sqrt(r1017475);
double r1017477 = sqrt(r1017473);
double r1017478 = r1017476 - r1017477;
double r1017479 = y;
double r1017480 = r1017479 + r1017474;
double r1017481 = sqrt(r1017480);
double r1017482 = sqrt(r1017479);
double r1017483 = r1017481 - r1017482;
double r1017484 = r1017478 + r1017483;
double r1017485 = z;
double r1017486 = r1017485 + r1017474;
double r1017487 = sqrt(r1017486);
double r1017488 = sqrt(r1017485);
double r1017489 = r1017487 - r1017488;
double r1017490 = r1017484 + r1017489;
double r1017491 = t;
double r1017492 = r1017491 + r1017474;
double r1017493 = sqrt(r1017492);
double r1017494 = sqrt(r1017491);
double r1017495 = r1017493 - r1017494;
double r1017496 = r1017490 + r1017495;
return r1017496;
}