newtondd(일반적 분할차분표)

function b=newtondd(x, y)

% x,y: 입력데이터

n=length(x); b=zeros(n,n); b(:,1)=y(:);

%분할차분표 작성

for i=2:n

    for j=1:n-i+1

        b(j,i)=(b(j+1,i-1)-b(j,i-1))/(x(i+j-1)-x(j));

    end

end