from future import division from sympy import from sympy import diff, Symbol, sin
x, k= symbols(‘x k ’, commutative=True) m = 3 iter = 50 f0= sin(x) fx= sin(x) sol = f0.evalf(subs={’x‘: 0}) real= f0.evalf(subs={’x’:m}) for i in range(1,iter): fx= diff(fx,x) k = fx.evalf(subs={’x’: 0})/math.factorial(i) sol = sol + (3i)k