In this study, a decomposition method for approximating the solution of systems of fractional integro-differential equations is implemented. The fractional derivative is considered in the Caputo sense. The approximate solutions are calculated in the form of a convergent series with easily computable components. Numerical results show that this approach is easy to implement and accurate when applied to integro-differential equations. © 2006 Elsevier Ltd.