function [e_val,e_vec]=power_el(A,v0,e)
  v0=v0/norm(v0,inf)
  v=A*v0
  l0=norm(v0,inf)
  l=norm(v,inf)
  v=v/norm(v,inf)
  while abs(l-l0)>e
    v0=v
    l0=l
    v=A*v0
    l=norm(v,inf)
    v=v/norm(v,inf)
  endwhile
  e_val=l*sign(v(1)/v0(1))
  e_vec=v
endfunction