A=[ 0 1/2 0 0 0 0; 0 0 1/3 0 1/2 0; 1/4 1/2 0 0 0 0; 1/4 0 0 0 1/2 0; 1/4 0 1/3 0 0 1; 1/4 0 1/3 1 0 0; ]; %Matrix A from lab p1=[1;0;0;0;0;0]; p2=[0;1;0;0;0;0]; p3=[0;0;1;0;0;0]; p4=[0;0;0;1;0;0]; p5=[0;0;0;0;1;0]; p6=[0;0;0;0;0;1]; for n=1:100 p1=A*p1; p2=A*p2; p3=A*p3; p4=A*p4; p5=A*p5; end