前辈们，帮忙看看助我一臂之力

张翔10

天字：一：号

所有的代码在这里了： function[C,L]=lagran(x,y);w=length(x);n=w-1;L=zeros(w,w);for k=1:n+1;    V=1;    for j=1:n+1;        if k~=j

我应该

V=conv(V,ploy(x(j))/(x(k)-x(j)));        end    end    L(k,:)=VendC=y*Lclear allclcclfx=-5:0.1:5;y=5./(1+x.^2);[C,L]=lagran(x,y);xx=-5:0.5:5;yy=ployval(C,xx);hold onplot(xx,yy,'b',x,y,'.')xp=-5:0.01:5;z=5./(1+xp.^2);plot(xp,z,'r')grid onlegend ('lagran','差值节点','原函数')