function rnorm = rMNorm(mu,cov,T)
% rMultiNorm generates a pxT array  of normal draws
%  from the p-dim N(mu,cov)
%
rnorm = repmat(mu,1,T) +chol(cov)'*randn(length(mu),T);