The global stability analysis for the mathematical model of an infectious disease is discussed here. The endemic equilibrium is shown to be globally stable by using a modification of the Volterra–Lyapunov matrix method. The basis of the method is the combination of Lyapunov functions and the Volterra–Lyapunov matrices. By reducing the dimensions of the matrices and under some conditions, we can easily show the global stability of the endemic equilibrium. To prove the stability based on Volterra–Lyapunov matrices, we use matrices with the symmetry properties (symmetric positive definite). The results developed in this paper can be applied in more complex systems with nonlinear incidence rates. Numerical simulations are presented to illustrate the analytical results.