Examples Internal Tools III (susgen)
c
c=== This routine generates a gaussian centered at 1000 with a sigma of
100
c=== This routine is called ONLY if you turn SUSGEN on
c=== User you can modify this routine to simulate your data if you wish!
real*4 h,v,t
10 call random(h)
h = 2000. * h
v=exp(-0.5*((h-1000.)/100.)**2)
call random(t)
if (t .gt. v) goto 10
h = h
Return
end