Univariate modified Fourier methods for second order boundary value problems

We develop and analyse 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 condition number and the availability of an optimal, diagonal preconditioner. This method offers a uniform convergence rate of , however we show that by the inclusion of an additional 2M basis functions, this figure can be increased to for any positive integer M.