function testnewtdd
	x = [1,2,3]
	y = [2,3,1]
	plot(x,y,'b*.');
   hold on;

	n = numel(x)
	c = newtdd(x,y,n)
	
	xx = linspace(x(1),x(n),50)
	yy = nest(n-1,c,xx,x);
	plot(xx,yy,'r.');
   hold off;
end