double f(double x, double y, double z, double t) {
double r591525 = x;
double r591526 = 1.0;
double r591527 = r591525 + r591526;
double r591528 = sqrt(r591527);
double r591529 = sqrt(r591525);
double r591530 = r591528 - r591529;
double r591531 = y;
double r591532 = r591531 + r591526;
double r591533 = sqrt(r591532);
double r591534 = sqrt(r591531);
double r591535 = r591533 - r591534;
double r591536 = r591530 + r591535;
double r591537 = z;
double r591538 = r591537 + r591526;
double r591539 = sqrt(r591538);
double r591540 = sqrt(r591537);
double r591541 = r591539 - r591540;
double r591542 = r591536 + r591541;
double r591543 = t;
double r591544 = r591543 + r591526;
double r591545 = sqrt(r591544);
double r591546 = sqrt(r591543);
double r591547 = r591545 - r591546;
double r591548 = r591542 + r591547;
return r591548;
}