I am writing a fortran 95 code for Brownian dynamics. So I need to insert random numbers at each time-step. To generate the random numbers or ‘noise’, which needs to be Gaussian white noise, I am calling a subroutine which employs the Box-Muller algorithm as follows –