% Eksempel, inverspotens metoden

clear all;   

A=[1 2 3 2; 2 12 13 11; -2 3 0 2; 4 5 7 2]  	

x=[1;0;0;0]  		  % x = startvektor

k=7			      % k= antall iterasjoner  

[V D]=eig(A)      

 a=0;          


 C=A-a*eye(4);
 

for j=1:k
  
  y=C\x; 
  
  [t , r] = max(abs(y)); 
  
  nu=a + 1/y(r), x=y/y(r)     % nu blir tilnaermet lik egnverdien etterhvert
			                  % x blir tilnaermet den tilh. egenvektoren
  
  error =  max(abs(nu*x-A*x)) % angir maal for  feilen

 end

