  It is given an algorithm to obtain generalized power asymptotic expansions of
the solutions of the Einstein equations arising for several homogeneous
cosmological models. This allows to investigate their behavior near the initial
singularity or for large times. An implementation of this algorithm in the CAS
system Maple V Release 4 is described and detailed calculations for three
equations are shown.
