plot int (1/sqrt(x^2-1)) dx, arccosh(x)