double f(double x, double y, double z, double t) {
double r570913 = x;
double r570914 = 1.0;
double r570915 = r570913 + r570914;
double r570916 = sqrt(r570915);
double r570917 = sqrt(r570913);
double r570918 = r570916 - r570917;
double r570919 = y;
double r570920 = r570919 + r570914;
double r570921 = sqrt(r570920);
double r570922 = sqrt(r570919);
double r570923 = r570921 - r570922;
double r570924 = r570918 + r570923;
double r570925 = z;
double r570926 = r570925 + r570914;
double r570927 = sqrt(r570926);
double r570928 = sqrt(r570925);
double r570929 = r570927 - r570928;
double r570930 = r570924 + r570929;
double r570931 = t;
double r570932 = r570931 + r570914;
double r570933 = sqrt(r570932);
double r570934 = sqrt(r570931);
double r570935 = r570933 - r570934;
double r570936 = r570930 + r570935;
return r570936;
}