function m=histogram_bins(y);
m=floor(2*iqr(y))/(length(y)^(1/3));
end

histogram_bins(Z)

ans =

    1.5791

histogram_bins(B)

ans =

   17.4452

 histogram(B,100,18) % dont pay attention to 100