% Eksempel, potensmetoden med skalering

clear all;

figure(3);

L=1; U=1;
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;


A=[0.5,1.5; 1, 1]   


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

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

plot([0,V(1,2)],[0,V(2,2)],'r'); 
%hold on;
			
x=[-0.5;1];

k=10

% plot([0 x(1)],[0 x(2)],'g'); hold on; pause(1)

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

  [t, r] = max(abs(y)); 	% r lagrer komponenten som gir maks.

  mu=y(r), x=y/mu    		% skalerer med mu 
  
  error=  max(abs(mu*x-A*x))  	% angir et maal for feilen
  
end