clear all
clf
A=[0 1; -1 0]
EA=eig(A)
x=[6;4]
for k=0:20   
B(:,k+1)=A^k*x;
end
B;
plot(B(1,:),B(2,:),'r-*')