This study presents numerical methods-fractional difference and Adomian decomposition-for solution of a dynamic system containing fractional derivative of order α, 0<α≤1. The fractional derivative is described in the Caputo sense. The Adomian decomposition method provides the solution in the form of a convergent power series with easily computable components. Then the diagonal Padé approximants are effectively used in the analysis to capture the essential behavior of the solution. The presented schemes are introduced in a general way so that they can be used in applied sciences. © 2008 Asian Network for Scientific Information.