
% After EigenPlay

x = -1+2*rand(2,1);
x0 = x;
f = figure; hold on;
xlim([-2.5,2.5]);
ylim([-2.5,2.5]);

for i = 1:500

		b = mod(round(100*rand(1)),2);
		b = 1;
		if b == 1
			A = a;
			%A = [1/2 3/4; -3/5 11/10];
		else
			r = sqrt(5);
			%A = [2/r -1/r; 1/r 2/r];
			A = [1/2 3/3.9; -3/5 11/10];
		end

		plot(x(1),x(2),'*','MarkerSize',5);
		x = A*x;
		drawnow
		disp(['A^',num2str(i),'(',num2str(x0(1)),',',num2str(x0(2)),')']);
		pause
end
		



