% Power iteration

function [v,lambda] = power_itr(A,v,N)

for k=1:N
    
    w = A*v;
    v = w/norm(w);
    lambda = v'*A*v;
    
end


