double f(double x, double y, double z, double t) {
double r450852 = x;
double r450853 = 1.0;
double r450854 = r450852 + r450853;
double r450855 = sqrt(r450854);
double r450856 = sqrt(r450852);
double r450857 = r450855 - r450856;
double r450858 = y;
double r450859 = r450858 + r450853;
double r450860 = sqrt(r450859);
double r450861 = sqrt(r450858);
double r450862 = r450860 - r450861;
double r450863 = r450857 + r450862;
double r450864 = z;
double r450865 = r450864 + r450853;
double r450866 = sqrt(r450865);
double r450867 = sqrt(r450864);
double r450868 = r450866 - r450867;
double r450869 = r450863 + r450868;
double r450870 = t;
double r450871 = r450870 + r450853;
double r450872 = sqrt(r450871);
double r450873 = sqrt(r450870);
double r450874 = r450872 - r450873;
double r450875 = r450869 + r450874;
return r450875;
}