clf
ax = 45 + 10*randn(1,500);
ay = 90 - ax + 5*randn(1,500);
bx = ax + 10;
by = ay + 10;
subplot(2,2,1)
scatter(ax,ay,'b.'); hold on
scatter(bx,by,'r.')
xlim([0, 100]);
ylim([0, 100]);
xlabel('dimension #1');
ylabel('dimension #2');

subplot(2,2,3)
[y x] = hist([ax; bx]',20);
bar(x,y,1)
xlim([0, 100]);
xlabel('dimension #1');

subplot(2,2,2)
[y x] = hist([ay; by]',20);
barh(x,y,1)
ylim([0, 100]);
ylabel('dimension #2');