function v = ANipoleval(x,y,t)
% Aitken-Neville algorithms for evaluation of interpolating polynomial
% x passed interpolation nodes, y the data values, and t the evaluation 
% point 
for i=1:length(y)
  for k=i-1:-1:1
      y(k) = y(k+1)+(y(k+1)-y(k))*(t-x(i))/(x(i)-x(k));
  end
end
v = y(1);