function [x,iter,vres]=potenze(P,eps,maxiter);
% Risolve (I-P')x=0 con il metodo delle potenze: M=I, N=P'
% Input: P stocastica, eps per la condizione di arresto, maxiter=numero max
% di iterazione
% Output: approssimazione di x, iterazioni iter effettuate, vettore dei
% residui
n=size(P,1);
% x=ones(n,1)/n;
x = rand(n,1); x=x/sum(x);
iter=0;
res=norm(P'*x-x,1);
vres=res;
while (res>eps & iter<maxiter)
iter=iter+1;
y=P'*x;
res=norm(y-x,1);
y=y/sum(y);
vres=[vres,res];
x=y;
end
