X(:,1) = [1 1 1 1]';

 A = [ 0 0 0.1 .5; 0.2 0 0 0; 0 .9 0 0;0 0 .9 0.85];

 for i =1:30
  X(:,i+1) = A*X(:,i);
 end