integrate 1/((x^2 + a^2)*Sqrt[x^2 + a^2 + b^2])