function y=lagrange(x0,y0,x);
a=length(x0);b=length(x);
for i=1:b
z=x(i);
s=0;
for j=1:a
p=1;
for k=1:a
if k~=j
p=p*(z-x0(k))./(x0(j)-x0(k));
end
end
s=p*y0(j)+s;
end
y(i)=s;
end
评论0