function [lambda, v] = power_method ( A, x, TOL, Nmax )

for i=1:Nmax
    v=x/norm(x);
    x=A*v;
    lambda=x'*v
    if norm(x-lambda*v)< TOL break;
    end
end

lambda
v