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