A Condition for Stability in an SIR Age Structured Disease Model with Decreasing Survival Rate

Similar documents
An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Smoking Tobacco Experiencing with Induced Death

On Ψ-Conditional Asymptotic Stability of First Order Non-Linear Matrix Lyapunov Systems

Stability Analysis of a delayed HIV/AIDS Epidemic Model with Saturated Incidence

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Impulsive Differential Equations. by using the Euler Method

Midterm exam 2, April 7, 2009 (solutions)

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Logistic equation of Human population growth (generalization to the case of reactive environment).

Microscopic Flow Characteristics Time Headway - Distribution

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

Elementary Differential Equations and Boundary Value Problems

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

Institute of Actuaries of India

Midterm Examination (100 pts)

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

Wave Equation (2 Week)

EXERCISE - 01 CHECK YOUR GRASP

Double Slits in Space and Time

On the Speed of Heat Wave. Mihály Makai

CSE 245: Computer Aided Circuit Simulation and Verification

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

3(8 ) (8 x x ) 3x x (8 )

EE 434 Lecture 22. Bipolar Device Models

Charging of capacitor through inductor and resistor

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

H is equal to the surface current J S

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

SUPERCRITICAL BRANCHING DIFFUSIONS IN RANDOM ENVIRONMENT

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

3.4 Repeated Roots; Reduction of Order

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

Available at Vol. 1, No. 1 (2006) pp

DE Dr. M. Sakalli

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

Transfer function and the Laplace transformation

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

Mundell-Fleming I: Setup

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

Stability and Bifurcation in a Neural Network Model with Two Delays

A SEIR Epidemic Model with Infectious Population Measurement

Control System Engineering (EE301T) Assignment: 2

Lecture 2: Current in RC circuit D.K.Pandey

B) 25y e. 5. Find the second partial f. 6. Find the second partials (including the mixed partials) of

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

arxiv: v1 [math.ap] 6 Jan 2016

Discussion 06 Solutions

Chapter 6 Differential Equations and Mathematical Modeling

Reliability Analysis of a Bridge and Parallel Series Networks with Critical and Non- Critical Human Errors: A Block Diagram Approach.

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

14.02 Principles of Macroeconomics Problem Set 5 Fall 2005

