Computational and Mathematical Methods in Medicine

Volume 2012 (2012), Article ID 826052, 8 pages

doi:10.1155/2012/826052

Research Article

Global Stability Analysis of SEIR Model with Holling Type II Incidence Function

Mohammad A. Safi¹ and Salisu M. Garba²

¹Department of Mathematics, The Hashemite University, Zarqa 13115, Jordan

²Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria 0002, South Africa

Received 27 June 2012; Revised 9 August 2012; Accepted 12 September 2012

Academic Editor: Jacek Waniewski

Copyright © 2012 Mohammad A. Safi and Salisu M. Garba. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A deterministic model for the transmission dynamics of a communicable disease is developed and rigorously analysed. The model, consisting of five mutually exclusive compartments representing the human dynamics, has a globally asymptotically stable disease-free equilibrium (DFE) whenever a certain epidemiological threshold, known as the basic reproduction number (), is less than unity; in such a case the endemic equilibrium does not exist. On the other hand, when the reproduction number is greater than unity, it is shown, using nonlinear Lyapunov function of Goh-Volterra type, in conjunction with the LaSalle's invariance principle, that the unique endemic equilibrium of the model is globally asymptotically stable under certain conditions. Furthermore, the disease is shown to be uniformly persistent whenever .

1. Introduction

Mathematical models have been widely used to gain insight into the spread and control of emerging and reemerging disease. The dynamics of these models is usually determined by a threshold quantity known as the basic reproduction number (denoted by ), which is defined as the number of secondary cases generated by an infected individual in a completely susceptible population [1–5]. Characteristically, when is less than unity, a small influx of infected individuals will not generate large outbreaks, and the disease dies out in time. On the other hand, when exceeds unity, the disease will persist. A basic epidemic model supports at least two equilibria (a disease-free equilibria and endemic equilibria); plays important role in the study of equilibria of a model. Several models found in the literature [2, 4–15] have been used to show that when crosses the threshold, , a transcritical bifurcation takes place. That is, asymptotic local stability is transferred from the disease-free state to the new (emerging) endemic (positive) equilibria. In some cases, it can be shown that the transfer of asymptotic stability is independent of initial conditions; that is, it is global; see, for instance, [6, 8–12, 16]. Establishing global properties of a dynamical system using Lyapunov function is generally a nontrivial problem. This is owing to the fact that there are no systematic methods for constructing Lyapunov function for infectious disease models with standard incidence rate [17]. The most successful approach to the problem is the direct Lyapunov method which involves the use of quadratic function of the form or by using nonlinear Lyapunov function of Goh-Volterra type of the form . However, other methods used in establishing global properties of some epidemic models include Dulac’s criterion to eliminate the existence of the periodic solution and prove the global stability by the Poincaré Bendixson theorem [18] and those reported in Kamgang and Sallet [19] and Qiao et al. [20].

Let , and denote the number of susceptible individuals, infectious individuals, and the total size of the population at time , respectively. Further, let be the average number of contacts that is sufficient to transmit infection (effective contact rate). Then, the force of infection, given by , represents the average number of contacts a susceptible individual makes with infectious individuals per unit time. If (i.e., the contact rate depends on the total population, ), then the incidence function is called mass action incidence. If (a constant), then the incidence function is called standard incidence [4]. These two functions are widely used in modeling the transmission dynamics of the human diseases [1, 21]. Another widely used incidence function is the Holling type II incidence function, given by , with , [22–28].

The nonlinear incidence function of type was first introduced by Capasso and Serio [22], in their study of cholera epidemic. The main justification for using such a functional form of the incidence function stems from the fact that the number of effective contacts between infective individuals and susceptible individuals may saturate at high infective levels due to crowding of infective individuals or due to the preventive measures (and behavioral changes) taken by the susceptible individuals in response to the severity of the disease [23, 25–27].

A number of mathematical models have been developed in the literature to gain insights into the transmission dynamics of diseases with subpopulation (compartments) [4, 6–14, 29–31]. The choice of which compartment to include in a model depends on the characteristic of the particular disease being modelled and the purpose of the model [4].

