double f(double x, double y, double z, double t) {
double r296851 = x;
double r296852 = 1.0;
double r296853 = r296851 + r296852;
double r296854 = sqrt(r296853);
double r296855 = sqrt(r296851);
double r296856 = r296854 - r296855;
double r296857 = y;
double r296858 = r296857 + r296852;
double r296859 = sqrt(r296858);
double r296860 = sqrt(r296857);
double r296861 = r296859 - r296860;
double r296862 = r296856 + r296861;
double r296863 = z;
double r296864 = r296863 + r296852;
double r296865 = sqrt(r296864);
double r296866 = sqrt(r296863);
double r296867 = r296865 - r296866;
double r296868 = r296862 + r296867;
double r296869 = t;
double r296870 = r296869 + r296852;
double r296871 = sqrt(r296870);
double r296872 = sqrt(r296869);
double r296873 = r296871 - r296872;
double r296874 = r296868 + r296873;
return r296874;
}