function [bw, numbins] = binwidth(data);
range = max(data) - min (data);
samples = size (data);
bw = range/(1+log(samples(1)));
numbins = range/bw;