% Création d'une matrice de transition quelconque de taille N
N=3;                    % dimesion de la matrice à créer
P=rand(N,N);            % création d'une matrice à lignes non normalisées
P=diag(1./sum(P,2))*P;  % division de chaque ligne par la somme sur la ligne
P                       % affichage de la matrice obtenue
sum(P,2)                % affichage de la somme de chaque ligne (doit donner 1)
eig(P)                  % spectre de la matrice (doit être de module 1)
