started working on getting the integration bounds.  All that's really left to do for this first part is figure out how to set up a call to a net exp(-
beta
phi) that is amenable to the multi-dimensional integration function I wrote