// saisie de la matrice de transition A=[0.5,0.5,0.1;0.4,0.3,0.2;0.1,0.2,0.7]; // saisie de la matrice identité I=[1,0,0;0,1,0;0,0,1]; W=A-I; // calcul du noyau de A-I X=kernel(W); // procédure de normalisation du vecteur norme=X(1)+X(2)+X(3); stationnaire=X/norme; // affichage des composantes de la distribution stationnaire afficher(stationnaire);