The objective of this paper is to develop a solution method for the differential-algebraic equation(DAE) derived from constrained muti-body dynamic systems. Mechanical systems are often modeled as bodies and joints. Differential equations of motion are formulated for bodies. Since the bodies are connected by joint, the differential variables must satisfy the kinematic constraint equations that come from the joints. Difficulties are arised due to drift of the differential variables off the constraint equations. An optimization method is adopted to correct the drift of the differential variables. To demonstrate the efficiency of the proposed method a slider-crank mechanism is analyzed dynamically. Identical results are obtained as these from the commercial program DADS. Dynamic analysis of a High Mobility Multi-purpose Wheeled. Vehicle(HMMWV) is carried out to show the practicalism of the proposed method.