Abstract
We develop a numerical method for computing approximations for the solutions of a system of third order boundary value problems associated with odd order obstacle problems. Such a problem arise in physical oceanography and can be studied in the framework of variational inequality theory. We study the convergence analysis of the present method and we show that it gives numerical results which are better than the other available results. Numerical example is presented to illustrate the applicability of the new method.