A Numerical Approach to Nonlinear Two-Point Boundary Value Problems for ODEs