def _F5(self, z):
"""
Calculate integral of function f5 from Hellstrom (1991)
"""
B = self._delta*self._beta2 + self._beta1*self._beta12/self._gamma
C = self._beta2*self._beta - B*self._gamma
S = - (self._beta2*self._gamma - self._beta*B)
denom = (self._beta**2 - self._gamma**2)
F5 = np.exp(self._beta*z) / denom \
* (C*np.cosh(self._gamma*z) + S*np.sinh(self._gamma*z))
return F5
评论列表
文章目录