x=np.linspace(-1,1,100)y=(x+1)*np.exp(2*x)y1=1+3*xy2=1+3*x+8*x**2plt.figure(figsize=[8,5])plt.plot(x,y,'k',linewidth=3,label='Function')plt.plot(x,y1,'b',linewidth=3,label='1st order approximation')plt.plot(x,y2,'r',linewidth=3,label='2nd order approximation')plt.legend()plt.xticks([-1,0,1])plt.show()