In this paper, we consider a method for finding certain eigenvalues of generalized eigenvalue problems in a given domain of the complex plane. We also discuss the relation between the presented method and the Lanczos method briefly. The presented method provides a good performance in parallel computations. A numerical example that was obtained on a PC cluster is included.