The Ludwig integration algorithm for triangular subregions

The well-known Ludwig algorithm for numerical treatment of integrals with oscillatory integrands using rectangular subregions is modified so as to make use of triangular subregions. Algebraic expressions are given suitable for programming into a numerical integration subroutine. All of the advantages of the Ludwig scheme are retained.