alpha_star function

Calculates parameter α\alpha^* of q(b)q^*(b) to optimize the evidence based lower bound (ELBO) in survregVB.fit and survregVB.frailty.fit.