The computer implementation of the recursive projection methods in flexible multibody dynamics is described. The organization of the computer program which is used to automatically construct and numerically solve the system of loosely coupled dynamic equations expressed in terms of the absolute and joint coordinates is discussed. The main processor computer program consists of three main modules: constraint module, mass module and force module. The constraint module is used for the numerical evaluation of the relationship between the absolute and joint accelerations. The mass module is used for the numerical evaluation of the system inertia matrix as well as the nonlinear Coriolis and centrifugal forces associated with the absolute and elastic coordinates. The force module is used to evaluate the generalized external and elastic forces associated with the absolute, joint and elastic coordinates.