A new recursive algorithm based on matrix form for Kane equations for calculation of closed form dynamic equations for simple kinematical chain is configured. It has computational efficiency 0(n3) and is used for motion simulation for space manipulators with rotation joints. This algorithm allows to consider links flexibility and usage of parallel calculations.