Operational Quadrature Methods for Wiener-Hopf Integral Equations

Abstract
We study the numerical solution of Wiener-Hopf integral equations by a class of quadrature methods which lead to discrete Wiener-Hopf equations, with quadrature weights constructed from the Fourier transform of the kernel (or rather, from the Laplace transforms of the kernel halves). As the analytical theory of Wiener-Hopf equations is likewise based on the Fourier transform of the kernel, this approach enables us to obtain results on solvability and stability and error estimates for the discretization. The discrete Wiener-Hopf equations are solved by using an approximate Wiener-Hopf factorization obtained with FFT. Numerical experiments with the Milne equation of radiative transfer are included.

This publication has 20 references indexed in Scilit: