different result compared to numpy

Hello there,

I am trying to use this repo to replace numpy but get different result.

I put range as the minimum of the input and the maximum of the input. But I found out that the result is missing some maximum value.

For example,

`test_case = np.array([1, 1, 2, 2, 3, 3, 10, 10])
freq, bins = np.histogram(test_case, range(np.min(test_case), np.max(test_case + 1)))
result = histogram1d(test_case, bins=9, range=(np.min(test_case), np.max(test_case)))`

Start by reproducing the reported example with the shown test_case and compare histogram1d against numpy.histogram using the stated range and nine bins. Trace how the maximum boundary is handled and confirm the corrected result includes the expected maximum values.

Written by the indexing model from the issue text.