Demonstration of Inverse Transform Sampling.
Relies on a numerical approximation to the probit function found here. Values are sampled uniformly from the interval [0, 1] and mapped through the probit function. If X is a random variable in Uniform(0, 1), then Y=probit(X) is a normally distributed random variable.