% Eksempel, inverspotens metoden

clear all;   

A=[0.5 1.5; 1  1]  	

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

k=7			      % k= antall iterasjoner  

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

 a=-0.2;          % "tipper" at denne verdien er naer -0.5


 C=A-a*eye(2);
 

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

