“…In recent years, much research has been focused on the numerical solution of systems of differentialalgebraic equations. Some numerical methods have been developed, such as implicit Runge-Kutta method [1], Padé approximation method [2,4,5,16], homotopy perturbation method [14], Adomain decomposition method [8,9] and variation iteration method [15]. Consider the following system of differential-algebraic equations of fractional order [18] In this paper, the homotopy perturbation method (HPM) is used to give a semi-analytical solution of the problem (1.1).…”