import sympy as sp x=sp.Symbol('x') div=sp.diff(x**2,x,1) sp.pprint(div) x0=10. print div.subs(x,x0)