In this paper the discrete eigenvalues of elliptic second order differential operators in $L^2(mathbb{R}^n)$, $n in mathbb{N}$, with singular $delta$- and $delta$-interactions are studied. We show the self-adjointness of these operators and derive equivalent formulations for the eigenvalue problems involving boundary integral operators. These formulations are suitable for the numerical computations of the discrete eigenvalues and the corresponding eigenfunctions by boundary element methods. We provide convergence results and show numerical examples.