4 C The next three subroutines are for Monte Carlo generation
5 C according to a given function FHB. One calls first HIFUN
6 C with assigned channel number I, low and up limits. Then to
7 C generate the distribution one can call HIRND(I) which gives
8 C you a random number generated according to the given function.
10 SUBROUTINE HIFUN(I,XMIN,XMAX,FHB)
13 FNORM=GAUSS1(FHB,XMIN,XMAX,0.001)
15 XX(I,J)=XMIN+(XMAX-XMIN)*(J-1)/200.0
17 RR(I,J)=GAUSS1(FHB,XMIN,XDD,0.001)/FNORM