Having chosen the acceptance probabilities, we have guaranteed that if the P(n)i has equilibrated, that is if it is equal to pi, it will remain equilibrated. Next we need to find the circumstances for convergence to equilibrium.
The dynamical equation can be written as P(n)i=∑jMijP(n−1)j with the matrix M given by Mij=δij[1−∑kTi→kAi→k]+Tj→iAj→i. Summing over i shows that ∑iMij=1, and since ∑kTi→k=1, and Ai→k≤1, the elements of the matrix satisfy Mij≥0. The matrix M is therefore a stochastic matrix.