* Nom du fichier: tmuon.kumac
*
* KUMAC pour montrer la distribution exponetielle
* MP 3/01
*                             Generer echantillon exponentiel
ne = 100000
sigma lam=2.19703
sigma x=rndm(array([ne]))
sigma x=-lam*log(x)
his/cre/1d  100 'Distribution temps de desintegration' 100 0. 25.
vec/hfill x 100
 *              Sauvegarde du graphique sur fichier ps
for/file 20 'tmuon1.eps
meta     20 -113
opt logy
opt stat 
zone 1 1
his/plot 100
close 20
wait
*               Differentes limites inferieures
his/cre/1d  101 'Distribution t sup  2.5' 90  2.5 25.
his/cre/1d  102 'Distribution t sup  5.0' 80  5.0 25.
his/cre/1d  103 'Distribution t sup  7.5' 70  7.5 25.
his/cre/1d  104 'Distribution t sup 10.0' 60 10.0 25.
vec/hfill x 101
vec/hfill x 102
vec/hfill x 103
vec/hfill x 104
for/file 20 'tmuon2.eps
meta     20 -113
zone 2 2
his/plot 101
his/plot 102
his/plot 103
his/plot 104
close 20
 
 
 
 
 
 
 
 
 opt liny