function [C,D]=newpoly(x,y) %牛頓插值
n=length(x);
D=zeros(n);
D(:,1)=y';
for j=2:n  %計(jì)算差商矩陣
    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  %構(gòu)造插值多項(xiàng)式
    C=conv(C,poly(x(k)));
    m=length(C);
    C(m)=C(m)+D(k,k);
end