Sigma_star function

Calculates parameter $\Sigma^*$ of $q^*(\beta)$ to optimize the evidence based lower bound (ELBO) in survregVB.fit.