def EV_flow(mean,std,Q):
EV = quad(lambda x: lognorm_pdf(x,mean,std)*x, 0, Q)
EV_2 = (1 - lognorm_cdf(Q,mean,std))*Q
EV_total = EV[0] + EV_2
return EV_total
微信
公众号
扫码关注公众号