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

*              Nombre d'occurance des 49 chiffres dans les 12 derniers tirages 
m  = 12
his/cre/1d 10 '                       ' 13 -0.5 12.5
his/cre/1d 11 'Distribution binominale' 13 -0.5 12.5
opt linx
set hcol 0
set/norm 10 1.00
*              Calculer la distribution binominale
sigma kf   = array(13,1&1&2&6&24&120&720&5040&40320&362880&3628800&39916800&479001600)
sigma nf   = array(13)*kf(13)
sigma nmkf = order(kf,-kf)
sigma k    = array(13,0#12)
sigma nmk  = array(13,12#0)
sigma bino = (0.143)**k * (1.-0.143)**nmk * nf/kf/nmkf
*              12  essais pour les 49 chiffres (source: Nice Matin) 
vec/cre nima(49) r 2 2 1 2 2 1 2 3 2 1 1 1 2 0 2 1 3 3 1 3 3 2 2 0 2 3 0 3 1 1 0 2 1 2 3 0 1 2 2 1 1 3 3 2 0 3 3 0 3 
zone 1 1
vec/hfi nima(1:49) 10
hi/plo 10
his/put/con 11 bino
hi/plo 11 'cs'
atitle 'No succes' 'No occurence' 
for/close 20