How to numerically compute the eigenvalues of the Laplacian in a triangular domain with Dirichlet boundary condition?

Consider an arbitrary triangle. Now impose the homogeneous Dirichlet boundary condition. How to compute the eigenvalues and eigenvectors of the Laplacian $$-\nabla^2 = - \frac{\partial^2}{\partial x^2} - \frac{\partial^2}{\partial y^2}$$ numerically?

• You should probably state what is intended by Dirichlet boundary conditions for the reader's sake. Commented Jul 23, 2014 at 21:28

where $\Delta$ is the triangle. In the first transformation we multiply by test functions from a chosen set $V$ which satisfy a homogeneous Dirichlet condition. In the second transformation we integrate both sides. In the third transformation we integrate by parts. In the last transformation we use the boundary condition of the test functions.
If we now choose to have $u$ approximated by continuous piecewise linear functions on a finite element mesh of the triangle which satisfy the boundary condition, and use test functions $v$ of the same type, then the last equation is an ordinary matrix eigenvalue problem.