A numerical scheme for solving variable order Caputo-Prabhakar fractional integro-differential equation

In this paper we apply the Chebyshev polynomials method for the numerical solution of a class of variable-order fractional integro-differential equations with initial conditions. Moreover, a class of variable-order fractional integro-differential equations with fractional derivative of Caputo-Prabhakar sense is considered. The main aim of the Chebyshev polynomials method is to derive four kinds of operational matrices of Chebyshev polynomials. With such operational matrices, an equation is transformed into the products of several dependent matrices, which can also be viewed as the system of linear equations after dispersing the variables. Finally, numerical examples have been presented to demonstrate the accuracy of the proposed method, and the results have  been compared with the exact solution.