function bs = binwidth(data);
range = max(data) - min (data);
samples = size (data);
bs = range/(1+log(samples(1)));