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

eig(A)

%% 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 ; 1]
% inverz hatvanymodszer
for i = 1:10
  z = inv(A) * w
  mu = 1 / norm(z,Inf)
  Mu(i) = mu
  w = z * mu
  W(:,i) = w
endfor

[sv, se] = eig(A)