func3p1.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 13 行

%func3p1泰勒级数逼近
y1=[1,0]
y2=[-1/6,0,1,0]
y3=[1/120,0,-1/6,0,1,0]
x=linspace(0,pi,200)
xx=linspace(0,pi/2,100)
yy0=sin(x)
yy1=polyval(y1,xx)
yy2=polyval(y2,xx)
yy3=polyval(y3,xx)
plot(x,yy0,'.y',xx,yy1,'-r',xx,yy2,'-b',xx,yy3,'-k')
%逼近效果请见图
%在此，以三次多项式逼近最为佳