double f(double x, double y, double z, double t) {
double r524786 = x;
double r524787 = 1.0;
double r524788 = r524786 + r524787;
double r524789 = sqrt(r524788);
double r524790 = sqrt(r524786);
double r524791 = r524789 - r524790;
double r524792 = y;
double r524793 = r524792 + r524787;
double r524794 = sqrt(r524793);
double r524795 = sqrt(r524792);
double r524796 = r524794 - r524795;
double r524797 = r524791 + r524796;
double r524798 = z;
double r524799 = r524798 + r524787;
double r524800 = sqrt(r524799);
double r524801 = sqrt(r524798);
double r524802 = r524800 - r524801;
double r524803 = r524797 + r524802;
double r524804 = t;
double r524805 = r524804 + r524787;
double r524806 = sqrt(r524805);
double r524807 = sqrt(r524804);
double r524808 = r524806 - r524807;
double r524809 = r524803 + r524808;
return r524809;
}