integrate erf( a / sqrt(x^2+b^2)) from x = 0 to infinity