newpoly.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 18 行

function[c,d]=newpoly(x,y)  
% 这里 x为n个节点的横坐标所组成的向量，y为纵坐标所组成的向量。
% c为所求的牛顿插值多项式的系数构成的向量。
n=length(x);
d=zeros(n,n);
d(:,1)=y';
for j=2:n
    for k=j:n
     d(k,j)=(d(k,j-1)-d(k-1,j-1))/(x(k)-x(k-j+1));
  end
end
c=d(n,n);
for k=(n-1):-1:1
   c=conv(c,poly(x(k)));
   m=length(c);
   c(m)=c(m)+d(k,k);
end