We investigate a numerical solution of the covariant Bethe-Salpeter equation in the Euclidean space for heavy
meson with gluon ladder in the Landau gauge and scalar confinement. A new approach is presented to solve the non-linear eigenvalue problem with suitable bases and fictitious eigenvalue parameters. We obtain unphysical states when the equation is solved for timelike spectra. We also present how to cover the singularity of a free quark propagator and Schwinger-Dyson equation when extrapolated to the timelike domain in the real plane.