Comparison of the observational data and the empirical orbit superposed model: the 3D velocity distributions v_r, v_φ, and v_θ in different spatial bins. We divided the system into 7 × 5 bins in r versus θ, but eliminating the first row with r < 8 kpc and only showing the odd number of columns here. In each sub-panel, the solid and dashed grey curves are the mean and 1σ uncertainty of the distribution constructed from observational data. The red, dark blue, and light blue curves were constructed by models with different radial distributions of the underlying DM mass. The model with a DM mass of log(ρ₀[M_⊙/kpc³]) = 6.2, r_s = 70 kpc, and γ = 1.0 (red) matches the data well. The velocity distribution of stars, especially v_r, strongly constrains the radial distribution of the underlying DM distributions.