integrate sqrt(1+4x^2/b^2) dx from -b to b