The density profiles at z = 0 of 142 dark matter halos with masses between 10⁸ and 10⁹   M_⊙. The dashed line indicates the limit of the resolution corresponding to the gravitational softening. The gray surface show the domain covered by the pseudo-isothermal profiles (Eq. (9)) adopted in our simulations, for the minimum and maximum values of central densities and radii (r_max). Only two halos among the 144 extracted strongly deviate from the NFW profile and are not seen in this plot. The reason is that at z = 0, they both undergo a major merging.