In the present paper, we study the first and second order of accuracy difference schemes for the approximate solution of the inverse problem for hyperbolic–parabolic equations with unknown time-independent source term. The unique solvability of constructed difference schemes and the stability estimates for their solutions are obtained. The proofs are based on the spectral representation of the self-adjoint positive definite operator in a Hilbert space.