  Using the supersymmetric formalism we compute exactly at finite $N$ the
expectation of the Wilson loop in the Gaussian Unitary Ensemble and derive an
exact formula for the spectral density at finite $N$. We obtain the same result
by a second method relying on enumerative combinatorics and show that it leads
to a novel proof of the Harer-Zagier series formula.
