
x=[-2; 1; 3];
fx=[-2; 11; 17];

coeff=newt_intp(x,fx)

pol_coeffs=diag(coeff);