function power_iteration(a, iterations, v = ones(size(a,1))) l = 1 for i = 1:iterations v = a * v j = findmax(abs.(v))[2] l = v[j] v /= l end (l,v) end