In this manuscript, a new method is introduced for solving multi-order fractional differential equations. By transforming the fractional differential equations into an optimization problem and using polynomial basis functions, we obtain the system of algebraic equation. Then, we solve the system of nonlinear algebraic equation and obtain the coefficients of polynomial expansion. Also, we show the convergence of the method. Some numerical examples are presented which illustrate the theoretical results and the performance of the method.