compute_posterior_intervals function

Compute Posterior Intervals