Spectral integration and two-point boundary value problems

A numerical method for two-point boundary value problems with constant coefficients is developed which is based on integral equations and the spectral integration matrix for Chebyshev nodes. The method is stable, achieves superalgebraic convergence, and requires O(N log N) operations, where N is the number of nodes in the discretization. Although stable spectral methods have been constructed in the past, they have generally been based on reformulating the recurrence relations obtained through spectral differentiation in an attempt to avoid the ill-conditioning introduced by that process.