runmcmc_cpall function

Estimate posterior distributions for the 0, 1, or 2 changepoint case.