Numeric Trap Volume Integration

Added the functions for numeric trap volume integration. This works for
a lattice of arbitrary dimension, but still need to integrate with the
3Body wave procedure.

DSB