With rapid development of methods for dynamic systems modeling, those with less computation effort are becoming increasingly attractive for different applications. This paper introduces a new form of Kane's equations expressed in the matrix notation. The proposed form can efficiently lead to equations of motion of multi-body dynamic systems particularly those exposed to large number of nonholonomic constraints. This approach can be used in a recursive manner resulting in governing equations with considerably less computational operations. In addition to classic equations of motion, an efficient matrix form of impulse Kane formulations is derived for systems exposed to impulsive forces.