In this paper, integration methods for stiff constrained mechanical systems are studied to analyze dynamic response of the stiff mechanical system efficiently. The stiff, non-stiff systems are identified by using eigenvalues of the Jacobian of the systems. To integrate both stiff system and non-stiff system efficiently a new switching method between the non-stiff differential equation solver and the stiff differential equation solver is presented.