Since this is the first Google result for “js gaussian random” in my experience, I feel an obligation to give an actual answer to that query.
The Box-Muller transform converts two independent uniform variates on (0, 1) into two standard Gaussian variates (mean 0, variance 1). This probably isn’t very performant because of the sqrt
, log
, and cos
calls, but this method is superior to the central limit theorem approaches (summing N uniform variates) because it doesn’t restrict the output to the bounded range (-N/2, N/2). It’s also really simple:
// Standard Normal variate using Box-Muller transform.
function gaussianRandom(mean=0, stdev=1) {
let u = 1 - Math.random(); // Converting [0,1) to (0,1]
let v = Math.random();
let z = Math.sqrt( -2.0 * Math.log( u ) ) * Math.cos( 2.0 * Math.PI * v );
// Transform to the desired mean and standard deviation:
return z * stdev + mean;
}