%A=[0,0,0,0,0,1/3,1/4;1/3,0,1/3,0,0,0,1/4;0,1,0,0,0,1/3,1/4;0,0,1/3,0,1/2,0,0;0,0,1/3,1/2,0,1/2,0;1/3,0,0,1/2,0,0,1/4;1/3,0,0,0,0,0,0]
R= transpose([0 1/3 0 0 0 1/3 1/3])
x=transpose([1/7 1/7 1/7 1/7 1/7 1/7 1/7])

while R.^k.*x ~= R.^(k+1).*x
k=k+1;
end
k
R.^k.*x
R.^(k+1).*x