function do_hist (data, quant, max, scale)
% help make pretty plots

if (nargin < 4),
  scale = 1;
end

[n,x] = hist(data, ((quant+1)/2):quant:(max+quant));
bar(x,n./scale);
set(gca, 'XLim', [0 max+1/2]);
