数值分析程序李庆扬版.docVIP

数值分析部分程序 Lagrange function y=Lagrange(x,y,xstar) n=length(y); s=0; for i=1:n z1=1;z2=1; for j=1:n if j~=i z1=z1.*(xstar-x(j)); z2=z2.*(x(i)-x(j)); end end z=z1./z2; s=s+z*y(i); end y=s; x=[0.5610,0.5628,0.56401,0.56521]; y=[0.83741,0.82659,0.82577,0.82495]; xstar=[0.5625,0.5635,0.5645]; y=Lagrange(x,y,xstar) Runge t=-5:0.1:5; ft=1./(1+t.*t); t1=-5:1:5; ft1=1./(1+t1.*t1); y1=Lagrange(t1,ft1,t); plot(t,ft,t,y1); QJS function b=QJS(A,x) n=length(A); b(1)=A(1); c(1)=A(1); for k=1:1:n-1; b(1+k)=b(k)*x+A(k+1); end s=b(n) for i=1:1:n-2; c(i+