function [XX,Niter] = ExempleGaussSeidel(N,B,X,epsilon,k); 

% [XX,Niter] = ExempleGaussSeidel(N,B,X,epsilon,k);
%
%Génère l'algorithme de Gauss-Seidel lorsque M est la 
%matrice de taille (N,N) à diagonale dominante constituée de 3, la 
%sur et la sous-diagonale étant constituée de -1. 
%Le second membre du système de Cramer est B. L'algorithme 
%est initié à X; le nombre maximal d'itérations autorisées 
%est k; Niter correspond au nombre d'itérations nécessaires 
%pour que l'erreur entre un itéré et son successeur passe sous
%le seuil de epsilon. 

M=diag(3*ones(N,1)) - diag(ones(N-1,1),1) - diag(ones(N-1,1),-1); 
[XX,Niter] = GaussSeidel1(M,B,X,epsilon,k); 