In this paper, we use Bernstein polynomial to approximate the analytical solution of a system of third-order nonlinear boundary problems associated with obstacle, unilateral and contact problems. A new approach is illustrated to provide the solution in the form of Bernstein polynomial, which dramatically reduces the size of work. An example is given to illustrate the efficiency of the method.