Abstract
In this paper, we develop a new cubic spline method for computing approximate solution of a system of fourth-order boundary value problems associated with obstacle, unilateral and contact problems. It is shown that the present method is of order two and gives approximations which are better than those produced by some other collocation, finite difference and spline methods. Numerical examples are presented to illustrate the applicability of the new method.