int 1/(a^2+x^2) dx = (1/a)*arctg (x/a)+C