A solution algorithm is developed for frictional contact problems solved using the boundary integral equation method when the loading is proportional. The algorithm takes advantage of certain characteristics of the boundary integral equation formulation of contact problems and of proportional loading. These characteristics lead to a solution process requiring that the incremental-iterative procedure commonly used in solving contact problems numerically be applied only to a reduced set of equations representing points within the potential contact zone.