Chaos, Solitons and Fractals
|
|
- Cynthia Nelson
- 5 years ago
- Views:
Transcription
1 Chaos, Solitons and Fractals 42 (2009) Contents lists available at ScienceDirect Chaos, Solitons and Fractals journal homepage: Solutions of the SIR models of epidemics using HAM Fadi Awawdeh *, A. Adawi, Z. Mustafa Department of Mathematics, Hashemite University, Jordan article info abstract Article history: Accepted 3 April 2009 In this paper, we investigate the accuracy of the Homopy Analysis Method (HAM) for solving the problem of the spread of a non-fatal disease in a population. The advantage of this method is that it provides a direct scheme for solving the problem, i.e., without the need for linearization, perturbation, massive computation and any transformation. Mathematical modeling of the problem leads to a system of nonlinear ODEs. MATLAB 7 is used to carry out the computations. Graphical results are presented and discussed quantitatively to illustrate the solution. Ó 2009 Elsevier Ltd. All rights reserved. 1. Introduction In the modeling transmission dynamics of a communicable disease, it is common to divide the population into disjoint classes (compartments) whose sizes change with time. The infection status of any individual in a population can be Susceptible, when the person is healthy and susceptible to the disease (dened by S), Exposed, when the person is in a latent period but n yet infectious (dened by E), Infected, when the individual carries the disease and is infectious (dened by I), or Removed, when the person has recovered and is at least temporarily immune or has died because of disease (dened by R). In some diseases such as HIV, there is no recovery. In her diseases, if an infected person recovered he/she may be susceptible again. A sequence of letters, such as SEIR, describes the movement of individuals between the classes: susceptible become latent, then infectious and finally recover with immunity. To model diseases which confer permanent immunity and which are endemic because of births of new susceptible, SIR or SEIR models with vital dynamics are suitable. Vital dynamics is needed to avoid explosion of the population size. Models of SEIRs or SIRs types are used to model diseases with temporary immunity and in cases where there is no immunity, models are named SIS or SEIS. The last S points the individual becoming susceptible again, after recovery. Such models may be appropriate for gonorrhea, for instance. Epidemic models has been widely used in different forms for studying epidemiological processes such as the spread of influenza [9] and SARS [2,8,11] and even for the spread of rumors [13,14]. Epidemic models are also applied to modeling of STI epidemics, but n all epidemic models are suitable for STIs since the sexual network plays an important role in spread of disease. Pair-formation models are a type of ordinary differential equation models that have sometimes been used to study STI transmission in populations. They incorporate the repeated contacts within partnerships which happen frequently in real sexual networks. They were first developed in 1988 by Dietz and Hadeler [3] to study STIs in monogamous partnerships. In this model if two susceptible individuals form a pair then they can be considered temporarily immune as long as they do n separate and have no contacts with her partners. This aspect influences transmission dynamics considerably, especially when the disease is first introduced, since the vast majority of existing pairs are susceptible. * Corresponding author. addresses: awawdeh@hu.edu.jo (F. Awawdeh), adawi@hu.edu.jo (A. Adawi), zmagablh@hu.edu.jo (Z. Mustafa) /$ - see front matter Ó 2009 Elsevier Ltd. All rights reserved. doi: /j.chaos
2 3048 F. Awawdeh et al. / Chaos, Solitons and Fractals 42 (2009) This is a simple model, due to Kermack and MacKendrick [5], of an epidemic of an infectious disease in a population. We assume the population consists of three types of individuals, whose numbers are dened by the letters S; I and R (which is why this is called an SIR model). All these of course are functions of time. SðtÞ is the number of susceptible, who do n have the disease but could get it. IðtÞ is the number of infectives, who have the disease and can transmit it to hers. RðtÞ is the number of removed, who cann get the disease or transmit it: either they have a natural immunity, or they have recovered from the disease and are immune from getting it again, or they have been placed in isolation, or they have died. The mathematical model does n distinguish between these possibilities. The following simple SIR model [1,4,10,12] is tested to show the efficiency of the HAM to solve such models. Assume there is a steady constant rate between susceptible and infectives and that a constant proportion of these constant result in transmission. Then in time dt; ds of the susceptible become infective, where ds ¼ bsidt; and b is a positive constant. If c > 0 is the rate at which current infectives become isolated, then di ¼ bsidt cidt; and the number of new isolates dr is given by dr ¼ cidt: If we let dt! 0, then the following nonlinear system of ODEs determines the progress of the disease: ds dt ¼ bsi; di ¼ bsi ci; dt ð2þ dr dt ¼ ci with initial conditions Sð0Þ ¼N S ; Ið0Þ ¼N I ; Rð0Þ ¼N R : In this paper, the HAM [6,7] is used to solve the epidemic models. An example is tested, and the obtained results suggest that newly improvement technique introduces a promising tool and powerful improvement for solving nonlinear equations. Unlike all previous techniques, the proposed technique provides us with a family of solution expressions in auxiliary parameter h. As a result, the convergence region and rate of solution depend upon the auxiliary parameter h and thus can be greatly enlarged by means of choosing a proper value of h. This provides a convenient way to adjust and control convergence region and rate of solution given by the technique. Liao [6] suggested the so-called h-curves to determine the value of h. 2. Basic idea of the HAM Consider N½yðtÞŠ ¼ 0; where N is any operator, yðtþ is unknown function and t the independent variable. Let y 0 ðtþ dene an initial guess of the exact solution yðtþ, h 0 an auxiliary parameter, HðtÞ 0 an auxiliary function, and L an auxiliary linear operator with the property L½yðtÞŠ ¼ 0 when yðtþ ¼0. Then using q 2½0; 1Š as an embedding parameter, we construct such a homopy ð1 qþl½/ðt; qþ y 0 ðtþš qhhðtþn½/ðt; qþš ¼ b H½/ðt; qþ; y 0 ðtþ; HðtÞ; h; qš: It should be emphasized that we have great freedom to choose the initial guess y 0 ðtþ, the auxiliary linear operator L, the non-zero auxiliary parameter h, and the auxiliary function HðtÞ. Enforcing the homopy (4) to be zero, i.e., bh½/ðt; qþ; y 0 ðtþ; HðtÞ; h; qš ¼0 ð1þ ð3þ ð4þ we have the so-called zero-order deformation equation ð1 qþl½/ðt; qþ y 0 ðtþš ¼ qhhðtþn½/ðt; qþš: When q ¼ 0, the zero-order deformation equation (5) becomes /ðt; 0Þ ¼y 0 ðtþ; and when q ¼ 1, since h 0 and HðtÞ 0, the zero-order deformation equation (5) is equivalent to /ðt; 1Þ ¼yðtÞ: ð5þ ð6þ ð7þ
3 F. Awawdeh et al. / Chaos, Solitons and Fractals 42 (2009) Thus, according to (6) and (7), as the embedding parameter q increases from 0 to 1, /ðt; qþ varies continuously from the initial approximation y 0 ðtþ to the exact solution yðtþ. Such a kind of continuous variation is called deformation in homopy. By Taylor s theorem, /ðt; qþ can be expanded in a power series of q as follows: /ðt; qþ ¼y 0 ðtþþ X1 y m ðtþq m where y m ðtþ ¼ 1 o m /ðt; qþ m! oq m : If the initial guess y 0 ðtþ, the auxiliary linear parameter L, the non-zero auxiliary parameter h, and the auxiliary function HðtÞ are properly chosen so that 1. The solution /ðt; qþ of the zero-order deformation equation (5) exists for all q 2½0; 1Š. 2. The deformation derivative om /ðt;qþ j oq m exists for m ¼ 1; 2; The power series (8) of /ðt; qþ converges at q ¼ 1. Then, we have under these assumptions the solution series /ðt; 1Þ ¼y 0 ðtþþ X1 y m ðtþ: For brevity, define the vector! y nðtþ ¼fy 0 ðtþ; y 1 ðtþ; y 2 ðtþ;...; y n ðtþg: According to the definition (9), the governing equation of y m ðtþ can be derived from the zero-order deformation equation (5). Differentiating the zero-order deformation equation (5) m times with respective to q and then dividing by m! and finally setting q ¼ 0, we have the so-called mth-order deformation equation where and L½y m ðtþ v m y m 1 ðtþš ¼ hhðtþr m ðy m 1 ðtþþ 1 R m ðy m 1 ðtþþ ¼ ðm 1Þ! v m ¼ 0; m 6 1 : 1; m > 1 o m 1 N½/ðt; qþš oq m 1 Ne that the high-order deformation equation (12) is governing by the linear operator L, and the term R m ðy m 1 ðtþþ can be expressed simply by (13) for any nonlinear operator N. According to the definition (13), the right-hand side of Eq. (12) is only dependent upon y m 1 ðtþ. Thus, we gain y 1 ðtþ; y 2 ðtþ;... by mean of solving the linear high-order deformation equation (12) one after the her in order. 3. Solution of the SIR model by HAM In order to explicitly construct approximate non-perturbative solutions of the system described by Eqs. (1) (3), HAM well addressed in [6,7] is employed. The advantage of this method is that it provides a direct scheme for solving the problem. To apply the HAM, we choose S 0 ðtþ ¼N S ; I 0 ðtþ ¼N I ; R 0 ðtþ ¼N R as initial approximations of SðtÞ; IðtÞ and RðtÞ. Let q 2½0; 1Š be the so-called embedding parameter. The HAM is based on a kind of continuos mappings SðtÞ!/ 1 ðt; qþ; IðtÞ!/ 2 ðt; qþ; RðtÞ!/ 3 ðt; qþ such that, as the embedding parameter q increases from 0 to 1, / i ðt; qþ varies from the initial approximation to the exact solution. To ensure this, choose such auxiliary linear operators as L i ½/ i ðt; qþš ¼ o/ iðt; qþ ; i ¼ 1; 2; 3 with the property L i ½C i Š¼0 ð8þ ð9þ ð10þ ð11þ ð12þ ð13þ
4 3050 F. Awawdeh et al. / Chaos, Solitons and Fractals 42 (2009) where C i are integral constants. We define the nonlinear operators N 1 ½/ i ðt; qþš ¼ o/ iðt; qþ þ b/ i ðt; qþ/ 2 ðt; qþ; N 2 ½/ i ðt; qþš ¼ o/ iðt; qþ b/ 1 ðt; qþ/ i ðt; qþ c/ i ðt; qþ; N 3 ½/ i ðt; qþš ¼ o/ iðt; qþ c/ 2 ðt; qþ: Let h i 0 and H i ðtþ 0 dene the so-called auxiliary parameter and auxiliary function, respectively. Using the embedding parameter q, we construct a family of equations ð1 qþl½/ 1 ðt; qþ S 0 ðtþš ¼ qh 1 H 1 ðtþn 1 ½/ 1 ðt; qþš; ð1 qþl½/ 2 ðt; qþ I 0 ðtþš ¼ qh 2 H 2 ðtþn 2 ½/ 2 ðt; qþš; ð1 qþl½/ 3 ðt; qþ R 0 ðtþš ¼ qh 3 H 3 ðtþn 3 ½/ 3 ðt; qþš; subject to the initial conditions / 1 ð0; qþ ¼S 0 ; / 2 ð0; qþ ¼I 0 ; / 3 ð0; qþ ¼R 0 : By Taylor s theorem, we expand / i ðt; qþ by a power series of the embedding parameter q as follows: / 1 ðt; qþ ¼S 0 ðtþþ X1 / 2 ðt; qþ ¼I 0 ðtþþ X1 / 3 ðt; qþ ¼R 0 ðtþþ X1 S m ðtþq m ; I m ðtþq m ; R m ðtþq m where S m ðtþ ¼ 1 o m / 1 ðt; qþ m! oq m ; I m ðtþ ¼ 1 o m / 2 ðt; qþ m! oq m ; R m ðtþ ¼ 1 o m / 3 ðt; qþ m! oq m : From the so-called mth-order deformation equations (12) and (13), we have L½S m ðtþ v m S m 1 ðtþš ¼ h 1 H 1 ðtþr m ðs m 1 ðtþþ; L½I m ðtþ v m I m 1 ðtþš ¼ h 2 H 2 ðtþr m ði m 1 ðtþþ; L½R m ðtþ v m R m 1 ðtþš ¼ h 3 H 3 ðtþr m ðr m 1 ðtþþ; S m ð0þ ¼0; I m ð0þ ¼0; R m ð0þ ¼0: ð14þ By the h-curves [6], it is reasonable to use h i ¼ 1. Using H i ðtþ ¼1, the mth-order deformation equation (14) for m P 1 becomes Z " # t Xm 1 S m ðtþ ¼v m S m 1 ðtþ ðsþþb S k ðsþi m 1 k ðsþ ds; I m ðtþ ¼v m I m 1 ðtþ R m ðtþ ¼v m R m 1 ðtþ S 0 m 1 0 Z t I 0 m 1 0 Z t 0 k¼0 " # Xm 1 ðsþ b S k ðsþi m 1 k ðsþþci m 1 ðsþ ds; k¼0 R 0 m 1 ðsþ ci m 1ðsÞ ds: 4. Numerical results and discussion For numerical results the following values, for parameters, are considered: N S N I N R b c
5 F. Awawdeh et al. / Chaos, Solitons and Fractals 42 (2009) Fig. 1. Pls of 20 terms approximations for SðtÞ; IðtÞ and RðtÞ versus time. Five and nine terms approximations for SðtÞ, IðtÞ and RðtÞ are calculated and presented below. Five terms approximations: S 5 ðtþ ¼ X5 I 5 ðtþ ¼ X5 R 5 ðtþ ¼ X5 Nine terms approximations: S 9 ðtþ ¼ X9 S m ðtþ ¼20 3t 0:045t 2 þ 0:028t 3 þ 0: t 4 0: t 5 ; I m ðtþ ¼15 þ 2:7t þ 0:018t 2 0:02817t 3 0: t 4 þ 0: t 5 ; R m ðtþ ¼10 þ 0:3t 0:054t 2 0:00012t 3 þ 0: t 4 0: t 5 : S m ðtþ ¼20 3t 0:045t 2 þ 0:028t 3 þ 0: t 4 0: t 5 0: t 6 þ 0: t 7 þ 0: t 8 0: t 9 ; I 9 ðtþ ¼ X9 I m ðtþ ¼15 þ 2:7t þ 0:018t 2 0:02817t 3 0: t 4 þ 0: t 5 þ 0: t 6 0: t 7 0: t 8 þ 0: t 9 ; R 9 ðtþ ¼ X9 R m ðtþ ¼10 þ 0:3t 0:054t 2 0:00012t 3 þ 0: t 4 0: t 5 þ 0: t 6 þ 0: t 7 0: t 8 0: t 9 : The HAM yields rapidly convergent series solution by using a few iterations. For the convergence of the HAM the reader is referred to [6]. Fig. 1 shows pls of 20 terms approximations of SðtÞ; IðtÞ and RðtÞ. Fig. 1 illustrates the case when we introduce a small number of infectives I 0 ¼ 1 into a susceptible population. An epidemic will occur and the number of infectives increases; the maximum infective population (climax of epidemic) I max ¼ 240:2564 will occur when S has decreased to the value 100. As time goes on (t!1) you travel along the curve to the right, eventually approaching I ¼ 0 and the disease died out. The epidemic will end as I! 0 with S approaching some positive value S 1 ¼ 3:4699. S 1 is the eventual population who were never infective. lim IðtÞ ¼lim I 20ðtÞ ¼0 and lim SðtÞ ¼lim S 20 ðtþ ¼S 1 ¼ 3:4699: t!1 t!1 t!1 t!1 5. Discussion and conclusion The analytical approximations to the solutions of the epidemic models are reliable and confirm the power and ability of the HAM as an easy device for computing the solution of nonlinear problems. The method avoids the difficulties and massive computational work that usually arise from parallel techniques and finite-difference method.
6 3052 F. Awawdeh et al. / Chaos, Solitons and Fractals 42 (2009) References [1] Biazar J. Solution of the epidemic model by Adomian decomposition method. Appl Math Comput 2006;173: [2] Bai Y, Jin Z. Prediction of SARS epidemic by BP neural networks with online prediction strategy. Chaos, Solitons & Fractals 2005;26: [3] Dietz K, Hadeler KP. Epidemiological models for sexually transmitted diseases. J Math Biol 1988;26:1 25. [4] Gao S, Teng Z, Xie D. Analysis of a delayed SIR epidemic model with pulse vaccination. Chaos, Solitons & Fractals 2009;40: [5] Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proc R Soc Lond A 1927;115: [6] Liao SJ. Beyond perturbation: introduction to the homopy analysis method. Boca Raton: Chapman & Hall/CRC Press; [7] Liao SJ. On the homopy analysis method for nonlinear problems. Appl Math Comput 2004;174: [8] Lipsitch M, Cohen T, Cooper B, Robins JM, Ma S, James L. Transmission dynamics and control of severe acute respiratory syndrome. Science 2003;300: [9] Liu Z, Lai YC, Ye N. Propagation and immunization of infection on general networks with bh homogeneous and heterogeneous components. Phys Rev 2003;E67: [10] Pang G, Chen L. A delayed SIRS epidemic model with pulse vaccination. Chaos, Solitons & Fractals 2007;34: [11] Riley S. Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions. Science 2003;300: [12] Satsuma J, Willox R, Ramani A, Grammaticos B, Carstea AS. Extending the SIR epidemic model. Physics A 2004: [13] Zanette DH. Critical behavior of propagation on small-world networks. Phys Rev 2001;E64:050901(R). [14] Zanette DH. Dynamics of rumor propagation on small-world networks. Phys Rev 2002;E65:
A new method for homoclinic solutions of ordinary differential equations
Chaos, Solitons and Fractals xxx (2007) xxx xxx www.elsevier.com/locate/chaos A new method for homoclinic solutions of ordinary differential equations F. Talay Akyildiz a, K. Vajravelu b, *, S.-J. Liao
More informationSeries solutions of non-linear Riccati differential equations with fractional order
Available online at www.sciencedirect.com Chaos, Solitons and Fractals 40 (2009) 1 9 www.elsevier.com/locate/chaos Series solutions of non-linear Riccati differential equations with fractional order Jie
More informationIntroduction to SEIR Models
Department of Epidemiology and Public Health Health Systems Research and Dynamical Modelling Unit Introduction to SEIR Models Nakul Chitnis Workshop on Mathematical Models of Climate Variability, Environmental
More informationA NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD
April, 4. Vol. 4, No. - 4 EAAS & ARF. All rights reserved ISSN35-869 A NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD Ahmed A. M. Hassan, S. H. Hoda Ibrahim, Amr M.
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx (2009) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A one-step optimal homotopy
More informationAnalytic solution of fractional integro-differential equations
Annals of the University of Craiova, Mathematics and Computer Science Series Volume 38(1), 211, Pages 1 1 ISSN: 1223-6934 Analytic solution of fractional integro-differential equations Fadi Awawdeh, E.A.
More informationAN AUTOMATIC SCHEME ON THE HOMOTOPY ANALYSIS METHOD FOR SOLVING NONLINEAR ALGEBRAIC EQUATIONS. Safwan Al-Shara
italian journal of pure and applied mathematics n. 37 2017 (5 14) 5 AN AUTOMATIC SCHEME ON THE HOMOTOPY ANALYSIS METHOD FOR SOLVING NONLINEAR ALGEBRAIC EQUATIONS Safwan Al-Shara Department of Mathematics
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 16 (2011) 2730 2736 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Homotopy analysis method
More informationAPPROXIMATING THE FORTH ORDER STRUM-LIOUVILLE EIGENVALUE PROBLEMS BY HOMOTOPY ANALYSIS METHOD
APPROXIMATING THE FORTH ORDER STRUM-LIOUVILLE EIGENVALUE PROBLEMS BY HOMOTOPY ANALYSIS METHOD * Nader Rafatimaleki Department of Mathematics, College of Science, Islamic Azad University, Tabriz Branch,
More informationAnalysis of Numerical and Exact solutions of certain SIR and SIS Epidemic models
Journal of Mathematical Modelling and Application 2011, Vol. 1, No. 4, 51-56 ISSN: 2178-2423 Analysis of Numerical and Exact solutions of certain SIR and SIS Epidemic models S O Maliki Department of Industrial
More informationNewton-homotopy analysis method for nonlinear equations
Applied Mathematics and Computation 188 (2007) 1794 1800 www.elsevier.com/locate/amc Newton-homotopy analysis method for nonlinear equations S. Abbasbandy a, *, Y. Tan b, S.J. Liao b a Department of Mathematics,
More informationEpidemics in Networks Part 2 Compartmental Disease Models
Epidemics in Networks Part 2 Compartmental Disease Models Joel C. Miller & Tom Hladish 18 20 July 2018 1 / 35 Introduction to Compartmental Models Dynamics R 0 Epidemic Probability Epidemic size Review
More informationSolving a modified nonlinear epidemiological model of computer viruses by homotopy analysis method
Mathematical Sciences (201) 12:211 222 https://doi.org/10.1007/s40096-01-0261-5 (045679().,-volV)(045679().,-volV) ORIGINAL PAPER Solving a modified nonlinear epidemiological model of computer viruses
More informationMathematical Analysis of Epidemiological Models: Introduction
Mathematical Analysis of Epidemiological Models: Introduction Jan Medlock Clemson University Department of Mathematical Sciences 8 February 2010 1. Introduction. The effectiveness of improved sanitation,
More informationOptimal Control of a Delayed SIRS Epidemic Model with Vaccination and Treatment
Acta Biotheor DOI.7/s44-5-9244- REGULAR ARTICLE Optimal Control of a Delayed SIRS Epidemic Model with Vaccination and Treatment Hassan Laarabi Abdelhadi Abta Khalid Hattaf Received: 2 February 24 / Accepted:
More informationFixed Point Analysis of Kermack Mckendrick SIR Model
Kalpa Publications in Computing Volume, 17, Pages 13 19 ICRISET17. International Conference on Research and Innovations in Science, Engineering &Technology. Selected Papers in Computing Fixed Point Analysis
More informationMULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS
MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS Hossein Jafari & M. A. Firoozjaee Young Researchers club, Islamic Azad University, Jouybar Branch, Jouybar, Iran
More informationA Note on the Spread of Infectious Diseases. in a Large Susceptible Population
International Mathematical Forum, Vol. 7, 2012, no. 50, 2481-2492 A Note on the Spread of Infectious Diseases in a Large Susceptible Population B. Barnes Department of Mathematics Kwame Nkrumah University
More informationMATHEMATICAL MODELS Vol. III - Mathematical Models in Epidemiology - M. G. Roberts, J. A. P. Heesterbeek
MATHEMATICAL MODELS I EPIDEMIOLOGY M. G. Roberts Institute of Information and Mathematical Sciences, Massey University, Auckland, ew Zealand J. A. P. Heesterbeek Faculty of Veterinary Medicine, Utrecht
More informationForecast and Control of Epidemics in a Globalized World
Forecast and Control of Epidemics in a Globalized World L. Hufnagel, D. Brockmann, T. Geisel Max-Planck-Institut für Strömungsforschung, Bunsenstrasse 10, 37073 Göttingen, Germany and Kavli Institute for
More informationMathematical Epidemiology Lecture 1. Matylda Jabłońska-Sabuka
Lecture 1 Lappeenranta University of Technology Wrocław, Fall 2013 What is? Basic terminology Epidemiology is the subject that studies the spread of diseases in populations, and primarily the human populations.
More informationAnalytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy Analysis Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(21) No.4,pp.414-421 Analytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy
More informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationFigure The Threshold Theorem of epidemiology
K/a Figure 3 6. Assurne that K 1/ a < K 2 and K 2 / ß < K 1 (a) Show that the equilibrium solution N 1 =0, N 2 =0 of (*) is unstable. (b) Show that the equilibrium solutions N 2 =0 and N 1 =0, N 2 =K 2
More informationModelling of the Hand-Foot-Mouth-Disease with the Carrier Population
Modelling of the Hand-Foot-Mouth-Disease with the Carrier Population Ruzhang Zhao, Lijun Yang Department of Mathematical Science, Tsinghua University, China. Corresponding author. Email: lyang@math.tsinghua.edu.cn,
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More information(mathematical epidemiology)
1. 30 (mathematical epidemiology) 2. 1927 10) * Anderson and May 1), Diekmann and Heesterbeek 3) 7) 14) NO. 538, APRIL 2008 1 S(t), I(t), R(t) (susceptibles ) (infectives ) (recovered/removed = βs(t)i(t)
More informationDynamic pair formation models
Application to sexual networks and STI 14 September 2011 Partnership duration Models for sexually transmitted infections Which frameworks? HIV/AIDS: SI framework chlamydia and gonorrhoea : SIS framework
More informationGlobal Analysis of an Epidemic Model with Nonmonotone Incidence Rate
Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate Dongmei Xiao Department of Mathematics, Shanghai Jiaotong University, Shanghai 00030, China E-mail: xiaodm@sjtu.edu.cn and Shigui Ruan
More informationMathematical Modeling and Analysis of Infectious Disease Dynamics
Mathematical Modeling and Analysis of Infectious Disease Dynamics V. A. Bokil Department of Mathematics Oregon State University Corvallis, OR MTH 323: Mathematical Modeling May 22, 2017 V. A. Bokil (OSU-Math)
More informationSpotlight on Modeling: The Possum Plague
70 Spotlight on Modeling: The Possum Plague Reference: Sections 2.6, 7.2 and 7.3. The ecological balance in New Zealand has been disturbed by the introduction of the Australian possum, a marsupial the
More informationChaos, Solitons and Fractals
Chaos, Solitons and Fractals 41 (2009) 962 969 Contents lists available at ScienceDirect Chaos, Solitons and Fractals journal homepage: www.elsevier.com/locate/chaos A fractional-order hyperchaotic system
More informationA comparison of delayed SIR and SEIR epidemic models
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 2, 181 190 181 A comparison of delayed SIR and SEIR epidemic models Abdelilah Kaddar a, Abdelhadi Abta b, Hamad Talibi Alaoui b a Université
More informationAn Analytical Scheme for Multi-order Fractional Differential Equations
Tamsui Oxford Journal of Mathematical Sciences 26(3) (2010) 305-320 Aletheia University An Analytical Scheme for Multi-order Fractional Differential Equations H. M. Jaradat Al Al Bayt University, Jordan
More informationNon-Linear Models Cont d: Infectious Diseases. Non-Linear Models Cont d: Infectious Diseases
Cont d: Infectious Diseases Infectious Diseases Can be classified into 2 broad categories: 1 those caused by viruses & bacteria (microparasitic diseases e.g. smallpox, measles), 2 those due to vectors
More informationAn analytic approach to solve multiple solutions of a strongly nonlinear problem
Applied Mathematics and Computation 169 (2005) 854 865 www.elsevier.com/locate/amc An analytic approach to solve multiple solutions of a strongly nonlinear problem Shuicai Li, Shi-Jun Liao * School of
More informationOn the convergence of the homotopy analysis method to solve the system of partial differential equations
Journal of Linear and Topological Algebra Vol. 04, No. 0, 015, 87-100 On the convergence of the homotopy analysis method to solve the system of partial differential equations A. Fallahzadeh a, M. A. Fariborzi
More informationPARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 4, Winter 211 PARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS Dedicated to Herb Freedman on the occasion of his seventieth birthday
More informationResearch Article On a New Reliable Algorithm
Hindawi Publishing Corporation International Journal of Differential Equations Volume 2009, Article ID 710250, 13 pages doi:10.1155/2009/710250 Research Article On a New Reliable Algorithm A. K. Alomari,
More informationSpread of Malicious Objects in Computer Network: A Fuzzy Approach
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 8, Issue 2 (December 213), pp. 684 7 Applications and Applied Mathematics: An International Journal (AAM) Spread of Malicious Objects
More informationStability of SEIR Model of Infectious Diseases with Human Immunity
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 1811 1819 Research India Publications http://www.ripublication.com/gjpam.htm Stability of SEIR Model of Infectious
More informationHOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (21), 89 98 HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS Hossein Jafari and M. A. Firoozjaee Abstract.
More informationUnderstanding the contribution of space on the spread of Influenza using an Individual-based model approach
Understanding the contribution of space on the spread of Influenza using an Individual-based model approach Shrupa Shah Joint PhD Candidate School of Mathematics and Statistics School of Population and
More informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
More informationA Mathematical Analysis on the Transmission Dynamics of Neisseria gonorrhoeae. Yk j N k j
North Carolina Journal of Mathematics and Statistics Volume 3, Pages 7 20 (Accepted June 23, 2017, published June 30, 2017 ISSN 2380-7539 A Mathematical Analysis on the Transmission Dynamics of Neisseria
More informationNumerical comparison of methods for solving linear differential equations of fractional order
Chaos, Solitons and Fractals 31 (27) 1248 1255 www.elsevier.com/locate/chaos Numerical comparison of methods for solving linear differential equations of fractional order Shaher Momani a, *, Zaid Odibat
More informationStability Analysis of a SIS Epidemic Model with Standard Incidence
tability Analysis of a I Epidemic Model with tandard Incidence Cruz Vargas-De-León Received 19 April 2011; Accepted 19 Octuber 2011 leoncruz82@yahoo.com.mx Abstract In this paper, we study the global properties
More informationME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included
ME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included sysid Mathematica 6.0.3, DynPac 11.01, 1ê13ê9 1. Introduction Description of the Model In this notebook, we include births and deaths in the
More informationEpidemics in Complex Networks and Phase Transitions
Master M2 Sciences de la Matière ENS de Lyon 2015-2016 Phase Transitions and Critical Phenomena Epidemics in Complex Networks and Phase Transitions Jordan Cambe January 13, 2016 Abstract Spreading phenomena
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 6 (2) 63 75 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns The scaled boundary FEM for nonlinear
More informationA New Mathematical Approach for. Rabies Endemy
Applied Mathematical Sciences, Vol. 8, 2014, no. 2, 59-67 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39525 A New Mathematical Approach for Rabies Endemy Elif Demirci Ankara University
More informationInternational Journal of Modern Mathematical Sciences, 2012, 3(2): International Journal of Modern Mathematical Sciences
Article International Journal of Modern Mathematical Sciences 2012 3(2): 63-76 International Journal of Modern Mathematical Sciences Journal homepage:wwwmodernscientificpresscom/journals/ijmmsaspx On Goursat
More informationThe death of an epidemic
LECTURE 2 Equilibrium Stability Analysis & Next Generation Method The death of an epidemic In SIR equations, let s divide equation for dx/dt by dz/ dt:!! dx/dz = - (β X Y/N)/(γY)!!! = - R 0 X/N Integrate
More informationDynamics of Disease Spread. in a Predator-Prey System
Advanced Studies in Biology, vol. 6, 2014, no. 4, 169-179 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/asb.2014.4845 Dynamics of Disease Spread in a Predator-Prey System Asrul Sani 1, Edi Cahyono
More informationA numerical analysis of a model of growth tumor
Applied Mathematics and Computation 67 (2005) 345 354 www.elsevier.com/locate/amc A numerical analysis of a model of growth tumor Andrés Barrea a, *, Cristina Turner b a CIEM, Universidad Nacional de Córdoba,
More informationComparison of homotopy analysis method and homotopy perturbation method through an evolution equation
Comparison of homotopy analysis method and homotopy perturbation method through an evolution equation Songxin Liang, David J. Jeffrey Department of Applied Mathematics, University of Western Ontario, London,
More informationThe Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
Applied Mathematics, 05, 6, 665-675 Published Online September 05 in SciRes http://wwwscirporg/journal/am http://dxdoiorg/046/am056048 The Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
More informationMathematical modelling and controlling the dynamics of infectious diseases
Mathematical modelling and controlling the dynamics of infectious diseases Musa Mammadov Centre for Informatics and Applied Optimisation Federation University Australia 25 August 2017, School of Science,
More informationJournal of Engineering Science and Technology Review 2 (1) (2009) Research Article
Journal of Engineering Science and Technology Review 2 (1) (2009) 118-122 Research Article JOURNAL OF Engineering Science and Technology Review www.jestr.org Thin film flow of non-newtonian fluids on a
More informationPreservation of local dynamics when applying central difference methods: application to SIR model
Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER
More informationHomotopy Perturbation Method for the Fisher s Equation and Its Generalized
ISSN 749-889 (print), 749-897 (online) International Journal of Nonlinear Science Vol.8(2009) No.4,pp.448-455 Homotopy Perturbation Method for the Fisher s Equation and Its Generalized M. Matinfar,M. Ghanbari
More informationHomotopy Analysis Transform Method for Time-fractional Schrödinger Equations
International Journal of Modern Mathematical Sciences, 2013, 7(1): 26-40 International Journal of Modern Mathematical Sciences Journal homepage:wwwmodernscientificpresscom/journals/ijmmsaspx ISSN:2166-286X
More informationObservation and model error effects on parameter estimates in susceptible-infected-recovered epidemiological models
Georgia State University ScholarWorks @ Georgia State University Public Health Faculty Publications School of Public Health 26 Observation and model error effects on parameter estimates in susceptible-infected-recovered
More informationSection 8.1 Def. and Examp. Systems
Section 8.1 Def. and Examp. Systems Key Terms: SIR Model of an epidemic o Nonlinear o Autonomous Vector functions o Derivative of vector functions Order of a DE system Planar systems Dimension of a system
More informationECS 289 / MAE 298, Lecture 7 April 22, Percolation and Epidemiology on Networks, Part 2 Searching on networks
ECS 289 / MAE 298, Lecture 7 April 22, 2014 Percolation and Epidemiology on Networks, Part 2 Searching on networks 28 project pitches turned in Announcements We are compiling them into one file to share
More informationHETEROGENEOUS MIXING IN EPIDEMIC MODELS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number 1, Spring 212 HETEROGENEOUS MIXING IN EPIDEMIC MODELS FRED BRAUER ABSTRACT. We extend the relation between the basic reproduction number and the
More informationSocial Influence in Online Social Networks. Epidemiological Models. Epidemic Process
Social Influence in Online Social Networks Toward Understanding Spatial Dependence on Epidemic Thresholds in Networks Dr. Zesheng Chen Viral marketing ( word-of-mouth ) Blog information cascading Rumor
More informationElectronic appendices are refereed with the text. However, no attempt has been made to impose a uniform editorial style on the electronic appendices.
This is an electronic appendix to the paper by Alun L. Lloyd 2001 Destabilization of epidemic models with the inclusion of realistic distributions of infectious periods. Proc. R. Soc. Lond. B 268, 985-993.
More informationNetwork Modelling for Sexually. Transmitted Diseases
Network Modelling for Sexually Transmitted Diseases Thitiya Theparod, MSc Bsc (Hons) Submitted for the degree of Doctor of Philosophy at Lancaster University 2015 Abstract The aim of this thesis is to
More informationChaos synchronization of nonlinear Bloch equations
Chaos, Solitons and Fractal7 (26) 357 361 www.elsevier.com/locate/chaos Chaos synchronization of nonlinear Bloch equations Ju H. Park * Robust Control and Nonlinear Dynamics Laboratory, Department of Electrical
More informationarxiv: v1 [math.ds] 16 Oct 2014
Application of Homotopy Perturbation Method to an Eco-epidemic Model arxiv:1410.4385v1 [math.ds] 16 Oct 014 P.. Bera (1, S. Sarwardi ( and Md. A. han (3 (1 Department of Physics, Dumkal College, Basantapur,
More informationResearch Article A Delayed Epidemic Model with Pulse Vaccination
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2008, Article ID 746951, 12 pages doi:10.1155/2008/746951 Research Article A Delayed Epidemic Model with Pulse Vaccination
More informationThree Disguises of 1 x = e λx
Three Disguises of 1 x = e λx Chathuri Karunarathna Mudiyanselage Rabi K.C. Winfried Just Department of Mathematics, Ohio University Mathematical Biology and Dynamical Systems Seminar Ohio University November
More information(Received 1 February 2012, accepted 29 June 2012) address: kamyar (K. Hosseini)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.14(2012) No.2,pp.201-210 Homotopy Analysis Method for a Fin with Temperature Dependent Internal Heat Generation
More informationChapter 4 An Introduction to Networks in Epidemic Modeling
Chapter 4 An Introduction to Networks in Epidemic Modeling Fred Brauer Abstract We use a stochastic branching process to describe the beginning of a disease outbreak. Unlike compartmental models, if the
More informationV. G. Gupta 1, Pramod Kumar 2. (Received 2 April 2012, accepted 10 March 2013)
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9(205 No.2,pp.3-20 Approimate Solutions of Fractional Linear and Nonlinear Differential Equations Using Laplace Homotopy
More informationA Mathematical Model for the Spatial Spread of HIV in a Heterogeneous Population
A Mathematical Model for the Spatial Spread of HIV in a Heterogeneous Population Titus G. Kassem * Department of Mathematics, University of Jos, Nigeria. Abstract Emmanuel J.D. Garba Department of Mathematics,
More informationOn epidemic models with nonlinear cross-diffusion
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 On epidemic models with nonlinear cross-diffusion S. Berres a and J. Gonzalez-Marin
More informationSynchronization of an uncertain unified chaotic system via adaptive control
Chaos, Solitons and Fractals 14 (22) 643 647 www.elsevier.com/locate/chaos Synchronization of an uncertain unified chaotic system via adaptive control Shihua Chen a, Jinhu L u b, * a School of Mathematical
More informationANALYTICAL APPROXIMATE SOLUTIONS OF THE ZAKHAROV-KUZNETSOV EQUATIONS
(c) Romanian RRP 66(No. Reports in 2) Physics, 296 306 Vol. 2014 66, No. 2, P. 296 306, 2014 ANALYTICAL APPROXIMATE SOLUTIONS OF THE ZAKHAROV-KUZNETSOV EQUATIONS A. JAFARIAN 1, P. GHADERI 2, ALIREZA K.
More informationSI j RS E-Epidemic Model With Multiple Groups of Infection In Computer Network. 1 Introduction. Bimal Kumar Mishra 1, Aditya Kumar Singh 2
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(2012) No.3,pp.357-362 SI j RS E-Epidemic Model With Multiple Groups of Infection In Computer Network Bimal Kumar
More informationDynamical models of HIV-AIDS e ect on population growth
Dynamical models of HV-ADS e ect on population growth David Gurarie May 11, 2005 Abstract We review some known dynamical models of epidemics, given by coupled systems of di erential equations, and propose
More informationGLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS
CANADIAN APPIED MATHEMATICS QUARTERY Volume 13, Number 4, Winter 2005 GOBA DYNAMICS OF A MATHEMATICA MODE OF TUBERCUOSIS HONGBIN GUO ABSTRACT. Mathematical analysis is carried out for a mathematical model
More informationOn the Spread of Epidemics in a Closed Heterogeneous Population
On the Spread of Epidemics in a Closed Heterogeneous Population Artem Novozhilov Applied Mathematics 1 Moscow State University of Railway Engineering (MIIT) the 3d Workshop on Mathematical Models and Numerical
More informationEpidemics and information spreading
Epidemics and information spreading Leonid E. Zhukov School of Data Analysis and Artificial Intelligence Department of Computer Science National Research University Higher School of Economics Social Network
More informationBehavior Stability in two SIR-Style. Models for HIV
Int. Journal of Math. Analysis, Vol. 4, 2010, no. 9, 427-434 Behavior Stability in two SIR-Style Models for HIV S. Seddighi Chaharborj 2,1, M. R. Abu Bakar 2, I. Fudziah 2 I. Noor Akma 2, A. H. Malik 2,
More informationComputers and Mathematics with Applications. A new application of He s variational iteration method for the solution of the one-phase Stefan problem
Computers and Mathematics with Applications 58 (29) 2489 2494 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A new
More informationA Time Since Recovery Model with Varying Rates of Loss of Immunity
Bull Math Biol (212) 74:281 2819 DOI 1.17/s11538-12-978-7 ORIGINAL ARTICLE A Time Since Recovery Model with Varying Rates of Loss of Immunity Subhra Bhattacharya Frederick R. Adler Received: 7 May 212
More informationDynamical behavior of a pest management model with impulsive effect and nonlinear incidence rate*
Volume 30, N. 2, pp. 381 398, 2011 Copyright 2011 SBMAC ISSN 0101-8205 www.scielo.br/cam Dynamical behavior of a pest management model with impulsive effect and nonlinear incidence rate* XIA WANG 1, ZHEN
More informationApplied Mathematics Letters
Applied athematics Letters 25 (212) 156 16 Contents lists available at SciVerse ScienceDirect Applied athematics Letters journal homepage: www.elsevier.com/locate/aml Globally stable endemicity for infectious
More informationBifurcation Analysis in Simple SIS Epidemic Model Involving Immigrations with Treatment
Appl. Math. Inf. Sci. Lett. 3, No. 3, 97-10 015) 97 Applied Mathematics & Information Sciences Letters An International Journal http://dx.doi.org/10.1785/amisl/03030 Bifurcation Analysis in Simple SIS
More informationHomotopy Analysis Method for Nonlinear Jaulent-Miodek Equation
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol. 5, No.,, pp. 8-88 Hootopy Analysis Method for Nonlinear Jaulent-Miodek Equation J. Biazar, M. Eslai Departent of Matheatics, Faculty
More informationAnalysis with Variable Step Size Strategy of Some SIR Epidemic Models
Karaelmas Fen ve Müh. Derg. 6(1):203-210, 2016 Karaelmas Fen ve Mühendislik Dergisi Journal home page: http://fbd.beun.edu.tr Research Article Analysis with Variable Step Size Strategy of Some SIR Epidemic
More informationLecture VI Introduction to complex networks. Santo Fortunato
Lecture VI Introduction to complex networks Santo Fortunato Plan of the course I. Networks: definitions, characteristics, basic concepts in graph theory II. III. IV. Real world networks: basic properties
More informationModels of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor. August 15, 2005
Models of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor August 15, 2005 1 Outline 1. Compartmental Thinking 2. Simple Epidemic (a) Epidemic Curve 1:
More informationHOMOTOPY ANALYSIS METHOD FOR SOLVING COUPLED RAMANI EQUATIONS
HOMOTOPY ANALYSIS METHOD FOR SOLVING COUPLED RAMANI EQUATIONS A. JAFARIAN 1, P. GHADERI 2, ALIREZA K. GOLMANKHANEH 3, D. BALEANU 4,5,6 1 Department of Mathematics, Uremia Branch, Islamic Azan University,
More informationImproving homotopy analysis method for system of nonlinear algebraic equations
Journal of Advanced Research in Applied Mathematics Vol., Issue. 4, 010, pp. -30 Online ISSN: 194-9649 Improving homotopy analysis method for system of nonlinear algebraic equations M.M. Hosseini, S.M.
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More information