// saisie de la distribution initiale en janvier
U=[0.5;0.3;0.2];
//saisie de la matrice de transition 3X3
A=[0.5,0.5,0.1;0.4,0.3,0.2;0.1,0.2,0.7];
//saisie des composantes de la distribution stationnaire
f=17/46; g=7/23; h=15/46; 
//calcul de la distribution à la fin du premier mois
V=A*U;
//initialisation du compteur d'itérations(nombre de mois écoulés)
i=1;
//boucle for sur la précision attendue des calculs 
for k=1:20
    //boucle while pour la condition d'arrêt
    //tant que nous sommes éloignés de la distribution stationnaire, on applique A 
while abs(V(1)-f)>10^(-k) 
    while abs(V(2)-g)>10^(-k)
        while abs(V(3)-h)>10^-k i=i+1; V=A^i*U; 
        end 
    end 
end
//Affichage du nombre de mois écoulés, de la précision et de la distribution obtenue
afficher("En "+string(i)+" itérations, le système est à 10^"+string(-k)+" près de la ..
distribution stationnaire.")
afficher(V(1)); afficher(V(2)); afficher(V(3));
end