% Eksempel, potensmetoden med skalering

clear all;


A=[2 1; 4  5]   


[V D]=eig(A)    	
	
			
x=[1;0];

k=6


for j=1:k
    
  
  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