* Nom du fichier: bino.kumac
*
* KUMAC pour demontrer la distribution binomiale
*                                      MP 03/01
*              Sauvegarde du graphique sur fichier ps
for/file 20 'bino.eps'
meta     20 -113

*              Compter le nombre de piles dans [m] essais
m  = 10
vec/cre np(1) r 0.
his/cre/1d 10 '                       ' 10 0.5 10.5
his/cre/1d 11 'Distribution binominale' 10 0.5 10.5
opt linx
set hcol 0
set/norm 10 1.00
*              Calculer la distribution binominale
sigma rf   = array([m],1#[m])
sigma rf   = prod(rf)
sigma nf   = array([m])*rf([m])*2**(-[m])
sigma nmrf = prod(array(10,1&1&2&3&4&5&6&7&8&9))
sigma nmrf = order(nmrf,-nmrf)
sigma bino = nf/rf/nmrf
*              Repeter les [m] essais [ne] fois
ne = 1000
zone 2 2
do i=1,[ne]
   sigma pm = array([m])
   sigma pm = rndm(pm)
   sigma np = vsum(pm gt 0.5)
   call hfill(10,np,0.,1.)
*              Graphique apres 1, 10, 100 et [ne] essais
   if ([i]=1).or.([i]=10).or.([i]=100).or.([i]=[ne])  then
      hi/plo 10
      his/put/con 11 bino
      hi/plo 11 'cs'
      atitle 'No piles' 'No essais'
      if ([i]=1)  then
       text 8. .8 'n=1' .3 
      elseif ([i]=10) then
       text 8. .3 'n=10' .3 
      elseif ([i]=100) then
       text 8. .2 'n=100' .3 
      else
       text 8. .2 'n=1000' .3 
      endif
   endif
enddo
for/close 20                                          
hi/del 0
ve/del *