  This paper introduces an algorithmic approach to the analysis of bifurcation
of limit cycles from the centers of nonlinear continuous differential systems
via the averaging method. We develop three algorithms to implement the
averaging method. The first algorithm allows to transform the considered
differential systems to the normal formal of averaging. Here, we restricted the
unperturbed term of the normal form of averaging to be identically zero. The
second algorithm is used to derive the computational formulae of the averaged
functions at any order. The third algorithm is based on the first two
algorithms that determines the exact expressions of the averaged functions for
the considered differential systems. The proposed approach is implemented in
Maple and its effectiveness is shown by several examples. Moreover, we report
some incorrect results in published papers on the averaging method.
