integrate x*sqrt(4x^2/(x^2+a^2)^2+1) from 0 to 1