Starting from the Мещерский equations, the higher order D'Alembert-Lagrange principle for mechanical system of variable mass is obtained, and different kinds of higher order differential equations of motion for holonomic mechanical system of variable mass are derived. The result shows that these equations extend and optimize the correlative theory of holonomic mechanics.