摘要

We develop a numerical method by using operational matrices of fractional order integrations and differentiations to obtain approximate solutions to a class of coupled systems of fractional order partial differential equations (FPDEs). We use shifted Legendre polynomials in two variables. With the help of the aforesaid matrices, we convert the system under consideration to a system of easily solvable algebraic equation of Sylvester type. During this process, we need no discretization of the data. We also provide error analysis and some test problems to demonstrate the established technique.