46. Let y = ln r. Then dy = dr, and so. = [ sin (ln r) cos (ln r)

Availability Analysis of Repairable Computer Systems and Stationarity Detection

CONTRIBUTION TO IMPULSIVE EQUATIONS

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

The Overlapping Generations growth model. of Blanchard and Weil

symmetric/hermitian matrices, and similarity transformations

2. The Laplace Transform

ANSWERS TO EVEN NUMBERED EXERCISES IN CHAPTER 11

cycle that does not cross any edges (including its own), then it has at least

4.3 Design of Sections for Flexure (Part II)

WEIBULL FUZZY PROBABILITY DISTRIBUTION FOR RELIABILITY OF CONCRETE STRUCTURES

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

1 Finite Automata and Regular Expressions

The Effect of an Unobservable Factor on Interest Rates in a Pure Exchange Economy

The transition:transversion rate ratio vs. the T-ratio.

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

CHAPTER 9 Compressible Flow

Nonlocal Symmetries and Exact Solutions for PIB Equation

On General Solutions of First-Order Nonlinear Matrix and Scalar Ordinary Differential Equations

ON PRICING CONTINGENT CLAIMS IN A TWO INTEREST RATES JUMP-DIFFUSION MODEL VIA MARKET COMPLETIONS

Asymptotic instability of nonlinear differential equations

Lecture 26: Leapers and Creepers

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

The Optimal Timing of Transition to New Environmental Technology in Economic Growth

The Firm s Value and Timing of Adopting E-Commerce Using Real Options Analysis under Uncertainty *

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

The Science of Monetary Policy

Chapter 17 Handout: Autocorrelation (Serial Correlation)

PhysicsAndMathsTutor.com

AN INTRODUCTION TO FOURIER ANALYSIS PROF. VEDAT TAVSANOĞLU

An Introduction to Malliavin calculus and its applications

Transcription:

A Condiion for abiliy in an I Ag rucurd Disas Modl wih Dcrasing urvival a A.K. upriana, Edy owono Dparmn of Mahmaics, Univrsias Padjadjaran, km Bandung-umng 45363, Indonsia fax: 6--7794696, mail: asupria@yahoo.com.au; aksupriana@bdg.cnrin.n.id Financial and Indusrial Mahmaics Group ITB, Bandung 43, Indonsia Absrac In his papr w prsn an I modl for disas ransmission wih an assumpion ha individuals in h undr-laying dmographic populaion xprinc a monoonically dcrasing survival ra. W show ha h rsuls in an analogous I modl for disas ransmission ar h spcial cas of h I modl in his papr. W found ha hr is a hrshold for h disas ransmission drmining h xisnc and h absnc of h ndmic quilibrium. W invsiga h sabiliy of his quilibrium via a Gronwall-lik inqualiy horm. Unlik in h I modl, h hrshold for h xisnc is no quivaln o h hrshold for h sabiliy of h quilibrium. W provid an addiional condiion which consisnly gnralizs h rsuls in h I modl.. Kywor : Disas Modling, I Modl, Thrshold Numbr, abiliy of an Equilibrium Poin I. Inroducion Ag srucur is among h imporan facors affcing h dynamics of a populaion in rlaion o h sprad of conagious disass. To sudy h ffc of ag srucur in h dynamics of conagious disass, a las hr ar wo approachs, firs by dvloping a populaion modl wih coninuous ag [,] and scond by dvloping a populaion wih ag groups [3]. A modl of I disas ransmission is sudid in [] and a modl I disas ransmission is sudid in [] by assuming coninuous ag. An I modl only fis o disass ha caus an infciv individual rmains infciv for lif. To incras ralism, in his papr w prsn a modl for an I disas ransmission by assuming coninuous ag. r w assum individuals in h undr-laying populaion xprinc a monoonically dcrasing survival ra as hir ag gos by. W also assum ha hr is a dnsiy-dpndn bu ag-indpndn birh ra. W show ha hr is an ndmic hrshold, blow which h disas will sop, and abov which h disas will say ndmic. Th rsuls in h I disas modl in [] gnraliz ino h I modl. II. Th Mahmaical Modl Th modl discussd hr is h gnralizaion of h modl in [] o includ an comparmn as an amp o incras h ralism of h modl. Throughou h papr w us h following noaions: N = Toal numbr of individuals in h populaion = Th numbr of suscpibl individuals in h populaion I = Th numbr of infciv individuals in h populaion = Th numbr of rcovr or immun individual in h populaion B = Th rcruimn ra or h birh ra β = Th ransmission probabiliy of h disas W assum ha h populaion N is dividd ino hr comparmns,, I, and, such ha N = I. To includ ag srucur, suppos ha hr xiss Q ( a ), a funcion of ag dscribing h fracion of human populaion who survivs o h ag of a or mor, such ha, Q = and Q ( a ) is a non-ngaiv and monoonically dcrasing for a. If i is assumd ha lif xpcancy is fini, hn Q a = L< and aq. a < Furhr, l also assum ha N,, I, and dnos, rspcivly, h numbrs of N,, I, and who surviv a im. Thn w hav

N N B Q ( a) =. inc h pr capia ra of infcion in h populaion a im is β I, hn h numbr of suscpibl a im is givn by a β I s B Q =. 3 If h ra of rcovry is γ hn h numbr of infciv a im is givn by βi γ a a I = I B Q βi a = I B. Q Furhrmor, considring ha = N I hn w hav = N BQ( a) βi a BQ( a) βi a I B Q βi a = B. Q I is clar ha N =, =, I =, and =. 6 nc, quaions (3), (4), (5), and (6) consiu an I ag srucurd disas modl. III. Th Exisnc of a Thrshold Numbr In his scion w will show ha hr is a hrshold numbr for h modl discussd abov. L us considr h following i sysm of quaions which has h sam bhavior wih h sysm (3) o (6) whnvr : N B Q a I a β =, 7 = B Q, 8 βi a I = B Q 9 βi a = BQ( a) Equaions (7) show ha h valu of N 4 5 is consan, hnc h quaions for h agsrucurd I modl rduc o hr quaions, (8) o. Th Equilibrium of h sysm is givn by (, I, ) wih I saisfying I a I a γ = B Q a. I is asy o s ha (, I, ) = ( N,,) is h disas-fr quilibrium. To find a non-rivial quilibrium (an ndmic quilibrium), w could obsrv h following. Ia βb Q( a) =. β I Th L of () is a funcion of I, say Ia a γ f ( I) B Q( a) β I This funcion is monoonically dcrass wih f( I ) I Ia B Q a I βi B aq f( I ) I Ia B Q I βi = I 3 4 5 Thrfor, a uniqu non-rivial valu of if and only if occurs B aq( a) >. 6 inc h L of (6) drmins h occurrnc of h non-rivial valu of I, hn i will b rfrd as a hrshold numbr of h modl. nc, an ndmic quilibrium (, I, ) ( N,,) occurs if and only if >. IV. Th abiliy of h Equilibria To invsiga h sabiliy of h quilibria w us h mhod in [] and us h lmma hrin.

LEMMA 4.. (BAUE, ). L f b a boundd non-ngaiv funcion which saisfis an sima of h form f f f( a) ( a), whr f is a non-ngaiv funcion wih f = and a is a non-ngaiv funcion wih a <. Thn f. = POOF. []. I is also showd in [] ha h lmma is sill ru if h inqualiy in h lmma is rplacd by f f sup f s ( a). 7 a s Th following lmma is h xnsion of Braur s lmma. LEMMA 4.. L f j, j =, b boundd nonngaiv funcions saisfying f f sup f s ( a), a s f f sup f s ( a), a s whr f j is non-ngaiv wih f j = and j ( a) is non-ngaiv wih j ( a) <. Thn li m f ( ), j j = =,. POOF. f sup{ f, f } sup sup{ f, f }sup{ ( a), ( a)} a s f sup{ f, f } sup sup{ f, f }sup{ ( a), ( a)} a s and hnc, sup{ f, f } sup{ f, f } sup a s sup{ f, f}sup{ ( a), ( a)} From Lmma 4. w conclud ha sup{ f, f} =, and his is suffic o show ha f j =, j =,. 4.. Th sabiliy of h disas-fr quilibrium. W invsiga h sabiliy of h disas-fr quilibrium for h cas of <. Considr h following inqualiis. βi a a β I nc w hav, I = I And aβ sup I 8 a s βi a B Q ( a)( ) I B Q ( a)( aβ sup I ) 9 a s = βi a B Q β a s B Q ( a)( a sup I ) Morovr, sinc I = and B Q a aβ = < hn using Lmma 4. w conclud ha I ( ) =. Nx, l us s h xprssion a γ B which, Q a aβ if = B, can b wrin in h Q a aβ form. nc, if < hn B Q a aβ < (Appndix ). Furhrmor, sinc = hn using Lmma 4. w conclud ha =. Consqunly, = ( N I) = N This shows ha h disas-fr quilibrium ( I, I ) = ( N, ) is globally sabl. V 4.. Th sabiliy of h ndmic quilibrium. Th ndmic quilibrium (, I, ) appars only if. L us s h > prurbaions of I and, rspcivly, by u and v. Dfin I = I u and subsiu his quaniy ino quaion (4) o obain h following calculaions.

I u = I β[ I u] a B Q u = I I βi βu a a B Q Ia = B Q a I u Ia β a B Q a Ia u = B Q ( a)( ) I Ia B Q( a)( ) u Ia β a BQ( a) Ia u = BQ( a)( ) I u Ia β a B Q Ia B Q ( a)( ) I nc, w hav Ia β a s B Q ( a) a sup u Ia ( ) u B Q I sup us BQ ( a ) β a Ia a s By dfining f = u, I a a γ β a = B Q a, and I a a γ = f B Q ( a)( ) I, hn w hav f f sup f s ( a). a s f = W s ha ha and i can b shown ( a) is non-ngaiv wih ( a) < (s Appndix ). Thn by Lmma 4. w hav f =, mans ha I = I. Nx, dfin = v and subsiu hs quaniis ino quaion (5) o obain h following calculaions: v = β[ I u] a B Q ( a)( ) v = βi βu a a B Q = B Q a Ia Ia β u s a B Q Ia = ( ) Ia BQ a u Ia β a B Q v B Q Ia v = BQ( a)( ) I u s a β a B Q Ia BQ( a)( ) Ia BQ ( a ) βa sup a s u nc, w hav Ia ( ) Ia sup a s β v B Q I us BQ a By dfining f = u, f = v Ia a = BQ βa, and ( I a ) a γ f = BQ a I, hn w hav f f sup a s f s ( a). W can show ha in Appndix 3 ha f = and ( a) is non-ngaiv wih ( a) <. Thn by Lmma 4. w hav f =, mans ha =. Finally, sinc N is a consan, =, and I ( ) I = hn (, I, ) is globally sabl.

V. Concluding marks In his papr w hav discussd an agsrucurd I disas modl wih a dcrasing survival ra. W found a hrshold numbr for h xisnc and uniqunss of an ndmic quilibrium, ha is, B aq( a). As is h cas of h I disas modl discuss in [], an ndmic quilibrium appars if > and disappars if <. In h I disas modl, h hrshold for h xisnc of h quilibrium is also h hrshold for h sabiliy of h quilibrium. owvr, in our cas in which hr is a rcovr comparmn, hr is an addiional condiion for h quilibrium o b sabl. r w found ha hr is a sabl ndmic quilibrium if > and B Q ( a ) a< is a sabl disas-fr quilibrium if, and hr < and B Q (. W noic ha a ) a< his condiion is consisn wih ha in [] if h rcovry ra γ =, sinc in his cas is quivaln o. nc, w conclud ha h I modl in [] is naurally nsd in h I modl discussd in his papr. VI. frncs [] F. Braur. A modl for an I disas in an ag-srucurd populaion. Discr and Coninuous Dynamical ysms ris B. (), 57-64. []. Busnbrg, M. Ianlli, and.. Thim. Global bhavior of an ag-srucurd pidmic modl, IAM J. Mah. Anal. (99), 65-8. [3].W. hco. An ag-srucurd modl for prussis ransmission. Mah. Biosc. 45 (997), 89-36. Appndix If < hn ( a) = BQ( a) βa <. I is sraighforward from h dfiniion ha a = B and Q a aβ γ = BQ( a) aβ. Appndix W claim ha f j = and j ( a) is non-ngaiv wih j ( a) <. I is clar ha, if X, Y {, V} wih X Y hn f = j XIYa X X X BQ( a)( ) I =. To winss ha, l us procd for j a < ( a) as follows. Dfin ax IVa a g( x) = Bβ Q( a) β x which is a dcrasing funcion of x. W s ha IVa aγ a g( IV) Bβ = Q( a). βiv inc I V is an quilibrium valu, hn w hav g( I V ) =. Furhrmor, β ( a). g = B aq ( a) = g IVa Considring ha is dcrasing funcion hn ( a) = g < g( I V ) =. Appndix 3 If < hn Ia a γ a = B Q βa is lss han on. Using h rsul in Appndix w hav Ia B Q βa BQ( a) βa <.