function p = aitken_neville(x,f,t)

% function p = aitken_neville(x,f,t)
% polynomial interpolant p(t) to data (x,f)

p = f; 
for k=1:length(f)-1;
   w = (x(1+k:end)-t) ./ (x(1+k:end)-x(1:end-k));
   p = w.*p(1:end-1) + (1-w).*p(2:end);
end;
