function [lambda,v]=invpotens(A,x0,mu)

% Finner egenverdi og tilhørende egenvektor som er nærmest innparameter mu

v = x0;

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

As = A-mu*eye(length(x0));

[L,U]=lu(As);

y = U\(L\v);
[tmp,ind]=max(abs(y));

for i=1:20
  
  y = U\(L\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;
  
end
lambda = 1/lambda+mu;

