function [counts, bins] = histg(x, chs) if nargin==1, chs = 'stem'; end x = x(:); bins = 0:255; [counts, bins] = hist(x, bins); switch chs case 'stem' h = stem(bins, counts, 'k'); set(h,'Marker','none') case 'bar' bar(bins, counts, 1, 'k'); case 'stairs' stairs([bins-0.5 bins(end)+0.5], [counts counts(end)], 'k') end y = round(max(counts)/40); axis([bins(1) bins(end) min(counts) max(counts)+ y])