double f(double x, double y, double z, double t) {
double r483850 = x;
double r483851 = 1.0;
double r483852 = r483850 + r483851;
double r483853 = sqrt(r483852);
double r483854 = sqrt(r483850);
double r483855 = r483853 - r483854;
double r483856 = y;
double r483857 = r483856 + r483851;
double r483858 = sqrt(r483857);
double r483859 = sqrt(r483856);
double r483860 = r483858 - r483859;
double r483861 = r483855 + r483860;
double r483862 = z;
double r483863 = r483862 + r483851;
double r483864 = sqrt(r483863);
double r483865 = sqrt(r483862);
double r483866 = r483864 - r483865;
double r483867 = r483861 + r483866;
double r483868 = t;
double r483869 = r483868 + r483851;
double r483870 = sqrt(r483869);
double r483871 = sqrt(r483868);
double r483872 = r483870 - r483871;
double r483873 = r483867 + r483872;
return r483873;
}