examp8_3.m

《高等应用数学问题MATLAB求解》的PDF版及习题解答

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,'--')