function [coe,DD]=newtonp(x,y)
%Newton 插值函數
%輸入參數 x= [x0 x1 ... xN]
%        y= [y0 y1 ... xN]
%輸出參數 n= Newton系數
N=length(x)-1;  %這里說明了給定的參數的個數就決定了數據擬合的階數
DD=zeros(N+1,N+1);
DD(1:N+1,1)=y';
for k=2:N+1
    for m=1:N+2-k
        DD(m,k)=(DD(m+1,k-1)-DD(m,k-1))/(x(m+k-1)-x(m));
    end
end
a=DD(1,:);
coe=a(N+1);
for k=N:-1:1
    coe=[coe a(k)]-[0 coe*x(k)];
end


%補兗：牛頓擬合的原理
%對于給定的數據系列{（x0,y0),(x1,y1),...(xN,yN)}可以被下面的多項式循環引用，滿足
%下面的N個數據系列{（x0,y0),(x1,y1),...(xN,yN)}的等式
%nN(x)=a0+a1(x-x0)+a2(x-x0)(x-x1)+...=nN-1(x)+aN(x-x0)(x-x1)...(x-xN-1)同時滿足
%n0(x)=a0
%由此得到
%a0=y0
%a1=(y1-a0)/(x1-x0)=(y1-y0)/(x1-x0)=Df0
%a2=(Df1-Df0)/(x2-x0)=D^2f0
%aN=D^Nf0
%