Fix bugs and gaps in multivariate distributions

* rand() for MultivariateDistribution's always emits columns
* Return full covariance matrix for Dirichlet and Multinomial
* Add cov() as alias for var() for MultivariateDistribution