double f(double x, double y, double z, double t) {
double r753789 = x;
double r753790 = 1.0;
double r753791 = r753789 + r753790;
double r753792 = sqrt(r753791);
double r753793 = sqrt(r753789);
double r753794 = r753792 - r753793;
double r753795 = y;
double r753796 = r753795 + r753790;
double r753797 = sqrt(r753796);
double r753798 = sqrt(r753795);
double r753799 = r753797 - r753798;
double r753800 = r753794 + r753799;
double r753801 = z;
double r753802 = r753801 + r753790;
double r753803 = sqrt(r753802);
double r753804 = sqrt(r753801);
double r753805 = r753803 - r753804;
double r753806 = r753800 + r753805;
double r753807 = t;
double r753808 = r753807 + r753790;
double r753809 = sqrt(r753808);
double r753810 = sqrt(r753807);
double r753811 = r753809 - r753810;
double r753812 = r753806 + r753811;
return r753812;
}