def integrate(fn=kroupa, bins=np.logspace(-2,2,500)):
xax = (bins[:-1]+bins[1:])/2.
integral = (bins[1:]-bins[:-1]) * (fn(bins[:-1])+fn(bins[1:])) / 2.
return xax,integral
微信
公众号
扫码关注公众号