Average mass in each Lagrangian shell 〈M〉_Lagr[M_⊙] (top) and the Lagrangian radii R_Lagr [pc] (bottom), both normalized by the cluster mass at the beginning of the simulation. We show the N-body and MOCCA results for the Lagrangian radii. The two plot show the shell representing 1%, 10%, 30%, 50%, 70%, and 90%. The differences between the Lagrangian radii in the N-body and MOCCA simulations arise from how escaping stars are treated. In the N -body simulations, stars with positive energy are not immediately removed from the system; instead, they remain until they cross a boundary set at ten times the half-mass radius. These stars are still included in the calculation of Lagrangian radii, which affects the results compared to MOCCA.