Meshless methods are a promising field of numerical methods recently introduced to computational electromagnetics. The potential of conformal and multi-scale modeling and the possibility of dynamic grid refinements are very attractive features that appear more naturally in meshless methods than in classical methods. The radial point interpolation method (RPIM) uses radial basis functions for the approximation of spatial derivatives. In this publication an eigenvalue solver is introduced for RPIM in electromagnetics. Eigenmodes are calculated on the example of a cylindrical resonant cavity. It is demonstrated that the computed resonance frequencies converge to the analytical values for increasingly fine spatial discretization. The computation of eigenmodes is an important tool to support research on a time-domain implementation of RPIM. It allows a characterization of the method's accuracy and to investigate stability issues caused by the possible occurrence of non-physical solutions