Three panels depicting cumulative mass profiles of stars at 708 Myr under different scenarios: BFSR (solid blue), observational data (solid black), and ensemble averages from multiple simulations with varied parameters (dashed red with 1 σ uncertainty in grey). Top and centre panels: stars above the completeness limit derived from a random realisation of the distribution of initial coordinates and velocities, as well as the IMF, respectively. Bottom panel: extended to all stars, averaging all 15 simulations. It also shows the Jacobi mass, M_J(r), as a function of the Jacobi radius (by inverting Eq. (3), dotted green).