This paper introduces a new method to obtain the spectral accuracy solutions to higher order differential equations and singularly perturbed boundary value problems (BVPs). Legendre polynomials (LPs) Pn (x) have been used ...