The authors have employed the multistage variational iteration method (MVIM) to solve a system of first order differential equations with quadratic nonlinearity of the form
subject to certain initial conditions. Earlier some authors have solved () by variational iteration method (VIM) on . The present technique enables the authors to go for larger domain. Comparison of MVIM is made with VIM and fourth order Runge-Kutta method and the outcome seems to be satisfactory. In the conclusion, the authors have mentioned several advantages of MVIM over other methods. They are hopeful of applying MVIM to problems with oscillatory solutions in their future work.