This routine is a utility routine that does not check the
input arguments for proper dimensions and type.
References
For further discussion, see Bayesian Statistics and Marketing
by Rossi, Allenby and McCulloch, Chapter 2.
http://faculty.chicagogsb.edu/peter.rossi/research/bsm.html