bayesNormalPrediction function

Obtain a Bayesian prediction interval