* Nom du fichier: gauss.kumac
*
* KUMAC pour montrer la distribution normale et sa cumulative
*                                      MP 1/95
*                             Histogrammes pour la presentation
his/cre/1d 500 'Distribution normale et cumulative' 500 -5. 5.
set/max 500 1.
set/min 500 0.
his/cre/1d 502 'FREQ, ERF et ERFC' 500 -5. 5.
set/max 502  2.
set/min 502 -1.
*                             Cumulative de la Gaussienne
sigma x=array(500,-5.#5.)
sigma f=freq(x)
*                             Distribution normale
sigma g=(1./sqrt(2.*3.1415927))*exp(-x**2/2.)
*                             Presentation
for/file 20 'gauss1.eps'
meta     20 -113
his/put/cont 500 f
his/plot 500
his/put/cont 500 g
his/plot 500 's'
close 20
wait
*                             Comparaison FREQ, ERF et ERFC
for/file 20 'gauss2.eps'
meta     20 -113
his/put/cont 502 f
his/plot 502
sigma f=erf(x)
his/put/cont 502 f
his/plot 502 's'
sigma f=erfc(x)
his/put/cont 502 f
his/plot 502 's'
for/close 20