Adding multivariate normal functions.
Continued implementation of heterogeneity in source strength.