In this Letter, approximate analytical solutions of systems of Fractional Differential Equations (FDEs) are derived by the Homotopy-Perturbation Method (HPM). The fractional derivatives are described in the Caputo sense. The solutions are obtained in the form of rapidly convergent infinite series with easily computable terms. Numerical results reveal that HPM is very effective and simple for obtaining approximate solutions of nonlinear systems of FDEs. (C) 2007 Elsevier B.V. All rights reserved.