Newton插值函数

%author:野渡无人
%create date:2007-12-07
function yi=newton(x,y,xi)
% x,y为已知节点向量值；yi为xi处所求插值点
d_xi=xi-x;
n=length(x);
L=zeros(n);
L(:,1)=y';
for i=2:n
    for j=1:n-i+1
        L(j,i)=(L(j+1,i-1)-L(j,i-1))/(x(i+j-1)-x(j));
    end
end
yi=L(1,1);
for i=2:n
    yi=yi+L(1,i)*prod(d_xi(1:i-1));
end