Boundary value problems for nonlinear integro-differential equations of the neutral type are investigated. A scheme for approximating the boundary value problem solution using cubic splines of defect two is proposed and substantiated. A model example illustrating the proposed approximation scheme is considered.