Generate random number between 0 and 1 with gaussian distributes
I wrote a blog post on how to generate random numbers with any given distribution:
http://ericlippert.com/2012/02/21/generating-random-non-uniform-data/
Summing up, the algorithm you want is:
- Work out the desired probability distribution function such that the area under a portion of the curve is equal to the probability of a value being randomly generated in that range.
- Integrate the probability distribution to determine the cumulative distribution.
- Invert the cumulative distribution to get the quantile function.
- Transform your uniformly-distributed-over-(0,1) random data by running it through the quantile function.
Of course if you already know the quantile function for your desired distribution then you don't need to do steps one through three.
You say you want a generator for normally distributed (gaussian) random numbers between 0 and 1.
First of all the normal distribution is not bounded...the function you show in your example generates normally distributed random numbers with a 0.0 mean and a standard deviation of 1.0
You can generate normally distributed random values of any mean and standard deviation by multiplying the value you get from this function times the desired standard deviation and then adding the desired mean...
The code is OK as is - the problem is a misunderstanding of the gaussian (normal) distribution which has a range of -inf to +inf...
about 2/3 of the time the value you get will be between +/- 1 standard deviation....about 95% of the time the value will be between +/1 3 times the standard deviation...