function [x,lambda]=powermethod(A,x0,k)
%calcola autovett/valore dominante
%facendo k passi del metodo delle potenze
x=x0;
% perche' non x=A^k*x0; ?
for i=1:k
   x=A*x;
   x=x*(1/norm(x));
%   lambda=(x'*A*x)/(x'*x);
%   norm(A*x-lambda*x)/norm(A*x)
end
lambda=(x'*A*x)/(x'*x);