This paper deals with implementation of the multistage Adomian decomposition method (MADM) to solve a class of nonlinear programming (NLP) problems, which are reformulated with a nonlinear system of fractional differential equations. The multistage strategy is used to investigate the relation between an equilibrium point of the fractional order dynamical system and an optimal solution of the NLP problem. The preference of the method lies in the fact that the multistage strategy gives this relation in an arbitrary longtime interval, while the Adomian decomposition method (ADM) gives the optimal solution just only in the neighborhood of the initial time. The numerical results taken by the fractional order MADM show that these results are compatible with the solution of NLP problem rather than the ADM. Furthermore, in some cases the fractional order MADM can perform more rapid convergency to the optimal solution of optimization problem than the integer order ones.