Numerical methods of high-order accuracy for nonlinear boundary value Problems