function yy=nalagr(x,y,xx)
%用途： Lagrange插值法數(shù)值求解
%格式： yy=nalagr(x,y,xx) x是節(jié)點(diǎn)向量，y是節(jié)點(diǎn)上的函數(shù)插值，xx是插值點(diǎn)
% （可以是多個(gè)），yy返回插值
m=length(x);n=length(y);
if m~=n, error('x!=y');end
s=0;
for i=1:n
    t=ones(1,length(xx));
    for j=1:n
        if j~=i,
            t=t.*(xx-x(j))/(x(i)-x(j));
        end
    end
    s=s+t*y(i);
end
yy=s;