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

*              Nombre de tirages sans sortie d'un num\'ero
his/cre/1d 10 '                       ' 50 200 400
his/cre/1d 11 'Distribution gaussienne' 50 200 400
opt stat
set hcol 0
*              Calculer la distribution gaussienne
sigma k   = array(50,202#398)
sigma gau = (49*4/sqrt(6.2832)/14)*exp(-0.5*(k-322)**2/14**2)
*              echantillon des 49 numeros (source: Nice Matin 'reussite') 
vec/cre nima(49) r 338 319 298 331 321 354 342 316 312 305 298 322 322 321 315 332 305 342 346 319 316 310 323 324 318 320 323 315 285 312 327 324 290 329 323 343 324 354 316 329 307 325 334 315 335 308 306 321 336
zone 1 1
vec/hfi nima(1:49) 10
hi/plo 10
his/put/con 11 gau
hi/plo 11 'cs'
atitle 'Reussite' 'No occurence' 
for/close 20