The Marsaglia polar method[1] is a pseudo-random number sampling method for generating a pair of independent standard normal random variables.[2] Standard normal random variables are frequently used in computer science, computational statistics, and in particular, in applications of the Monte Carlo method. The polar method works by choosing random points (x, y) in the square −1 < x < 1, −1 < y < 1