Scripted HTML version of Foils prepared 12 November 1996

Foil 3 Generation of Gaussian Distributions

From CPS615-End of Basic Overview of Random Numbers and First Part of Monte Carlo Integration Delivered Lectures of CPS615 Basic Simulation Track for Computational Science -- 22 October 96. by Geoffrey C. Fox *
Secs 83.5
"A Small Miracle" asserts that:
If x1 x2 are uniformly distributed in [0,1] -- Then:
  • g1 = (-2lnx1)1/2 cos 2px2
  • g2 = (-2lnx1)1/2 sin 2px2
are Gaussianly distributed
with mean = 0 and standard deviation = 1
while g1 and g2 are independent.
Proof: Consider
Integral:
with g1 g2 going to Polar coordinates (r=radius, angle)
and then transform to x1and x2 by
(-2lnx1)1/2 = radius i.e. x1=exp(-r2/2) and
2px2 = angle



© Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Fri Aug 15 1997