% Esta rutina genera el espacio de estados S la distribucion inicial mu y
% la matriz de transición P del Modelo de Ehrenfest. Para ser usado junto
% con la rutina generalchain_live.m. Ejecutar primero esta rutina y luego
% general_chain_live.

N=10; % 100 bolitas
S=0:N; % El espacio de estados es 0, ..., N.
mu=zeros(1,N); % Empezamos 
mu(1)=1;       % con todas las bolitas en la urna B. Ojo que mu(1) es la probabilidad de 0 bolitas en la urna A.  
P=zeros(N+1,N+1); % Definimos la matriz de transición
for i=2:N
    P(i,i+1)=.5*(N-(i-1))/N;
    P(i,i-1)=.5*(i-1)/N;
    P(i,i)=.5;
end
P(1,1)=.5;
P(1,2)=.5;
P(N+1,N+1)=.5;
P(N+1,N)=.5;