In this paper, we present a numerical method for solving a class of fractional differential equations (FDEs). Based on Bernstein Polynomials (BPs) basis, new matrices are utilized to reduce the multi-term orders fractional differential equation to a system of algebraic equations. Convergence analysis is shown by several theorems. Illustrative examples are included to demonstrate the validity and applicability of this method.
