double f(double x, double y, double z, double t) {
double r570604 = x;
double r570605 = 1.0;
double r570606 = r570604 + r570605;
double r570607 = sqrt(r570606);
double r570608 = sqrt(r570604);
double r570609 = r570607 - r570608;
double r570610 = y;
double r570611 = r570610 + r570605;
double r570612 = sqrt(r570611);
double r570613 = sqrt(r570610);
double r570614 = r570612 - r570613;
double r570615 = r570609 + r570614;
double r570616 = z;
double r570617 = r570616 + r570605;
double r570618 = sqrt(r570617);
double r570619 = sqrt(r570616);
double r570620 = r570618 - r570619;
double r570621 = r570615 + r570620;
double r570622 = t;
double r570623 = r570622 + r570605;
double r570624 = sqrt(r570623);
double r570625 = sqrt(r570622);
double r570626 = r570624 - r570625;
double r570627 = r570621 + r570626;
return r570627;
}