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)