On Derivation of Motion Equations for Systems with Non-Holonomic High-Order Program Constraints