x=linspace(0,1,6);
y=rand(size(x));

DD=DDtable(x,y);
c=diag(DD);

plotnewtoninterp(x,y,c);