projected density profile and the circularized effective radius [$\Re\equiv\sqrt{\sae\sbe}$]{} (where [$\sae$]{} and [$\sbe$]{} are the major and minor semi-axis of the effective isodensity ellipse). We fit (over the radial range $0.1 \, \lsim \, R/\Re \, \lsim \, 10$) the circularized projected density profiles of the end-products with the $R^{1/m}$ Sersic (1968) law: $$\label{eqser}
I(R)=\Ie \,
\exp\left\{-b(m)\left[ \left( \frac{R}{\Re} \right)^{1/m} -1 \right]\right\},$$ where $\Ie\equiv I(\Re)$ and $b(m)\simeq 2m-1/3+4/405m$ (Ciotti & Bertin 1999). In the fitting procedure $m$ is the only free parameter, because $\Re$ and $\Ie$ are determined by their measured values obtained by particle count. In addition, we measure the central velocity dispersion $\sgz$, obtained by averaging the projected velocity dispersion over the circularized surface density profile within an aperture of $\Re/8$. Some of these structural parameters are reported in Table 2 for the five simulations described above, as well as for three additional simulations, which start from different initial conditions (see Section \[secdis\]).

Results {#secres}
=======

In Newtonian gravity, collisionless systems reach virialization through violent relaxation in few dynamical times, as predicted by the theory (Lynden-Bell 1967) and confirmed by numerical simulations (e.g. van Albada 1982). On the other hand, due to the non linearity of the theory 