Does np.divide not work here? And why not use `np.cumsum(pmf * binnums, dtype=float) / cpmfl`?