* Nom du fichier: pile.kumac
*
* KUMAC pour jouer au pile et face
*                                      MP 12/94
*              Sauvegarde du graphique sur fichier ps
for/file 20 'pile.eps'
meta     20 -113
*              [n] essais de pile ou face
n  = 50000
*              Definitions du vecteur a contenir les v.a.
vec/cre rr([n]) r
*              Definitions du vecteur a contenir 1. partout
vec/cre tt([n]) r [n]*1.
*              Creer un vecteur rr de [n] pseudo-aleatoires 0 < rr < 1.
sigma !noprint
sigma rr = rndm(rr)
*              Creer un vecteur de [n] pseudo-aleatoires 0 (face) ou 1 (pile)
sigma pf = sumv(rr gt 0.5)
*              Calculer la frequence relative de pile/face jusqu'a l'element i
sigma fp = pf/sumv(tt)
sigma ff = tt-fp
*              Dessiner le contenu des vecteurs fp et ff des frequences
his/cre/1d 1 '    ' [n] 0. [n]
his/cre/1d 2 '    ' [n] 0. [n]
max 1 1.
min 1 0.
opt logx
set hcol 0
his/put/cont 1 fp
his/put/cont 2 ff
his/plot 1
set hcol 2
his/plot 2 's'
atitle 'No essais' 'Frequence'
for/close 20
return