double f(double x, double y, double z, double t) {
double r431642 = x;
double r431643 = 1.0;
double r431644 = r431642 + r431643;
double r431645 = sqrt(r431644);
double r431646 = sqrt(r431642);
double r431647 = r431645 - r431646;
double r431648 = y;
double r431649 = r431648 + r431643;
double r431650 = sqrt(r431649);
double r431651 = sqrt(r431648);
double r431652 = r431650 - r431651;
double r431653 = r431647 + r431652;
double r431654 = z;
double r431655 = r431654 + r431643;
double r431656 = sqrt(r431655);
double r431657 = sqrt(r431654);
double r431658 = r431656 - r431657;
double r431659 = r431653 + r431658;
double r431660 = t;
double r431661 = r431660 + r431643;
double r431662 = sqrt(r431661);
double r431663 = sqrt(r431660);
double r431664 = r431662 - r431663;
double r431665 = r431659 + r431664;
return r431665;
}