credible_plot function

Plot credible intervals for a, R, and b