INTRODUCTION. The robot dynamics modeling and simulation problem has been studied for the last three decades intensively. In particular, the forward dynamics problem of a robot is a very relevant issue, which there is still to say about in terms of efficient computation algorithms, that can be also simple to understand, to develop and to implement, above all for practical robots, robots with many links and/or with flexible links (Featherstone, 1987), (Featherstone and Orin, 2000), (Sciavicco and Siciliano, 2000), (Khalil and Dombre, 2002). Indeed, in these cases the methods based on the classical Lagrange formulation give rise to an analytical model with numerous terms that may be difficult to use. The methods based on the Newton-Euler...