In this paper, a fractional order differential system for modeling human T-cell lymphotropic virus I (HTLV-I) infection of CD4+ T-cells is studied and its approximate solution is presented using a multi-step generalized differential transform method. The method is only a simple modification of the generalized differential transform method, in which it is treated as an algorithm in a sequence of small intervals (i.e. time step) for finding accurate approximate solutions to the corresponding systems. The solutions obtained are also presented graphically. © 2011 Elsevier Ltd. All rights reserved.