double f(double x, double y, double z, double t) {
double r560162 = x;
double r560163 = 1.0;
double r560164 = r560162 + r560163;
double r560165 = sqrt(r560164);
double r560166 = sqrt(r560162);
double r560167 = r560165 - r560166;
double r560168 = y;
double r560169 = r560168 + r560163;
double r560170 = sqrt(r560169);
double r560171 = sqrt(r560168);
double r560172 = r560170 - r560171;
double r560173 = r560167 + r560172;
double r560174 = z;
double r560175 = r560174 + r560163;
double r560176 = sqrt(r560175);
double r560177 = sqrt(r560174);
double r560178 = r560176 - r560177;
double r560179 = r560173 + r560178;
double r560180 = t;
double r560181 = r560180 + r560163;
double r560182 = sqrt(r560181);
double r560183 = sqrt(r560180);
double r560184 = r560182 - r560183;
double r560185 = r560179 + r560184;
return r560185;
}