double f(double x, double y, double z, double t) {
double r33141022 = x;
double r33141023 = 1.0;
double r33141024 = r33141022 + r33141023;
double r33141025 = sqrt(r33141024);
double r33141026 = sqrt(r33141022);
double r33141027 = r33141025 - r33141026;
double r33141028 = y;
double r33141029 = r33141028 + r33141023;
double r33141030 = sqrt(r33141029);
double r33141031 = sqrt(r33141028);
double r33141032 = r33141030 - r33141031;
double r33141033 = r33141027 + r33141032;
double r33141034 = z;
double r33141035 = r33141034 + r33141023;
double r33141036 = sqrt(r33141035);
double r33141037 = sqrt(r33141034);
double r33141038 = r33141036 - r33141037;
double r33141039 = r33141033 + r33141038;
double r33141040 = t;
double r33141041 = r33141040 + r33141023;
double r33141042 = sqrt(r33141041);
double r33141043 = sqrt(r33141040);
double r33141044 = r33141042 - r33141043;
double r33141045 = r33141039 + r33141044;
return r33141045;
}