Abstract
The numerical solution of second order differential equations is studied in the limit as h → 0. The form of the error is derived in a number of singular cases where the simple formulation breaks down.