n this paper, we propose a numerical algorithm for solving system of fractional differ-
ential equations by using the homotopy analysis transform method. The homotopy analysis trans-
form method is the combined form of the homotopy analysis method and Laplace transform
method. The solutions of our modeled equations are calculated in the form of convergent power
series with easily computable components. The numerical results shows that the approach is easy
to implement and accurate when applied to various fractional differential equations