function LXX = LagrangeNewton(X,Y,XX); 

%function LXXN = LagrangeNewton(X,Y,XX); 
%
%Evalue le polynôme de Lagrange de degré N (interpolant aux entrées de X
%les entr�es de Y, X et Y étant deux vecteurs de même longueur N+1) en les
%entrées sucessives d'un vecteur ligne XX. Le polynôme de Lagrange est ici
%calculé suivant la méthode de Newton. Ce programme utilise comme sous
%programme la routine DiffDiv (TP5, exercice 1, question 2). 

A=DiffDiv(X,Y); 
M=length(A);
MM=length(XX); 
LXX = A(1)*ones(1,MM);  
for i=2:M
    LXX=LXX.*(XX-X(M-i+1))+A(i); 
end
