We consider expectations of the form $E[trh_1 (X^N_1) \dots trh_r(X^N_r)]$, where $X^N_i$ are self-adjoint polynomials in various independent classical random matrices and $h_i$ are smooth test function and obtain a large $N$ expansion of these quantities, building on the framework of polynomial approximation and Bernstein-type inequalities recently developed by Chen, Garza-Vargas, Tropp, and van Handel. As applications of the above, we prove the higher-order asymptotic vanishing of cumulants for smooth linear statistics, establish a Central Limit Theorem, and demonstrate the existence of formal asymptotic expansions for the free energy and observables of matrix integrals with smooth potentials. This talk is based on joint work with Benoît Collins