limits = np.arange(mn, mx + d / 2, d)