function D=Newtons_Interp(x,y);
x=x(:); y=y(:);
n=length(x);
D=zeros(n);
D(:,1)=y;
for i=1:n-1;
	D(1:end-i,1+i)=diff(D(1:end-i+1,i))./(x(i+1:end)-x(1:end-i));
end;
D=D(1,:);
