N = 10; % desired number of samples of each variable
mu = [10; 10]; % vector of means
cov = [3 1; 1 3]; % covariance matrix
samples = mvnrnd(mu, cov, N);