function P = lagrange (x,y); 

%  P = lagrange (x,y); 
% 
% Calcule (via la méthode d'Aitken) le polynôme d'interpolation de Lagrange
% interpolant les valeurs y aux points x; les vecteurs x et y sont deux
% vecteurs lignes de même longueur N>= 2; les entrées de x sont des nombres
% complexes distincts deux à deux; les entrées de y sont des nombres
% complexes. Le polynôme P en sortie (de degré N-1) est ici exprimé comme
% la liste de ses N coefficients (dans l'ordre des monômes décroissants).

N = length(x);
if N == 1
    P= y(1);  
else
    P1 = lagrange (x(2:N),y(2:N));
    P2 = lagrange (x(1:N-1),y(1:N-1));
    P11 = [P1 0] - x(1)*[0 P1]; 
    P22 = [P2 0] - x(N)*[0 P2];
    P=(P11 - P22)/(x(N)-x(1));  
end
end 

    
    
