This paper presents a novel algebraic method for the construction of Lyapunov functions to study global stability of the disease-free equilibrium points of deterministic epidemic ordinary differential equation models with mass action and standard incidence functions. The method is named as Jacobian-Determinant method. In our method, a direct algebraic procedure that also relies only on determinant of the Jacobian matrix of the infected subsystem is developed to determine a threshold quantity, R'_R0′">0 akin to the basic reproduction number, R_R0">0 of such class of models. The developed technique is applied on a wide variety of models to construct Lyapunov functions to study the global stability of the infection-free critical points. Further, implementation of our method reveals that the threshold quantity is the same as (or the square) of the basic reproduction numbers as obtained using the next-generation matrix method. It is further observed that even for models that do not use the standard or mass action incidence, the threshold quantity is still related to the basic reproduction numbers as obtained with the next-generation matrix method.

link to the article: https://www.sciencedirect.com/science/article/pii/S0960077923003041