An analytical approach to determining the coefficients in Lyapunov direct method: With application to an age-structured epidemiological model

Abstract

Lyapunov direct method has been employed as a powerful tool to show the global stability of an equilibrium in differential equations. One widely cited Lyapunov function for population dynamic models takes the form L(t)=∑i=1nciLi(t) with Li(t)=xi(t)−xi∗−xi∗ln(xi(t)/xi∗) (if xi∗>0) and Li(t)=xi(t) (if xi∗=0), a combination of functions involving different state variables xi(t) of the model system with (x1∗,x2∗,…,xn∗) being an equilibrium. However, two challenges hinder the efficient applications of Lyapunov direct method: (a) determining the coefficients ci; and (b) rearranging the time derivative L′(t) along the trajectories of the system to show it is negative (semi-)definite. This study is to propose an easy-to-follow analytical approach to tackle these two challenges, which will be illustrated through an application to an epidemiological model with vaccination age. Furthermore, the Lyapunov functional for the endemic equilibrium can be reformulated to investigate the global stability for the disease-free equilibrium and a family of Lyapunov functionals can be proposed for the same purpose. It is expected that the approach can be further applied to other age-structured models and be extended to analyze more complicated models with other heterogeneous factors.