function X=gausssamp(mu,C,n)
% X=gausssamp(mu,C,n)
% Sample from Gaussian distribution.
X=randn(n,length(mu))*chol(C)+repmat(mu,[n 1]);
