Fixing Multiple Boson Routines and Adding Bose-Fermi Routines

Still having problems getting correct answers with multiple bosons.
Seems to have to do with the way two-body observables are calculated.
Still needs to be debugged. Started adding bose-fermi functions.