Sampling from a Gaussian distribution

What subroutine should I use to sample from a Gaussian distribution?

You can use one of the following:

 CALL FLNRRN (RGAUSS)

returns one number normally distributed

 CALL FLNRR2 (RGAUS1, RGAUS2)

returns two such numbers, independent of each other