clear all
clf
A=[.5 -.6;.75 1.1]
eig(A)
x=[6;4]
for k=0:40   
B(:,k+1)=A^k*x;
end
B;
plot(B(1,:),B(2,:),'r-*')
% hold on
% x=[200000;800000]
% for k=0:20
% B(:,k+1)=A^k*x;
% end
% B
% plot(B(1,:),B(2,:),'b-*')