The classical SIR or SIRS model assumed that the disease incubation is negligible so that each susceptible individual (in the class) once infected becomes infectious (and move to class) and later recovers (to move to class) where they acquire permanent or temporary immunity [32]. However, more general models than SIR or SIRS models assumed that susceptible individuals, once infected, first go through the latent period (in the class) before becoming infectious; the resulting models are of SEIR or SEIRS type, depending on whether recovered individuals acquired permanent or temporary immunity [33, 34].

The model considered in this study is based on SEIR (: susceptible, : exposed, : infected, : recovered) where recovered individuals acquire permanent immunity, so that they will not become infected again. This is owing to the fact that for many viral diseases such as measles, smallpox, rubella, HIV/AIDS, and influenza recovered individuals confer lifelong immunity [29, 33]. Huang and Takeuchi [29] study the classical SIR, SIS, SEIR, and SEI models with time delay and a general incidence function. Safi et. al. [13] consider the effect of periodic fluctuations on the transmission dynamics of a communicable disease using SEIRS model, subject to quarantine and isolation; the authors show that adding periodicity to the autonomous model does not alter the threshold dynamics of the model with respect to the control of the disease in the population. Zhang and Ma [35] considered the global dynamics of an SEIR model with saturating contact rate. Korobeinikov [36] established global asymptotic dynamics of SEIR and SIR models with several parallel infectious stages. Li et al. [33] analyzed the global dynamics of an SEIR model with vertical transmission and bilinear incidence. Li and Jin [34] consider global stability of an SEIR epidemic model with infectious force in latent, infected, and immune period. It should be stated, however, that the aforementioned three studies [33, 34, 36] considered mass action (bilinear) incidence to model the infection. This study focuses on the mathematical modeling of the transmission dynamics of an arbitrary disease with educated (counsel) and uneducated infectious stages.

This study has three important differences from those reported in [33, 34, 36]. The first is established essential qualitative features of an SEIR model with Holling type II incidence function. The second is that a public health education and counselling are offered to infected individuals (). Recent studies further reinforced the widely held belief that one key strategy for preventing and controlling the spread of communicable disease such as HIV (especially in resource-poor nations) is to provide HIV-related public health education and counselling (such as sexual education and awareness of the risk and life-threatening consequences of HIV/AIDS) which would, hopefully, lead to reduction in risky sexual behavior and safer lifestyle within the community [7]. In addition to the above extensions, rigorous qualitative analysis will be provided for the resulting SEIR model. In particular, the paper gives special emphasis on the global asymptotic stability of the disease-free equilibrium and endemic equilibrium.

The paper is organized as follows. The model is formulated and analysed for its basic dynamical features in Section 2. The stability analysis of the model is presented in Sections 3 and 4.

2. Model Formulation

The total population at time , denoted by , is subdivided into five compartments of susceptible , exposed (those who have been infected but are not yet infectious) , uneducated infected individuals , educated infected individuals , and recovered individuals, so that

The susceptible population is increased by the recruitment of individuals into the population (assumed susceptible), at a rate . Susceptible individuals may acquire infection, following effective contact with infected individuals (in the or class) at a rate , where

In (2), is the effective contact rate (contact capable of leading to infection), while the modification parameter accounts for the assumed reduction in disease transmission by educated infected individuals in comparison to uneducated infected individuals in the class. The population of susceptible individuals is further decreased by natural death (at a rate ). Thus, the rate of change of the susceptible population is given by

The population of exposed individuals is generated by the infection of susceptible individuals (at the rate ). This population is decreased by development of disease symptoms (at a rate ) and natural death (at a rate ), so that

The population of uneducated infected individuals is generated at the rate . It is decreased by natural recovery (at a rate ), education (at a rate ), natural death (at the rate ), and disease-induced death (at a rate ). This gives

The population of educated infected individuals is generated by the education of infected individuals (at the rate ). This population is decreased by recovery (at a rate ), natural death (at the rate ), and disease-induced death (at a rate ). It is assumed that the disease-induced mortality rate of educated infected individuals is low in comparison with uneducated infected individuals. Hence, the rate of change of this population is given by

Finally, the population of recovered individuals is generated by the recovery of uneducated and educated infected individuals (at the rates and , resp.). It is decreased by natural death (at the rate ), so that

