2dhist.py

import numpy as np
import matplotlib.pyplot as plt

nbins = 25

delr=[0, 10]; delv = [0,10]
# np.meshgrid might be better here instead of this x,y stuff

rngex = np.linspace(min(delr),max(delr),nbins+1)
rngey = np.linspace(min(delv),max(delv),nbins+1)

x, y = np.meshgrid(rngex[:-1],rngey[:-1])
x = x.flatten()
y = y.flatten()
print(f"{x}\n{y}")

plt.hist2d(x,y,bins=[rngex,rngey])
plt.show()