clear all
clf
A=[.95 .03;.05 .97]
x=[600000;400000]
for k=0:100
B(:,k+1)=A^k*x;
end
B
plot(B(1,:),B(2,:),'k-*')
% hold on
% x=[200000;800000]
% for k=0:20
% B(:,k+1)=A^k*x;
% end
% B
% plot(B(1,:),B(2,:),'b-*')