1、数值分析第二次上机题目：在区间-1,1上分别取 n=10、20 用两组等距节点对龙格函数 做多251)(xf项式插值及三次样条插值，对每个 n 值，分别画出插值函数及 的图形x程序：function y1=lagrange(x0,y0,x1)n=length(x0);syms x;for k=1:nl(k)=x/x;for p=1:nif p=kl(k)=l(k)*(x-x0(p)/(x0(k)-x0(p);endendendz=0;for k=1:nz=z+l(k)*y0(k);endy1=subs(z,x,x1);程序：n=10 时n=10;x0=-1:2/n:1;y0=1./(1+25*
2、x0.2);x=-1:.001:1;y1=lagrange(x0,y0,x);y=1./(1+25*x.2);plot(x,y,x,y1,-.,x0,y0,p);legend(Runge Function,插值函数, 插值节点);title(n=10 时的 Lagrange 插值的龙格现象);xlabel(x);ylabel(y);结果：程序：n=20 时n=20;x0=-1:2/n:1;y0=1./(1+25*x0.2);x=-1:.001:1;y1=lagrange(x0,y0,x);y=1./(1+25*x.2);plot(x,y,x,y1,-.,x0,y0,p);legend(Runge Function,插值函数, 插值节点);title(n=20 时的 Lagrange 插值的龙格现象);xlabel(x);ylabel(y);结果：
