function P=Flagrange(X,Y)//X est le vecteur qui contient le x_i et Y celui qui contient les y_i;
//P est le polynôme d'interpolation de Lagrange
n=length(X);// n+1 est le nombre de points. n sera donc le degré du polynôme P 
x=poly(0,"x");P=0;//x est l'indéterminée du polynôme P 
for i=1:n, L=1; 
for j=[1:i-1,i+1:n] L=L*(x-X(j))/(X(i)-X(j));end //on a effectué les différences divisées
P=P+L*Y(i); 
end 
disp(P)// on affiche P
endfunction
//on prend 3 points (1;6) ,(-3;-8) et (8;-4) [voir l'exemple traité en début de document]
X=[1,-3,8]
Y=[6,-8,-4]
P=Flagrange(X,Y)