Epidemiology and Chaos Jonathan Wills Math 053: Chaos! Prof. Alex Barnett Fall Abstract:

Size: px
Start display at page:

Download "Epidemiology and Chaos Jonathan Wills Math 053: Chaos! Prof. Alex Barnett Fall Abstract:"

Transcription

1 Abstract: Epidemiology and Chaos Jonathan Wills Math 053: Chaos! Prof. Alex Barnett Fall 2009 Mathematical epidemiological models for the spread of disease through a population are used to predict the prevalence of a disease or to study the impacts of treatment or prevention measures. Initial conditions for these models are measured from statistical data collected from a population since these initial conditions can never be exact, the presence of chaos in mathematical models has serious implications for the accuracy of the models as well as how epidemiologists interpret their findings. This paper confirms the chaotic behavior of a model for dengue fever by investigating sensitive depence, chaotic attractors, and Lyapunov exponents under a variety of initial conditions. Mathematical Epidemiological Models: Mathematical models for the spread of disease through a population have a variety of uses. A good mathematical model can be used to predict the prevalence of a disease within a population. Different parameters in the model reflect the average lifespan of a population, the survival rate for a disease, the infectiousness of a disease, the rate of loss of immunity over time, the rate of recovery from the disease, or a number of other factors that vary by disease. Varying these constants allows epidemiologists to approximate the effects of prevention measures such as vaccinations or treatment measures. Models also allow epidemiologists to predict whether a disease or virus will die out over time under certain conditions. From the models, mathematicians can calculate reproduction numbers and other measures of how much the disease will spread at a given time, or threshold values necessary for certain conditions such as eradication to be achieved (Hethcote). Finding which conditions will achieve eradication lets epidemiologists know where it will be most effective to concentrate their efforts. More advance models differentiate between age groups or seasonal variations in disease prevalence caused by climate changes or school terms. The standard epidemiological model for disease prevalence is the MSEIR model, a set of differential equations that categorizes percentages of the population, as shown in the figure below (Hethcote). The M portion of the population represents those infants who have temporary maternal passive immunity from the disease. Over time, those in M will move in to the S group, those susceptible to the disease. If a susceptible individual has contact with an infectious individual and the contact is sufficient to contract the disease, then the susceptible one will move into E, the exposed group. Those in E will sp a certain amount of time specific to the disease before the start to show symptoms and become infectious, at which point they will move into the I group. If they survive the disease, they recover and move into the R group, and will remain immune to the disease for an amount of time specific to each disease. In the MSEIR model, it is assumed for the sake of mathematical simplicity that the birth rate is the same as the death rate; the total population stays constant and there is 1

2 constant input into the M and S categories. More complicated mathematical models include populations that vary over time. MSEIR represents a general case that can be adapted to better reflect many specific diseases. For some cases, the time between exposure and display of symptoms or the ability of mothers to transfer antibodies to their fetuses may be negligible, in which case the E or M group would be omitted from the model. Different parameters in each model will further reflect the specific properties of a single disease. Limitations: While these models are powerful tools for predicting the behavior of diseases within a population, it is important to remember that they are drastic simplifications of reality. Disease dynamics in real life are not deterministic human behavior can affect outcomes for any set of initial conditions. These models are merely approximations of what should happen for sufficiently large populations. They assume constant total population (although they account for births and deaths at an equal rate) and they assume that the population will mix homogeneously (that each individual will come into contact with the same constant number of other individuals). It is also impossible to obtain perfect information about the constants that fit into the model to characterize the disease. Statistical information about real diseases is gathered, but one must always account for unreported cases of the disease or instances of the disease that have been affected by other variables. More complicated models require more and more parameters to fully model the spread of a disease; when choosing how to study a disease, a mathematician or epidemiologist must balance the accuracy of the mathematical model against the feasibility of obtaining all the necessary parameters to fully characterize the a disease. Chaos and sensitive depence further complicate this dynamic. Epidemiologists can only approximate initial conditions, and the possibility that two very similar sets of initial conditions could have results that diverge at an exponential rate has serious ramifications for the accuracy of their models. The production of vaccines takes a significant amount of time, and which diseases or strains of a virus the mathematical models predict will be most prevalent in the future and will be most affected by mass vaccinations govern their production. 2

3 Chaos in Dengue Fever: I investigated chaos in an adaptation of the SIRS model (shown to the right, and taken from a study by Aguiar, Kooi, and Stollenwerk) for two coexisting strains of dengue fever. Dengue is essentially a human disease, but is passed through a population by transfer of blood, most often by mosquitoes. In many epidemiological models for dengue fever, it is assumed that mosquitoes are mixing homogeneously with the human population (Derouich et al). Note also that in this case, the creators of the model deemed the time and number of infants with maternal passive immunity and the length of time between exposure and display of symptoms to be negligible; hence the M and E groups have been omitted. In this model, those infected with one strain or another of dengue will temporarily have cross immunity to the other strain. Those with the first strain cannot get the second until a certain period of time after recovering from the first, and vice versa. Chaos has been found in this model, as seen in some regions of the bifurcation diagram to the right, obtained by varying the parameter (Aguiar). Mathematical analysis was done using Matlab. Total population was set at 100 so that results would reflect percentages of the population rather than absolute numbers. The other parameters that approximate those for dengue fever are as follows: µ = 1/65, = 52, = 2, =2, and the parameter, a measure of infectiousness that can reflect what proportion of the population is vaccinated, will be varied. I first kept constant at 0.9 to search for sensitive depence. After finding an orbit that appeared to be chaotic, I mapped two trajectories with very similar initial conditions (uo = [ ; ; ; ; ; ; ; ; ; ] and uo2 = uo + 1e-4 * [1; zeros(8,1); -1]). The initial distance between the two sets of starting values was e-04. These trajectories diverged at an exponential rate, as is apparent in the figure below. 3

