Indeed, I've used randn() wherever a normally distributed random number is needed. This is computed using the Box-Muller method.