

			M = [ 	.95 0 .05;
				 	.05 .8 .15;
					0 .2 .8]; 


			K = 3;
			M = -10+20*rand(K);
			for k = 1:K; M(:,k) = M(:,k)/norm(M(:,k)/1); end;

			

			v = rand(K,1);
			v = v/norm(v,1);

			for i = 1:500;

					w = v;
					v = M*v;
					v = v/norm(v,1);
					
					if norm(v - w) < .01 || rank([v,w]) == 1
							disp(['Convergence after ',num2str(i),' iterations.']);
							return
					end

					
			end