4 To calculate the rate of divergence, I mapped time against the log of the distance between the two trajectories (results below). There is a clear upward slope followed by a plateau at which the trajectories cannot diverge any more. I measured the slope of the region before the plateau from time 10 until the first time that the distance between the two trajectories was greater than 1. This slope gave me an approximation for the leading Lyapunov exponent. To make my approximation more accurate, I repeated the calculations for 15 different sets of initial values and averaged the results. This yielded a weakly positive average leading Lyapunov exponent of h = when = 0.9. To further investigate the effects of the parameter, I repeated this estimate of the leading Lyapunov exponent for values of between 0 and 3. Each reported leading Lyapunov exponent is the average of the calculated Lyapunov exponent for 15 different initial conditions. My results are shown below on the right. On the left for comparison are four Lyapunov exponents for the same model in the smaller range 0 < < 1, calculated by Aguiar, Kooi, and Stollenwerk with very different methods. My calculations seem to confirm the small peak of 0.4 < h < 0.6 for close to 0.8, but does not capture the sharp drop-off of h as decreases below 0.5 or the downward spikes in the 0.5 < < 1 range, most likely because I only calculated h in 0.1 intervals of. My calculations also seem to confirm the chaos found by Aguiar et al. in the range 2 < < 3, but my model also finds some chaos in the 1.5 < < 2 range, which may be erroneous. 4

5 My results, shown to the right above, give an approximation for the Lyapunov exponent for a given, but suffer from some limitations. First, I was only able to average 15 different initial conditions for each, and I calculated the Lyapunov exponent only for values that are multiples of 0.1; thus, the model has limited predictive capability for intermediate values of between those calculated. Furthermore, there is no guarantee that the chosen orbits will remain chaotic when is changed. The Lyapunov exponents remain weakly positive for the majority of values, but the calculated Lyapunov exponent could be influenced positively or negatively by my methods for calculations of slope. It should be used only to draw general conclusions about the positive Lyapunov exponents of the model. Having identified chaos in the model, I wanted to characterize it further. Below is the chaotic attractor for the first initial conditions used in the calculations for Lyapunov exponent. It is shown on a graph of total percent of susceptible individuals in the population versus total percentage of infectious individuals in the population. Percent recovered in the population is omitted because it can be determined by subtracting infectious and susceptible populations from the total population. Graphed in three dimensions, the strange attractor lies on a plane for which the total population adds up to 100%; therefore, graphing the chaotic attractor in two dimensions still retains all information. Conclusions: My findings confirmed and expanded upon the results found by other researchers of this model with other methods. My findings cannot conclusively confirm chaos in the dynamics of actual dengue fever, ruled by a system involving complex and changing interactions and behaviors that are not necessarily determined by initial conditions. 5

6 Rather it confirms the presence of chaos for many values of in this particular model for the disease, calling into question the predictive capabilities of the model and raising the possibility of sensitive depence in the actual disease dynamics to the extent that the disease reflects the model. A better understanding of the chaos in this and other models for epidemics can allow epidemiologists to calculate threshold values for asymptotic stability to approach eradication of a disease, and to better understand the limitations the models have and the effects their treatment programs might have. 6

7 Sources: Aguiar, Maíra, et al. Multi-Strain Deterministic Chaos in Dengue Epidemiology, A Challenge for Computational Mathematics. Numerical Analysis and Applied Mathematics International Conference. Vol 2. Eds. Simos, T. et al. American Institute of Physics, Web. 2 Dec Derouich, M, et al. A Model of dengue fever. PubMed: Biomed Eng Online. Vol 2, Web 2 Dec < Feng, Zhilan, and Jorge Velasco-Hernández. Competitive Exclusion in a vector-host model for the dengue fever. Journal of Mathematical Biology. Vol 35: pp , Web. 2 Dec Grenfell, B. T. Chance and Chaos in Measles Dynamics. Journal of the Royal Statistical Society 54.2 (1992): JSTOR. Web. 10 Nov Hethcote, Herbert. The Mathematics of Infectious Diseases. SIAM Review 42.4 (2000): JSTOR. Web. 10 Nov Hethcote, Herbert. The Basic Epidemiology Models The Mathematical Understanding of Infectious Disease Dynamics. Eds. Stefan Ma and Yingcun Xia. University of Iowa, Iowa City: Web. 1 Dec < Kleczkowski, A. Nonlinear Dynamics and Chaos in Infectious Diseases. University of Cambridge. Web. 10 Nov fqmat/atelier/ddd/1091.doc+seirs+epidemiology+chaos&cd=6&hl=en&ct=clnk &gl=us&client=safari Li, Michael, James Muldowney, and P. Van Den Driessche. Global Stability of SEIRS Model in Epidemiology. Canadian Applied Mathematics Quarterly 7.4 (1999): < Schaffer, W. M. Chaos in Ecology and Epidemiology. Chaos in Biological Systems. Ed. Degn, H., A. V. Holden, and L. F. Olsen. New York: Plenium Press,

