A = [2 3 ; 1 4]
v = [1 ; 0]
clear('V','La')

%% hatvanymodszer
for i = 1:8
  z = A*v
  la = norm(z,Inf)
  La(i) = la
  v = z / la
  V(:,i) = v
endfor

clear('Mu')
w = [1 ; 0]
% inverz hatvanymodszer
for i = 1:5
  z = inv(A) * w
  mu = 1 / norm(z,Inf)
  Mu(i) = mu
  w = z * mu
  W(:,i) = w
endfor

[sv, se] = eig(A)