^d, where <b>d</b> is the spatial dimension.  The function includes an optional
      argument <b>bpt</b>, which is an array allowing the user to specify the lower and upper 
      endpoints of integration.  However, in the experience of the installer, this extra argument
      does not seem to have the correct effect, resulting in quadrature rules that are not
      correctly shifted to the new domain.  A note has been sent to the author asking for
      clarification.  In the meantime, the user can comput