csqrt
double complex csqrt(double complex);
w = clog(z + csqrt(z + 1) * csqrt(z - 1));
z2 = csqrt(zz);