function [lambda,v]=potens(A,x0)

% Finner største egenverdi lambda og tilhørende egenvektor v til A, 
% ved hjelp av potensmetoden. 

v = x0;
ind = 2;    % Velger phi(x)=x(ind);
lambda_max = 5;

% For Aitken akselerasjon
r1=0; 
r2=0; 

% Plotting av to kjente egenvektorer+første gjetting x0
%vex1 = [-1;1];
%vex2 = [1;1];
%figure(1)
%quiver(0,0,vex1(1),vex1(2),'r-','linewidth',2)
%hold on
%quiver(0,0,vex2(1),vex2(2),'r-','linewidth',2)
%hold on
%quiver(0,0,x0(1),x0(2),'b--')
%hold on
%axis([-1.2 1.2 0 1.5])
%pause

for i=1:1000
  
  y = A*v;
  r = y(ind)/v(ind);
  v = y/y(ind);
  lambda = r;
  
  % Aitken
  if i < 3
    lambda = r;
  elseif i>2
    lambda = r-(r-r1)^2/(r-2*r1+r2); 
  end
  r2=r1;
  r1=r;
  
  %quiver(0,0,v(1),v(2),'b--')
  %hold on
  %pause
  
  if abs(lambda-lambda_max) < 10^(-8)
    break
  end
end
%hold off
iterations=i

