The Jeffreys model of heat conduction is a system of two partial differential equations of mixed hyperbolic and parabolic character. The analysis of an initial-boundary value problem for this system is given. Existence and uniqueness of a weak solution of the problem under very weak regularity assumptions on the data is proved. A finite difference approximation of this problem is discussed as well. Stability and convergence of the discrete problem are proved.