  We discuss spectral properties of the Laplacian with multiple ($N$) point
interactions in two-dimensional bounded regions. A mathematically sound
formulation for the problem is given within the framework of the self-adjoint
extension of a symmetric (Hermitian) operator in functional analysis. The
eigenvalues of this system are obtained as the poles of a transition matrix
which has size $N$. Closely examining a generic behavior of the eigenvalues of
the transition matrix as a function of the energy, we deduce the general
condition under which point interactions have a substantial effect on
statistical properties of the spectrum.
