To provide mathematically rigorous eigenvalue bounds for the Steklov eigenvalue problem, an enhanced version of the eigenvalue estimation algorithm developed by the third author is proposed, which removes the requirements of the positive definiteness of bilinear forms in the formulation of eigenvalue problems. In practical eigenvalue estimation, the Crouzeix--Raviart finite element method (FEM) along with quantitative error estimation is adopted. Numerical experiments for eigenvalue problems defined on a square domain and an L-shaped domain are provided to validate the precision of computed eigenvalue bounds.

15 pages, 5 figures, 3 tables