def METHOD_NAME(N=10000, rseed=0):
    rng = np.random.default_rng(rseed)
    X = rng.standard_normal(N)
    dx, bins = knuth_bin_width(X, return_bins=True)
    assert_allclose(len(bins), 58)
    dx2 = knuth_bin_width(X)
    assert dx == dx2
    with pytest.raises(ValueError):
        knuth_bin_width(rng.random((2, 10)))