integrate (x^2 + 1)^(1/2) dx from 4000pi to 4002pi