<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2013.410A2011</article-id><article-id pub-id-type="publisher-id">AM-38059</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Global Stability Analysis of a Delayed SEIQR Epidemic Model with Quarantine and Latent
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>iantian</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yakui</surname><given-names>Xue</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mathematics, North University of China, Taiyuan, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>xyk5152@163.com(YX)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>30</day><month>09</month><year>2013</year></pub-date><volume>04</volume><issue>10</issue><fpage>109</fpage><lpage>117</lpage><history><date date-type="received"><day>June</day>	<month>12,</month>	<year>2013</year></date><date date-type="rev-recd"><day>July</day>	<month>12,</month>	<year>2013</year>	</date><date date-type="accepted"><day>July</day>	<month>19,</month>	<year>2013</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
     
   In this paper, we study a kind of the delayed SEIQR infectious disease model with the quarantine and latent, and get the threshold value which determines the global dynamics and the outcome of the disease. The model has a disease-free equilibrium which is unstable when the basic reproduction number is greater than unity. At the same time, it has a unique endemic equilibrium when the basic reproduction number is greater than unity. According to the mathematical dynamics analysis, we show that disease-free equilibrium and endemic equilibrium are locally asymptotically stable by using Hurwitz criterion and they are globally asymptotically stable by using suitable Lyapunov functions for any <inline-formula><inline-graphic xlink:href="dit_284ee28f-a80d-4a2a-ac8f-72b348be4d8e.png" xlink:type="simple"/></inline-formula> Besides, the SEIQR model with nonlinear incidence rate is studied, and the <inline-formula><inline-graphic xlink:href="dit_87230881-e679-44b3-9cfa-accae72773ee.png" xlink:type="simple"/></inline-formula> that the basic reproduction number is a unity can be found out. Finally, numerical simulations are performed to illustrate and verify the conclusions that will be useful for us to control the spread of infectious diseases. Meanwhile, the <inline-formula><inline-graphic xlink:href="dit_87690537-83a9-4dc1-a0bc-f7297e666e78.png" xlink:type="simple"/></inline-formula> will effect changing trends of <inline-formula><inline-graphic xlink:href="dit_8de48f15-bf27-42a6-8a3d-72b187d46c78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="dit_91f47ffe-e8db-40cf-b541-6e08b379ce30.png" xlink:type="simple"/></inline-formula> in system (1), which is obvious in simulations. Here, we take <inline-formula><inline-graphic xlink:href="dit_c2f2d155-7e4a-40f2-bea4-2743cda6d27a.png" xlink:type="simple"/></inline-formula> as an example to explain that.  
   
    
 
