In this paper, the variational iteration method is used to solve a system of third-order boundary value problem associated with obstacle, unilateral and contact problems. Numerical solution obtained by the method is of high accuracy. The numerical example compared with those considered by other authors shows that the method is more efficient. [ABSTRACT FROM AUTHOR]