 clf;   % On efface les figures
 clear; % On nettoie la memoire
 n=input('Nombre total de boules (ex. 10) = ? '); 
 k=input('Nombre initial de boules dans la premiere urne (ex. 2) = ? '); 
 m=input('Nombre de termes dans le thm ergodique (ex. 300) = ? ');
 X(1)=k; % X(i)= nombre de particules dans la premiere urne a l'étape i
 for i=1:m-1
  % Notez que la Bernoulli de parametre p est une binomiale B(1,p)
  % Ici, on simule une bernoulli de parametre p=1-X(i)/n sur -1,+1 a partir 
  % d'une bernoulli de meme parametre sur 0,1.
  u=2*rbinom(1,1,1-(X(i)/n))-1;
  X(i+1)=X(i)+u;
 end
 histo(X,n+1,1);
 hold on;
 Y=dbinom(0:n,n,1/2);
 plot([0:n],Y*m,'r');
 title('Mesures invariantes empirique et theorique');
 legend('Mesure invariante simulee','Loi binomiale')
