In this study, we have studied the multi-order differential equations. The model we have followed agrees with initial value problem which, in its turn, has a group of linear ordinary differential equations. This paper's aim is programming a Variational Iteration Method (VIM) using Wolfram's Mathematica. Variational Iteration Method offers a study that introduce approximated solutions of the multi-order ordinary differential equations. Several examples of different order have been resolved.