import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

xedges = [0, 1, 2, 3]
yedges = [0, 1, 2, 3, 4]
x = np.array([0, 0.1, 0.2, 1., 1.1, 2., 2.1])
y = np.array([0, 0.1, 0.2, 1., 1.1, 2., 3.3])

out = np.histogram2d(x, y, bins=(xedges, yedges))

plt.scatter(x, y)
plt.grid()