Operational matrices to solve nonlinear Riccati differential equations of arbitrary order

In this paper, an effective numerical method to achieve the numerical solution of nonlinear Riccati differential equations of arbitrary (integer and fractional) order has been developed. For this purpose, the fractional order of the Chebyshev functions (FCFs) based on the classical Chebyshev polynomials of the first kind have been introduced, that can be used to obtain the solution of these equations. Also, the operational matrices of fractional derivative and product for the FCFs have been constructed. The obtained results illustrated demonstrate that the suggested approaches are applicable and valid.