integrate x^2/(sqrt(x^2+a^2)*(sqrt(x^2+a^2)-b))