Initialement, les N boules sont dans l'urne A.
A chaque étape k, il y a un nombre N_A de boules dans l'urne A.
La matrice ligne L_k correspond au nombre de boules dans l'urne A à chaque étape k : L_k=l_{1,j} avec l_{1,j}=1 si j=N_A et 0 sinon avec j compris entre 0 et N.
La matrice de passage M d'une étape à la suivante sera soit M_{AB} soit M_{BA} selon que la boule choisie au hasard se situe dans l'urne A ou l'urne B.
M_{AB} diminue le nombre de boules de A, sauf s'il n'y en a pas, auquel cas elle l'augmente. M_{BA} fait le contraire.
On génère donc un nombre aléatoire p. Si p<=N_A/N alors on utilise la matrice de passage M_{AB} et sinon on utilise M_{BA}.