A simulation of the long-term evolution of the synthetic Flora (top), Koronis (middle) and Eos (bottom) families in the proper semi-major axis a_p vs. the pole latitude β plane. Left: objects larger than D > 30 km, which almost do not evolve in β. Right: objects with D ≤ 30 km, with the initial conditions denoted by empty circles and an evolved state at 1 Gyr denoted by full circles. The sizes of symbols correspond to the actual diameters D. The initial conditions for Flora correspond to an isotropic size-independent velocity field with α = 3.25 and v_esc = 95 m s^-1, and a uniform distribution of poles (i.e. sinβ). We increase the number of objects 10 times compared to the observed members of the Flora (Koronis and Eos as well) family to improve statistics. We retain their size distribution, of course. The objects in Flora family are discarded from these plots when they left the family region (eccentricity e_p = 0.1to0.18, inclination sinI_p = 0.05to0.13), because they are affected by strong mean-motion or secular resonances (ν₆ in this case). Thermal parameters were set as follows: the bulk density ρ_bulk = 2500 kg m^-3, the surface density ρ_surf = 1500 kg m^-3, the thermal conductivity K = 0.001 W m^-1 K^-1, the thermal capacity C_t = 680 J kg^-1, the Bond albedo A = 0.1, and the infrared emissivity ϵ = 0.9. The time step for the orbital integration is dt = 91 days and dt_spin = 10³ yr for the (parallel) spin integration. The parameters for Koronis and Eos are chosen similarly, only for Koronis do we use v_esc = 100 m s^-1, and v_esc = 225 m s^-1 and ρ_surf = 2500 kg m^-3 for Eos.