data = randn(1,7500); %// example data
n = 21; %// number of bins
[y, x] = hist(data, 21);
y = y/(x(2)-x(1))/numel(data); %// normalize so that total area is 1
plot(x, y)