* Nom du fichier: coulomb.kumac
*
* KUMAC pour simuler la diffusion multiple de Coulomb
*                                      MP 2/95
*                             Histogrammes pour la presentation
mess 'CPU Time 1' $CPTIME 
ne = 3500
si y=array([ne])
si th=array([ne])
zone 1 1
fun2 501 fcoulomb.f(x,y) 20 -1. 1. 20 -1. 1. SURF1
his/cre/2d 500 'Distribution y vs. Theta ' 20 -1. 1. 20 -1. 1.
his/cre/proy 500
his/cre/prox 500
his/cre/slix 500 10
his/cre/sliy 500 10
mess 'CPU Time 2' $CPTIME 
do i=1,[ne]
 call hrndm2(501,y([i]),th([i]),1.)
 call hfill(500,y([i]),th([i]),1.)
enddo
mess 'CPU Time 3' $CPTIME 
*                             Sauvegarde de l'echantillon
vec/write y  yotst.txt
vec/write th totst.txt
mess 'CPU Time 4' $CPTIME 'Results written to file'
*                             Sauvegarde du graphique sur fichier eps
for/file 20 'coulomb1.eps'
meta     20 -113
*              Options de presentation graphique
opt stat
*                             Presentation bidimensionelle
his/plot 500                  ***** DEMO OPTIONS *****
atit 'y' 'theta'
for/close 20                                                                   
wait
*                             Distributions marginales
for/file 20 'coulomb2.eps'
meta     20 -113
zone 1 2
his/plot 500.prox
atit 'y' 'evenements'
his/plot 500.proy             
atit 'theta' 'evenements'
for/close 20
wait       
*                             Distributions conditionelles
for/file 20 'coulomb3.eps'
meta     20 -113
zone 2 5
his/plot 500.slix.1
his/plot 500.slix.2
his/plot 500.slix.3
his/plot 500.slix.4
his/plot 500.slix.5
his/plot 500.slix.6
his/plot 500.slix.7
his/plot 500.slix.8
his/plot 500.slix.9
his/plot 500.slix.10
for/close 20
wait
*                             Distributions conditionelles
for/file 20 'coulomb4.eps'
meta     20 -113
his/plot 500.sliy.1
his/plot 500.sliy.2
his/plot 500.sliy.3
his/plot 500.sliy.4
his/plot 500.sliy.5
his/plot 500.sliy.6
his/plot 500.sliy.7
his/plot 500.sliy.8
his/plot 500.sliy.9
his/plot 500.sliy.10
for/close 20
wait