  The magnetic field integral equation for axially symmetric cavities with
perfectly conducting surfaces is discretized according to a high-order
convergent Fourier--Nystr\"om scheme. The resulting solver is used to determine
eigenwavenumbers and normalized magnetic eigenfields to very high accuracy in
the entire computational domain.
