x=sp.symbols('x')
G=sp.symbols('G')
c = lambda x: G/(x**2+G**2)
u1=sp.integrate(c(x),x)
G=sp.symbols('G',positive=True)
c = lambda x: G/(x**2+G**2)
u2=sp.integrate(c(x),x)
G=sp.symbols('G',real=True)
c = lambda x: G/(x**2+G**2)
u3=sp.integrate(c(x),x)
print u1,'\n',u2,'\n',u3
This is what I get
-I*log(-I*G + x)/2 + I*log(I*G + x)/2
atan(x/G)
0