Solving a multi-order fractional differential equation using adomian decomposition

An algorithm has been developed to convert the multi-order fractional differential equation: D(*)(alpha)y(t) = f(t,y(t),D(*)(beta 1)y(t),...,D-*(beta)y(t)), y((k))(0) = ck, k = 0,...m, where m < alpha <= m+ 1, 0 < beta(1) < beta(2) < (. . .) < beta(n)< alpha and D-*(alpha) denotes Caputo fractional derivative of order a into a system of fractional differential equations. Further Adomian decomposition method has been employed to solve the system of fractional differential equations. Some illustrative examples are presented. (c) 2006 Elsevier Inc. All rights reserved.