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