8 Appix: Matlab Code First code is for plotting divergent trajectories, plotting chaotic attractor in 3d and 2d, and plotting graph for calculations of Lyapunov exponents: %Searching for Chaos in the MSEIR model for infectious diseases %Jonathan Wills, Fall 2009 %Dengue fever is a multi-strain system that can be modeled by a permutation %of the SIR model: here are the differential eqns. %Constants: N = 100; %puts values in percentages of total population M = (1/65); %years^(-1) %average life span/death rate Y = 52; %years^(-1) rate of recovery from infection, In- >Rn, where n=1 or n=2 B = 2*Y; %infection rate from either of groups I1 or I2 A = 2; %years^(-1) %rate of loss of cross immunity from R1->S1 %the constant Q will be varied. %the rate of infection of S1 by I12 is Q*B if 0 %Differential Equations: S -B/N*S*(I1+Q*I21)-B/N*S*(I2+Q*I12)+M*(N-S); %susceptible I1 B/N*S*(I1+Q*I21)-(Y+M)*I1; %infected with strain 1 but not strain 2 I2 B/N*S*(I2+Q*I12)-(Y+M)*I2; R1 Y*I1-(A+M)*R1; %recovered from strain 1, temporarily cross immune R2 Y*I2-(A+M)*R2; S1 -B/N*S1*(I2+Q*I12)+A*R1-M*S1; %immune to strain 1 but susceptible to strain 2 S2 -B/N*S2*(I1+Q*I21)+A*R2-M*S2; %and vice versa I12 B/N*S1*(I2+Q*I12)-(Y+M)*I12; %rezinfected from S1 with strain 2 I21 B/N*S2*(I1+Q*I21)-(Y+M)*I21; R Y*(I12+I21)-M*R; %recovered from both strains %now comes the differential equations part %eps = 0.1; %F [S; I1; I2; R1; R2; S1; S2; I12; I21; R]; %commented out the above functions because it was simpler to write them all %out as one function F below. Q = 0.9; %varied parameter eps = 0.1; F [-B/N*u(1)*(u(2)+Q*u(9))-B/N*u(1)*(u(3)+Q*u(8))+M*(N-u(1)); 8

9 B/N*u(1)*(u(2)+Q*u(9))-(Y+M)*u(2); B/N*u(1)*(u(3)+Q*u(8))-(Y+M)*u(3); Y*u(2)-(A+M)*u(4); Y*u(3)-(A+M)*u(5); -B/N*u(6)*(u(3)+Q*u(8))+A*u(4)- M*u(6); -B/N*u(7)*(u(2)+Q*u(9))+A*u(5)-M*u(7); B/N*u(6)*(u(3)+Q*u(8))- (Y+M)*u(8); B/N*u(7)*(u(2)+Q*u(9))-(Y+M)*u(9); Y*(u(8)+u(9))-M*u(10)]; uo = [ ; ; ; ; ; ; ; ; ; ]; %uo = [27;16;14;12;10;8;6;4;2;1]; %uo = [27;16;14;22;0;12;6;0;2;1]; %uo = [ ;0;0.0119;0.0048;0.0538; ;9.0077;0.0088;0; ]; %uo = [97;1;1;0;0;0;0;0;0;1]; s1 = ode45(f,[0 500],uo,odeset('abstol',1e-9,'nonnegative',1:10)); ts1 = 1:0.01:500; us1 = deval(s1,ts1)'; eps = 0.1; F [-B/N*u(1)*(u(2)+Q*u(9))-B/N*u(1)*(u(3)+Q*u(8))+M*(N-u(1)); B/N*u(1)*(u(2)+Q*u(9))-(Y+M)*u(2); B/N*u(1)*(u(3)+Q*u(8))-(Y+M)*u(3); Y*u(2)-(A+M)*u(4); Y*u(3)-(A+M)*u(5); -B/N*u(6)*(u(3)+Q*u(8))+A*u(4)- M*u(6); -B/N*u(7)*(u(2)+Q*u(9))+A*u(5)-M*u(7); B/N*u(6)*(u(3)+Q*u(8))- (Y+M)*u(8); B/N*u(7)*(u(2)+Q*u(9))-(Y+M)*u(9); Y*(u(8)+u(9))-M*u(10)]; uo2 = uo + 1e-4 * [1; zeros(8,1); -1]; %uo2 = [ ;0;0.0119;0.0048;0.0538; ;9.0077;0.0088;0; ]; s2 = ode45(f,[0 500],uo2,odeset('abstol',1e-9,'nonnegative',1:10)); ts2 = 1:0.01:500; us2 = deval(s2,ts2)'; a=1:numel(ts1); b=1:numel(ts2); plot3 (us1(a,1)+us1(a,6)+us1(a,7), us1(a,2)+us1(a,3)+us1(a,8)+us1(a,9), us1(a,4)+us1(a,5)+us1(a,10)); hold on; plot3 (us2(b,1)+us2(b,6)+us2(b,7), us2(b,2)+us2(b,3)+us2(b,8)+us2(b,9), us2(b,4)+us2(b,5)+us2(b,10), 'r'); axis vis3d a =1:1:numel(ts1); b =1:1:numel(ts2); plot3 (us1(a,1)+us1(a,6)+us1(a,7), us1(a,2)+us1(a,3)+us1(a,8)+us1(a,9), us1(a,4)+us1(a,5)+us1(a,10)); plot3 (us2(b,1)+us2(b,6)+us2(b,7), us2(b,2)+us2(b,3)+us2(b,8)+us2(b,9), us2(b,4)+us2(b,5)+us2(b,10), 'r'); xlabel('susceptible'); ylabel('infectious'); zlabel('recovered'); hold off; hold on; plot (us1(a,1)+us1(a,6)+us1(a,7), us1(a,2)+us1(a,3)+us1(a,8)+us1(a,9)); xlabel('susceptible'); ylabel('infectious'); hold off; plot(ts1,us1);title('1e-9 acc +ve'); axis([ e3 1e3]) hold on; 9

