* Nom du fichier: gauss.kumac
*

*                             Histogrammes pour la presentation
his/cre/1d 500 'Distribution N=10' 500 -5. 5.
his/cre/1d 501 'Distribution N=20' 500 -5. 5.
his/cre/1d 502 'Distribution N=100' 1000 -5. 15.
set/max 500 1.
set/min 500 0.
set/max 501 1.
set/min 501 0.
set/max 502 1.
set/min 502 0.
sigma x1=array(500,-5.#5.)
sigma x2=array(1000,-5.#15.)


*                             Distribution normale
sigma mu=0.05*10
sigma si=sqrt(mu)
sigma g1=(1./sqrt(2.*3.1415927)/si)*exp(-(x1-mu)**2/2./(si**2))
sigma mu=0.05*20
sigma si=sqrt(mu)
sigma g2=(1./sqrt(2.*3.1415927)/si)*exp(-(x1-mu)**2/2./(si**2))
sigma mu=0.05*100
sigma si=sqrt(mu)
sigma g3=(1./sqrt(2.*3.1415927)/si)*exp(-(x2-mu)**2/2./(si**2))


*                             Distribution binomiale
sigma pp=0.05
sigma ratio=gamma(11)/gamma(x1+1)/gamma(11-x1)
sigma b1=ratio*(pp**x1)*((1-pp)**(10-x1))
sigma ratio=gamma(21)/gamma(x1+1)/gamma(21-x1)
sigma b2=ratio*(pp**x1)*((1-pp)**(20-x1))
sigma ratio=100**x2/gamma(x2+1)
sigma b3=ratio*(pp**x2)*((1-pp)**(100-x2))


*                             Distribution poisson
sigma mu=0.05*10
sigma p1=exp(-mu)*mu**x1/gamma(x1+1)
sigma mu=0.05*20
sigma p2=exp(-mu)*mu**x1/gamma(x1+1)
sigma mu=0.05*100
sigma p3=exp(-mu)*mu**x2/gamma(x2+1)



*                             dessin
*for/file 20 'gauss1.eps'
*meta     20 -113
his/put/cont 500 g1
his/plot 500 
his/put/cont 500 p1
his/plot 500 s
his/put/cont 500 b1
his/plot 500 s
wait
his/put/cont 501 g2
his/plot 501 
his/put/cont 501 p2
his/plot 501 s
his/put/cont 501 b2
his/plot 501 s
wait
his/put/cont 502 g3
his/plot 502 
his/put/cont 502 p3
his/plot 502 s
his/put/cont 502 b3
his/plot 502 s
*close 20