double f(double x, double y, double z, double t) {
double r23180 = x;
double r23181 = 1.0;
double r23182 = r23180 + r23181;
double r23183 = sqrt(r23182);
double r23184 = sqrt(r23180);
double r23185 = r23183 - r23184;
double r23186 = y;
double r23187 = r23186 + r23181;
double r23188 = sqrt(r23187);
double r23189 = sqrt(r23186);
double r23190 = r23188 - r23189;
double r23191 = r23185 + r23190;
double r23192 = z;
double r23193 = r23192 + r23181;
double r23194 = sqrt(r23193);
double r23195 = sqrt(r23192);
double r23196 = r23194 - r23195;
double r23197 = r23191 + r23196;
double r23198 = t;
double r23199 = r23198 + r23181;
double r23200 = sqrt(r23199);
double r23201 = sqrt(r23198);
double r23202 = r23200 - r23201;
double r23203 = r23197 + r23202;
return r23203;
}