Second, fourth and sixth order methods are developed and analysed for the numerical solution of non-linear second order boundary value problems.
The methods arise from a two-step recurrence relation involving exponential terms, these being replaced by Padé approximants .
The methods are tested on two problems from the literature