  The time evolution of a low-energy two-dimensional Gaussian wave packet in
ABC-stacked $n$-layer graphene (ABC-NLG) is investigated. Expectation values of
the position $(x,y)$ of center-of-mass and the total probability densities of
the wave packet are calculated analytically using the Green's function method.
These results are confirmed using an alternative numerical method based on the
split-operator technique within the Dirac approach for ABC-NLG, which
additionally allows to include external fields and potentials. The main
features of the zitterbewegung (trembling motion) of wave packets in graphene
are demonstrated and are found to depend not only on the wave packet width and
initial pseudospin polarization, but also on the number of layers. Moreover,
the analytical and numerical methods proposed here allow to investigate wave
packet dynamics in graphene systems with an arbitrary number of layers and
arbitrary potential landscapes.
