integrate sqrt(x-1/x+1)*dx