function j = simula(P,N,i)
%
% j = simula(P,N)
% j = simula(P,N,i)
%
% Simula la catena di Markov di matrice di transizione P, 
% per N transizioni
n = length(P);
if (nargin == 2)
  i = randdisc(ones(1,n)/n);
end
for k = 1:N
  j = randdisc(P(i,:));
  i = j;
end
