<?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.2014.521313</article-id><article-id pub-id-type="publisher-id">AM-52113</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Stability of Highly Pathogenic Avian Influenza Epidemic Model with Saturated Contact Rate
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>huqin</surname><given-names>Che</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</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></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Likang</surname><given-names>Ma</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Jinxi Industries Group Co., Ltd., Taiyuan, China</addr-line></aff><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>2641511342@qq.com(HC)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>01</day><month>12</month><year>2014</year></pub-date><volume>05</volume><issue>21</issue><fpage>3365</fpage><lpage>3371</lpage><history><date date-type="received"><day>29</day>	<month>September</month>	<year>2014</year></date><date date-type="rev-recd"><day>28</day>	<month>October</month>	<year>2014</year>	</date><date date-type="accepted"><day>18</day>	<month>November</month>	<year>2014</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 present a highly pathogenic Avian influenza epidemic model with saturated contact rate. According to study of the dynamics, we calculated the basic reproduction number of the model. Through the analysis of this model, we have the following conclusion: if R
  <sub>0</sub> ≤ 1, there is only one disease-free equilibrium which is globally stable, the disease will die; if R
  <sub>0</sub> &gt; 1, there is only one endemic equilibrium which is globally stable, disease will be popular.
 
</p></abstract><kwd-group><kwd>Avian Influenza</kwd><kwd> The Basic Reproduction Number</kwd><kwd> Stability</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Avian influenza virus belongs to the influenza A virus. According to the difference of the pathogenic Avian influenza virus in chicken and turkey, we divided it into three levels: high, medium, low/non pathogenic. Because of the bird flu virus hemagglutinin structure characteristics, general infected birds, when the virus genetic reassortment during replication, causing structural changes. The Avian influenza virus which acquires the ability to infect people, can make adult infect the Avian influenza disease. The highly pathogenic Avian influenza has a high death rate, which is about 100 percent for birds and more than 70 percent for humans [<xref ref-type="bibr" rid="scirp.52113-ref1">1</xref>] .</p><p>Up to now, we found that Avian influenza virus subtypes which can directly infect human are: H5N1, H7N1, H7N2, H7N3, H7N7, H9N2, H7N9, subtype. Among them, the new subtype Avian influenza virus H7N9 subtype, which was first discovered in 2013 March, and the high pathogenic Avian influenza H5N1 subtype, are particularly noteworthy. The Avian influenza virus not only caused human casualties, but also hit the poultry industry.</p><p>Avian influenza virus (H5N1 subtype Avian influenza virus) had 12 human infections in Hongkong, including 6 deaths. According to the WHO’s statistics, 393 people around the world in more than 15 countries and regions were infected, including 248 deaths; the mortality rate was 0.63. Since 2003, 31 people have been infected with Avian influenza in China, including 21 deaths. Although people infection probability is not high, the mortality rate is very high. By March 2013, the world has reported a total of human infection of highly pathogenic H5N1 Avian influenza in 622 cases, including 371 deaths. The distribution of cases in 15 countries, including China, is found in 45 cases, 30 cases of death. Most of human infections with H5N1 Avian influenza are young people and children. In March 2013, human infection with H7N9 Avian influenza was first found in China. By May 1, 2013, Shanghai, Anhui, Jiangsu, Zhejiang, Beijing, Henan, Shandong, Jiangxi, Hunan, Fujian and other 10 provinces (city) have reported 127 confirmed cases, including 26 death cases. The cases in the majority are with the elderly, men more than women [<xref ref-type="bibr" rid="scirp.52113-ref2">2</xref>] .</p><p>The present study found that people are infected with bird flu infection for the virus’s birds, but the transmission is still to be defined. Research suggests that the main way of human infection with the H5N1 Avian influenza is in close contact with poultry, high-risk behavior, including the slaughter, feathers and processing of infected poultry. In a few cases, when children frequently play in the backyard poultry regional, the exposure poultry manure was also thought to be a source of infection. At present, most of the researched evidences indicate the presence of bird-human transmission; there may be an environment-human (birds excrement pollution) transmission, and a few non sustained human-human transmission of H5N1. Now H7N9 Avian influenza patients are infected through direct contact with poultry or the items, environment which are contaminated by feces. Human infection with H7N9 Avian influenza is still sporadic, although there have been cases of aggregation of individual families; but it is not currently found that the virus has sustained human-human transmission capacity.</p><p>At present, some authors have researched some Avian influenza model [<xref ref-type="bibr" rid="scirp.52113-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.52113-ref4">4</xref>] , they had constructed a mathe- matical model which interprets the spread of Avian influenza from the bird world to the human world. Literature [<xref ref-type="bibr" rid="scirp.52113-ref5">5</xref>] has introduced a piecewise treatment function. When the number of the infective had not exceeded the maxi- mum treatment capacity, the treatment rate was proportional to the number of the infective. When the number of the infective had exceeded the maximum treatment capacity, it took maximum saturation treatment value. Lit- erature [<xref ref-type="bibr" rid="scirp.52113-ref6">6</xref>] -[<xref ref-type="bibr" rid="scirp.52113-ref8">8</xref>] have studied of SIR model with saturated treatment rate. In the literature [<xref ref-type="bibr" rid="scirp.52113-ref8">8</xref>] , the saturated</p><p>treatment rate is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x6.png" xlink:type="simple"/></inline-formula>, where r is the cure rate and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x7.png" xlink:type="simple"/></inline-formula> is the parameters of infection which is due to delayed</p><p>treatment. The conclusion had indicated that in the prevention and treatment of Avian flu drugs under the condition of limited, culling of infected poultry was the most effective way to control the spread of Avian flu in humans. Literature [<xref ref-type="bibr" rid="scirp.52113-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.52113-ref10">10</xref>] have proposed the model of highly pathogenic Avian influenza.</p><p>In this paper, the authors propose a Avian influenza model with saturated contact rate [<xref ref-type="bibr" rid="scirp.52113-ref11">11</xref>] ,</p><disp-formula id="scirp.52113-formula84"><label>(1.1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-7402380x8.png"  xlink:type="simple"/></disp-formula><p>In system (1.1), the human is divided into three compartments: Susceptible (S), infected (I), recovery (R). The birds are divided into susceptible poultry (X) and infected poultry (Y). The parameters c and b are respectively the natural birth rate of Avian and human. d and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x9.png" xlink:type="simple"/></inline-formula> are respectively the natural mortality of poultry and human. m and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x10.png" xlink:type="simple"/></inline-formula> are respectively the poultry and human mortality due to illness. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x11.png" xlink:type="simple"/></inline-formula>stands for infectious rate of susceptible poultry to infected poultry. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x12.png" xlink:type="simple"/></inline-formula>stands for infected poultry of the infection rate of susceptible individuals. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x13.png" xlink:type="simple"/></inline-formula>is the recovery rate that infects individuals through treatment. When Y is small, the contact ratio, infected poultry and susceptible poultry, is approximatively proportional to the Y; With the increase of Y, the con-</p><p>tact rate gradually reaches saturation. When Y is very large, it is close to a constant<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x14.png" xlink:type="simple"/></inline-formula>. The same way to ex-</p><p>plain<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x15.png" xlink:type="simple"/></inline-formula>, that is to say, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x16.png" xlink:type="simple"/></inline-formula>is a parameter, which is effects of infectious diseases, when the contact rate of</p><p>the disease is saturated. The transmission process of the disease is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p></sec><sec id="s2"><title>2. The Existence of the Equilibrium Point</title><p>The first four equations of system (1.1) do not contain R, by the method of van den Driessche and Watmough Diekmann [<xref ref-type="bibr" rid="scirp.52113-ref12">12</xref>] , we only need to study the following system:</p><disp-formula id="scirp.52113-formula85"><label>(2.1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-7402380x17.png"  xlink:type="simple"/></disp-formula><p>It’s easy to get a disease-free equilibrium of the system (2.1) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x18.png" xlink:type="simple"/></inline-formula>is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x19.png" xlink:type="simple"/></inline-formula>, by the method of van den Driessche and Watmough, we can get the basic reproductive number of the system (2.1)</p><disp-formula id="scirp.52113-formula86"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x20.png"  xlink:type="simple"/></disp-formula><p>By the positive of the endemic equilibrium point, we can get that if R<sub>0</sub> &gt; 1, there is a unique endemic equilibrium<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x21.png" xlink:type="simple"/></inline-formula>, which satisfied:</p><disp-formula id="scirp.52113-formula87"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x22.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.52113-formula88"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x23.png"  xlink:type="simple"/></disp-formula><p>So, we can get the following theorem:</p><p>Theorem 2.1 If R<sub>0</sub> ≤ 1, the system (2.1) only exists the disease-free equilibrium<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x24.png" xlink:type="simple"/></inline-formula>; when R<sub>0</sub> &gt; 1, there exists only one endemic equilibrium</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x25.png" xlink:type="simple"/></inline-formula>.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Dissemination flow diagram</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-7402380x26.png"/></fig></sec><sec id="s3"><title>3. The Stability of the Disease-Free Equilibrium</title><p>The Jacobian matrix of system (2.1) is</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x27.png" xlink:type="simple"/></inline-formula>,</p><p>which in the disease-free equilibrium E<sub>0</sub>, is</p><disp-formula id="scirp.52113-formula89"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x28.png"  xlink:type="simple"/></disp-formula><p>The characteristic equation of the Jacobian matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x29.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.52113-formula90"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x30.png"  xlink:type="simple"/></disp-formula><p>Here, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x31.png" xlink:type="simple"/></inline-formula>, λ denotes the indeterminate of the polynomial. If and only if R<sub>0</sub> ≤</p><p>1, all roots of this characteristic equation have negative real parts. It implies that E<sub>0</sub> is locally asymptotically stable. Then we can get the following theorem:</p><p>Theorem 3.1 If R<sub>0</sub> ≤ 1, the disease-free equilibrium E<sub>0</sub> is locally asymptotically stable; if R<sub>0</sub> &gt; 1, the disease- free equilibrium E<sub>0</sub> is unstable.</p><p>Next, we discuss the global stability of E<sub>0</sub>, considering the Liapunov function [<xref ref-type="bibr" rid="scirp.52113-ref13">13</xref>] <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x32.png" xlink:type="simple"/></inline-formula>, then</p><disp-formula id="scirp.52113-formula91"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x33.png"  xlink:type="simple"/></disp-formula><p>when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x34.png" xlink:type="simple"/></inline-formula>, we can get<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x35.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x36.png" xlink:type="simple"/></inline-formula> has no other closed trajectory in addition to E<sub>0</sub>. Then we can get the following theorem:</p><p>Theorem 3.2 If<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x37.png" xlink:type="simple"/></inline-formula>, the disease-free equilibrium E<sub>0</sub> is globally asymptotically stable.</p></sec><sec id="s4"><title>4. The Stability of the Endemic Equilibrium</title><p>The Jacobian matrix of system (2.1), which in the endemic equilibrium E<sub>+</sub>, is</p><disp-formula id="scirp.52113-formula92"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x38.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.52113-formula93"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x39.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.52113-formula94"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x40.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.52113-formula95"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x41.png"  xlink:type="simple"/></disp-formula><p>Thus, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x42.png" xlink:type="simple"/></inline-formula>evaluated is stable if and only if so are A and C. Then the characteristic equation of the matrix A is</p><disp-formula id="scirp.52113-formula96"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x43.png"  xlink:type="simple"/></disp-formula><p>where,</p><disp-formula id="scirp.52113-formula97"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x44.png"  xlink:type="simple"/></disp-formula><p>Obviously, when R<sub>0</sub> &gt; 1, h<sub>1</sub> &gt; 0, h<sub>2</sub> &gt; 0. By the Harwize criterion, it implies that the characteristic roots of matrix A have negative real parts. then, we analyze the characteristics equation of the matrix C</p><disp-formula id="scirp.52113-formula98"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x45.png"  xlink:type="simple"/></disp-formula><p>Easily, if and only if R<sub>0</sub> &gt; 1, all roots of this characteristic equation have negative real parts. So, we can get the following theorem:</p><p>Theorem 4.1 If R<sub>0</sub> &gt; 1, the endemic equilibrium E<sub>+</sub> is locally asymptotically stable.</p><p>Next, we discuss the global stability of E<sub>+</sub>, considering the Liapunov function</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x46.png" xlink:type="simple"/></inline-formula>, then derivated W<sub>2</sub> about system (2.1),</p><disp-formula id="scirp.52113-formula99"><graphic  xlink:href="http://html.scirp.org/file/7-7402380x47.png"  xlink:type="simple"/></disp-formula><p>By the relationship of arithmetic mean and geometric mean, we can know<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x48.png" xlink:type="simple"/></inline-formula>. That is to say,</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x49.png" xlink:type="simple"/></inline-formula>. If and only if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x50.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-7402380x51.png" xlink:type="simple"/></inline-formula>. Thus, by LaSalle Invariance Principle [<xref ref-type="bibr" rid="scirp.52113-ref14">14</xref>] , we can get the following theorem:</p><p>Theorem 4.2 If R<sub>0</sub> &gt; 1, the endemic equilibrium E<sub>+</sub> is globally asymptotically stable.</p></sec><sec id="s5"><title>5. The Numerical Simulations</title><p>In this section, the model has been analyzed numerically by using MAT-LAB 7.5.0, to confirm the results we have obtained. For this purpose we have used the following set of parameter values:</p><p>When c = 2; β = 0.01; d = 0.03; ω = 0.014; m = 0.97; b = 1; α = 0.069; ε = 0.63; γ = 0.301; δ = 0.05, then R<sub>0</sub> = 0.9333. These parameters values satisfy the condition R<sub>0</sub> &lt; 1. Thus, the two images of <xref ref-type="fig" rid="fig2">Figure 2</xref> show that if R<sub>0</sub> &lt; 1, the disease-free equilibrium E<sub>0</sub> is global asymptotically stable.</p><p>When c = 2; β = 0.01; d = 0.03; ω = 0.02; m = 0.97; b = 1; α = 0.069; ε = 0.63; γ = 0.301; δ = 0.05 then R<sub>0</sub> = 1.3333; These parameters values satisfy the condition R<sub>0</sub> &gt; 1. Thus, the two image of <xref ref-type="fig" rid="fig3">Figure 3</xref> shows that if R<sub>0</sub> &gt; 1, the unique endemic equilibrium E<sup>*</sup> is global asymptotically stable.</p></sec><sec id="s6"><title>6. Conclusion</title><p>In this paper we have discussed the stability of highly pathogenic Avian influenza epidemic model with saturated contact rate. We get the reproduction number R<sub>0</sub>; it is the threshold which is endemic or not. If R<sub>0</sub> ≤ 1,</p><fig-group id="fig2"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> c = 2; β = 0.01; d = 0:03; ω = 0.014; m = 0.97; b = 1; α = 0.069; ε = 0.63; γ = 0.301; δ = 0.05; R<sub>0</sub> = 0.9333.</title></caption><fig id ="fig2_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-7402380x52.png"/></fig><fig id ="fig2_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-7402380x53.png"/></fig></fig-group><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> c = 2; β = 0.01; d = 0.03; ω = 0.02; m = 0.97; b = 1; α = 0.069; ε = 0.63; γ = 0.301; δ = 0.05; R<sub>0</sub> = 1.3333.</title></caption><fig id ="fig3_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-7402380x54.png"/></fig><fig id ="fig3_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-7402380x55.png"/></fig></fig-group><p>the disease-free equilibrium is global asymptotic stability which implies that the disease will be extinct. If R<sub>0</sub> &gt; 1, the endemic equilibrium E<sub>+</sub> is globally asymptotically stable, that implies that the disease will sustain and lead to epidemic disease eventually. In order to prevent epidemics, through the analysis of the model, we can get: reducing the activities of the poultry market; culling or burying the infected poultry. That the infected poultry feces get disinfected or burned can prevent the pollution of the surrounding water. The government must strictly manage cooked poultry meat, and meat sources. Traditional farming poultry from backyard changes the concentration of scale breeding, slaughter processing and scientific transport, raising poultry breeding, circulation biological safety level, thereby reduces the population of live poultry or poultry exposure opportunity. Early discovery, early diagnosis of Avian influenza are the implementation of case isolation and diagnosis and treatment timely, effectively and reasonably.</p></sec><sec id="s7"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.52113-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Iwami, S., Takeuchi, Y. and Liu, X.N. (2007) Avian-Human Influenza Epidemic Model. Mathematical Biosciences, 207, 1-25. http://dx.doi.org/10.1016/j.mbs.2006.08.001</mixed-citation></ref><ref id="scirp.52113-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Feng, Z.J. (2013) Human Infect with Avian Influenza H7N9 Treatment Program. Encyclopedia of Medical Network, 4.  
http://www.baikemy.com/jiankangkepu/10339395895809</mixed-citation></ref><ref id="scirp.52113-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Iwami, S., Takeuchi, Y., Korobeinikov, A. and Liu, X. (2008) Prevention of Avian Influenza Epidemic: What Policy Should We Choose. Journal of Theoretical Biology, 252, 732-741. http://dx.doi.org/10.1016/j.mbs.2006.08.001</mixed-citation></ref><ref id="scirp.52113-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Rao, D.M., Chernyakhovsky, A. and Rao, V. (2009) Modeling and Analysis of Global Epidemiology of Avian Influenza. Environmental Modeling &amp; Software, 24, 124-134. http://dx.doi.org/10.1016/j.envsoft.2008.06.011</mixed-citation></ref><ref id="scirp.52113-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Wang, W. (2006) Backward Bifurcation of an Epidemic Treatment. Mathermatical Biosciences, 201, 58-71.  
http://dx.doi.org/10.1016/j.mbs.2005.12.022</mixed-citation></ref><ref id="scirp.52113-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, X. and Liu, X. (2008) Backford Bifurcation of an Epidemic Model with Saturated Treatment Function. Journal of Mathematical Analysis and Applications, 348, 433-443. http://dx.doi.org/10.1016/j.jmaa.2008.07.042</mixed-citation></ref><ref id="scirp.52113-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Capasso, V. and Serio, G. (1978) A Generalization of the Kermack-Mckendrick Deterministic Epidemic Model. Mathematical Biosciences, 42, 43-61. http://dx.doi.org/10.1016/0025-5564(78)90006-8</mixed-citation></ref><ref id="scirp.52113-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Guo, S.M., Guo, L.N. and Li, X.Z. (2010) Analysis of an Avian Influenza Epidemic Model with Saturation Treatment. Mathematics in Practice and Theory, 40, 134-137.</mixed-citation></ref><ref id="scirp.52113-ref9"><label>9</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Li</surname><given-names> Y.-L. </given-names></name>,<etal>et al</etal>. (<year>2009</year>)<article-title>Study on SI Transmission Model of Highly Pathogenic Avian Influenza</article-title><source> Journal of Anhui Agricultural Sciences</source><volume> 37</volume>,<fpage> 13603</fpage>-<lpage>13605</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.52113-ref10"><label>10</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Liang</surname><given-names> R.H. </given-names></name>,<etal>et al</etal>. (<year>2007</year>)<article-title>The Construction and Application of SEID Model for the Disseminate Mechanism of Highly Pathogenic Avian Influenza (HPAI)</article-title><source> Journal of Xinyang Normal University: Natural Science Edition</source><volume> 20</volume>,<fpage> 262</fpage>-<lpage>265</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.52113-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Lasalle, J.P. (1986) The Stability and Control of Discrete Process. Springer-Verlag, New York.  
http://dx.doi.org/10.1007/978-1-4612-1076-4</mixed-citation></ref><ref id="scirp.52113-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">van den Driessche, P. and Watmough, J. (2007) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 196, 1679-1684.</mixed-citation></ref><ref id="scirp.52113-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Ma, Z.-E. and Zhou, Y.-C. (2001) The Theory of Ordinary Differential Equation Qualitative and Stability. Science Press, Beijing.</mixed-citation></ref><ref id="scirp.52113-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Lasalle, J.P. (1976) The Stability of Dynamical Systems. Regional Conference Series in Applied Mathematics Philadelphia, SIAM, 21, 418-420.</mixed-citation></ref></ref-list></back></article>