We investigate some numerical methods for computing approximate solutions of a system of second order boundary value problems associated with obstacle, unilateral and contact problems. We show that cubic spline method gives approximations which are better than that computed by higher order spline and finite difference techniques.