© 2025 the author(s), published by De Gruyter, Berlin/Boston.Let { ? n = { ?1,n, ?, ? dn, n}}n be a sequence of finite multisets of real numbers such that dn ? ? as n ? ?, and let f: ? ? ?d ? ? be a Lebesgue measurable function defined on a domain ? with 0 < ?d (?) < ?, where ?d is the Lebesgue measure in ?d. We say that { ?n}n has an asymptotic distribution described by f, and we write { ?n}n ? f, if (?) lim n? ? 1/dn ? =dni=1 F (?i, n) = 1/?d (?) F (f (x)) dx for every continuous function F with bounded support. If ?n is the spectrum of a matrix An, we say that { An }n has an asymptotic spectral distribution described by f and we write { An }n ? ? f. In the case where d = 1, ? is a bounded interval, ?n ? f(?) for all n, and f satisfies suitable conditions, Bogoya, Böttcher, Grudsky, and Maximenko proved that the asymptotic distribution (?) implies the uniform convergence to 0 of the difference between the properly sorted vector [ ?1,n, ?, ?d n, n ] and the vector of samples [ f (x 1, n), ?, f (xd n, n) ], i.e., (??) lim n? ? maxi = 1, ?, dn|f (xi, n) - ?? n (i), n|= 0, where x1, n, ?, xd n, n is a uniform grid in ? and ? n is the sorting permutation. We extend this result to the case where d ? 1 and ? is a Peano-Jordan measurable set (i.e., a bounded set with ? d(??) = 0). We also formulate and prove a uniform convergence result analogous to (??) in the more general case where the function f takes values in the space of k × k matrices. Our derivations are based on the concept of monotone rearrangement (quantile function) as well as on matrix analysis arguments stemming from the theory of generalized locally Toeplitz sequences and the observation that any finite multiset of numbers ?n = { ?1, n, ?, ?d n, n } can always be interpreted as the spectrum of a matrix An = diag (?1, n, ?, ?d n, n). The theoretical results are illustrated through numerical experiments, and a reinterpretation of them in terms of vague convergence of probability measures is hinted.