The operational matrix of fractional derivative of the fractional-order Chebyshev functions and its applications

In this paper, we introduce a family of fractional-order Chebyshev functions based on the classical Chebyshev polynomials. We calculate and derive the operational matrix of derivative of fractional order $gamma$ in the Caputo sense using the fractional-order Chebyshev functions. This matrix yields to low computational cost of numerical solution of fractional order differential equations to the solution of a system of algebraic equations. Several numerical examples are given to illustrate the accuracy of our method. The results obtained, are in full agreement with the analytical solutions and numerical results presented by some previous works.