We use a direct numerical integration of the Vlasov equation in spherical symmetry with a background gravitational potential to determine the evolution of a collection of particles in different models of a galactic halo. Such a collection is assumed to represent a dark matter inhomogeneity which reaches a stationary state determined by the virialization of the system. We describe some features of the stationary states and, by using several halo models, obtain distinctive signatures for the evolution of the inhomogeneities in each of the models.