function [sajate,sajatv]=invhatvany(A,v0,e)
  lnull=1/norm(v0,inf)
  v=inv(A)*v0
  l=1/norm(v,inf)
  v=v/norm(v,inf)
  while abs(l-lnull)>e
    v0=v
    lnull=l
    v=inv(A)*v0
    l=1/norm(v,inf)
    v=v/norm(v,inf)
  endwhile
  sajate=l*sign(v(1)/v0(1))
  sajatv=v
endfunction