n=1000;
gamma=0.75;
mu=[10 12 8 22 4];
mu_mat=ones([n 1])*mu;
scale=10;
C=gamma*scale*randn(5)^10+(1-gamma)*eye(5);
ndist=[randn([1 n]); randn([1 n]); randn([1 n]); randn([1 n]); randn([1 n])]';
gdist=ndist*C;
offdist=gdist+mu_mat;


scatter(offdist(:,1),offdist(:,2))


				% 
				% C =
				% 
				%    1.0e+03 *
				% 
				%     0.1476    0.4185    0.1509   -0.0997    0.0641
				%     0.5691    0.4719    0.3385   -0.0129    0.4526
				%    -0.3655    0.0890   -0.4012   -0.2623    0.0393
				%     1.2967   -0.3477    0.7049    0.5604    0.9309
				%    -1.6184   -0.9829   -1.2949   -0.2985   -0.7401
				% 

