def cumint(fn=kroupa, bins=np.logspace(-2,2,500)): xax,integral = integrate(fn,bins) return integral.cumsum() / integral.sum()