In this paper, the method is developed to differential–algebraic equations systems. More effective method is presented and illustrated by numerical example. The solution for a differential–algebraic equation can be expanded up to arbitrary order using MAPLE computer algebra systems. First we calculate power series of the given equations system then transform it into Padé series form, which give an arbitrary order for solving differential–algebraic equation numerically.