Abstract
A new operational matrix of fractional integration of arbitrary order for generalized Laguerre polynomials is derived. The fractional integration is described in the Riemann-Liouville sense. This operational matrix is applied together with generalized Laguerre tau method for solving general linear multiterm fractional differential equations (FDEs). The method has the advantage of obtaining the solution in terms of the generalized Laguerre parameter. In addition, only a small dimension of generalized Laguerre operational matrix is needed to obtain a satisfactory result. Illustrative examples reveal that the proposed method is very effective and convenient for linear multiterm FDEs on a semi-infinite interval.