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)