Abstract
In the present work, a numerical technique for solving a general form of nonlinear fractional order integro-differential equations (GNFIDEs) with linear functional arguments using Chebyshev series is presented. The recommended equation with its linear functional argument produces a general form of delay, proportional delay, and advanced non-linear arbitrary order Fredholm–Volterra integro-differential equations. Spectral collocation method is extended to study this problem as a matrix discretization scheme, where the fractional derivatives are characterized in the Caputo sense. The collocation method transforms the given equation and conditions to an algebraic nonlinear system of equations with unknown Chebyshev coefficients. Additionally, we present a general form of the operational matrix for derivatives. The introduced operational matrix of derivatives includes arbitrary order derivatives and the operational matrix of ordinary derivative as a special case. To the best of authors’ knowledge, there is no other work discussing this point. Numerical test examples are given, and the achieved results show that the recommended method is very effective and convenient.