n this paper, we propose a numerical algorithm for solving system of fractional differential equations by using the homotopy analysis transform method. The homotopy analysis transform 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