double f(double x, double y, double z, double t) {
double r522736 = x;
double r522737 = 1.0;
double r522738 = r522736 + r522737;
double r522739 = sqrt(r522738);
double r522740 = sqrt(r522736);
double r522741 = r522739 - r522740;
double r522742 = y;
double r522743 = r522742 + r522737;
double r522744 = sqrt(r522743);
double r522745 = sqrt(r522742);
double r522746 = r522744 - r522745;
double r522747 = r522741 + r522746;
double r522748 = z;
double r522749 = r522748 + r522737;
double r522750 = sqrt(r522749);
double r522751 = sqrt(r522748);
double r522752 = r522750 - r522751;
double r522753 = r522747 + r522752;
double r522754 = t;
double r522755 = r522754 + r522737;
double r522756 = sqrt(r522755);
double r522757 = sqrt(r522754);
double r522758 = r522756 - r522757;
double r522759 = r522753 + r522758;
return r522759;
}