Abstract: The direct differentiation method is used for design sensitivity analysis of multibody system dynamics which is described by differential/algebraic equations with holonomic constraints. The results can be used to optimization of multibody system dynamics with generic objective functions and constraints. An example of a planar manipulator with two links is analyzed to test the given methods.