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