randn()

I would suggest to add a function randn() which generates pseudo-random numbers according to the normal distribution.

In several cases it happened to me that I needed to generate a randomly perturbed initial condition to test the sensitivity of the solution by performing then extended multi-precision computations.

Thanks in advance!