clf
N=60;
M=20;
S=30;
X=zeros(2,N);
B=rand(2,2)
A=[ 1.1 0; 0 .7];
A=B*A*inv(B)
for k=1:M
x=[S*(rand-1/2);S*(rand-1/2)];
 for j = 1:N
    X(:,j)=A^(j-1)*x;     
 end   
plot(X(1,:),X(2,:),'g-*')
hold on
plot(X(1,1),X(2,1),'b*')  
plot(X(1,N),X(2,N),'r*')
%xlabel({'Owls',''});
%ylabel({'Rats',''});
%axis([-16 16 -16 16]);
end