</p></abstract><kwd-group><kwd>SEIQR Model; Lyapunov Function; Delay; Global Stability; Nonlinear Incidence Rate; Simulations</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Many people have been paying attention to the study of some epidemics, and have accumulated a lot of experience. By establishing reasonable mathematical models, they put forward the measures which controlled the spread of epidemics effectively. And many scholars researched specific diseases and considered the diseases with incubation period, recovery time, quarantine and so on [1-6]. So many epidemics were controlled. Generally speaking, when epidemics spread, there are many kinds of delays, which include immunity period delay [7-9], infectious period delay, incubation period delay. In [<xref ref-type="bibr" rid="scirp.38059-ref10">10</xref>], Enatsu et al. studied stability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates, at the same time, they proved disease-free equilibrium was globally asymptotically stable and endemic equilibrium was permanent under certain conditions. At the same time, global stability of an SIR (where S, I, R denote the number of susceptible individuals, infectious individuals, recovery individuals) epidemic model with constant infectious period was studied by Zhang et al.</p><p>[<xref ref-type="bibr" rid="scirp.38059-ref11">11</xref>], they showed the endemic equilibrium was globally asymptotically stable with appropriate Lyapunov functions. And in [<xref ref-type="bibr" rid="scirp.38059-ref12">12</xref>], Gao et al. discussed pulse vaccination of an SEIR (E denote the number of exposed individuals) epidemic model with delay and bilinear incidence. Meanwhile, impulsive vaccination of SEIR epidemic model with time delay and nonlinear incidence rate was researched by Zhao et al. [<xref ref-type="bibr" rid="scirp.38059-ref13">13</xref>], and showed the pulse system that was similar to the pulse system with bilinear incidence rate. Besides, on the basis of [<xref ref-type="bibr" rid="scirp.38059-ref13">13</xref>], Xu and Ma introduced the saturated incidence rate. Meanwhile, they showed disease-free equilibrium and endemic equilibrium were globally asymptotically stable under certain condition in [<xref ref-type="bibr" rid="scirp.38059-ref14">14</xref>]. However, in addition to the bilinear incidence rate, nonlinear incidence rate and saturated incidence rate, there were some scholars who studied the non-monotone incidence. For example, an SIRS epidemic model with pulse vaccination and non-monotonic incidence rate was discussed by Zhang et al. [<xref ref-type="bibr" rid="scirp.38059-ref15">15</xref>], and they proved the disease-free equilibrium and endemic equilibrium were asymptotically stable under certain conditions. Besides, some scholars studied a delayed SEIQR (Q denote the number of quarantined individuals) epidemic model with pulse vaccination and the quarantine measure, and they showed that the disease-free equilibrium of the system was globally attractive and endemic equilibrium was permanent under certain conditions. In this paper, we study a delayed SEIQR epidemic model without pulse on the basis of [14,16].</p><p>The organization of this paper is as follows: In Section 2, SIQR epidemic model and its basic reproduction number and existence of equilibrium are given. In Section 3, the local stability of endemic equilibrium and disease-free equilibrium is showed by using Hurwitz criterion. By using suitable Lyapunov functions and LaSalle’s invariance principle, we prove the disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than unity and the endemic equilibrium is globally asymptotically stable when the basic reproduction number is greater than unity. At the same time, the system with the nonlinear incidence rate is discussed in Section 3. In Section 4, presents the numerical simulations of the system followed by a conclusion in Section 3. At last, a brief discussion is given in Section 5 to conclude this work.</p></sec><sec id="s2"><title>2. Establishment of the Model</title><p>We establish the following SEIQR epidemic model, Here <img src="11-7401649\ccfe691a-4cc2-47ff-9483-a5ac6520a417.jpg" /> represents the number of individuals who are susceptible to disease, that is, who are not yet infected at time t. <img src="11-7401649\d94c885c-2c56-4721-a5e1-56b07b6bd7cc.jpg" />is the number of individuals who are infected but hardly infectious. So we think they can’t infect other people, but they need to be quarantined. <img src="11-7401649\0d314b84-7fe1-4992-87b3-d899d59f8e97.jpg" />represents the number of infected individuals who are infectious and are able to spread the disease by contact with susceptible individuals. <img src="11-7401649\c76f2ef6-d7f2-4431-96e3-3864c8fad303.jpg" />is the number of infectious individuals who are quarantined at time t. <img src="11-7401649\a9bbb720-2148-435a-97ca-d513c56689e0.jpg" />represents the number of recovered individuals at time t.</p><disp-formula id="scirp.38059-formula24380"><label>(1)</label><graphic position="anchor" xlink:href="11-7401649\ad7e7758-23a0-413c-af27-efe0c03562bc.jpg"  xlink:type="simple"/></disp-formula><p>The initial conditions for system (1) are</p><p><img src="11-7401649\2458d94e-7aef-46ca-8145-ee06e5f2d6a6.jpg" /></p><p>And the feasible region of the model with the initial conditions above is</p><p><img src="11-7401649\53d7dddb-bbc8-4ef2-9d95-b10c68ddddfd.jpg" /></p><p>Here, we presume that</p><p><img src="11-7401649\e4ed4d09-d5b3-4ef0-8a89-00f796a3bc4a.jpg" /></p><p>It is easy to show that <img src="11-7401649\fe821662-0f8d-4a96-a437-944ad85ab788.jpg" /> is positively invariant with respect to system (1).</p><p>Where all the parameters are positive constants, <img src="11-7401649\17bfb228-b34d-4b4f-bd19-44caeba8ef63.jpg" />is the recruitment rate of the susceptible population, <img src="11-7401649\26ef8a7f-083a-4c59-87a7-5fad5942d6ea.jpg" />, <img src="11-7401649\3eba8089-f88c-459d-ae65-fb57336ec375.jpg" /><img src="11-7401649\2f518b5c-3aba-4abd-8d53-d098419959de.jpg" /><img src="11-7401649\8f91334b-f531-47dc-b3e8-4abc853dbaa6.jpg" /><img src="11-7401649\e548fb04-2bf6-4068-914c-6f052f03b7b4.jpg" />are the natural death rate of the susceptible, exposed, infectious, quarantine and recovered respectively, <img src="11-7401649\78b90c45-db9f-45e1-8e61-869e7e94e57b.jpg" />is the disease transmission coefficient, <img src="11-7401649\07732207-8e3c-4535-b44a-c9a470a2b574.jpg" />is the death rate due to disease without quarantine, <img src="11-7401649\d0fb3168-3726-4bc0-a7af-15fa68d7def2.jpg" />is the death rate due to disease after quarantine, <img src="11-7401649\93aa6644-7641-4e53-9227-2716ece4fd8f.jpg" />is the recovery rate after quarantine, <img src="11-7401649\5d91f490-e4c7-4357-b79e-b4cf8fda3604.jpg" />is the recovery rate without quarantine, <img src="11-7401649\6a2c078c-2551-4432-851f-87447ee1b8ce.jpg" />, <img src="11-7401649\8f2a2b13-0da9-4d0c-8b7d-f99a7d01162f.jpg" />are quarantine rate of<img src="11-7401649\94edc719-19a6-42a0-b10c-887c55391fcb.jpg" />, <img src="11-7401649\8b4b432d-1639-4422-903c-985c54b006c3.jpg" />respectively, <img src="11-7401649\8cae9df3-27f0-40f2-af78-271ae48f23de.jpg" />is the recovery rate of <img src="11-7401649\e418ec8b-6770-4b12-85f0-be1068fee4bd.jpg" /> and <img src="11-7401649\43156ac4-36ec-4b9a-90b8-0439f01b6ab4.jpg" /> is the latent period of the epidemic.</p><p>Because the variables R and Q do not appear in the first three equations in system (1), we further simplify system (1) and then obtain the following model</p><disp-formula id="scirp.38059-formula24381"><label>(2)</label><graphic position="anchor" xlink:href="11-7401649\fec55a2a-ac15-4fcb-ab41-7bbfedb6b0e5.jpg"  xlink:type="simple"/></disp-formula><p>In this paper, we are concerned with system (2).</p><p>The initial conditions for system (2) are</p><p><img src="11-7401649\05167d4d-cdc4-480f-bc0a-26c71f665cde.jpg" /></p><p>And the feasible region of the model with the initial condition above is</p><p><img src="11-7401649\da85236b-2505-4fbe-9a63-a6565634b6bb.jpg" /></p><p>Here, we presume that</p><p><img src="11-7401649\2e9d93e6-6ea9-4944-92a9-3d40f46101a6.jpg" /></p><p>It is easy to show that <img src="11-7401649\8864edef-45d0-4e22-973e-fe37ffaa8dd7.jpg" /> is positively invariant with respect to system (2).</p><p>According to the practical significance of the epidemic model, system (2) always has a disease-free equilibrium</p><p><img src="11-7401649\063241da-79a1-409f-a7de-670201b88c30.jpg" /></p><p>Denote the basic reproduction number of system (2)</p><p><img src="11-7401649\d30f0528-3c04-4624-b427-ff4e64eee559.jpg" /></p><p>Define <img src="11-7401649\ad5efcd6-10e7-40fc-84d2-1ca9a9728d9c.jpg" /> If the basic reproductive number <img src="11-7401649\44a7f47f-531c-41a7-8aba-5f6eeb6b5249.jpg" /> system (2) has an unique endemic equilibrium</p><p><img src="11-7401649\6c0a03d4-d9d9-45e3-b400-ed62222cf1a6.jpg" /></p></sec><sec id="s3"><title>3. The Stability of Equilibrium</title><p>In this section, we discuss the local stability of endemic equilibrium and disease-free equilibrium of system (2) by analyzing the corresponding characteristic equations respectively. By defining reasonable Lyapunov functions, we resolve the global dynamics of equilibriums without requiring any extra conditions. In addition, system (2) with nonlinear incidence is studied.</p><sec id="s3_1"><title>3.1. Stability of Disease-Free Equilibrium</title><p>Theorem 3.1.1. If<img src="11-7401649\1d9631c3-771f-4796-94b6-88c35cbd9aa7.jpg" />, the disease-free equilibrium <img src="11-7401649\46ad65d1-18f5-4410-9f7f-bc5d14859511.jpg" /> of system (2) is locally asymptotically stable for any <img src="11-7401649\65f6cbfb-4381-4859-9f37-576270e07e4d.jpg" /> in<img src="11-7401649\3cf00690-e194-464b-af51-ad89248effea.jpg" />. If<img src="11-7401649\1bbfa61b-2f1b-4da7-bc63-139d22eb6eb3.jpg" />, it is unstable for any <img src="11-7401649\4b7e373f-bc9c-4d5d-9055-354f7da21b86.jpg" /> in<img src="11-7401649\556390a6-9ed0-4839-b802-c64d1718162a.jpg" />.</p><p>Proof. The characteristic matrix at the disease-free equilibrium <img src="11-7401649\8831313f-2200-40d5-a7ca-2890162bd13a.jpg" /></p><p><img src="11-7401649\b965c1b8-dd84-4d0c-b076-f37feb907dc0.jpg" /></p><p>When <img src="11-7401649\39b937d2-94b3-4e79-961d-b9094de31119.jpg" /> the characteristic equation at the disease-free equilibrium <img src="11-7401649\41e2af59-4296-4b1b-9f6d-a4a74f787214.jpg" /> of system (2) takes the form</p><p><img src="11-7401649\d352463a-8823-4a67-b530-7df86e6c9b62.jpg" /></p><p>Clearly, system (2) always has two negative real roots</p><p><img src="11-7401649\dfa3bc1a-c5bf-41d5-aac5-c94c2f9cdb66.jpg" /></p><p>All other roots are given by the roots of equation</p><p><img src="11-7401649\8f7b0393-d677-460f-a15e-0f4d94da1a18.jpg" /></p><p><img src="11-7401649\055ba22d-8154-475b-9bb5-bc4cad2551dc.jpg" /></p><p>Assume <img src="11-7401649\efbf02c0-72f9-4401-82ad-44e7b66d0f11.jpg" /> <img src="11-7401649\36bd2e9c-97e0-4bd9-b8b3-19de275b3572.jpg" /></p><p>That is, <img src="11-7401649\9447e515-bf13-4d08-a84a-d5099aa63768.jpg" /></p><p>Because <img src="11-7401649\665f0c5a-7758-4b78-9b96-75ac49efcee4.jpg" /> <img src="11-7401649\c1d94a0e-1c3e-4815-9b9a-d742c8a40c33.jpg" /> which is contradictory. So <img src="11-7401649\5b2c7cb1-a83c-4737-81de-4017a1711a40.jpg" /> Therefore the disease-free equilibrium <img src="11-7401649\1aa595ce-56b3-493e-8fef-7f6c5c6a71bd.jpg" /> of system (2) is locally asymptotically stable.</p><p>If <img src="11-7401649\1b9ed7b2-6dab-4576-b785-f49eb6662a0c.jpg" /> let</p><p><img src="11-7401649\0bfe04c4-1165-4241-bd8d-ed7a8d498451.jpg" /></p><p>so there is a positive real root at least. The disease-free equilibrium <img src="11-7401649\9e20560f-ad73-4ff7-82be-e5f0b5cda7ff.jpg" /> of system (2) is unstable.</p><p>When <img src="11-7401649\c0c6503e-4cc4-4a02-a9e5-4765c2aa7d90.jpg" /> it is easy for us to prove the disease-free equilibrium <img src="11-7401649\698dce3c-a9c3-458d-8371-b8fedc89a12a.jpg" /> of system (2) is locally asymptotically stable.</p><p>Theorem 3.1.2. If<img src="11-7401649\d52c513e-6c46-4833-ab69-ec02597e7133.jpg" />, the disease-free equilibrium <img src="11-7401649\243aad63-bbee-46be-86bd-ca3a808dafb9.jpg" /> of system (2) is globally asymptotically stable for any <img src="11-7401649\f8c10384-5967-44d2-b951-54a4f8e1d169.jpg" /> in<img src="11-7401649\f0697541-bf5a-409f-ab26-36878ecd2adf.jpg" />.</p><p>Proof. For <img src="11-7401649\e9f24e84-505f-458a-8fc1-0084b24c4eb7.jpg" /> define a differentiable Lyapunov function</p><p><img src="11-7401649\12229c60-51c8-443b-993e-f1df19fc8dc5.jpg" /></p><p><img src="11-7401649\0d288897-cbde-4f1e-a9fe-4e3ad97eb1ba.jpg" /></p><p>Obviously, <img src="11-7401649\cf62ecbd-f03c-48cd-b072-e18e2ca69f9b.jpg" /></p><p>Calculating the derivative of <img src="11-7401649\1408ada6-a1f4-4cab-a623-c29d4c90d84e.jpg" /> along positive solutions of system (2), it follows that</p><p><img src="11-7401649\fd8b42c2-a73c-4f10-a8c1-1d7b11249b1c.jpg" /></p><p>According to the feasible region, <img src="11-7401649\900b62c0-a3c9-469d-99a0-93d199398f2a.jpg" /></p><p>So</p><p><img src="11-7401649\024628a8-6d33-4f7a-86f3-cef85e633191.jpg" /></p><p>That is,</p><p><img src="11-7401649\854a0cea-6219-45cb-81c8-9a691883251f.jpg" /></p><p>And when <img src="11-7401649\a99438ab-a848-42d5-b931-cfd11ac89fb4.jpg" /></p><p>While <img src="11-7401649\d61cde4d-0f59-49dd-9848-e20b3da2a575.jpg" /> if and only if,</p><p><img src="11-7401649\24770707-2e86-49b1-b04b-caf2386a3a82.jpg" /></p><p>For all t, it is easy to show that <img src="11-7401649\730919a8-cd3e-4ea0-af81-f8b4ccc608e1.jpg" /> is the largest invariant subset of the set <img src="11-7401649\55cc4fa0-8296-4f45-9f7a-cf6218beef2b.jpg" /> Because of LaSalle’s invariance principle, disease-free equilibrium <img src="11-7401649\5223b724-0ced-4c4b-99ce-2bd925b798b6.jpg" /> of system (2) is globally asymptotically stable. This completes the proof.</p></sec><sec id="s3_2"><title>3.2. The Stability of Endemic Equilibrium</title><p>Theorem 3.2.1. For any<img src="11-7401649\6d8c52ff-3740-4e0e-9282-6c1366cee3fb.jpg" />, if <img src="11-7401649\3eac56f6-7127-470d-bdeb-00eb2b36621b.jpg" /> the endemic equilibrium <img src="11-7401649\8f1a01aa-f39a-42cc-bfeb-d25d4e628646.jpg" /> of system (2) is locally asymptotically stable in <img src="11-7401649\4070ebb1-d8a0-4ca9-8460-bf543f788235.jpg" /></p><p>Proof. The characteristic matrix at the endemic equilibrium <img src="11-7401649\6b12f57c-3cc2-4f82-a80d-d219d2c7269f.jpg" /></p><p><img src="11-7401649\00a02116-52e0-471f-adb8-bca6746e3e6b.jpg" /></p><p>Order <img src="11-7401649\c12d4fbd-33d2-4bf3-be3b-be66913bb9f5.jpg" /></p><p>The characteristic equation at the endemic equilibrium <img src="11-7401649\93b1ac5c-049d-494a-b044-5db828ca2f22.jpg" /> is</p><p><img src="11-7401649\d3f543c9-08f5-4746-b520-fcdaee5581fb.jpg" /></p><p>Clearly, system (2) always has a negative real root</p><p><img src="11-7401649\9d51092d-ae45-4b35-bdd6-971df07ad446.jpg" /></p><p>When <img src="11-7401649\c45fd59b-cc61-42ce-91ef-39315e095e0a.jpg" /> all other roots are given by the roots of equation</p><p><img src="11-7401649\4fbbbb51-4bb3-4080-bfeb-696473cb050a.jpg" /></p><p><img src="11-7401649\f13d63c8-24cf-4e93-a5ff-895cf6d767b2.jpg" />so</p><p><img src="11-7401649\cb30b368-f402-40e2-a6ef-5067be79110c.jpg" /></p><p>So according to Hurwitz criterion, the endemic equilibrium <img src="11-7401649\5870dbe6-4150-422c-9f93-56d397e9a17b.jpg" /> of system (2) is locally asymptotically stable.</p><p>When <img src="11-7401649\9a9f68d4-477c-404a-ad47-8126c7a2e979.jpg" /> all other roots are given by the roots of equation</p><p><img src="11-7401649\edea6fc4-64d9-4e06-9303-e8226e94f5fa.jpg" /></p><p>Simplify, we can get</p><p><img src="11-7401649\5f0d61a7-62b6-44f5-aae3-b72464fb11df.jpg" /></p><p>Let</p><p><img src="11-7401649\42f114fd-87de-4ee6-8747-b4e621b22e48.jpg" /></p><p>Then</p><disp-formula id="scirp.38059-formula24382"><label>(3)</label><graphic position="anchor" xlink:href="11-7401649\e0f10f37-9a7a-4567-99a4-7ad5e1687abc.jpg"  xlink:type="simple"/></disp-formula><p>Let <img src="11-7401649\7d62429a-e0d2-4096-87c7-64e2c980fae5.jpg" /> is the root of Equation (3), on substituting to <img src="11-7401649\ec83f6e7-6083-4702-9a43-de9bdffe7890.jpg" /> Equation (3), we derive that</p><p><img src="11-7401649\f46b8b46-f15d-4a0e-9222-9f332a4c791a.jpg" /></p><p>Separating real and imaginary parts, it follows that</p><p><img src="11-7401649\d33e085d-dda8-4bae-b586-08aa23ca41a8.jpg" /></p><p><img src="11-7401649\71392714-5657-43c6-a688-70de1b20af2c.jpg" /></p><p>Then we can get</p><p><img src="11-7401649\a2678c9e-dc64-41cc-ae14-7063a08d6036.jpg" /></p><p>Order</p><p><img src="11-7401649\62496d6a-1958-4926-aca8-91453a500018.jpg" /></p><p>Letting <img src="11-7401649\a8105de3-930f-4812-9079-be47e0aeecbb.jpg" /> then Equation (4) becomes</p><p><img src="11-7401649\b72aff0c-9ced-4876-abc7-71edf5dba373.jpg" /></p><p>Here</p><p><img src="11-7401649\05ab6826-aae7-4b90-8465-692c4b7d10d4.jpg" /></p><p>Application of the conclusions of [<xref ref-type="bibr" rid="scirp.38059-ref17">17</xref>], we can know that positive <img src="11-7401649\2278f0a2-5ce9-4fe5-8573-3f6b8d7eb739.jpg" /> doesn’t exist. Hence <img src="11-7401649\0714fe84-7546-4294-a2a0-d53470255563.jpg" /> also doesn’t exist. There are not pure imaginary roots in system (2). Therefore all the roots have negative real component. So endemic equilibrium <img src="11-7401649\5caa6b00-6b6c-49b1-a128-59d9324ce6e8.jpg" /> of system (2) is locally asymptotically stable.</p><p>Theorem 3.2.2. If <img src="11-7401649\bd5173ff-0d88-491e-8f02-d1082491e851.jpg" /> when <img src="11-7401649\c6e3a437-4c54-45ec-9470-3ff1f9f10022.jpg" /> the endemic equilibrium <img src="11-7401649\f5785905-8c82-44af-8515-c4673acca70c.jpg" /> of system (2) is globally asymptotically stable in <img src="11-7401649\635f9230-fb4e-465a-92d9-71504a7e95ee.jpg" /></p><p>Proof. Define a differentiable Lyapunov function</p><p><img src="11-7401649\b3d9d2ef-c723-49c3-8349-b838e105bc10.jpg" /></p><p><img src="11-7401649\bd070d44-c71e-481e-8240-ae12918f7c41.jpg" /></p><p><img src="11-7401649\23d8b412-2f32-4f7b-ba19-b016b1163dd5.jpg" />both of them are real numbers. The function is positive definite. Calculating the derivative of <img src="11-7401649\af321b23-6385-4516-a266-5b5f108d677a.jpg" /> along positive solutions of system (2), it follows that</p><p><img src="11-7401649\51ca57fa-2dfa-4c7e-a657-377212c436bd.jpg" /></p><p>On substituting <img src="11-7401649\5c46d583-1d20-4d14-aa09-29dad8d64f4b.jpg" /> we have</p><p><img src="11-7401649\12285d15-4f63-4a20-aad0-0a9a3cc97a13.jpg" /></p><p>Let <img src="11-7401649\9044f7f4-0a69-4198-8044-4bd4605a3c4d.jpg" /></p><p>So <img src="11-7401649\f22f8690-ea3a-43ee-a95c-81666faccceb.jpg" /> In addition, when <img src="11-7401649\795ff82c-e643-4d03-b3b6-c80494c2dc79.jpg" /> if and only if <img src="11-7401649\cf116f9f-6402-4272-ae81-36d63c8cc9e7.jpg" /></p><p>It is easy to show that <img src="11-7401649\1abb1082-5459-4d2b-9cf5-3849d262f11b.jpg" /> is the largest invariant subset of the set <img src="11-7401649\f2b05c77-0989-47f5-be99-bd91ce90de24.jpg" /> Because of LaSalle’s invariance principle, the endemic equilibrium <img src="11-7401649\e214771d-bd7f-41cd-9046-12f0ab4a3c8c.jpg" /> of system (2) is globally asymptotically stable when<img src="11-7401649\c1440326-df9e-4e2a-8a60-a7a90890c7be.jpg" />. This completes the proof.</p><p>Theorem 3.2.3. If <img src="11-7401649\2d4626fb-6843-4ca5-b05a-87881f2b66d9.jpg" /> when <img src="11-7401649\2b8cbeb0-e7d0-444c-bd67-d89aba1a4eee.jpg" /> the endemic equilibrium <img src="11-7401649\7bd3d706-836c-490e-9864-729d932f83cc.jpg" /> of system (2) is globally asymptotically stable in <img src="11-7401649\cd1875ee-7b93-4656-ace8-4c896c0258c0.jpg" /></p><p>Proof. For <img src="11-7401649\d72e7194-b6f8-44a6-a12b-6bb8eca5da34.jpg" /> define a differentiable Lyapunov function</p><p><img src="11-7401649\47482c8d-69d1-478e-b32d-7dd4a9b0ca8d.jpg" /></p><p>Order</p><p><img src="11-7401649\a94f981c-6839-439d-b221-b6e44044e7d3.jpg" /></p><p><img src="11-7401649\b9fdb3f7-f7a1-42b6-97ff-b23a4f14f163.jpg" />both of them are real numbers. Let</p><p><img src="11-7401649\09e7a6a0-653a-4019-a0de-0dd1c79f14ce.jpg" /></p><p>Then the derivative of <img src="11-7401649\c40bf7ae-3b7c-4ef4-ba74-53bb3c573e66.jpg" /> along the solution of system (2) satisfies</p><p><img src="11-7401649\715c47cd-7352-43b1-ac65-4be1827f7923.jpg" /></p><p><img src="11-7401649\40022877-c1c7-4070-af11-bdfd477aca97.jpg" /></p><p>Then</p><p><img src="11-7401649\ecbe4c0f-e2a3-4a9b-b9ca-9c1d947d0364.jpg" /></p><p>Simplify, we can get</p><p><img src="11-7401649\48a3a6b2-48bf-4cca-a6c8-4bbc06683233.jpg" /></p><p>Order <img src="11-7401649\eade008d-acee-4b4d-8f46-90a608fde11a.jpg" /> <img src="11-7401649\d5783a83-fb65-418f-a87d-b816467ae30b.jpg" /> Besides, when <img src="11-7401649\f95a3c8a-8025-4734-bd2b-8c19407c866d.jpg" /> if and only if <img src="11-7401649\d22f8d03-9a6e-42b3-a843-de5e282c67b1.jpg" /></p><p>It is easy to show that <img src="11-7401649\2a96c2b8-38e0-4e0b-9813-9e71919360a7.jpg" /> is the largest invariant subset of the set <img src="11-7401649\74fdeed4-2aff-45ec-8a1d-bac64924d57f.jpg" /> Because of LaSalle’s invariance principle, the endemic equilibrium <img src="11-7401649\b4927799-50c8-4c91-aa40-7d510dbb2ebd.jpg" /> of system (2) is globally asymptotically stable when<img src="11-7401649\f9a9054a-1489-4b1b-9339-398b559c2082.jpg" />. This completes the proof.</p></sec><sec id="s3_3"><title>3.3. The SEIQR Epidemic Model with Nonlinear Incidence Rate</title><p>Zhao et al. studied delay SEIR epidemic model with the nonlinear incidence rate like <img src="11-7401649\cbe1a1e7-f928-44f8-add9-8543b969bcd5.jpg" /> in the case of pulse. In this paper, the model without pulse is discussed.</p><p><img src="11-7401649\74a9bd42-b4c4-4826-b7ed-0b6bba932d65.jpg" /></p><p>It is easy to show disease-free equilibrium is globally asymptotically stable, endemic equilibrium is locally asymptotically stable. The ways we use are similar to that in system (1), here they are omitted.</p></sec></sec><sec id="s4"><title>4. The Numerical Simulations</title><p>In this section, we study system (1) numerically. According to the different datas that can reflect the actual situation, we get the different simulation images to prove our conclusions obviously (Figures 1-9).</p><p>Here, according to the different actual situations, while take different parameters, we can get different simulation diagrams of the disease-free equilibrium. At the same time, we find out the disease will die out after much more time when <img src="11-7401649\8fbbd2c6-ccf9-4e90-ac9c-4f43828bab75.jpg" /> increases. For example,</p><p><img src="11-7401649\8d647d90-0311-40b6-aeb8-60a9a49c7c3a.jpg" /></p><p>Here <img src="11-7401649\84b82105-f25a-47d6-b33c-6da9f220bea8.jpg" /> see <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>Here <img src="11-7401649\151a28b8-5c84-45ab-a23c-fbb5345c5b7c.jpg" /> see <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>When take different<img src="11-7401649\db49f878-71a6-459a-8126-1d511b39ee26.jpg" />, we can get different simulation images. In other words, the <img src="11-7401649\5a35c7c4-a6df-4227-97b8-c5b833831fa7.jpg" /> increases when <img src="11-7401649\06c53f03-c169-42bc-a81d-87d64ebaddd1.jpg" /> increases, which is obvious in Figures 3 and 4. And then it is easy for us to find that how <img src="11-7401649\9ba70585-c551-49b7-8f71-d8d49d26ae8c.jpg" /> effects changing trends of<img src="11-7401649\952771de-3d3b-4f6e-bb50-5209bddd70f7.jpg" />, <img src="11-7401649\f536f3fe-15db-4f03-b483-381c939eef07.jpg" />, <img src="11-7401649\cd22b686-5f6e-40be-ac5b-1f101e0557fc.jpg" />, <img src="11-7401649\5ddf28f3-67b0-4524-be51-0fdd738c0467.jpg" />,<img src="11-7401649\abfb9feb-f402-4687-ab5e-c35f8accd22b.jpg" />.</p><p><img src="11-7401649\dbad6973-d356-405b-86eb-a128171828c7.jpg" /></p><p>Here <img src="11-7401649\24dc6211-15d4-4b3d-959b-c986b8c81e14.jpg" /> see <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p><img src="11-7401649\5c1ee199-1c1a-4120-a6f4-5d6800ae9327.jpg" /></p><p>Here <img src="11-7401649\cf380f89-002f-4795-b431-0f58c2c644a8.jpg" /> see <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>At last, if the basic reproduction number is much larger and we will get new diagrams. For example, let</p><p><img src="11-7401649\dc2025b4-af57-4747-ad76-ebc4039b784b.jpg" /></p><p>Here <img src="11-7401649\376dc71a-dbaf-429e-bc80-75268463b694.jpg" /> see <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p><p>At the same time, the changing trends of <img src="11-7401649\70cd772c-577b-4e26-a173-cfcf955e20d3.jpg" /> and <img src="11-7401649\0e81192a-674a-4fc1-bc74-c3e1ed94ae1d.jpg" /> are shown in Figures 6 and 7. And the time which <img src="11-7401649\61c012b0-8928-42c8-a992-9b4cb6a06eda.jpg" /> comes peak will become large as <img src="11-7401649\f3b741b8-2a6b-4e56-b11c-1245a75bffad.jpg" /> increases, the <img src="11-7401649\41ff4348-8ba3-4492-a2ae-f26f72e6c26b.jpg" /> will decrease. For example, <img src="11-7401649\2ec8eb9d-9402-431b-8b54-1f23c9168fe1.jpg" />, see <xref ref-type="fig" rid="fig8">Figure 8</xref>. In addition, when <img src="11-7401649\f70f28fa-c3b6-4c5f-bd29-8d9db0a18416.jpg" /> changes, <img src="11-7401649\0ce23d1f-d5dc-4dca-9c53-12454c5547bb.jpg" />will change. And we can find out that when <img src="11-7401649\281e7338-98f5-4fff-954a-1d0454403c4b.jpg" /> That is, disease will be endemic disease while <img src="11-7401649\011f1ef5-9303-408b-9704-25a2aa036787.jpg" /> see <xref ref-type="fig" rid="fig9">Figure 9</xref>.</p></sec><sec id="s5"><title>5. Discussions</title><p>In this paper, a kind of a delayed SEIQR epidemic model with the quarantine and latent is studied. Using Hurwitz criterion, the local stability of the disease-free equilibrium and endemic equilibrium of system (2) is proved. For any time delay <img src="11-7401649\811c20ad-810e-40b2-8e3a-c16b38d8d04d.jpg" /> we prove the disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than unity and the endemic equilibrium is globally asymptotically stable when the basic reproduction number is greater than unity by means of suitable Lyapunov functions and LaSalle’s invariance principle. So the delay is harmless to system (2). From the biological point of view, the delay here has no influence on the transmission of diseases. However, in [<xref ref-type="bibr" rid="scirp.38059-ref16">16</xref>], the disease-free equilibrium is periodic and globally attractive. At the same time, the disease will be endemic after some period of time. Above all, we consider that <img src="11-7401649\c5661855-a556-45da-98a8-7539141faaf4.jpg" /> is quarantined and can recover in this model, which will effect changing trends of<img src="11-7401649\ab2edfb3-04fb-427d-8bd8-4cea5900457f.jpg" />, <img src="11-7401649\53fa6a70-67f5-4b2b-92ab-47bdf3381c4d.jpg" />, <img src="11-7401649\955fbefc-24ca-4e38-96f4-765e6000e35b.jpg" />, <img src="11-7401649\1808134d-a852-47b0-b0c1-558e013dbf4a.jpg" />,<img src="11-7401649\6e828431-45b3-434a-9566-dbb42390dcbe.jpg" />. Here, we take <img src="11-7401649\347569ee-bb82-4f34-b9ad-f012f3b95fac.jpg" /> as an example to explain that. Meanwhile, the simulation image which <img src="11-7401649\ea8d0b0d-329d-4a7a-ab22-cb7860abaca8.jpg" /> changes as <img src="11-7401649\13f1596c-2cf3-4a00-b433-3dbae63b2420.jpg" /> can be obtained and we can find out <img src="11-7401649\bd940e30-f296-434a-95e2-5d1a811d8f60.jpg" /> which the basic reproduction number is a unity. Those are useful for us to control epidemics. At last, the conclusions above are verified by numerical simulations.</p></sec><sec id="s6"><title>6. Acknowledgements</title><p>This research was supported by the National Science Foundation of China (10471040) and the National Sciences Foundation of Shanxi Province (2009011005-1).</p></sec><sec id="s7"><title>REFERENCES</title></sec><sec id="s8"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.38059-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">X. B. Liu and L. J. Yang, “Stability Analysis of an SEIQV Epidemic Model with Saturated Incidence Rate,” Nonlinear Analysis: Real World Applications, Vol. 13, No. 6, 2012, pp. 2671-2979.</mixed-citation></ref><ref id="scirp.38059-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">M. Y. Li, J. R. Graef, L. C. Wang and J. Karsai, “Global Dynamics of a SEIR Model with Varying Total Population Size,” Mathematical Biosciences, Vol. 160, No. 2, 1999, pp. 191-213.http://dx.doi.org/10.1016/S0025-5564(99)00030-9</mixed-citation></ref><ref id="scirp.38059-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">J. Zhang and Z. E. Ma, “Global Dynamics of an SEIR Epidemic Model with Saturating Contact Rate,” Mathematical Biosciences, Vol. 185, No. 1, 2003, pp. 15-32.http://dx.doi.org/10.1016/S0025-5564(03)00087-7</mixed-citation></ref><ref id="scirp.38059-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">A. Onofrio, “Stability Properties of Pulse Vaccination Strategy in SEIR Epidemic Model,” Mathematical Biosciences, Vol. 179, No. 1, 2002, pp. 57-72.http://dx.doi.org/10.1016/S0025-5564(02)00095-0</mixed-citation></ref><ref id="scirp.38059-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">G. H. 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, 2005, pp. 1177-1184.http://dx.doi.org/10.1016/j.chaos.2004.11.062</mixed-citation></ref><ref id="scirp.38059-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">B. K. Mishra and N. Jha, “SEIQRS Model for the Transmission of Malicious Objects in Computer Network,” Applied Mathematical Modelling, Vol. 34, No. 3, 2010, pp. 710-715. http://dx.doi.org/10.1016/j.apm.2009.06.011</mixed-citation></ref><ref id="scirp.38059-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Y. N. Kyrychkoa and K. B. Blyussb, “Global Properties of a Delayed SIR Model with Temporary Immunity and Nonlinear Incidence Rate,” Nonlinear Analysis: Real World Applications, Vol. 6, No. 3, 2005, pp. 495-507.http://dx.doi.org/10.1016/j.nonrwa.2004.10.001</mixed-citation></ref><ref id="scirp.38059-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">K. L. Cooke and P. van den Driessche, “Analysis of an SEIRS Epidemic Model with Two Delays,” Journal of Mathematical Biology, Vol. 35, No. 2, 1996, pp. 240-260.http://dx.doi.org/10.1007/s002850050051</mixed-citation></ref><ref id="scirp.38059-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">R. Xu, Z. E. Ma and Z. P. Wang, “Global Stability of a Delayed SIRS Epidemic Model with Saturation Incidence and Temporary Immunity,” Computers and Mathematics with Applications, Vol. 59, No. 9, 2010, pp. 3211-3221.http://dx.doi.org/10.1016/j.camwa.2010.03.009</mixed-citation></ref><ref id="scirp.38059-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Y. Enatsu, E. Messina, Y. Muroya, Y. Nakata, E. Russo and A. Vecchio, “Stability Analysis of Delayed SIR Epidemic Models with a Class of Nonlinear Incidence Rates,” Applied Mathematics and Computation, Vol. 218, No. 9, 2012, pp. 5327-5336.http://dx.doi.org/10.1016/j.amc.2011.11.016</mixed-citation></ref><ref id="scirp.38059-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">F. P. Zhang, Z. Z. Li and F. Zhang, “Global Stability of an SIR Epidemic Model with Constant Infectious Period,” Applied Mathematics and Computation, Vol. 199, No. 1, 2008, pp. 285-291.http://dx.doi.org/10.1016/j.amc.2007.09.053</mixed-citation></ref><ref id="scirp.38059-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">S. J. Gao, L. S. Chen and Z. D. Teng, “Pulse Vaccination of an SEIR Epidemic Model with Time Delay,” Nonlinear Analysis: Real World Applications, Vol. 9, No. 2, 2008, pp. 599-607.http://dx.doi.org/10.1016/j.nonrwa.2006.12.004</mixed-citation></ref><ref id="scirp.38059-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Z. Zhao, L. S. Chen and X. Y. Song, “Impulsive Vaccination of SEIR Epidemic Model with Time Delay and Nonlinear Incidence Rate,” Mathematics and Computers in Simulation, Vol. 79, No. 3, 2008, pp. 500-510.http://dx.doi.org/10.1016/j.matcom.2008.02.007</mixed-citation></ref><ref id="scirp.38059-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">R. Xu and Z. E. Ma, “Global Stability of a Delayed SEIRS Epidemic Model with Saturation Incidence Rate,” Nonlinear Dynamics, Vol. 61, No. 1-2, 2010, pp. 229239. http://dx.doi.org/10.1007/s11071-009-9644-3</mixed-citation></ref><ref id="scirp.38059-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">X. B. Zhang, H. F. Huo, H. Xiang and X. Y. Meng, “Two Profitless Delays for the SEIRS Epidemic Disease Model with Nonlinear Incidence and Pulse Vaccination,” Applied Mathematics and Computation, Vol. 186, No. 1, 2007, pp. 516-529.http://dx.doi.org/10.1016/j.amc.2006.07.124</mixed-citation></ref><ref id="scirp.38059-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Y. Z. Pei, S. Y. Liu, S. J. Gao, S. P. Li and C. G. Li, “A Delayed SEIQR Epidemic Model with Pulse Vaccination and the Quarantine Measure,” Computers and Mathematics with Applications, Vol. 58, No. 1, 2009, pp. 135-145.</mixed-citation></ref><ref id="scirp.38059-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">H. M. Wei, X. Z. Li and M. Martcheva, “An Epidemic Model of a Vector-Borne Disease with Direct Transmission and Time Delay,” Journal of Mathematical Analysis and Applications, Vol. 342, No. 2, 2008, pp. 895-908.</mixed-citation></ref></ref-list></back></article>