func3p1.asv

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

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,'-',xx,yy1,'-','MarkerEdgeColor','g')%,yy2,xx,'-','MarkerEdgeColor','g',yy3,xx,'-','MarkerEdgeColor','y')