Currently, we only have 1D integration. For completeness we should have 1D-6D. Note: this project may be useful: https://github.com/jacobwilliams/quadrature-fortran See also #33