function [c]=nwtcoeff(y)
% calcola le differenze finite di y
% se y ha n elementi, c ha n-1 elementi
[nn,mm]=size(y);
if ((nn ~=1)&(mm ~=1)) error('y deve essere un vettore');end;
n=nn+mm-1;
c(1)=y(1);
q=1;
m=n;
while(m>2)
    q=q+1;
    for k=2:m
        y(k-1)=y(k)-y(k-1);
    end
c(q)=y(1);
m=m-1;
% debug
% y(1:m)
end

end
