Univariate modified Fourier methods for second order boundary value problems (Q2391027)

The author develops and analyses a new spectral-Galerkin method for the numerical solution of linear, second order differential equations with homogeneous Neumann boundary conditions. The basis functions for this method are the eigenfunctions of the Laplace operator subject to these boundary conditions. Due to this property this method has a number of beneficial features, including an \(\mathcal{O}(N^{2})\) condition number and the availability of an optimal, diagonal preconditioner. This method offers a uniform convergence rate of \(\mathcal{O}(N^{-3})\), however the author shows that by the inclusion of additional \(2M\) basis functions, this figure can be increased to \(\mathcal{O}(N^{-2M-3})\) for any positive integer \(M\).