Abstract
The Goodman-Lance method of adjoints for solving two point boundary value problems of systems of non-linear ordinary differential equations is extended to handle the most general boundary conditions and specialized to several important subcases.