A method of dynamic analysis of mechanical systems considering probabilistic properties is proposed in this paper. Probabilistic properties that result from manufacturing tolerances can be represented by means and standard deviations (or variances). The probabilistic characteristics of dynamic responses of constrained multi-body systems are obtained by two ways : the proposed analytical approach and the Monte Carlo simulation. The formerpaper, necessitates sensitivity information to calculate the standard deviations. In this a direct differentiation method is employed to find the sensitivities of constrained multi-body systems. To verify the accuracy of the proposed method, numerical examples are solved and the results obtained by using the proposed method are compared to those obtained by Monte Carlo simulation.