The theoretical analysis and computational implementation of factorization-based methods for the numerical solution of linear boundary value problems for ordinary differential equations are presented.

Polynomial interpolation to analytic functions can be very accurate, depending on the distribution of the interpolation nodes. However, in equispaced nodes and the like, besides being badly ...