function LXXN = NevilleAitken(a,b,N,f); 

% LXXN=NevilleAitken (a,b,N,f); 
% 
%Etant donné un segment [a,b], un entier strictement positif N, et une
%fonction de [a,b] dans R ou C déclarée préalablement 'inline', cette 
%routine calcule (suivant la méthode récursive de Neville-Aitken) le
%polynôme de Lagrange de degré N interpolant les valeurs de f aux N+1
%points équidistribués entre a et b (extrémités comprises).  

if N>1 
    LXXNaux1=NevilleAitken(a+(b-a)/N,b,N-1,f); 
    LXXNaux2=NevilleAitken(a,b-(b-a)/N,N-1,f);
    LXXNaux3=[LXXNaux1 0] - a*[0 LXXNaux1];
    LXXNaux4=[LXXNaux2 0] - b*[0 LXXNaux2];
    LXXN=(LXXNaux3-LXXNaux4)/(b-a); 
else
    LXXN=[f(b)-f(a) b*f(a)-a*f(b)]/(b-a); 
end





