double f(double x, double y, double z, double t) {
double r510109 = x;
double r510110 = 1.0;
double r510111 = r510109 + r510110;
double r510112 = sqrt(r510111);
double r510113 = sqrt(r510109);
double r510114 = r510112 - r510113;
double r510115 = y;
double r510116 = r510115 + r510110;
double r510117 = sqrt(r510116);
double r510118 = sqrt(r510115);
double r510119 = r510117 - r510118;
double r510120 = r510114 + r510119;
double r510121 = z;
double r510122 = r510121 + r510110;
double r510123 = sqrt(r510122);
double r510124 = sqrt(r510121);
double r510125 = r510123 - r510124;
double r510126 = r510120 + r510125;
double r510127 = t;
double r510128 = r510127 + r510110;
double r510129 = sqrt(r510128);
double r510130 = sqrt(r510127);
double r510131 = r510129 - r510130;
double r510132 = r510126 + r510131;
return r510132;
}