% polinomio soluzione 1+x+x^3
x = [0 1 -3 2]';
y = [1 3 -29 11]';
dd = differenze_divise(x, y) ;
dd

% polinomio interpolante
coeffs = interp_poly( x, y ) ;
coeffs