A previously published integration algorithm applicable to the numerical computation of integrals with rapidly oscillating integrands is generalized. The previous algorithm involved quadratic approximation of the phase function which was assumed to be real. The present generalization concerns approximation of a complex phase function by a polynomial of arbitrary degree. As before, the integrand is then written without approximation as a slowly varying function multiplied by the polynomial phase exponential and the slowly varying factor is approximated by a finite sum of Chebyshev polynomials. The integral is thus expressed as a sum of constituent integrals which are computed recursively via LU decomposition applied to a system of linear equations with a banded coefficients matrix. Examples are presented comparing various degree phase approximants.
Read full abstract