Thus, the model for the transmission dynamics of an infectious disease in the presence of educated (counsel) infected individuals is given by the following nonlinear system of differential equations:

The model (8) extends some SEIR models in the literature such as those in [33, 34, 36] by

(i)replacing the mass action incidence function with Holling type II incidence function, (ii)splitting the compartment of infected individuals into educated (counsel) and uneducated (noncounsel) infected individuals thereby allowing time-varying infection rate. The epidemiological implication of this assumption is that educated infected individuals transmit the disease at a reduced rate () in comparison to the uneducated infected individuals (due to behavioral changes). For instance, In Zambia, the decline in HIV incidence since early 1990s is attributed to behavioral changes [37]. Public health education campaigns have also been successfully implemented in numerous countries and communities, such as Uganda, Thailand, Zambia, and the US gay community [38, 39]. Between 1991 and 1998, HIV prevalence dramatically declined in Uganda from 21 to 9.8% (with a corresponding reduction in nonregular sexual partners by 65% coupled with greater levels of awareness about HIV/AIDS [38]. The Ugandan programme fostered community mobilization towards change in risky behavior, without increasing stigma [40].

Further, unlike in the aforementioned modelling studies, detailed rigorous mathematical analysis of the model (8) will be provided.

2.1. Basic Properties

Since the model (8) monitors human populations, all its associated parameters are nonnegative. Further, the following nonnegativity result holds.

Theorem 1. The variables of the model (8) are nonnegative for all the time. In other words, solutions of the model system (8) with positive initial data will remain positive for all time .

Proof. Let . Thus, . It follows from the first equation of the system (8) that

which can be rewritten as

Hence,

so that

Similarly, it can be shown that , , and , for all time .

The previous result can also be established using the method in of [41, Appendix ].

We claim the following result.

Lemma 2. The closed set

is positively invariant.

Proof. Adding all the equations of the model (8) gives,

Since , it follows that if . Thus, a standard comparison theorem [42] can be used to show that . In particular, if . Thus, the region is positively invariant. Further, if , then either the solution enters infinite time or approaches asymptotically. Hence, the region attracts all solutions in .

Since the region is positively invariant, it is sufficient to consider the dynamics of the flow generated by the model (8) in , where the usual existence, uniqueness, continuation results hold for the system [10].

3. Local Stability of Disease-Free Equilibrium (DFE)

The DFE of the model (8) is given by

The local stability of will be explored using the next generation operator method [43, 44]. Using the notation in [44], the nonnegative matrix, , of the new infection terms and the -matrix, , of the transition terms associated with the model (8) are given, respectively, by

It follows that the control reproduction number [4, 21], denoted by , where is the spectral radius, is given by

where

Using in [44, Theorem 2], the following result is established.

Lemma 3. The DFE of the model (8), given by (15), is locally asymptotically stable (LAS) if and unstable if .

The quantity measures the average number of new infections generated by a single infected individual in a population. Lemma 3 implies that the disease can be eliminated from the community (when ) if the initial sizes of the subpopulations of the model are in the basin of attraction of the DFE (). To ensure that disease elimination is independent of the initial sizes of sub-populations, it is necessary to show that the DFE is globally asymptotically stable (GAS) if . This is explored below.

3.1. Global Stability of DFE

Theorem 4. The DFE of the model (8), given by (15), is GAS in whenever .

Proof. Consider the following Lyapunov function:

with Lyapunov derivative (where a dot represents differentiation with respect to time) given by

Since all the parameters and variables of the model (8) are nonnegative (Theorem 1), it follows that for with if and only if . Hence, is a Lyapunov function on . Therefore, the largest compact invariant subset of the set where is the singleton . Thus, it follows, by the LaSalle's invariance principle [18, 45], that

Since and (from (21)), it follows that, for sufficiently small , there exist constants and such that for all and for all . Hence, it follows from the last equation of the model (8) that, for ,

Thus, by comparison theorem [42],

so that, by letting ,

Similarly (by using and ), it can be shown that

Thus, it follows from (24) and (25) that

Hence,

Similarly, it can be shown that

Thus, by combining (21), (27), and (28), it follows that every solution of the equations of the model (8), with initial conditions in , approaches as (for ).

The previous result shows that the disease can be eliminated from the community if the associated reproduction number of the model is less than unity.

4. Existence and Stability for Endemic Equilibrium Point

In this section, the possible existence and stability of endemic (positive) equilibria of the model (8) (i.e., equilibria where at least one of the infected components of the model is nonzero) will be explored.

4.1. Existence of Endemic Equilibrium Point (EEP)

First of all, the persistence of the disease in the population will be investigated below. The model system (8) is said to be uniformly persistent if there exists a constant such that any solution satisfies ([31, 46])

provided that (the interior of the region ).

Theorem 5. System (8) is uniformly persistent in if and only if .

Proof. The theorem can be proved by using the approach used to prove Proposition 3.3 of [32], by applying a uniform persistence result in [46] and noting that the DFE of the model (8) is unstable whenever (Lemma 3).

When , it follows (from Theorem (3)) that model (8) is uniformly persistent; by using in [47, Theorem 2.8.6] and in [42, Theorem D.3] it follows that the model (8) has at least one endemic equilibrium in . Thus the following result is established.

Lemma 6. The model has at least one endemic equilibrium, given by , whenever .

The uniqueness of the endemic equilibrium will be investigated at the end of the next subsection.

4.2. Global Stability of Endemic Equilibrium for Special Case

Here, the global asymptotic stability property of the endemic equilibrium of the model (8) is given for a special case when educated (counsel) infected individuals do not transmit infection (). The model (8), with , then reduces to:

where

The reproduction number of the model (30) is given by

We claim the following result.

Theorem 7. The endemic equilibrium of the reduced model, given by (30), is GAS in if . Thus it is unique.

Proof. Consider the reduced model, given by (30). Let , so that the associated endemic equilibrium exists. Further, consider the following nonlinear Lyapunov function:

with Lyapunov derivative

It can be shown from (30) that, at endemic steady state,

Using the relations (35) in (34) gives

Adding and subtracting and in (36) gives

Rewriting the last four terms of (37) in terms of (by using the relations in (35)) gives

Since is an increasing function and , it follows that

finally, since the arithmetic mean exceeds the geometric mean, then

Further, since all the model parameters are nonnegative, it follows that for with if and only if , . Hence, is a Lyapunov function on . Therefore, the largest compact invariant subset of the set where is the singleton . Thus, it follows, by the laSalle's invariance Principle [18, 45], that , and as . Since , it follows that, for sufficiently small , there exist constants such that for all . It follows from the fourth equation of (30) that, for ,

Thus, by comparison theorem [42],

Hence, by letting , we have

Similarly, by using , it can be shown that

Thus, it follows from (44) and (45) that

Hence, . In a similar way, it can be shown that . Thus, every solution to the equations of the reduced model, with initial condition in , approaches the endemic equilibrium of the reduced system (30) as for . This shows that the endemic equilibrium is unique.

5. Conclusions

In this paper, an SEIR epidemic model with a nonlinear (Holling type II) incidence rate is designed and analysed. Some of the theoretical and epidemiological findings of the study are as follows.

(i)The model (8) has a locally stable disease-free equilibrium whenever the associated reproduction number is less than unity. (ii)The DFE of the model (8) is shown to be globally asymptotically stable when . (iii)The model (8) is uniformly persistent in if and only if . (iv)The endemic equilibrium of the reduced model (8), with , is shown to be globally asymptotically stable, when it exists.

Acknowledgments

S. Garba acknowledges, with thanks, the support on part of the University of Pretoria Research and Development Project (RDP). The authors are grateful to the two anonymous reviewers and the handling editor for their constructive comments, which have enhanced the paper.

References

R. M. Anderson and R. M. May, Eds., Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, New York, NY, USA, 1991.

E. Beretta and D. Breda, “An SEIR epidemic model with constant latency time and infectious period,” Mathematical Biosciences and Engineering, vol. 8, no. 4, pp. 931–952, 2011.

K. Dietz, “Transmission and control of arbovirus disease,” in Epidemiology, K. L. Cooke, Ed., p. 104, SIAM, Philadelphia, Pa, USA, 1975.

H. W. Hethcote, “Mathematics of infectious diseases,” SIAM Review, vol. 42, no. 4, pp. 599–653, 2000.

C. Rebelo, A. Margheri, and N. Bacaër, “Persistence in seasonally forced epidemiological models,” Journal of Mathematical Biology, vol. 64, no. 6, pp. 933–949, 2012.

Z. Feng and J. X. Velasco-Hernández, “Competitive exclusion in a vector-host model for the dengue fever,” Journal of Mathematical Biology, vol. 35, no. 5, pp. 523–544, 1997.

S. M. Garba and A. B. Gumel, “Mathematical recipe for HIV elimination in Nigeria,” Journal of the Nigerian Mathematical Society, vol. 29, pp. 1–66, 2010.

S. M. Garba, A. B. Gumel, and J. M. S. Lubuma, “Dynamically-consistent non-standard finite difference method for an epidemic model,” Mathematical and Computer Modelling, vol. 53, no. 1-2, pp. 131–150, 2011.

A. B. Gumel, “Global dynamics of a two-strain avian influenza model,” International Journal of Computer Mathematics, vol. 86, no. 1, pp. 85–108, 2009.

H. W. Hethcote and H. R. Thieme, “Stability of the endemic equilibrium in epidemic models with subpopulations,” Mathematical Biosciences, vol. 75, no. 2, pp. 205–227, 1985.

M. A. Safi and A. B. Gumel, “Mathematical analysis of a disease transmission model with quarantine, isolation and an imperfect vaccine,” Computers and Mathematics with Applications, vol. 61, no. 10, pp. 3044–3070, 2011.

M. A. Safi and A. B. Gumel, “Qualitative study of a quarantine/isolation model with multiple disease stages,” Applied Mathematics and Computation, vol. 218, no. 5, pp. 1941–1961, 2011.

M. A. Safi, M. Imran, and A. B. Gumel, “Threshold dynamics of a non-autonomous SEIRS model with quarantine and isolation,” Theory in Biosciences, vol. 131, no. 1, pp. 19–30, 2012.

O. Sharomi, C. N. Podder, A. B. Gumel, E. H. Elbasha, and J. Watmough, “Role of incidence function in vaccine-induced backward bifurcation in some HIV models,” Mathematical Biosciences, vol. 210, no. 2, pp. 436–463, 2007.

T. Zhang and Z. Teng, “On a nonautonomous SEIRS model in epidemiology,” Bulletin of Mathematical Biology, vol. 69, no. 8, pp. 2537–2559, 2007.

L. D. Wang and J. Q. Li, “Global stability of an epidemic model with nonlinear incidence rate and differential infectivity,” Applied Mathematics and Computation, vol. 161, no. 3, pp. 769–778, 2005.

C. Vargas-De-León, “On the global stability of SIS, SIR and SIRS epidemic models with standard incidence,” Chaos, Solitons & Fractals, vol. 44, no. 12, pp. 1106–1110, 2011.

J. K. Hale, Ordinary di Erential Equations, John Wiley & Sons, New York, NY, USA, 1969.

J. C. Kamgang and G. Sallet, “Computation of threshold conditions for epidemiological models and global stability of the disease-free equilibrium (DFE),” Mathematical Biosciences, vol. 213, no. 1, pp. 1–12, 2008.

M. Qiao, A. Liu, and U. Foryś, “Qualitative analysis of the SICR epidemic model with impulsive vaccinations,” Mathematical Methods in the Applied Sciences. In press

R. M. Anderson and R. M. Verlag, Population Biology of Infectious Diseases, Springer, New York, NY, USA, 1982.

V. Capasso and G. Serio, “A generalization of the Kermack-McKendrick deterministic epidemic model,” Mathematical Biosciences, vol. 42, no. 1-2, pp. 43–61, 1978.

J. Hou and Z. Teng, “Continuous and impulsive vaccination of SEIR epidemic models with saturation incidence rates,” Mathematics and Computers in Simulation, vol. 79, no. 10, pp. 3038–3054, 2009.

J. Zu, W. Wang, and B. Zu, “Evolutionary dynamics of prey-predator systems with holling type II,” Mathematical Biosciences and Engineering, vol. 4, no. 2, pp. 221–237, 2007.

W. M. Liu, S. A. Levin, and Y. Iwasa, “Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models,” Journal of Mathematical Biology, vol. 23, no. 2, pp. 187–204, 1985.

S. Ruan and W. Wang, “Dynamical behavior of an epidemic model with a nonlinear incidence rate,” Journal of Differential Equations, vol. 188, no. 1, pp. 135–163, 2003.

M. A. Safi and A. B. Gumel, “The effect of incidence functions on the dynamics of a quarantine/isolation model with time delay,” Nonlinear Analysis: Real World Applications, vol. 12, no. 1, pp. 215–235, 2011.

J. Zu, W. Wang, and B. Zu, “Evolutionary dynamics of prey-predator systems with Holling type II functional response,” Mathematical Biosciences and Engineering, vol. 4, no. 4, article 755, 2007.

G. Huang and Y. Takeuchi, “Global analysis on delay epidemiological dynamic models with nonlinear incidence,” Journal of Mathematical Biology, vol. 63, no. 1, pp. 125–139, 2011.

H. R. Thieme, “Global asymptotic stability in epidemic models,” in Equadiff, H. W. Knobloch and K. Schmidt, Eds., vol. 1017 of Lecture notes in Mathematics, pp. 608–615, Springer, Heidelberg, Germany, 1983.

H. R. Thieme, “Epidemic and demographic interaction in the spread of potentially fatal diseases in growing populations,” Mathematical Biosciences, vol. 111, no. 1, pp. 99–130, 1992.

M. Y. Li, J. R. Graef, L. Wang, and J. Karsai, “Global dynamics of a SEIR model with varying total population size,” Mathematical Biosciences, vol. 160, no. 2, pp. 191–213, 1999.

M. Y. Li, H. L. Smith, and L. Wang, “Global dynamics of an seir epidemic model with vertical transmission,” SIAM Journal on Applied Mathematics, vol. 62, no. 1, pp. 58–69, 2001.

G. Li and Z. Jin, “Global stability of a SEIR epidemic model with infectious force in latent, infected and immune period,” Chaos, Solitons and Fractals, vol. 25, no. 5, pp. 1177–1184, 2005.

J. Zhang and Z. Ma, “Global dynamics of an SEIR epidemic model with saturating contact rate,” Mathematical Biosciences, vol. 185, no. 1, pp. 15–32, 2003.

A. Korobeinikov, “Global properties of SIR and SEIR epidemic models with multiple parallel infectious stages,” Bulletin of Mathematical Biology, vol. 71, no. 1, pp. 75–83, 2009.

K. Fylkesnes, R. M. Musonda, M. Sichone, Z. Ndhlovu, F. Tembo, and M. Monze, “Declining HIV prevalence and risk behaviours in Zambia: evidence from surveillance and population-based surveys,” AIDS, vol. 15, no. 7, pp. 907–916, 2001.

D. Low-Beer and R. L. Stoneburner, “Behaviour and communication change in reducing HIV: is Uganda unique?” African Journal of AIDS Research, vol. 2, no. 1, pp. 9–21, 2003.

N. Hussaini, M. Winter, and A. B. Gumel, “Qualitative assessment of the role of public health education program on HIV transmission dynamics,” Mathematical Medicine and Biology, vol. 28, no. 3, pp. 245–270, 2011.

E. C. Green, D. T. Halperin, V. Nantulya, and J. A. Hogle, “Uganda's HIV prevention success: the role of sexual behavior change and the national response,” AIDS and Behavior, vol. 10, no. 4, pp. 335–346, 2006.

H. R. Thieme, Mathematics in Population, Princeton University Press, 2003.

H. L. Smith and P. Waltman, The Theory of the Chemostat, Cambridge University Press, 1995.

O. Diekmann, J. A. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.

P. Van Den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.

J. P. LaSalle, The Stability of Dynamical Systems, CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1976.

H. I. Freedman, S. Ruan, and M. Tang, “Uniform persistence and flows near a closed positively invariant set,” Journal of Dynamics and Differential Equations, vol. 6, no. 4, pp. 583–600, 1994.

N. P. Bhatia and G. P. Szeg, Dynamical Systems: Stability Theory and Applications, vol. 35 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 1967.