An adaptive quadrature numerical integration routine based upon a Newton-Cotes eight panel formula. Sub-interval sizes automatically adjust to meet accuracy requirements, thereby greatly enhancing the routines efficiency. User supplies integrand function as subroutine. Returns indicator of solution reliability. 16K.
