An exponential convergence rate has been obtained for approximate eigenvalues as well as for the eigenfunctions. The result is shown when the boundary of the underlying domain is smooth and also when the boundary contains corners. The computational results confirm the error estimates. The results shown in the document are valid for multiple or grouped eigenvalues. The method works for problems not attached as well.