def _f3(self, z): """ Calculate function f3 from Hellstrom (1991) """ f3 = np.exp(self._beta*z)*(np.cosh(self._gamma*z) + self._delta*np.sinh(self._gamma*z)) return f3