A sixth order method for the integration of two point boundary value problems

A finite difference method for obtaining sixth order accurate approximation to the solution of the two-point linear boundary value problem is given. The convergence of the method is proved. Numerical results for a typical problem are tabulated and in each case the observed error is compared with its theoretical estimate. The numerical results are compared with those obtained from an earlier method of the author and the method of Noumerov.