if num_bins is None:
            num_bins = int(np.ceil(4 * np.sqrt(np.sqrt(len(x)))))