examp8_3.m

来自「东北大学 薛定谔的高等应用数学问题MATLAB求解的源代码」· M 代码 · 共 6 行

x0=-1+2*[0:10]/10; y0=1./(1+25*x0.^2);
x=-1:.01:1; y=lagrange(x0,y0,x);  % Lagrange 插值
ya=1./(1+25*x.^2); plot(x,ya,x,y,':')

y1=interp1(x0,y0,x,'cubic'); y2=interp1(x0,y0,x,'spline');
plot(x,ya,x,y1,':',x,y2,'--')