In [1]:
from sympy import *
init_printing(use_unicode=False, wrap_line=False)
x = Symbol('x')
a = Symbol('a')

In [10]:
%matplotlib notebook

In [3]:
simplify( integrate( a / (a**2 + x**2)**(3/2) ,   a) )


Out[3]:
$$- \frac{1.0}{x^{1.0} \left(\frac{a^{2}}{x^{2}} + 1\right)^{0.5}}$$

In [12]:
def E(x):
    return  x*(1/abs(x) - 1/sqrt( x**2 + 1 ) )

In [15]:
plot(E(x) , (x,-3,3))


Out[15]:
<sympy.plotting.plot.Plot at 0x7fa09e94b4a8>

In [ ]: