double f(double x, double y, double z, double t) {
double r445200 = x;
double r445201 = 1.0;
double r445202 = r445200 + r445201;
double r445203 = sqrt(r445202);
double r445204 = sqrt(r445200);
double r445205 = r445203 - r445204;
double r445206 = y;
double r445207 = r445206 + r445201;
double r445208 = sqrt(r445207);
double r445209 = sqrt(r445206);
double r445210 = r445208 - r445209;
double r445211 = r445205 + r445210;
double r445212 = z;
double r445213 = r445212 + r445201;
double r445214 = sqrt(r445213);
double r445215 = sqrt(r445212);
double r445216 = r445214 - r445215;
double r445217 = r445211 + r445216;
double r445218 = t;
double r445219 = r445218 + r445201;
double r445220 = sqrt(r445219);
double r445221 = sqrt(r445218);
double r445222 = r445220 - r445221;
double r445223 = r445217 + r445222;
return r445223;
}