function pv=pvns(P)
% Calcola il vettore invariante di probabilita' della matrice stocastica P
% mediante risoluzione del sistema nonsingolare pig' A= [0...0 1]
% dove A e' la matrice ottenuta sostituendo l'ultima colonna di I-P
% con [1,...,1]'
n=size(P,1);
A=eye(n)-P;
A(:,n)=1;
v=zeros(n,1); v(n)=1;
pv=A'\v;


    