Frictionless contact problems are the simplest and classical contact problems, and the contact conditions of sticking, slipping, and separation mode all can be ascribed to complementary problems. Consequently, a smooth system of equations approach for the design and analysis of complementarity problems for frictionless contacts is presented. A compute program based on boundary element technique is given and applied to two practical contact examples. The validity and accuracy of the proposed method are demonstrated.