double f(double x, double y, double z, double t) {
double r23155790 = x;
double r23155791 = 1.0;
double r23155792 = r23155790 + r23155791;
double r23155793 = sqrt(r23155792);
double r23155794 = sqrt(r23155790);
double r23155795 = r23155793 - r23155794;
double r23155796 = y;
double r23155797 = r23155796 + r23155791;
double r23155798 = sqrt(r23155797);
double r23155799 = sqrt(r23155796);
double r23155800 = r23155798 - r23155799;
double r23155801 = r23155795 + r23155800;
double r23155802 = z;
double r23155803 = r23155802 + r23155791;
double r23155804 = sqrt(r23155803);
double r23155805 = sqrt(r23155802);
double r23155806 = r23155804 - r23155805;
double r23155807 = r23155801 + r23155806;
double r23155808 = t;
double r23155809 = r23155808 + r23155791;
double r23155810 = sqrt(r23155809);
double r23155811 = sqrt(r23155808);
double r23155812 = r23155810 - r23155811;
double r23155813 = r23155807 + r23155812;
return r23155813;
}