% tester igjen potensmetoden uten skalering

clear all;

figure(2);

L=10; U=10;
axis([-L U -L U]);
axis manual;
A=[-L U; 0 0];
B=[0 0; -L U];
plot(A(1,:),A(2,:));
hold on;
plot(B(1,:),B(2,:));
grid;
hold on;

A=[1,2; -1, -1]   %A har ingen reelle egenverdier

[V D]=eig(A)    % Egenverdiene til A er i og -i, 
		         % saa A har ikke en strengt dom. egenverdi
                 
x=[-2 ;4]  		% x = startvektor

k=8			    % k= antall iterasjoner

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

        % x-ene  tar  bare fire verdier som det byttes imellom. 

x=[1 ;-5]  		% x = ny startvektor



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