eta_star function

Calculates parameter $\eta^*$ of $q^*(\sigma^2_{\gamma})$ to optimize the evidence based lower bound (ELBO) in survregVB.frailty.fit.