LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
|
|
- Dennis Charles
- 6 years ago
- Views:
Transcription
1 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 processes converge to the trajectories of a two dimensional dynamical system. If the reproductive number, R, is greater than one then the limiting system has an interior fixed point so that the epidemic processes started with a nontrivial proportion of initially infected are endemic. We show global stability of the interior fixed point of the limiting dynamical system. 1. Introduction The SIRS epidemic models were introduced in 1933 by Kermack and McKendrick 5) to describe endemic infections. In any time interval t, t + t) a susceptible individual is equally likely to come into contact with any of the currently infected individuals and the average number of such contacts is proportional to a fixed fraction of currently infected individuals, θ 1 I t /N. Susceptible individuals are equally susceptible to the disease and the probability of infection upon a contact with an infected individual is θ 2. In the same time interval, any infected individual may recover and become temporarily immune at rate ρ and an immune individual may lose his immunity and become susceptible at rate 1. That is, the unit of time is taken to be equal to the average immune period. The population is assumed closed so that there is no inflow or outflow of individuals. Let N IN be the population size parameter. For t, let S t = the number of susceptible individuals at time t, I t = the number of infected individuals at time t, R t = the number of recovered immune) individuals at time t, and let s t = S t /N, i t = I t /N and r t = R t /N. Since S t + I t + R t N, the pair S t, I t ) or s t, i t ) completely describes the state of the system at any time t. We say that a transition occurs whenever an individual changes state. Under the assumptions, the process S t, I t ) is a pure jump Markov process Key words and phrases. Stochastic epidemic model, SIRS epidemic processes, Law of large numbers, Endemic infection, Threshold phenomena in epidemic processes. 1
2 2 R. G. DOLGOARSHINNYKH S θ 1 I ρ R Figure 1. SIRS model. The letters S, I, and R refer to the possible states an individual can assume in turn susceptible, infected and recovered immune). The individual transition rates are shown above the arrows. with the following instantaneous transition rates q{s, I) S 1, I + 1)} = θsi/n, q{s, I) S, I 1)} = ρi, q{s, I) S + 1, I)} = R, where θ = θ 1 θ 2. It will often be more convenient to work with the scaled process γt N = γ t = s t, i t ). This is a continuous time Markov Chain with state space K N = {s, i) : Ns, Ni) Z 2 +, s + i 1} and transition rates 1) q{s, i) s N 1, i + N 1 )} = Nθsi, q{s, i) s, i N 1 )} = Nρi, q{s, i) s + N 1, i)} = Nr. The transition rules 1) imply that for small h > Es t+h F t ) = s t + r t h θi t s t h + oh), Ei t+h F t ) = i t + θi t s t h ρi t h + oh), where F t contains all events that happen before or at time t. Letting h suggests that the random SIRS paths at least over short periods of time) may be approximated in mean by the solutions of the system of ordinary differential equations 2) ds t dt = r t θi t s t, di t dt = θi ts t ρi t.
3 i LLN FOR THE SIRS EPIDEMIC 3 1 Mean field: theta=3,rho= s Figure 2. Mean field. Trajectories of 2) for θ > ρ with started at initial points with i >. We will show in sec. 3 that the random SIRS paths, γ N, tend to solutions of system 2) as the populations size N tends to infinity. This is a first order approximation to the random paths, or the strong law of large numbers for the problem. Here and later we will often refer to the random SIRS paths by γ N or γ and to the solutions of the mean path ODE by γ. Trajectories of the dynamical system 2) lie in K = {s, i) : s, i, s + i 1}. The system has fixed points, that is points such that the rates of change in s and i at these points are all equal to. Solving for such points, we get two possibilities; either s = ρ θ s and i = θ ρ θ1 + ρ) = 1 ρ/θ 1 + ρ i or s = 1 and i =. When ρ < θ the fixed point γ = s, i ) is in the interior of the triangle K; when ρ θ the only fixed point of the dynamical system in K is s = 1, i =. We will show in sec. 5 that the solutions of 2) started at initial points with i > converge to γ when ρ < θ and to 1, ) if ρ θ as t tends to infinity. Once we establish that the random paths γ N converge to solutions γ of 2) we can conclude that if the rate of infection θ is less than or equal to the rate of recovery ρ, the infection tends to become rapidly extinct by following the mean path to the point where s = 1 and i =. On the other hand, if ρ is less than θ then with a significant fraction of the population initially infected we expect the epidemic to follow the mean path to the fixed point in the interior of K and become endemic; that is, the fraction of infected individuals in the population remains significant for a long time. In the case when ρ < θ, the state of the system turns out to be well approximated by the endemic level, γ, for a long time.
4 i i i 4 R. G. DOLGOARSHINNYKH N=1 N= population fractions s time N=4 N= population fractions s time N=25 N= population fractions s time Figure 3. Initial infection spread. The simulations are for population sizes N = 1, 4, 25 and for parameters θ = 2, ρ = 1. For this choice of parameters the fixed point is γ = s, i ) = 1/2, 1/4). All three simulations were started at s, i ) = 1/4, 1/2) N. The left plots show the dynamics of the infection spread in terms of fractions, s and i, of susceptibles and infecteds for the three population sizes over time period of length 5. As N increases the random paths are seen to follow the limiting deterministic paths more closely. In the plots on the right side values of s and i are plotted over time. We see that the fluctuations around the fixed point γ decrease with N.
5 LLN FOR THE SIRS EPIDEMIC 5 The existence of a threshold at the level ρ = θ, that is, qualitatively different behavior of the model for ρ < θ and ρ θ, is not unexpected. The reproduction number, R, is defined as the average) number of secondary cases that are produced by one infected individual in wholly susceptible population. It has been established that large epidemics can occur when R > 1, and the infection quickly becomes extinct when R < 1 see e.g. 6)). For the SIRS model with transitions given by 1) the reproduction number R = θ 1/ρ = θ/ρ. Therefore, the epidemics die out quickly when R < 1 or θ < ρ and but not when R > 1, that is when θ > ρ. The case R = 1 for the SIRS model also leads to extinction of the infection after a short period of time. 2. SIRS Processes as Time Changed Poisson Processes Before we proceed with the proof of our main theorem we construct a version of the random SIRS processes γ N by applying random time changes to a collection of Poisson processes see e.g. 3) chap 4, 4. In particular, let Y 1, Y 2, Y 3 be independent rate one Poisson processes and let y k ) = N 1 Y k N ) for k = 1, 2, 3 and N IN. If the initial population fractions γ = s, i ) are nonrandom then there is a version of γ satisfying ) ) s t = s y 1 θs u i u du + y 3 r u du 3) ) i t = i + y 1 θs u i u du y 2 In what follows it will be useful to define ) ρi u du. 4) β 1 γ) = θsi, β 2 γ) = ρi, β 3 γ) = r = 1 s i. Note that for any t > and k = 1, 2, 3 β k γ u ) du θ + ρ + 1)t < and therefore the Y k Nt) s are well defined for all t [, ). For random initial conditions we can construct the processes γ N satisfying 3) by first generating initial conditions independently of Y 1, Y 2, Y Proof of the Law of Large Numbers Fix T >. Let γ = γ N be a sequence of SIRS processes, indexed by the population size parameter N, satisfying 3) for t T. Theorem 1. If lim N γ = γ and γ is a solution of 2) on [, T ] with an initial condition γ then lim sup γ γ = a.s. N t T
6 6 R. G. DOLGOARSHINNYKH Proof. Since γ satisfies 3) for t T, s t s t s s + ỹ1 β 1 t γ u ) du) + ỹ3 β 3 γ u ) du) + β 1 γ u ) du β 1 t γ u ) du + β 3 γ u ) du β 3 γ u ) du where ỹ k s) = y k s) s are centered and scaled Poisson processes and β k s. are defined by 4). Let ɛ k t) = ỹ k βk γ u ) du) From the strong law of large numbers for a Poisson process, we have for v < and k = 1, 2, 3 lim sup ỹ k u) = a.s. N u v Since β k γ u ) θ + ρ + 1) it follows that lim sup ɛ k t) = N t T for k = 1, 2, 3. Because the functions β k are Lipschitz continuous as functions of γ, there exists a constant M > such that s t s t s s + ɛ 1 t) + ɛ 3 t) + a.s. M γ u γ u du. The same argument leads to a similar bound on i t ī t and combining the two we get ) γ t γ t 2 γ γ + ɛ 1 t) + ɛ 2 t) + ɛ 3 t) + M γ u γ u du ) 2 γ γ + sup k sup ɛ k t) + t T M γ u γ u du Hence by Gronwall s inequality stated below for proof see e.g. 4), chap 3) γ t γ t 2 γ γ + sup k and the conclusion of the theorem follows. Gronwall s Inequality. sup ɛ k t) ) e 2MT t T Let f be an integrable function on [, T ]. If M > and then ft) ɛ + M fs) ds, for t T, ft) ɛe Mt, for t T..
7 LLN FOR THE SIRS EPIDEMIC 7 4. Construction of Lyapunov Function Let fx) be a continuous function of x defined on an open set containing. Let x : [, ) R d solve ẋ t = fx t ). A function V x) defined in a neighborhood of x = is a Lyapunov function if it has continuous partial derivatives, V x) and the trajectory derivative of V satisfies V x t ) = grad V ) fx t ), where the dot denotes scalar multiplication and grad V = V/ x 1,..., V/ x d ). Constructing a Lyapunov function allows us to show global asymptotic stability of trajectories of 2) started at initial points such that i > ; that is, if γ is a solution of 2) such that i > then lim t γ t = γ. To construct a Lyapunov function for the mean path vector field we need to find a function that is positive on the domain and has a negative derivative along the directions of the vector field. We will find it convenient to switch to i, r) coordinates. In these coordinates our vector field becomes di dt = θ1 i t r t )i t ρi t dr dt = ρi t r t. Let x t = i t i, y t = r t r. Showing that, ) is globally asymptotically stable in the x, y) coordinates is equivalent to showing that γ is globally asymptotically stable in the old coordinates. In x, y) coordinates, dx dt = θ 1 x t + i ) y t + r ) ) x t + i ) ρx t + i ) dy dt = ρx t + i ) y t + r ) and simplifying we get dx dt = θx2 t θx t y t θ ρ 1 + ρ x t θ ρ 1 + ρ y t 5) dy dt = ρx t y t. Following the prescription that can be found in e.g. 2) we seek a Lyapunov function of the form vx, y) = F x) + Gy). Then dvx t, y t ) dt = df dx θx 2 t θx t y t θ ρ 1 + ρ x t θ ρ ) 1 + ρ y t + dg dy ρx t y t ).
8 8 R. G. DOLGOARSHINNYKH We further restrict our search to functions v such that the cross terms disappear, that is omitting index t to simplify expressions) df dx θxy θ ρ ) 1 + ρ y + dg ρx ; dy implying that df dx y df dx This suggests the possibility θx + θ ρ 1 + ρ θx + θ ρ 1 + ρ df dx = which when integrated gives and 6) F x) = x i log ) df dx x + i x x x + i, = dg dy ρx ) 1 ρx = dg dy 1 y = dg dy s y. dg dy = y s, 1 + x ), Gy) = y2 i 2s vx, y) = x i log 1 + x ) + y2. i 2s It is not hard to see that vx, y) on the domain of interest, namely on {x, y) x > i, y r, x + y s } and dv dt = df dx θx 2 θ ρ ) 1 + ρ x + dg dy y) = x θxx + i ) y 2 /s = θ x 2 + y 2 /ρ). x + i 5. Stability of Interior Fixed Points In this section we examine stability of the interior fixed point, γ, for the limiting deterministic dynamical system when ρ < θ. We show that for any initial condition with i > and hence x = i i > i lim γ t = γ. t Recall that x t = i t i and y t = r t r and it is enough to show that the system 5) is globally asymptotically stable. In this section, we continue to work in x, y) coordinates.
9 LLN FOR THE SIRS EPIDEMIC 9 Lemma 1. Let x t, y t ) be the solution of 5) started at x, y ), where x > i. There exists x > i such that x t > x for all t. Proof. Since vx, ) + as x i there exists an x < x such that vx, ) > vx, y ) and x > i. Suppose that the solution of 2) started at x, y ) intersects level x = x at some point x, y ). Let v be the Lyapunov function defined by 6). Since v is nonincreasing along the trajectories of 5), but by construction vx, y ) vx, y ), vx, y ) vx, ) > vx, y ) and we arrive at a contradiction. Therefore x t > x for all t. Lemma 2. Any solution of 5) started at a point x, y ) such that x > i is asymptotically stable. Proof. The linearization of 5) near the point,) is dx dt = θ ρ 1 + ρ x θ ρ 1 + ρ y := A x, y) dy dt = ρx y and A has eigenvalues λ 1,2 = θ 1 + ρ ± 1 ) + θ ρ T ρθ ρ), 1 + ρ with negative real parts, since θ > ρ >. Therefore, system 5) is Lyapunov stable, that is for any ε > there exists a δ > such that a solution started in a δ-neighborhood, U δ, of the stable point never leaves the ε-neighborhood of the stable point, U ε, see e.g. 1) chap. 3. Let ε > and δ be as above and x be as in Lemma 1. Suppose the trajectory started at x, y ) never enters U δ. Then the trajectory will lie in K U δ, where K = {x, y) x > x, y r, x + y s }. But dv dt = θ x2 + y 2 /ρ) < C in K U δ for some constant C >. So that vx t, y t ) = vx, y ) + v t dt vx, y ) Ct. If t, v will eventually have to become negative and we arrive at a contradiction. Therefore the trajectory falls into U δ eventually and hence x t, y t ) U ɛ for all t after that happens. Since ɛ was arbitrary the trajectory tends to the stable point as t goes to infinity.
10 1 R. G. DOLGOARSHINNYKH References [1] V. Arnold 1992) Ordinary Differential Equations. Springer, New York [2] Barabashin 197) Introduction to the Theory of Stability. Wolters- Noordhoff, Netherlands [3] S. Ethier and T. Kurtz 1985) Markov Processes. Characterization and Convergence. Wiley & Sons, NY [4] P. Hartman 1964) Ordinary Differential Equations. Wiley & Sons, NY [5] W. Kermack and A. McKendrick1933) Contributions to the mathematical theory of epidemics, III - Further studies of the problem of endemicity. Proceedings of the Royal Society of Edinburgh. Section A. Mathematics [6] D. Mollison, ed. 1995) Epidemic Models: Their Structure and Relation to Data. Cambridge, Cambridge
SAMPLE PATH LARGE DEVIATIONS FOR SIRS EPIDEMIC PROCESSES. Key words: SIRS epidemic models, sample path large deviations.
SAMPLE PATH LARGE DEVIATIONS FOR SIRS EPIDEMIC PROCESSES R. G. DOLGOARSHINNYKH Abstract. We study infection spread in continuous time SIRS epidemic models. When infection is supercritical, the proportions
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationIntroduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium
Introduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium point and the model reproduction number Prove the stability
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 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 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 informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
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 information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationMathematics for Economists
Mathematics for Economists Victor Filipe Sao Paulo School of Economics FGV Metric Spaces: Basic Definitions Victor Filipe (EESP/FGV) Mathematics for Economists Jan.-Feb. 2017 1 / 34 Definitions and Examples
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
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 informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 0
MATH 56A: STOCHASTIC PROCESSES CHAPTER 0 0. Chapter 0 I reviewed basic properties of linear differential equations in one variable. I still need to do the theory for several variables. 0.1. linear differential
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
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 informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationTime-varying Consumption Tax, Productive Government Spending, and Aggregate Instability.
Time-varying Consumption Tax, Productive Government Spending, and Aggregate Instability. Literature Schmitt-Grohe and Uribe (JPE 1997): Ramsey model with endogenous labor income tax + balanced budget (fiscal)
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 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 informationStochastic 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 informationEconomics 204 Fall 2011 Problem Set 2 Suggested Solutions
Economics 24 Fall 211 Problem Set 2 Suggested Solutions 1. Determine whether the following sets are open, closed, both or neither under the topology induced by the usual metric. (Hint: think about limit
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
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 informationMarkov-modulated interactions in SIR epidemics
Markov-modulated interactions in SIR epidemics E. Almaraz 1, A. Gómez-Corral 2 (1)Departamento de Estadística e Investigación Operativa, Facultad de Ciencias Matemáticas (UCM), (2)Instituto de Ciencias
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Department of Mathematics University of Manitoba Winnipeg Julien Arino@umanitoba.ca 19 May 2012 1 Introduction 2 Stochastic processes 3 The SIS model
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
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 informationAn asymptotic ratio characterization of input-to-state stability
1 An asymptotic ratio characterization of input-to-state stability Daniel Liberzon and Hyungbo Shim Abstract For continuous-time nonlinear systems with inputs, we introduce the notion of an asymptotic
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
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 informationTOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS
TOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS ALEX HUMMELS Abstract. This paper proves a theorem that gives conditions for the topological equivalence of linear ordinary differential
More informationNonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points
Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and
More informationMATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:
MATH 51H Section 4 October 16, 2015 1 Continuity Recall what it means for a function between metric spaces to be continuous: Definition. Let (X, d X ), (Y, d Y ) be metric spaces. A function f : X Y is
More informationEcon Lecture 3. Outline. 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n
Econ 204 2011 Lecture 3 Outline 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n 1 Metric Spaces and Metrics Generalize distance and length notions
More information4 Inverse function theorem
Tel Aviv University, 2014/15 Analysis-III,IV 71 4 Inverse function theorem 4a What is the problem................ 71 4b Simple observations before the theorem..... 72 4c The theorem.....................
More informationMath 361: Homework 1 Solutions
January 3, 4 Math 36: Homework Solutions. We say that two norms and on a vector space V are equivalent or comparable if the topology they define on V are the same, i.e., for any sequence of vectors {x
More informationConvergence of Feller Processes
Chapter 15 Convergence of Feller Processes This chapter looks at the convergence of sequences of Feller processes to a iting process. Section 15.1 lays some ground work concerning weak convergence of processes
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationMA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total
MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationSolutions to Assignment-11
Solutions to Assignment-. (a) For any sequence { } of positive numbers, show that lim inf + lim inf an lim sup an lim sup +. Show from this, that the root test is stronger than the ratio test. That is,
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 informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationAustralian Journal of Basic and Applied Sciences
AENSI Journals Australian Journal of Basic and Applied Sciences ISSN:1991-8178 Journal home page: www.ajbasweb.com A SIR Transmission Model of Political Figure Fever 1 Benny Yong and 2 Nor Azah Samat 1
More informationSolving the Poisson Disorder Problem
Advances in Finance and Stochastics: Essays in Honour of Dieter Sondermann, Springer-Verlag, 22, (295-32) Research Report No. 49, 2, Dept. Theoret. Statist. Aarhus Solving the Poisson Disorder Problem
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 informationGLOBAL 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 informationStatistics 253/317 Introduction to Probability Models. Winter Midterm Exam Monday, Feb 10, 2014
Statistics 253/317 Introduction to Probability Models Winter 2014 - Midterm Exam Monday, Feb 10, 2014 Student Name (print): (a) Do not sit directly next to another student. (b) This is a closed-book, closed-note
More informationPartially Observed Stochastic Epidemics
Institut de Mathématiques Ecole Polytechnique Fédérale de Lausanne victor.panaretos@epfl.ch http://smat.epfl.ch Stochastic Epidemics Stochastic Epidemic = stochastic model for evolution of epidemic { stochastic
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 informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationMA651 Topology. Lecture 10. Metric Spaces.
MA65 Topology. Lecture 0. Metric Spaces. This text is based on the following books: Topology by James Dugundgji Fundamental concepts of topology by Peter O Neil Linear Algebra and Analysis by Marc Zamansky
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week 2 November 6 November Deadline to hand in the homeworks: your exercise class on week 9 November 13 November Exercises (1) Let X be the following space of piecewise
More informationStochastic Differential Equations in Population Dynamics
Stochastic Differential Equations in Population Dynamics Numerical Analysis, Stability and Theoretical Perspectives Bhaskar Ramasubramanian Abstract Population dynamics in the presence of noise in the
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 informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationLimit Cycles and Dynamics of Rumor Models. A thesis. presented to. the faculty of the Department of Mathematics. East Tennessee State University
Limit Cycles and Dynamics of Rumor Models A thesis presented to the faculty of the Department of Mathematics East Tennessee State University In partial fulfillment of the requirements for the degree Master
More informationEXISTENCE AND UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH LINEAR PROGRAMS EMBEDDED
EXISTENCE AND UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH LINEAR PROGRAMS EMBEDDED STUART M. HARWOOD AND PAUL I. BARTON Key words. linear programs, ordinary differential equations, embedded
More information6. Age structure. for a, t IR +, subject to the boundary condition. (6.3) p(0; t) = and to the initial condition
6. Age structure In this section we introduce a dependence of the force of infection upon the chronological age of individuals participating in the epidemic. Age has been recognized as an important factor
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 9 217 Linearization of an autonomous system We consider the system (1) x = f(x) near a fixed point x. As usual f C 1. Without loss of generality we assume x
More informationSTABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS
LABORATOIRE INORMATIQUE, SINAUX ET SYSTÈMES DE SOPHIA ANTIPOLIS UMR 6070 STABILITY O PLANAR NONLINEAR SWITCHED SYSTEMS Ugo Boscain, régoire Charlot Projet TOpModel Rapport de recherche ISRN I3S/RR 2004-07
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
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 informationMarkov Chains and Pandemics
Markov Chains and Pandemics Caleb Dedmore and Brad Smith December 8, 2016 Page 1 of 16 Abstract Markov Chain Theory is a powerful tool used in statistical analysis to make predictions about future events
More informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System
More informationModeling with First-Order Equations
Modeling with First-Order Equations MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Radioactive Decay Radioactive decay takes place continuously. The number
More informationApproximating diffusions by piecewise constant parameters
Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional
More informationLie Symmetries Analysis for SIR Model of Epidemiology
Applied Mathematical Sciences, Vol. 7, 2013, no. 92, 4595-4604 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.36348 Lie Symmetries Analysis for SIR Model of Epidemiology A. Ouhadan 1,
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationExpected Values Estimated via Mean-Field Approximation are 1/N-Accurate
Expected Values Estimated via Mean-Field Approximation are 1/-Accurate icolas Gast To cite this version: icolas Gast. Expected Values Estimated via Mean-Field Approximation are 1/-Accurate. ACM SIGMETRICS
More informationA.Piunovskiy. University of Liverpool Fluid Approximation to Controlled Markov. Chains with Local Transitions. A.Piunovskiy.
University of Liverpool piunov@liv.ac.uk The Markov Decision Process under consideration is defined by the following elements X = {0, 1, 2,...} is the state space; A is the action space (Borel); p(z x,
More informationExistence and Uniqueness
Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect
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 informationIstván Faragó, Róbert Horváth. Numerical Analysis of Evolution Equations, Innsbruck, Austria October, 2014.
s István, Róbert 1 Eötvös Loránd University, MTA-ELTE NumNet Research Group 2 University of Technology Budapest, MTA-ELTE NumNet Research Group Analysis of Evolution Equations, Innsbruck, Austria 14-17
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationSIR Epidemic Model with total Population size
Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 7, Number 1 (2016), pp. 33-39 International Research Publication House http://www.irphouse.com SIR Epidemic Model with total Population
More informationNonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points
Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationMATH 217A HOMEWORK. P (A i A j ). First, the basis case. We make a union disjoint as follows: P (A B) = P (A) + P (A c B)
MATH 217A HOMEWOK EIN PEASE 1. (Chap. 1, Problem 2. (a Let (, Σ, P be a probability space and {A i, 1 i n} Σ, n 2. Prove that P A i n P (A i P (A i A j + P (A i A j A k... + ( 1 n 1 P A i n P (A i P (A
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for COMPUTATIONAL BIOLOGY A COURSE CODES: FFR 110, FIM740GU, PhD Time: Place: Teachers: Allowed material: Not allowed: June 8, 2018, at 08 30 12 30 Johanneberg Kristian
More informationd(x n, x) d(x n, x nk ) + d(x nk, x) where we chose any fixed k > N
Problem 1. Let f : A R R have the property that for every x A, there exists ɛ > 0 such that f(t) > ɛ if t (x ɛ, x + ɛ) A. If the set A is compact, prove there exists c > 0 such that f(x) > c for all x
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informationA Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and Cure of Infected Cells in Eclipse Stage
Applied Mathematical Sciences, Vol. 1, 216, no. 43, 2121-213 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.216.63128 A Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and
More informationAngle contraction between geodesics
Angle contraction between geodesics arxiv:0902.0315v1 [math.ds] 2 Feb 2009 Nikolai A. Krylov and Edwin L. Rogers Abstract We consider here a generalization of a well known discrete dynamical system produced
More informationOscillation Criteria for Certain nth Order Differential Equations with Deviating Arguments
Journal of Mathematical Analysis Applications 6, 601 6 001) doi:10.1006/jmaa.001.7571, available online at http://www.idealibrary.com on Oscillation Criteria for Certain nth Order Differential Equations
More informationProblem set 3: Solutions Math 207A, Fall where a, b are positive constants, and determine their linearized stability.
Problem set 3: s Math 207A, Fall 2014 1. Fin the fixe points of the Hénon map x n+1 = a x 2 n +by n, y n+1 = x n, where a, b are positive constants, an etermine their linearize stability. The fixe points
More informationON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS
PORTUGALIAE MATHEMATICA Vol. 55 Fasc. 4 1998 ON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS C. Sonoc Abstract: A sufficient condition for uniqueness of solutions of ordinary
More informationLocally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem
56 Chapter 7 Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem Recall that C(X) is not a normed linear space when X is not compact. On the other hand we could use semi
More informationDie-out Probability in SIS Epidemic Processes on Networks
Die-out Probability in SIS Epidemic Processes on etworks Qiang Liu and Piet Van Mieghem Abstract An accurate approximate formula of the die-out probability in a SIS epidemic process on a network is proposed.
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
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 informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Martin Rasmussen, Peter Giesl A Note on Almost Periodic Variational Equations Preprint Nr. 13/28 14. März 28 Institut für Mathematik, Universitätsstraße,
More information