double f(double x, double y, double z, double t) {
double r20850797 = x;
double r20850798 = 1.0;
double r20850799 = r20850797 + r20850798;
double r20850800 = sqrt(r20850799);
double r20850801 = sqrt(r20850797);
double r20850802 = r20850800 - r20850801;
double r20850803 = y;
double r20850804 = r20850803 + r20850798;
double r20850805 = sqrt(r20850804);
double r20850806 = sqrt(r20850803);
double r20850807 = r20850805 - r20850806;
double r20850808 = r20850802 + r20850807;
double r20850809 = z;
double r20850810 = r20850809 + r20850798;
double r20850811 = sqrt(r20850810);
double r20850812 = sqrt(r20850809);
double r20850813 = r20850811 - r20850812;
double r20850814 = r20850808 + r20850813;
double r20850815 = t;
double r20850816 = r20850815 + r20850798;
double r20850817 = sqrt(r20850816);
double r20850818 = sqrt(r20850815);
double r20850819 = r20850817 - r20850818;
double r20850820 = r20850814 + r20850819;
return r20850820;
}