function [lambda,v]=potens(A,x0)
%
% [lambda,v]=potens(A,x0)
%
% Finner største egenverdi lambda og tilhørende egenvektor v til A, 
% ved hjelp av potensmetoden. 
%

r1=0; r2=0; % For Aitken akselrasjon (se s. 363)
v = x0; 

ind = 1;    % Velger phi(x)=x(ind);

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