For the eigenvalue problem of the Laplacian, the conforming and non-conforming finite element methods are utilized to provide explicit lower and upper bounds for the eigenvalues.
Particularly, by applying Liu's method, the lower eigenvalue bounds can be easily obtained through the Crouzeix-Rarviart finite element method; such a method even works for non-convex domains.
The lower and upper eigenvalue bounds for both the unit square domain and an L-shaped domain are given in the notebook file.
The code is prepared as a Jupyter notebook. If you have Jupyter installed on your computer, just download the zipped file and open it in Jupyter.