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 [ ]: