% Eksempel, inverspotens metoden

clear all;   

A=[8 0 12; 1 -2 1; 0 3 0]  	

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

k=7			      % k= antall iterasjoner  

[V D]=eig(A)      

 a=-1.4;          


 C=A-a*eye(3);
 

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

