In this research, a numerical algorithm has been developed, which can be applied to the large deformation and large displacement contact problems between two deformable bodies. The contact conditions expressed in terms of the rate of angle change have been proposed considering the change in geometric shape and rate of contact force. A set of linear simultaneous equations is constructed by adding the geometric shape change and contact conditions to the original stiffness matrix. A new method to determine time increment has been proposed based on Euler method, in which the condition to prevent the contact bodies from penetrating and overrunning each other has been taken into consideration. Practical application to contact problem is extrusion in which bodies are sliding along the contact boundary.