*
* KUMAC pour fitter une fonction periodique
*                                      MP 4/95
*                                Creation echantillon
sigma n=25
sigma alpha=2.
sigma  t=array(n,0.#2.*pi)
sigma  a=alpha*sin(t)
sigma da=array(n)*0.1*alpha
sigma r1=rndm(array(n))
sigma r2=rndm(array(n))
sigma  a=sin(2.*pi*r1)*sqrt(-2.*log(r2))*da + a
his/cre/1d  100 'Echantillon fonction periodique' 25 0. 6.2832
his/put/con 100 a
his/put/err 100 da
*                                Fitter la fonction
opt nstat
opt fit
set fit 111111
vec/cre par(2) r 2.1 0.999
vec/cre stp(2) r 0.1 0.1
vec/cre erp(2) r 0   0
appl hminuit fin
migrad
scan 2 100 0. 5.
exit
fin
his/fit 100 'sinus.f' 'm' 2 par