Local Stability Analysis of an SVIR Epidemic Model

In this paper, we present an SVIR epidemic model with deadly deseases. Initially the basic formulation of the model is presented. Two equilibrium point exists for the system; disease free and endemic equilibrium point. The local stability of the disease free and endemic equilibrium exists when the basic reproduction number less or greater than unity, respectively. If the value of R0 less than or equal one then the desease free equilibrium point is locally asymptotically stable, and if its exceeds, the endemic equilibrium point is locally asymptotically stable. The numerical results are presented for illustration.


INTRODUCTION
Health problems are closely associated with a disease.It is motivating humans to study health as well as illness.One of the goals of a nation to achieve prosperity in the future is the creation of a healthy living environment.One of the indicators of a healthy environment involves a high level of quality health services.Quality health care issues need to be addressed systematically, effectively and efficiently.Analysis of the spread of a disease is one step to address the problem of quality health services.
The spread of disease due to virus or bacteria that enter into the body will lead to human health disorders and will affect the socio-economic development of the community.Efforts to prevent a spreading disease can work optimally if several stages such as research, development of various diagnostic tools, drugs, and new vaccines, has achieved.
Mathematical models are one of a useful tool for reviewing or analyse the patterns of a disease spread.According to [2] and [4] the spread of infectious diseases can be described mathematically through a model, such as SIR and SIRS.The most basic procedure in modelling the spread of disease is by using the compartment model.In this case, the population is divided into three different classes namely, Susceptible, Infected, and Recovered which then shortened to SIR.Vaccination can be considered as the addition of a class naturally into the epidemic base model for some types of diseases.
Along with the development of knowledge in the field of health, control of an epidemic disease can be done by vaccination action.Then the SIR model develops with the addition of a vaccinated population class.Some researchers have proposed models of the dynamics of vaccination.Among others is SVI type, which has been developed by [5], [6] and [7].Several studies have stated that the threshold

Local Stability Analysis of an SVIR Epidemic Model
Joko Hariantoianto 21 to determine the occurrence of endemic disease seen from the value of Basic Reproduction Number ( 0 ).The disease distribution model will be locally stable if  0 ≤ 1 and the disease will be present if  0 > 1. Paper [8] discusses the SVIR epidemic model for non-lethal diseases that does not involve the individual's mortality rate due to the disease.From these studies, we can build some further assumptions related to the SVIR epidemic model that developed in public life.
In this paper, we analyse the SVIR epidemic model for a deadly disease.We hope the model can be useful to analyse the disease outbreaks so that we can take optimal precautions.

METHODS
This research used literature study with literature sources from some reputable journal.The simulation of model SVIR used mathematical software, i.e.Maple 2016.

Formulation Model
The following are assumptions used in modelling: 1. Closed population (no migration) 2. Births occur in every sub populations and enter in vulnerable subpopulations.
3. Natural death occurs on each sub populations at the same rate.4. Illness can cause death.5. Individuals who have cured cannot return to susceptive class (permanent cure).6.Short incubation period.
7. Vaccination has given to susceptible subpopulations.8.We do not distinguish the rate of recovery for a child has been infected and for the adult.9. Vaccination will reach the level of immunity over time and finally enters the sub populations that heals.10.Individuals who have vaccinated become infected if they lose immunity.11.Contact between infected and susceptive individuals can lead to disease transmission, and so are contact between infected and those have vaccinated.
We investigate the basic model formulation by dividing the total population into four compartments, such as S(t) for susceptible, V(t) for vaccinated, I(t) for infected and R(t) for recovered individuals.Let  is the rate of contact that is sufficient to transmit the disease.We also assume a constant recovery rate .The rate at which the susceptible population is vaccinated is .We assume that there can be disease related death and define  to be the rate of disease related death, while  is the rate of natural death that is not related to the disease.We assume that all newborns enter the susceptible class at the constant rate of .Let  1 be the average rate for susceptible individuals to obtain immunity and moved into recovered population.We do not distinguish the natural immunity and vaccine-induced immunity here because vaccine-induced immunity can also last for a long term.We assume that before obtaining immunity the vaccinees still have possibility of infection with a disease transmission rate  1 while contacting with infected individuals. 1 may be assume to be less than  because the vaccinating individuals may have some partial immunity during the process or they may recognize the transmission characters of the disease and hence decrease the effective contacts with infected individuals.Figure 1 shows the details of the population transfer diagram of the SVIR model.where (0) > 0, (0) > 0, (0) > 0, (0) = 0,  +  +  +  = 1, ∀  ≥ 0 and all parameters are positive.The parameters with their description are presented in Table 1.The transmission rate for vaccination individuals to be infected before gaining immunity  1

Local Stability Analysis of an SVIR Epidemic Model
The average rate for vaccination individuals to obtain immunity and move into recovered population The death rate due to disease  Recovery rate

. Equilibrium Point
Next, we will study system dynamics in equation (1).Since the last equation does not depend on other equations, we simply study the following system.

Local Stability Analysis of an SVIR Epidemic Model
Joko Hariantoianto 23 Furthermore, we performed a SVIR model stability analysis around the equilibrium point of the system (2).As a first step, we define system equilibrium points (disease free points and endemic points).
In general, according to [1] and [3] there exist two equilibrium points of epidemic model, i.e. the disease free equilibrium point and the endemic equilibrium point.The two equilibrium points are reviewed based on the existence of the disease in a population with a continuous time.The disease-free equilibrium point is the point at which the disease is unlikely to spread in an area because the infected population is equal to zero (I=0) for  → ∞.While the endemic equilibrium point is the point at which the disease must spread (I>0) for  → ∞, in a defined closed area.Theorem 1 There are two equilibrium points in the system (2), namely: a) The disease-free equilibrium point, that is, b) The endemic equilibrium point, that is, where  * is the positive root of  1  2 +  2  +  3 (1 − ) = 0, for C > 1, and