10 plot(ts2,us2); plot(ts1, log(sqrt(sum((us1-us2).^2,2)))); Second code is an adaptation of the first one, for plotting a graph of Lyapunov exponent against the parameter : %Searching for Chaos in the MSEIR model for infectious diseases %Jonathan Wills, Fall 2009 %Dengue fever is a multi-strain system that can be modeled by a permutation %of the SIR model: here are the differential eqns. %Constants: N = 100; %puts values in percentages of total population M = (1/65); %years^(-1) %average life span/death rate Y = 52; %years^(-1) rate of recovery from infection, In- >Rn, where n=1 or n=2 B = 2*Y; %infection rate from either of groups I1 or I2 A = 2; %years^(-1) %rate of loss of cross immunity from R1->S1 Q = 0.9; steps = 15; %increase to get more lyapunov exponents to average for a single Q lyap = zeros(10,50); %different starting values for each calculation of the lyapunov exponent for a single Q lyap(:,1)=uo; h = zeros(1,steps); %list of lyapunov exponenets for a single Q to average h2 = zeros(1,10); %list of average lyapunov exponents to plot hold on; for Q2 = 1:30; %vary this to get more accurate graphs of Q vs. lyapunov exponenet Q = Q2/10; eps = 0.1; F [-B/N*u(1)*(u(2)+Q*u(9))-B/N*u(1)*(u(3)+Q*u(8))+M*(N-u(1)); B/N*u(1)*(u(2)+Q*u(9))-(Y+M)*u(2); B/N*u(1)*(u(3)+Q*u(8))-(Y+M)*u(3); Y*u(2)-(A+M)*u(4); Y*u(3)-(A+M)*u(5); -B/N*u(6)*(u(3)+Q*u(8))+A*u(4)- M*u(6); -B/N*u(7)*(u(2)+Q*u(9))+A*u(5)-M*u(7); B/N*u(6)*(u(3)+Q*u(8))- (Y+M)*u(8); B/N*u(7)*(u(2)+Q*u(9))-(Y+M)*u(9); Y*(u(8)+u(9))-M*u(10)]; % below are a variety of initial conditions i used throughout the % investigations %uo = [27;16;14;12;10;8;6;4;2;1]; %uo = [ ; ; ; ; ; ; ; ; ; ]; uo = [ ; ; ; ; ; ; ; ; ; ]; 10

11 for z=1:steps; uo = lyap(:,z); uo2 = uo + 1e-4 * [1; zeros(8,1); -1]; %uo2 = [ ;0;0.0119;0.0048;0.0538; ;9.0077;0.0088;0; ]; %[ts,us] = ode45(f,[0 500],uo); T = 500; s1 = ode45(f,[0 T],uo,odeset('abstol',1e-9,'nonnegative',1:10)); ts = 1:0.5:T; us1 = deval(s1,ts)'; s2 = ode45(f,[0 T],uo2,odeset('abstol',1e-9,'nonnegative',1:10)); us2 = deval(s2,ts)'; disp(q2); %indicator of progress if 0 %plot for each calculation of lyap exp plot(ts, log(abs(us1(:,1)-us2(:,1)))); xlabel('time'); ylabel('log(distance between two trajectories)'); dist = sqrt(sum((us1-us2).^2,2)); ff = find(dist>1); fg = numel(ff); if fg > 0 gg = ff(1); if fg <1 gg = 250; %an approximation for slope when lyapunov exponent is less than zero (any stopping point should work because there will be no plateau) lyap(:,z+1) = us1(gg,:); h(z)=(1/(gg-10)*((log(dist(gg)))-log(dist(10)))); avgh = sum(h)/numel(h); %averages h for a single Q h2(q2) = avgh; %records the average h for this value of Q, makes a list plot(0.1:0.1:3, h2); %plots Q vs. h 11

Introduction to SEIR Models

Introduction 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 information

A Note on the Spread of Infectious Diseases. in a Large Susceptible Population

A 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 information

Mathematical Analysis of Epidemiological Models: Introduction

Mathematical 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 information

