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

 A = [0 0 3.5; 0.15 0 0; 0 .5 .9];

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