compute_mu_s function

Compute posterior mean and precision of s