function XX=Jacobi(M,B,X,k); 

% XX=Jacobi(M,B,X,k); 
% Génère l'algorithme it�ratif de Jacobi à partir d'une
% matrice carrée M (réelle ou complexe) inversible dont tous les termes
% diagonaux sont non nuls. Ce programme fournit la solution XX du système de
% Cramer M*XX=B. L'algorithme est initié à X et itéré k fois. Cet algorithme
% converge en particulier si M est à diagonale dominante.

D=diag(diag(M));
E=D-M; 
XX=X; 
for i=1:k
    XX=D^(-1)*E*XX+D^(-1)*B; 
end
