  In this paper we study systems of autonomous algebraic ODEs in several
differential indeterminates. We develop a notion of algebraic dimension of such
systems by considering them as algebraic systems. Afterwards we apply
differential elimination and analyze the behavior of the dimension in the
resulting Thomas decomposition. For such systems of algebraic dimension one, we
show that all formal Puiseux series solutions can be approximated up to an
arbitrary order by convergent solutions. We show that the existence of Puiseux
series and algebraic solutions can be decided algorithmically. Moreover, we
present a symbolic algorithm to compute all algebraic solutions. The output can
either be represented by triangular systems or by their minimal polynomials.
