We consider a finite difference method of order four for nonlinear two-point boundary value problems. In linear case the finite difference schemes lead to a tridiagonal linear system. Numerical experiments support the theoretical results.