Basic Reproduction Number (𝑹 𝟎 )
The equilibrium point stability analysis on the SVIR model depends on the value of basic reproduction number  0 (the number of susceptible individuals who are then infected if they interact with patients in an entirely vulnerable population).A selection of this number is by observing the state of the endemic equilibrium point.
Consider Eq. ( 6), an endemic equilibrium point only applicable to positive roots ( * > 0) when C > 1.In conclusion, the value of a positive endemic equilibrium point lies in the value of C, so that it can be defined the value basic reproduction number, that is,  0 =  ( + )( +  + ) +  1  ( + )( +  1 )( +  + )

LOCAL STABILITY OF EQUILIBRIUM POINT
We will determine the stability of the equilibrium point of the SVIR model around the disease-free equilibrium and the endemic equilibrium point.Equations ( 3), (4), and (5) are nonlinear equations, so by using Taylor series at each equilibrium point, there will be a system of linear differential equations.Theorem 2 Defined: If  0 ≤ 1, then there is exist one (unique) equilibrium of system (1), that is, the diseasefree equilibrium point  0 .Furthermore, if  0 > 1, then there are two equilibrium points of system (1), that is, the disease-free equilibrium points  0 dan  * .ii.
If  0 ≤ 1, then the equilibrium point  0 is locally asymptotically stable and if  0 > 1, then the equilibrium point  0 is unstable.iii.
If  0 > 1, then the endemic equilibrium  * is locally asymptotically stable.

Proof:
(i) We have discussed the existence of the equilibrium point of disease-free in Theorem 1.In the discussion, we obtain a point of equilibrium when I = 0 and is singular if  0 ≤ 1.The equilibrium point is the disease-free equilibrium point.
(ii) The Jacobian Matrix of system (2) is So the Jacobian matrix around the point  0 = ( 0 ,  0 ,  0 ) = (  + ,  (+)( 1 +) , 0) is The characteristic equation is It is clear if  0 < 1 then all the eigenvalues of ( 0 ) are negative, hence the point  0 is locally asymptotically stable.Meanwhile, if  0 > 1, then there is a positive eigenvalue, hence the point  0 is unstable.
(iii) It is clear that the endemic point  * exist when  0 > 1, so the Jacobian matrix around the point  * = ( * ,  * ,  * ) is We can modify the element  11 ,  22 ,  33 , from matrix above into  Based on the results, we will interpret the results biologically as follows,  If  0 ≤ 1, then for  → ∞ and (, , ) that close enough to  0 = ( 0 ,  0 ,  0 ), the solution of system (1) will move to  0 = ( 0 ,  0 ,  0 ).It means that if  0 ≤ 1, then for the number of susceptible, vaccinated and infected individuals close to  0 = ( 0 ,  0 ,  0 ), the disease will not plague and tends to disappear indefinitely time.This condition is then called asymptotically stable around the equilibrium point  0 = ( 0 ,  0 ,  0 ), and we interpret as a disease-free equilibrium point. If  0 > 1, then for  → ∞ and (, , ) that close enough to  * = ( * ,  * ,  * ), the solution of system (1) will move to  * = ( * ,  * ,  * ).This means if  0 > 1, then for the number of susceptible, vaccinated and infected individuals close to  * = ( * ,  * ,  * ), the disease will plague but does not reach extinction in an infinite time.This condition is then called asymptotically stable around the equilibrium point  * = ( * ,  * ,  * ), and we interpret as an endemic equilibrium point.

NUMERICAL RESULTS
In this section, we investigate the numerical solution of the system (1) by using the Runge-Kutta order four scheme.The state variables are chosen with same initials conditions.The numerical results are shown in Figure 1 and 2. Figure 1, illustrate the fact, when the basic reproduction number less than unity, when the value of basic reproduction number exceeds than unity, Figure 2

CONCLUSION
From the results of endemic equilibrium point analysis, we define the basic reproduction number parameter, that is,  0 =  ( + )( +  + ) +  1  ( + )( +  1 )( +  + ) Furthermore,  0 is a necessary condition for the existence of the two points of equilibrium as well as its local stability.
Results of local stability analysis indicates if  0 ≤ 1, then the disease free equilibrium point  0 is locally asymptotically stable.It means if the terms  0 ≤ 1 fulfilled, then in a long time, there would have no spread of disease in susceptible and vaccinated subpopulations, or in other words, the epidemic will stop.Conversely, if  0 > 1 then the endemic equilibrium point E  is locally asymptotically stable.It means, in a long time, the disease will always exist in the population with the proportion of each subpopulation is equal to  * ,  * , dan  * .

Table 1
. Parameter Description Parameter Description  Natural death rate (and equally, the birth rate)  The transmission rate of disease when susceptible individuals contact with infected individuals  The rate of migration of susceptible individuals to the vaccination process  1 According to the Routh-Hurwitz criterion, all the eigenvalues of the above characteristic equations have a negative real part.So, the equilibrium point  * = ( * ,  * ,  * ) is locally asymptotically stable.