There is no built in function in MATLAB to do this.
My suggestion is:
x = rand(1000,1);
[n,xout] = hist(x);
width = 0.8*(xout(2)-xout(1));
base = 0;
clr = jet(length(n));
hold on
for k=1:length(n)
fill([xout(k)-width/2 xout(k)+width/2 xout(k)+width/2 xout(k)-width/2],...
[base base n(k) n(k)],clr(k,:))
end
hold off
see also:
