The multivariate normal distribution is easy to implement in Mathematica.
This article presents self-contained code which allows one to specify any
multidimensional normal distributions. Its use is illustrated with examples
and applications of the bivariate and trivariate normals. Random number
generation is also covered.