function DV=DiffDiv(X,Y); 

% DV = DiffDiv(X,Y); 
% Calcule (dans l'ordre inverse) la liste des N+1 différences divisées des
% entrées du vecteur ligne Y (dans cet ordre) par les entrées du vecteur
% ligne X. Les vecteurs X et Y sont des vecteurs ligne de même longueur N+1,
% dont les entrées sont des nombres complexes. Les entrées de X sont
% distinctes (TP5, exercice 1, question 1). 

N=size(X,2)-1; 
U=zeros(N+1,N+2); 
U(:,1)= conj(X'); 
U(:,2)= conj(Y'); 
DV = [Y(1)]; 
for j=2:N+1
    for i=1:N+2-j 
        U(i,j+1)= (U(i,j)-U(i+1,j))/(U(i,1) - U(i+j-1,1)); 
    end
    DV = [U(1,j+1) DV]; 
end

    