def METHOD_NAME():
    rng = np.random.default_rng(0)
    X = rng.standard_normal(100)
    counts, bins = histogram(X, bins="knuth")
    assert_allclose(counts, [2, 1, 13, 19, 15, 18, 14, 10, 8])
    # fmt: off
    assert_allclose(bins, [-2.32503077, -1.84420596, -1.36338114, -0.88255632, -0.4017315,
                           0.07909331, 0.55991813, 1.04074295, 1.52156777, 2.00239258])
    # fmt: on