summarize_uncertainty function

Get Monte Carlo confidence intervals