DifferentialEvolution.py 文件源码

def NPV(rho,g,eff,H,q_process,delta_t,p,Q,a,b):
  integrand = lambda t: np.exp(-r*t)*power(rho,g,eff,H,q_process)*delta_t*p
  PV = quad(integrand, 0, np.inf)
  NPVal = PV[0] - cost(a, b, Q)
  return NPVal

# check the function, it does not look right!!!