a = [2 4; 4 3];
//a = [4 2; 2 3];
m = 10;
t = 0:%pi*2/m:2*%pi;
z = [cos(t); sin(t)];
scf(1);
clf();
n = 10;
for i = 1:n,
	clf();
	x = z;
	w = a*x;
	z = zeros(w);
	for j = 1:(m+1),
		z(:,j) = w(:,j)/norm(w(:,j));
	end;
	plot(x(1,:),x(2,:),'ro');
	plot(w(1,:),w(2,:),'rx');
	for j = 1:(m+1),
		plot([x(1,j),w(1,j)],[x(2,j),w(2,j)],'r-');
	end;
	pause;
end;