The paper investigates an algorithm for the numerical solution of a parametric eigenvalue problem for the Helmholtz equation on the plane specially tailored for the accurate mathematical modeling of lasing modes of microring lasers. The original problem is reduced to a nonlinear eigenvalue problem for a system of Muller boundary integral equations. For the numerical solution of the obtained problem, we use a trigonometric Galerkin method, prove its convergence, and derive error estimates in the eigenvalue and eigenfunction approximation. Previous numerical experiments have shown that the method converges exponentially. In the current paper, we prove that if the generalized eigenfunctions are analytic, then the approximate eigenvalues and eigenfunctions exponentially converge to the exact ones as the number of basis functions increases. To demonstrate the practical effectiveness of the algorithm, we find geometrical characteristics of microring lasers that provide a significant increase in the directivity of lasing emission, while maintaining low lasing thresholds.