我们要研究前向差分和中心差分的差分逼近中的误差,将误差制表h=[1.E-3 1.E-4 1.E-5 1.E-6 1.E-7 1.E-8 1.E-9 1.E-10 1.E-11 1.E-12 1.E-13]
并绘制对数对数图。关于如何做到这一点的任何提示?这是我们的中锋和前锋的区别。
centdiff=(subs(f, x+h))/(2*h) - (subs(f, x-h))/(2*h)
framdiff=(subs(f, x+h) - f)/h
而我们的功能:
f=60*x-(x.^2+x+0.1).^6./(x+1).^6-10*x.*exp(-x);