 Solution

Put $D=|⟨L_N,f⟩-⟨{\hat{L}}_N,f⟩|$. Apply the Cauchy-Schwarz inequality to the average of the nonnegative eigenvalue differences, and then the Hoffman–Wielandt inequality:Taking nonnegative square roots gives the spectral Lipschitz bound from Frobenius distanceThe factor $N^{-1/2}$ is essential: the empirical spectral measure has total mass $1$, not $N$.