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 = 1;    % Velger phi(x)=x(ind);

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

for i=1:100
  
  y = A*v;
  r = y(ind)/v(ind);
  v = y/y(ind);
  
  % Aitken
  if i < 3
    lambda = r
  elseif i>2
    lambda = r-(r-r1)^2/(r-2*r1+r2)
  end
  r2=r1;
  r1=r;
  v
  pause
  
end


