A fast numerical method for fractional partial differential equations

In this paper, we use operational matrices of Chebyshev polynomials to solve fractional partial differential equations (FPDEs). We approximate the second partial derivative of the solution of linear FPDEs by operational matrices of shifted Chebyshev polynomials. We apply the operational matrix of integration and fractional integration to obtain approximations of (fractional) partial derivatives of the solution and the approximation of the solution. Then we substitute the operational matrix approximations in the FPDEs to obtain a system of linear algebraic equations. Finally, solving this system, we obtain the approximate solution. Numerical experiments show an exponential rate of convergence and hence the efficiency and effectiveness of the method.