GaussianJointDistCommands.m

% generate the data
D = [randn(1000,2)+1.5; randn(1000,2)-1.5];
%% There are lots of ways to visualize it
[H,C] = hist3(D,[20 20]); 
imagesc(D);
hist3(D,[20 20]);
imagesc(H);
[I,J]=meshgrid(C{1},C{2}); 
plot3(D(:,1),D(:,2),zeros(2000,1)-0.002,'r*');
hold on
surf(I,J,H/2000);
alpha(0.2);