% tester potensmetoden uten skalering

 clear all;
 
 figure(1);

 hold on;

 grid;


A=[0.5,1.5; 1, 1]

[V D]=eig(A)    % Egenverdiene til A er -0.5 og 2

x=[-10;10];  		% x = startvektor

k=6;			      % k= antall iterasjoner


%plotter del av egenrom utspent av V(:,2)

plot([-100*V(1,2),50*V(1,2)],[-100*V(2,2),50*V(2,2)],'r'); hold on;

for j=1:k
  plot([0 x(1)],[0 x(2)], 'g');
  hold on;pause(1)
  x=A*x;
end

x=[5;-5];

for j=1:k
  plot([0 x(1)],[0 x(2)], 'b');
  hold on;pause(1)
  x=A*x;
end