Thursday. Threshold and Sensitivity Analysis

Thursday. 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 information

Stability of SEIR Model of Infectious Diseases with Human Immunity

Stability 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 information

Qualitative Analysis of a Discrete SIR Epidemic Model

Qualitative 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

A Time Since Recovery Model with Varying Rates of Loss of Immunity

A 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 information

MODELING THE SPREAD OF DENGUE FEVER BY USING SIR MODEL. Hor Ming An, PM. Dr. Yudariah Mohammad Yusof

MODELING THE SPREAD OF DENGUE FEVER BY USING SIR MODEL. Hor Ming An, PM. Dr. Yudariah Mohammad Yusof MODELING THE SPREAD OF DENGUE FEVER BY USING SIR MODEL Hor Ming An, PM. Dr. Yudariah Mohammad Yusof Abstract The establishment and spread of dengue fever is a complex phenomenon with many factors that

More information

The Existence and Stability Analysis of the Equilibria in Dengue Disease Infection Model

The Existence and Stability Analysis of the Equilibria in Dengue Disease Infection Model Journal of Physics: Conference Series PAPER OPEN ACCESS The Existence and Stability Analysis of the Equilibria in Dengue Disease Infection Model Related content - Anomalous ion conduction from toroidal

More information

Spatial Heterogeneity in Epidemic Models

