We extend the operational matrices technique to design a spectral solution of nonlinear fractional differential equations (FDEs). The derivative is considered in the Caputo sense. The coupled system of two FDEs is considered, subjected to more generalized integral type conditions. The basis of our approach is the most simple orthogonal polynomials. Several new matrices are derived that have strong applications in the development of computational scheme. The scheme presented in this article is able to convert nonlinear coupled system of FDEs to an equivalent S-lvester type algebraic equation. The solution of the algebraic structure is constructed by converting the system into a complex Schur form. After conversion, the solution of the resultant triangular system is obtained and transformed back to construct the solution of algebraic structure. The solution of the matrix equation is used to construct the solution of the related nonlinear system of FDEs. The convergence of the proposed method is investigated analytically and verified experimentally through a wide variety of test problems.
Keywords: approximation; iterative methods; numerical simulation.