y= sqrt(x^2+1) - log (1/x + sqrt(x +1/x)