double f(double x, double y, double z, double t) {
double r614937 = x;
double r614938 = 1.0;
double r614939 = r614937 + r614938;
double r614940 = sqrt(r614939);
double r614941 = sqrt(r614937);
double r614942 = r614940 - r614941;
double r614943 = y;
double r614944 = r614943 + r614938;
double r614945 = sqrt(r614944);
double r614946 = sqrt(r614943);
double r614947 = r614945 - r614946;
double r614948 = r614942 + r614947;
double r614949 = z;
double r614950 = r614949 + r614938;
double r614951 = sqrt(r614950);
double r614952 = sqrt(r614949);
double r614953 = r614951 - r614952;
double r614954 = r614948 + r614953;
double r614955 = t;
double r614956 = r614955 + r614938;
double r614957 = sqrt(r614956);
double r614958 = sqrt(r614955);
double r614959 = r614957 - r614958;
double r614960 = r614954 + r614959;
return r614960;
}