def _f2(self, z): """ Calculate function f2 from Hellstrom (1991) """ f2 = np.exp(self._beta*z)*self._beta12/self._gamma \ * np.sinh(self._gamma*z) return f2