This paper presents a fractional-order nonlinear system of delay integro-differential equations which is used to model biological species living together. A numerical scheme based on the Chelyshkov polynomials is implemented to solve this fractional order system of integro-differential equations. The main advantages of the presented method is that it reduces under consideration problem to a system of nonlinear algebraic equations. Some test problems are considered to confirm validity and accuracy of the presented method. Moreover, the obtained numerical results are compared with those existing in the literature.Numerical treatment of fractional-order nonlinear system are derived in Sec. 4. In Sec. 5, an efficient numerical method is introduced to solve fractional-order nonlinear system of delay integro-differential equations (3). To confirm efficiency and accuracy of the presented method, we give some examples and numerical results in Sec. 6. Finally, a conclusion is drawn in Sec. 7.
Preliminary Remarks on Fractional CalculusFractional order calculus is a branch of calculus which deals with integration and differentiation operators of non-integer order. Although the fractional order operators enable to model a wider class of problems, there not exists a unique definition of fractional derivative. Among the several formulations of the generalized derivative, the Riemann-Liouville and Liouville-Caputo definition is the most commonly used, which can be described as follows (see e.g. [7,13]).