clear P
N=10;
A=[.95 .03;.05 .97]
x=[600000;400000]
A*x
plot(x(1),x(2),'r*')
hold on
for j=1:N;
    P(:,j)=A^j*x;
plot(P(1,j), P(2,j),'b*')
end
P
