建立如下拉格朗日插值函数： function y=lagrange(x0,y0,x); n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j)); end end
s=p*y0(k)+s; end y(i)=s; end
画图程序如下： x=[-5:1:5]; y=1./(1+x.^2); x0=[-5:0.001:5]; y0=lagrange(x,y,x0); y1=1./(1+x0.^2); plot(x0,y0,'r') hold on
--精品--
plot(x0,y1,'g')
注：画出的图形为n =10的图形
得到图形如下：
n=10的图像
--精品--
因篇幅问题不能全部显示，请点此查看更多更全内容