function T=AitkenInterpolation(x,y,z)
% AITKENINTERPOLATION interpolation using the Aitken Neville scheme
%   T=AitkenInterpolation(x,y,z) uses the ponts (x,y) to generate 
%   the Aitken-Neville Scheme for the scalar interpolation value z. 

n=length(x); T=zeros(n);
for i=1:n
  T(i,1)=y(i);
  for j=1:i-1
    T(i,j+1)=((x(i)-z)*T(i-1,j)+(z-x(i-j))*T(i,j))/(x(i)-x(i-j));
  end;
end;