The new Legendre–Gauss collocation method is easy to be implemented for
nonlinear problems, and possesses the spectral accuracy. Thereby, the smoother
the exact solution, themore accurate the numerical results. However, the existing
numerical methods for ordinary differential equations do not possess such
property.