In this paper, we use asymptotic techniques and the finite differences method to study the spectrum of the differential operator arising in the exponential stabilization of the following hyperbolic partial differential equation      with boundary conditions      We build a numerical scheme and investigate the eigenvalues location as a function of the positive feedback parameters a and b. We present also a qualitative study of the position of spectrum with respect to the vertical asymptote.