Created package newprodforecast. Got Stewart's code for Bayesian version.