integral (sqrt(1 + 4/a^4*x^2*e^(-2x^2/a^2)) )