n this study, an effective numerical method for solving fractional differential equa- tions using Chebyshev cardinal functions is presented. The fractional derivative is described in the Caputo sense. An operational matrix of fractional order integration is derived and is utilized to reduce the fractional differential equations to system of algebraic equations. In addition, illustrative examples are presented to demonstrate the efficiency and accuracy of the proposed method