% newton_table_test.m

x = [1, 3/2, 0, 2];
y = [3, 13/4, 3, 5/3];

F = newton_table(x, y)
a = newton_poly(x, y)



% >> newton_table_test
% F =
%     3.0000         0         0         0
%     3.2500    0.5000         0         0
%     3.0000    0.1667    0.3333         0
%     1.6667   -0.6667   -1.6667   -2.0000
% a =
%     3.0000    0.5000    0.3333   -2.0000
% >>