function y=neville(x,xi,fi)

% Neville's iterated interpolation: See section 3.1.
% First input is the point where P(x) should be evaluated.
% Second input is the vector of nodes (sample x-coordinates).
% Third input is either the corresponding vector of y-coordinates,
%  or the function f.

n=length(xi)-1;
Q=neville_table(x,xi,fi);

y=Q(n+1,n+1);