Abstract
In this paper, we study shifted Jacobi polynomials and develop a simple but highly accurate scheme for the numerical solution of coupled system of fractional differential equations. We derive some operational matrices of integration and differentiation of fractional order. By the application of these matrices we provide a theoretical treatment to approximate the solutions of the corresponding system. We use Matlab to perform necessary operations. The applicability of the technique is shown with some examples and the results are displayed graphically.