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