The cubic spline approximation to a two-point boundary value problem for the differential equation y″ + f(x)y′ + g(x)y = r(x) is shown to reduce to the solution of a three-term recurrence relationship. For the special case when f(x) is a constant, the approximation is shown to be simply related to a finite-difference representation and to have a local truncation error of order