Spatial Heterogeneity in Epidemic Models J. theor. Biol. (1996) 179, 1 11 Spatial Heterogeneity in Epidemic Models ALUN L. LLOYD AND ROBERT M. MAY University of Oxford, Department of Zoology, South Parks Road, Oxford OX1 3PS, U.K. (Received on

More information

Name Student ID. Good luck and impress us with your toolkit of ecological knowledge and concepts!

Name Student ID. Good luck and impress us with your toolkit of ecological knowledge and concepts! Page 1 BIOLOGY 150 Final Exam Winter Quarter 2000 Before starting be sure to put your name and student number on the top of each page. MINUS 3 POINTS IF YOU DO NOT WRITE YOUR NAME ON EACH PAGE! You have

More information

Analysis of Numerical and Exact solutions of certain SIR and SIS Epidemic models

Analysis 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 information

Project 1 Modeling of Epidemics

Project 1 Modeling of Epidemics 532 Chapter 7 Nonlinear Differential Equations and tability ection 7.5 Nonlinear systems, unlike linear systems, sometimes have periodic solutions, or limit cycles, that attract other nearby solutions.

More information

STABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL

STABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL VFAST Transactions on Mathematics http://vfast.org/index.php/vtm@ 2013 ISSN: 2309-0022 Volume 1, Number 1, May-June, 2013 pp. 16 20 STABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL Roman Ullah 1, Gul

More information

Mathematical Modeling and Analysis of Infectious Disease Dynamics

Mathematical 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 information

The Spreading of Epidemics in Complex Networks

The Spreading of Epidemics in Complex Networks The Spreading of Epidemics in Complex Networks Xiangyu Song PHY 563 Term Paper, Department of Physics, UIUC May 8, 2017 Abstract The spreading of epidemics in complex networks has been extensively studied

More information

2 One-dimensional models in discrete time

2 One-dimensional models in discrete time 2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations

More information

Fixed Point Analysis of Kermack Mckendrick SIR Model

Fixed 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 information

Global Stability of a Computer Virus Model with Cure and Vertical Transmission

Global 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 information

Spotlight on Modeling: The Possum Plague

Spotlight 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 information

A survey of mathematical models of Dengue fever

A survey of mathematical models of Dengue fever Georgia Southern University Digital Commons@Georgia Southern Electronic Theses and Dissertations Graduate Studies, Jack N. Averitt College of Spring 2015 A survey of mathematical models of Dengue fever

More information

Accepted Manuscript. Backward Bifurcations in Dengue Transmission Dynamics. S.M. Garba, A.B. Gumel, M.R. Abu Bakar

Accepted Manuscript. Backward Bifurcations in Dengue Transmission Dynamics. S.M. Garba, A.B. Gumel, M.R. Abu Bakar Accepted Manuscript Backward Bifurcations in Dengue Transmission Dynamics S.M. Garba, A.B. Gumel, M.R. Abu Bakar PII: S0025-5564(08)00073-4 DOI: 10.1016/j.mbs.2008.05.002 Reference: MBS 6860 To appear

More information

The SIR Disease Model Trajectories and MatLab

The SIR Disease Model Trajectories and MatLab The SIR Disease Model Trajectories and MatLab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 2013 Outline Reviewing the SIR

More information

Applications in Biology

Applications in Biology 11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety

More information

Simple Mathematical Model for Malaria Transmission

Simple Mathematical Model for Malaria Transmission Journal of Advances in Mathematics and Computer Science 25(6): 1-24, 217; Article no.jamcs.37843 ISSN: 2456-9968 (Past name: British Journal of Mathematics & Computer Science, Past ISSN: 2231-851) Simple

More information

A New Mathematical Approach for. Rabies Endemy

A 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 information

The effect of emigration and immigration on the dynamics of a discrete-generation population

The effect of emigration and immigration on the dynamics of a discrete-generation population J. Biosci., Vol. 20. Number 3, June 1995, pp 397 407. Printed in India. The effect of emigration and immigration on the dynamics of a discrete-generation population G D RUXTON Biomathematics and Statistics

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

Analysis of SIR Mathematical Model for Malaria disease with the inclusion of Infected Immigrants

Analysis of SIR Mathematical Model for Malaria disease with the inclusion of Infected Immigrants IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 14, Issue 5 Ver. I (Sep - Oct 218), PP 1-21 www.iosrjournals.org Analysis of SIR Mathematical Model for Malaria disease

More information

Dynamics of Disease Spread. in a Predator-Prey System

Dynamics 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 information

Figure The Threshold Theorem of epidemiology

Figure 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 information

Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations

Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger Project I: Predator-Prey Equations The Lotka-Volterra Predator-Prey Model is given by: du dv = αu βuv = ρβuv

More information

MATHEMATICAL MODELS Vol. III - Mathematical Models in Epidemiology - M. G. Roberts, J. A. P. Heesterbeek

MATHEMATICAL 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 information

Global Analysis of an SEIRS Model with Saturating Contact Rate 1

Global 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

Models 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 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 information

Modelling of the Hand-Foot-Mouth-Disease with the Carrier Population

Modelling 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 information

The Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time

The 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 information

On state-space reduction in multi-strain pathogen models, with an application to antigenic drift in influenza A

On state-space reduction in multi-strain pathogen models, with an application to antigenic drift in influenza A On state-space reduction in multi-strain pathogen models, with an application to antigenic drift in influenza A Sergey Kryazhimskiy, Ulf Dieckmann, Simon A. Levin, Jonathan Dushoff Supporting information

More information

Research Article A Delayed Epidemic Model with Pulse Vaccination

Research 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 information

THE STABILITY AND HOPF BIFURCATION OF THE DENGUE FEVER MODEL WITH TIME DELAY 1

THE STABILITY AND HOPF BIFURCATION OF THE DENGUE FEVER MODEL WITH TIME DELAY 1 italian journal of pure and applied mathematics n. 37 2017 (139 156) 139 THE STABILITY AND HOPF BIFURCATION OF THE DENGUE FEVER MODEL WITH TIME DELAY 1 Jinlan Guan 2 Basic Courses Department Guangdong

More information

GLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS

GLOBAL 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 information

Behavior Stability in two SIR-Style. Models for HIV

Behavior 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 information

Dynamical models of HIV-AIDS e ect on population growth

Dynamical 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 information

Delay SIR Model with Nonlinear Incident Rate and Varying Total Population

Delay 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 information

MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS

MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS Journal of Biological Systems, Vol. 11, No. 3 2003 325 335 c World Scientific Publishing Company MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS

More information

Network Models in Epidemiology

Network Models in Epidemiology Network Models in Epidemiology Considering Discrete and Continuous Dynamics arxiv:1511.01062v1 [q-bio.pe] 19 Oct 2015 Edward Rusu University of Washington email epr24@uw.edu 24 April 2014 Abstract Discrete

More information

Global Stability of SEIRS Models in Epidemiology

Global Stability of SEIRS Models in Epidemiology Global Stability of SRS Models in pidemiology M. Y. Li, J. S. Muldowney, and P. van den Driessche Department of Mathematics and Statistics Mississippi State University, Mississippi State, MS 39762 Department

More information

A comparison of existing measles models

A comparison of existing measles models A comparison of existing measles models by Clifford Allotey A thesis submitted to the Faculty of Graduate Studies of The University of Manitoba in partial fulfillment of the requirements of the degree

More information

Modeling the Spread of Epidemic Cholera: an Age-Structured Model

Modeling the Spread of Epidemic Cholera: an Age-Structured Model Modeling the Spread of Epidemic Cholera: an Age-Structured Model Alen Agheksanterian Matthias K. Gobbert November 20, 2007 Abstract Occasional outbreaks of cholera epidemics across the world demonstrate

More information

Epidemics in Complex Networks and Phase Transitions

Epidemics 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 information

Stability Analysis of a SIS Epidemic Model with Standard Incidence

Stability 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 information

Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate

Global 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 information

SI j RS E-Epidemic Model With Multiple Groups of Infection In Computer Network. 1 Introduction. Bimal Kumar Mishra 1, Aditya Kumar Singh 2

SI 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 information

arxiv: v2 [q-bio.pe] 3 Oct 2018

arxiv: v2 [q-bio.pe] 3 Oct 2018 Journal of Mathematical Biology manuscript No. (will be inserted by the editor Global stability properties of renewal epidemic models Michael T. Meehan Daniel G. Cocks Johannes Müller Emma S. McBryde arxiv:177.3489v2

More information

Modeling the Existence of Basic Offspring Number on Basic Reproductive Ratio of Dengue without Vertical Transmission

Modeling the Existence of Basic Offspring Number on Basic Reproductive Ratio of Dengue without Vertical Transmission International Journal on Recent and Innovation Trends in Computing and Communication ISSN: 232-869 Modeling the Existence of Basic Offspring Number on Basic Reproductive Ratio of Dengue without Vertical

More information

A comparison of delayed SIR and SEIR epidemic models

A 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 information

LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC

LAW 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 information

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission P. van den Driessche a,1 and James Watmough b,2, a Department of Mathematics and Statistics, University

More information

Mathematical modelling and controlling the dynamics of infectious diseases

Mathematical 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 information

PARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS

PARAMETER 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 information

Electronic appendices are refereed with the text. However, no attempt has been made to impose a uniform editorial style on the electronic appendices.

Electronic 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 information

The E ect of Occasional Smokers on the Dynamics of a Smoking Model

The E ect of Occasional Smokers on the Dynamics of a Smoking Model International Mathematical Forum, Vol. 9, 2014, no. 25, 1207-1222 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.46120 The E ect of Occasional Smokers on the Dynamics of a Smoking Model

More information

6.2 Brief review of fundamental concepts about chaotic systems

6.2 Brief review of fundamental concepts about chaotic systems 6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification

More information

... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré

... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out

More information

HETEROGENEOUS MIXING IN EPIDEMIC MODELS

HETEROGENEOUS 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 information

Non-Linear Models Cont d: Infectious Diseases. Non-Linear Models Cont d: Infectious Diseases

Non-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 information

Stochastic modelling of epidemic spread

Stochastic modelling of epidemic spread Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca

More information

The dynamics of disease transmission in a Prey Predator System with harvesting of prey

The dynamics of disease transmission in a Prey Predator System with harvesting of prey ISSN: 78 Volume, Issue, April The dynamics of disease transmission in a Prey Predator System with harvesting of prey, Kul Bhushan Agnihotri* Department of Applied Sciences and Humanties Shaheed Bhagat

More information

Introduction to Stochastic SIR Model

Introduction to Stochastic SIR Model Introduction to Stochastic R Model Chiu- Yu Yang (Alex), Yi Yang R model is used to model the infection of diseases. It is short for Susceptible- Infected- Recovered. It is important to address that R

More information

(mathematical epidemiology)

(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 information

AN ABSTRACT OF THE THESIS OF. Margaret-Rose W. Leung for the degree of Honors Baccalaureate of Science in Mathematics

AN ABSTRACT OF THE THESIS OF. Margaret-Rose W. Leung for the degree of Honors Baccalaureate of Science in Mathematics AN ABSTRACT OF THE THESIS OF Margaret-Rose W. Leung for the degree of Honors Baccalaureate of Science in Mathematics presented on June 5, 2012. Title: A Vector Host Model for Coinfection by Barley/Cereal

More information

Mathematical Model of Dengue Disease Transmission Dynamics with Control Measures

Mathematical Model of Dengue Disease Transmission Dynamics with Control Measures Journal of Advances in Mathematics and Computer Science 23(3): 1-12, 2017; Article no.jamcs.33955 Previously known as British Journal of Mathematics & Computer Science ISSN: 2231-0851 Mathematical Model

More information

GLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT

GLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY

More information

Mathematical Model of Dengue Disease Transmission with Severe DHF Compartment

Mathematical Model of Dengue Disease Transmission with Severe DHF Compartment BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 30(2) (2007), 143 157 Mathematical Model of Dengue Disease Transmission with Severe

More information

ME 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 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 information

Chapter 2 Chaos theory and its relationship to complexity

Chapter 2 Chaos theory and its relationship to complexity Chapter 2 Chaos theory and its relationship to complexity David Kernick This chapter introduces chaos theory and the concept of non-linearity. It highlights the importance of reiteration and the system

More information

A NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD

A 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 information

Epidemics in Two Competing Species

Epidemics in Two Competing Species Epidemics in Two Competing Species Litao Han 1 School of Information, Renmin University of China, Beijing, 100872 P. R. China Andrea Pugliese 2 Department of Mathematics, University of Trento, Trento,

More information

Hepatitis C Mathematical Model

Hepatitis C Mathematical Model Hepatitis C Mathematical Model Syed Ali Raza May 18, 2012 1 Introduction Hepatitis C is an infectious disease that really harms the liver. It is caused by the hepatitis C virus. The infection leads to

More information

Stability Analysis of an SVIR Epidemic Model with Non-linear Saturated Incidence Rate

Stability Analysis of an SVIR Epidemic Model with Non-linear Saturated Incidence Rate Applied Mathematical Sciences, Vol. 9, 215, no. 23, 1145-1158 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.41164 Stability Analysis of an SVIR Epidemic Model with Non-linear Saturated

More information

Asynchronous and Synchronous Dispersals in Spatially Discrete Population Models

Asynchronous and Synchronous Dispersals in Spatially Discrete Population Models SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 7, No. 2, pp. 284 310 c 2008 Society for Industrial and Applied Mathematics Asynchronous and Synchronous Dispersals in Spatially Discrete Population Models Abdul-Aziz

More information

Seasonally forced disease dynamics explored as switching between attractors

Seasonally forced disease dynamics explored as switching between attractors Physica D 148 (2001) 317 335 Seasonally forced disease dynamics explored as switching between attractors Matt J. Keeling, Pejman Rohani, Bryan T. Grenfell Zoology Department, University of Cambridge, Downing

More information

Development and Validation of. Statistical and Deterministic Models. Used to Predict Dengue Fever in. Mexico

Development and Validation of. Statistical and Deterministic Models. Used to Predict Dengue Fever in. Mexico Development and Validation of Statistical and Deterministic Models Used to Predict Dengue Fever in Mexico A thesis presented by Aditi Hota to the Applied Mathematics Department in partial fulfillment of

More information

Chaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad

Chaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad Chaos Dr. Dylan McNamara people.uncw.edu/mcnamarad Discovery of chaos Discovered in early 1960 s by Edward N. Lorenz (in a 3-D continuous-time model) Popularized in 1976 by Sir Robert M. May as an example

More information

GLOBAL STABILITY OF A VACCINATION MODEL WITH IMMIGRATION

GLOBAL STABILITY OF A VACCINATION MODEL WITH IMMIGRATION Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 92, pp. 1 10. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY

More information

NOTES ON CHAOS. Chaotic dynamics stem from deterministic mechanisms, but they look very similar to random fluctuations in appearance.

NOTES ON CHAOS. Chaotic dynamics stem from deterministic mechanisms, but they look very similar to random fluctuations in appearance. NOTES ON CHAOS. SOME CONCEPTS Definition: The simplest and most intuitive definition of chaos is the extreme sensitivity of system dynamics to its initial conditions (Hastings et al. 99). Chaotic dynamics

More information

Math 310: Applied Differential Equations Homework 2 Prof. Ricciardi October 8, DUE: October 25, 2010

Math 310: Applied Differential Equations Homework 2 Prof. Ricciardi October 8, DUE: October 25, 2010 Math 310: Applied Differential Equations Homework 2 Prof. Ricciardi October 8, 2010 DUE: October 25, 2010 1. Complete Laboratory 5, numbers 4 and 7 only. 2. Find a synchronous solution of the form A cos(ωt)+b

More information

Viable Control of an Epidemiological Model

Viable Control of an Epidemiological Model Viable Control of an Epidemiological Model arxiv:1510.01055v1 [math.oc] 5 Oct 2015 Michel De Lara Lilian Sofia Sepulveda Salcedo January 5, 2018 Abstract In mathematical epidemiology, epidemic control

More information

Control of Epidemics by Vaccination

Control of Epidemics by Vaccination Control of Epidemics by Vaccination Erik Verriest, Florent Delmotte, and Magnus Egerstedt {erik.verriest,florent,magnus}@ece.gatech.edu School of Electrical and Computer Engineering Georgia Institute of

More information

Available online at Commun. Math. Biol. Neurosci. 2015, 2015:29 ISSN:

Available online at   Commun. Math. Biol. Neurosci. 2015, 2015:29 ISSN: Available online at http://scik.org Commun. Math. Biol. Neurosci. 215, 215:29 ISSN: 252-2541 AGE-STRUCTURED MATHEMATICAL MODEL FOR HIV/AIDS IN A TWO-DIMENSIONAL HETEROGENEOUS POPULATION PRATIBHA RANI 1,

More information

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0

More information

Mathematical Biology - Lecture 1 - general formulation

Mathematical Biology - Lecture 1 - general formulation Mathematical Biology - Lecture 1 - general formulation course description Learning Outcomes This course is aimed to be accessible both to masters students of biology who have a good understanding of the

More information

ECS 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, 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 information

Chaotic motion. Phys 750 Lecture 9

Chaotic motion. Phys 750 Lecture 9 Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

More information

Transmission Dynamics of an Influenza Model with Vaccination and Antiviral Treatment

Transmission Dynamics of an Influenza Model with Vaccination and Antiviral Treatment Bulletin of Mathematical Biology (2010) 72: 1 33 DOI 10.1007/s11538-009-9435-5 ORIGINAL ARTICLE Transmission Dynamics of an Influenza Model with Vaccination and Antiviral Treatment Zhipeng Qiu a,, Zhilan

More information

Section 8.1 Def. and Examp. Systems

Section 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 information

1 Disease Spread Model

1 Disease Spread Model Technical Appendix for The Impact of Mass Gatherings and Holiday Traveling on the Course of an Influenza Pandemic: A Computational Model Pengyi Shi, Pinar Keskinocak, Julie L Swann, Bruce Y Lee December

More information

A FRACTIONAL ORDER SEIR MODEL WITH DENSITY DEPENDENT DEATH RATE

A FRACTIONAL ORDER SEIR MODEL WITH DENSITY DEPENDENT DEATH RATE Hacettepe Journal of Mathematics and Statistics Volume 4(2) (211), 287 295 A FRACTIONAL ORDER SEIR MODEL WITH DENSITY DEPENDENT DEATH RATE Elif Demirci, Arzu Unal and Nuri Özalp Received 21:6 :21 : Accepted

More information

A mathematical model for malaria involving differential susceptibility, exposedness and infectivity of human host

A mathematical model for malaria involving differential susceptibility, exposedness and infectivity of human host A mathematical model for malaria involving differential susceptibility exposedness and infectivity of human host A. DUCROT 1 B. SOME 2 S. B. SIRIMA 3 and P. ZONGO 12 May 23 2008 1 INRIA-Anubis Sud-Ouest

More information

Analysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION

Analysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION Analysis of bacterial population growth using extended logistic Growth model with distributed delay Tahani Ali Omer Department of Mathematics and Statistics University of Missouri-ansas City ansas City,

More information