The dynamic equations of motion of the constrained multibody mechanical system are mixed differential-algebraic equations (DAEs). The numerical solution of the DAE systems solved using ordinary-differential equation (ODE) solvers may suffer from constraint drift phenomenon. To solve this problem, Baumgarte proposed a constraint stabilization method in which a position and velocity terms were added in the second derivative of the constraint equation. Baumgarte’s method is a proportional-derivative (PD) type controller design. In this paper, an controller is included to form a proportional-integral-derivative (PID) controller so that the steady state error of the numerical integration can be reduced. Stability analysis methods in the digital control theory are used to find out the correct choice of the coefficients for the PID controller.