A new method for constructing the Lyapunov function and determining the stability of discrete time systems with a computer on the basis of the similarity transformation theory is presented. The system matrix of the system under discussion is applied directly instead of solving the discrete Lyapunov matrix equation.