N=evstr(x_dialog('Nombre de boules dans A ?','4'));
T =evstr(x_dialog('Nombre de tirages ?','25'));
X=grand(2,N,'unf',0,1);//positions initiales dans l'urne A
X(1,:)=X(1,:)-1
X(2,:)=2*(X(2,:)-0.5)
Y=zeros(1,T)
clf
f=gcf();f.pixmap='on';
f.figure_size= [610,469]*1.4;
f.figure_position =[10,10];
xset('window',0)
xsetech([0,0,1,0.5],[0,1,0,1])
  xset("color",4)
  xfrect(0,1,1,1)
   xset("color",32)
  xfrect(1,1,1,1)
  plot2d(0,0)
plot2d(X(1,: ),X(2,: ),style=-4)
for i=1:T// les tirages
  c=grand(1,1,'uin',1,N);
  if (X(1,c)<0) then
    X(1,c)=X(1,c)+1;
  else
    X(1,c)=X(1,c)-1;
  end
  Y(1,i)=sum(X(1,:)<0);
  clf
    xsetech([0,0,1,0.5],[0,1,0,1])
   xset("color",4)
  xfrect(-1,1,1,2)
   xset("color",32)
  xfrect(0,1,1,2)
   plot2d(X(1,: ),X(2,: ),-4,"012",axesflag=0)
   xsetech([0,0.5,1,0.5])
 plot2d([0,1],[N,N-1],5)
 plot2d([1:T],Y(1,: ),5)
 xtitle("Evolution du nombre de boules dans A");
  show_pixmap()//on peut accelerer considérablement en n'affichant pas chaque tirage
end


