integrate (sqrt(x^2+1)/x)