clear
hold off
DATA=[0 1 2 3 4 5; 4 6.1 9.8 35.2 59.4 41.7];
n=length(DATA(1,:));
A=divdifs(DATA')
X=0:.1:6;
m=length(X);
Z=A(n).*ones(1,m);
for i=n-1:-1:1
    Z=A(i)+(X-DATA(1,i)).*Z;
end
plot(X,Z,'k',DATA(1,:),DATA(2,:),'r*')
axis([0,5,0,70])
grid on