4-8.m

来自「MATLAB科学计算与工程分析源代码源程序3」· M 代码 · 共 14 行

 x=[-5:0.1:5];
y=1./(1+x.^2);
x=[-5:1:5];
y=1./(1+x.^2);
x0=[-5:0.1:5];
y0=lagrange(x,y,x0);
y1=1./(1+x0.^2);
y2=interp1(x,y,x0);

plot(x0,y0,'o');
hold on
plot(x0,y1,'--');
hold on
plot(x0,y2,'*')