function [b]=nwtpol(x,x1,h,y)
% calcola il polinomio di interpolazione di Newton 
% a partire dal vettore delle differenze finte y
% dal punto di inizio x1 e dal passo h
a=(x-x1)/h;
[nn,mm]=size(y);
if ((nn ~=1)&(mm ~=1)) error('y deve essere un vettore');end;
n=nn+mm-1;
b=y(n);
while (n>1)
   n=n-1;
   b=y(n) + b*(a-n+1)/n; 
end
end