Numerical method for nonlinear differential equation with arbitrary BC