generate N points following a a normal distribution fortran

Does anyone know how could easily generate N random numbers following a normal distribution, with a mean mu and a standard deviation sigma, in fortran90? Or even the logic process to produce the N values?