This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform
method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed
for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential
equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components.
Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to
implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for
solving many linear and nonlinear fractional differential equations