The Box–Muller transform method (also mentioned int he article) is kind of fun. You generate two Gaussians at once. You rejection sample from the unit square into the largest contained circle (generates a point about 3/4 of the time). You then roughly log-transform the distance along the radius (sending (0,0) to infinity) and the two points are now independent Gaussians.