Dynamical models of prototype gravastars were constructed in order to study their stability. The models are the Visser-Wiltshire three-layer gravastars, in which an infinitely thin spherical shell of stiff fluid divides the whole spacetime into two regions, where the internal region is de Sitter, and the external is Schwarzschild. It is found that in some cases the models represent the "bounded excursion" stable gravastars, in which the thin shell is oscillating between two finite radii, while in other cases they collapse until the formation of black holes. In the phase space, the region for the ``bounded excursion" gravastars is very small in comparison to that of black holes, but not empty. Therefore, although the existence of gravastars cannot be excluded from such dynamical models, our results do indicate that, even if gravastars indeed exist, they do not exclude the existence of black holes.

This paper presents a review of the fractal approach for describing the large scale distribution of galaxies. We start by presenting a brief, but general, introduction to fractals, which emphasizes their empirical side and applications rather than their mathematical side. Then we discuss the standard correlation function analysis of galaxy catalogues and many observational facts that brought increasing doubts about the reliability of this method, paying special attention to the standard analysis implicit assumption of an eventual homogeneity of the distribution of galaxies. Some new statistical concepts for analysing this distribution is presented, and without the implicit assumption of homogeneity they bring support to the hypothesis that the distribution of galaxies does form a fractal system. The Pietronero-Wertz's single fractal (hierarchical) model is presented and discussed, together with the implications of this new approach for understanding galaxy clustering.