“…This has further generated further and more intense interest in finding numerical methods which accurately and efficiently solve nonlinear fractional differential equations. Some of the methods include Adomian Decomposition Method (ADM) [7][8][9][10][11][12][13], the Variational Iteration Method (VIM) [13][14][15][16], Homotopy Analysis Method (HAM) [17][18][19], Homotopy perturbation Method (HPM) [13,20]. Many physical problems are governed by a system of differential-Algebraic Equations (DAE's), and the solution of these equations has been a subject of many investigations in recent years [7][8][9][14][15][16]18,19,21,22].…”