clear
%initialisation
nbr_pas =20;    %nombre de pas de temps
nbr_puces =2;
chain(1)=nbr_puces;        %?tat initial

%parcours de la chaine
for i=1:nbr_pas
    %r?cup?ration de l'?tat actuel
    chien1=chain(i);
    chien2=nbr_puces-chien1;
    
    %calcul des probabilit?s de transition (proportionnelle au nombre de puce)
    chien1_vers_2=chien1/nbr_puces;
    chien2_vers_1=1-chien1_vers_2;
    
    %passage ? l'?tat suivant
    if rand()<chien1_vers_2     %une puce passe du chien 1 au chien 2
        chain(i+1)=chien1-1;
    else                        %une puce passe du chien 2 au chien 1
        chain(i+1)=chien1+1;
    end
end

%Analyse des r?sultats
plot(chain)
axis([0 nbr_pas 0 nbr_puces])
saveas(gcf,'Markov1','pdf')