The influence of vaccination on the control of JE with a standard incidence rate of mosquitoes, pigs and humans.

In this article, a nonlinear mathematical model used for the impact of vaccination on the control of infectious disease, Japanese encephalitis with a standard incidence rate of mosquitoes, pigs and humans has been planned and analyzed. During the modeling process, it is expected that the disease spreads only due to get in touch with the susceptible and infected class only. It is also assumed that due to the effect of vaccination, the total human population forms a separate class and avoids contact with the infection. The dynamical behaviors of the system have been explored by using the stability theory of differential equations and numerical simulations. The local and global stability of the system for both equilibrium states under certain conditions has been studied. We have set up a threshold condition in the language of the vaccine-induced reproduction number R ( α 1 ) , which is fewer than unity, the disease dies in the absence of the infected population, otherwise, the infection remains in the population. Furthermore, it is found that vaccine coverage has a substantial effect on the basic reproduction number. Also, by continuous efforts and effectiveness of vaccine coverage, the disease can be eradicated. It is also found a more sensitive parameter for the transmission of Japanese encephalitis virus by using sensitivity analysis. In addition, numerical results are used to investigate the effect of some parameters happening the control of JE infection, for justification of analytical results. © Korean Society for Informatics and Computational Applied Mathematics 2020.

Introduction

Japanese encephalitis (JE) is an illness of the brain caused through the Japanese encephalitis virus (JEV) and is transmitted from the pig population to the human population through bloodsucking insects Culex mosquitoes, Centers for Disease Control, (CDC). The JEV is maintained in a natural enzootic cycle between reservoir population (especially pigs), mosquitoes and humans. The transmission of JEV, for the most part, happens in the country in horticultural zones, where floodwater system and rice development are related [1-4]. In the areas where JEV is endemic the annual incident range from 1 to 10 cases per 10,000 person per year affected children are under 16 years of age [5]. The time duration of infection to illness in infected humans is about 5–15 days and initial symptoms are fever, headache and vomiting [3, 6–8]. Most humans who are infected with only a few symptoms or have no symptoms and JE can be more deadly in advanced stages, National Vector Borne Disease Control Program, (NVBDCP).

The first human case of JEV was registered in Japan in 1935 and in India in 1955 [3, 9]. Unfortunately, there is no specific treatment available that can cure Japanese encephalitis, including the management of some basic supportive medical care and complications in medicine. Those who experience a more severe form of the disease will have to undergo treatment at the hospital. The finding of JE disease ought to be performed by utilizing a JEV-explicit immunoglobulin M (IgM) counteracting agent catch ELISA for serum and cerebrospinal liquid (CSF). JEV-explicit IgM can be estimated in the CSF of most patients by four days after the commencement of side effects and in serum by seven days after inception. The prevention of JE is mainly based on the control of mosquitoes and the vaccination of the human population. To control the mosquito population, some methods like using insecticides, reducing the number of breeding sites on your property, by killing the mosquito larvae where the insect’s breed and bed nets are used in endemic regions. Mathematical modeling of infectious disease has a long history for understanding the effects of vaccination. Through the speedy rate of vaccine development, medical science has saved millions of lives. Smallpox, influenza, polio and lots of supplementary diseases are conquering examples in this regard [10]. In the present life, it is bizarre and uncommon for a youngster not to get an antibody. An individual can get accessible immunizations for a malady common in his general vicinity. Some sickness immunizations are as of now created and anybody can take the antibody if the specific malady is in danger for it. Currently, several groups of vaccines available in the world like inactivated Vero cell culture-derived JE vaccine, SA 14-14-2 live attenuated JE vaccine, Ixiaro JE vaccine, etc. But Ixiaro is the only licensed JE vaccine available in the United States [11]. It was avowed in march 2009 for the use of the human population aged more than 17 years. Later on, May 2013, avowed for the use of children whose aged lies between 2 months to 16 years [1, 6, 12]. The essential inoculation plans for Ixiaro is two portions intramuscularly administrated on 0–28 days. For youngsters from 2 months to 2 years, each dose is 0.25 ml, for grownups and kids for more than 3 years, each dose is 0.5 ml. Prior, In India, the mouse brain-derived JE antibody was utilized, yet the productivity was insufficient to meet the national prerequisite [13, 14]. The SA 14-14-2 live attenuated JE vaccine has been utilized in China for over 20 years with an effectiveness of 80% for single-dose and 97.5% for a two-doses routine with 1 year interim among youngsters under 15 years old [11, 14]. Thus, the effectiveness of the JE vaccine is more than 90%, but the time duration of protection with the JE vaccine is not clear, but its effectiveness appears to decrease over time [1, 6].

Japanese encephalitis spreading fast in all states of India except Meghalaya, Himachal Pradesh, Jammu and Kashmir, Dadra, Daman, Punjab, Diu, Nagar Haveli, Gujarat, Lakshadweep and Sikkim [1]. This disease is seasonal, most human cases reported in north India in May–October. The highest rate of this disease reported from the states of Uttar Pradesh, Bihar, Andhra Pradesh, Assam, Kerala, Goa, Haryana, Karnataka, Tamil Nadu and West Bengal [5, 9]. In modeling of such infectious disease, many authors have already worked on epidemic models to minimize the infection of JE and maximize the rate of recovered population [15]. In recent years, mathematical models have been used to understand the interplay between various traditional measures such as the mobility of epidemic outbreaks [16, 17] and vaccination [5, 18–20] quarantine, behavior change [21].

When the immunization is accessible, a general inquiry emerges; How to allocate and implement this immunization. Unquestionably, we can’t inoculate all people to control the sickness from the human population on account of social and moral issues, significant expenses, and so forth. In this manner, the legitimate conveyance and usage systems are critical to controlling JEV using vaccines. The adequacy of such a vaccination assignment procedure can be resolved through the investigation of a mathematical model. Right now, we need to illuminate this significant issue and would like to give a valuable rule to the arrangement producer.

To appropriately actualize vaccination battles in the human population, the recognized and insightful idea may be to inoculate individuals having a spot with explicit social affairs or territories vulnerable to contamination. On account of mosquito dependent infectious diseases, the vaccination can be an important tool to help control the spread of infection in such infectious disease and encephalitis is one of them. It is pointed out that in the above model, vaccination has been studied without considering the effective role of the mosquito population which depends on the human population-related factor. In this work, we construct and analyze a modified mathematical model of De et al. [9], by incorporating the effect of vaccination on the spread of JE, which is the most effective control strategy against JE. In addition, we use a realistic standard incidence rate for direct contact between susceptible and infectious. We might want to express that we don’t consider any structure in the model analyzed down right now, these references can be utilized to improve the present work.

Kar et al. [19] have been studied the SIR epidemic mathematical model to minimize the infection rate and to maximize the recovered rate of the human population by using vaccination program and optimal control strategy to SIR model. Singh et al. [22] have found that the development rate of the mosquito population relies upon the increasing density of ecological discharge by the human population. Also, the growth rate of the mosquito populations can be controlled by controlling the increasing density of ecological discharge. Iannelli et al. [23] build an SIS mathematical model to contemplate the effect of inoculation by expecting the immunization age, super-infection, immaculate immunization and furthermore talked about the asymptotic behavior of the equilibriums state. Raul Prealta et al. [24], developed an SVIR mathematical model to study the effect of the vaccine by considering the immune loss rate, which depends on the vaccine-age under the condition that during vaccination individuals have the same immunity.

In the present investigation, modifications to the SVIR model for JEV-affected humans and SI model for JEV-affected pigs [9] are considered by reducing the vaccinated compartment of a pig. The main objective is to analyze the effect of vaccination on the human population for the prevention of Japanese encephalitis with a standard incidence rate of mosquitoes and pigs. In the case of mosquitoes there is nothing to explain because it is obvious that vaccination, medicine is not possible for that compartment. Hence we consider susceptible, vaccinated, infected and recovered (SVIR) compartment for humans and susceptible, infected (SI) compartment of mosquitoes and pigs (cf. the Fig. 1). Disease transmission of the susceptible human population depends on the standard incidence rate of infected mosquitoes. Similarly, JEV spreads among the susceptible mosquito population when it comes in contact with infected human and pig populations. In this way, JE transmission for the susceptible mosquito population should be dependent on the rate of standard incidence on infected pig and human populations.

Fig. 1

Schematic diagram representing dynamics for the transmission of JE disease in humans, pigs and mosquitoes

The organization of the manuscript is as follows: In the next Sect. 2, notations and assumptions along with full explanation are provided. In Sect. 3, the mathematical formulation of the JE endemic model with a schematic diagram of the transmission of JEV puts in place and derive the model equations. We provide the region of biologically feasible solutions in Sect. 4, together with their positivity and boundedness. Dynamical behavior (existence of equilibrium states, basic reproduction number and stability) of the framework has been discussed in Sect. 5. In Sect. 6, vaccine-induced reproduction number and model parameters sensitivity investigation has been discussed to study the effects of controllable parameters of the JE model on the basic reproduction number and infected human population in Sect. 7 while the simulations and discussions are presented in Sect. 8 and conclusions in Sect. 9.

Notations and assumptions

In our proposed model, it is assumed that the JEV starts its life cycle from the pig population and then spreads into the human population through Culex mosquitoes. It is also assumed that all the newborn’s pigs and mosquitoes are susceptible, but newborn humans are not susceptible. In addition, JEV transmitted only in susceptible individuals. Many variables and constants(parameters) symbols have utilized in our proposed mathematical epidemiological model, including four compartments for humans, two compartments for pigs and mosquitoes, which are given in Tables 1 and 2.

Table 1

Model variables and their descriptions

Symbol	Descriptions
S_h(t)	Susceptible human population (people who are in danger of getting tainted by disease) at any time t
I_h(t)	Infected human population (the people who can transmit disease) at any time t
V_h(t)	Vaccinated human population (the human who was vaccinated and are currently immune, however not totally) at any time t
R_h(t)	Recovered human population (the human who is now immune but not completely) at any time t
S_m(t)	Susceptible mosquito population (able to compress the disease) at any time t
I_m(t)	Infected mosquito population (equipped for transmitting the disease to the human population) at any time t
S_p(t)	Susceptible pig population (pigs who are able to contract the JE disease) at any time t
I_p(t)	Infected pig population (pigs who are equipped for transmitting the JE illness to susceptible mosquitoes) at any time t

Table 2

Model parameters, its description and their values used for numerical calculation

Symbol	Description	Value	Source
α_1	The rate at which humans are vaccinated	0.50	[9]
λ_h	The rate at which humans are recruited	1.14	Assumed
λ_m	The rate at which susceptible mosquitoes are generated	1.51	Assumed
λ_p	The rate at which susceptible pigs are recruited	0.405	Assumed
μ_p	Natural mortality of pigs	0.01	Assumed
μ_h	Natural mortality of humans	0.001538	Assumed
μ_m	Natural mortality of mosquitoes	0.03	Assumed
β_1	The rate at which vaccinated human goes to susceptible class of human	0.001	[9]
B	Average biting rate on susceptible human by infected mosquitoes per day	0.7	[9]
C	Average biting rate on susceptible pig by infected mosquitoes per day	0.9	[9]
λ	Disease-related death rate of infected human population	0.0222	[9]
τ_1,τ_2	The rate at which infected human becomes susceptible without and with medicine respectively	0.5	[9]
1/η_h	Mean viremic period of JEV (in days)	4	[1]
δ	Effectiveness of medicine on the human population	0.10	[9]
β_mh	JEV transmission potential rate from infected mosquitoes to susceptible humans	0.2	[9]
β_pm	JEV transmission potential rate from an infected pigs to susceptible mosquitoes	0.3	[9]
β_mp	JEV transmission potential rate from an infected mosquitoes to the susceptible pigs	0.3	[9]
α_2	The rate at which infected pigs go to the susceptible class of pigs by medicine	0.5	[9]
u_1	Drug control of pig population	0.5	[9]
N_h	The absolute human population at time t	2500	Assumed
N_m	The absolute mosquito population at time t	12,500	Assumed
N_p	The absolute pig population at time t	50	Assumed
f_1	Rate at which infected human goes to susceptible class at time t	0.02	Assumed

Model formulation

In our present investigation, we consider susceptible, vaccinated, infected and recovered (SVIR) model for human population, susceptible, infected (SI) models for mosquito and pig population. The model sub-divides the total human population into sub-population of susceptible humans , vaccinated humans , infected humans and recovered humans , so that . The total mosquito population is divided into sub-population of susceptible mosquitoes and infected mosquitoes , so that . Also, the total pig population is sub-divided into susceptible pigs and infected pigs such that . In SVIR model, humans are recruited at vaccinated birth rate and unvaccinated birth rate . It is assumed that the JE spreads among the due to contribution of the [5]. Hence, the transmission of JE disease to susceptible humans should be a standard incidence rate of an infected mosquito population and it is in the form of . Similarly, the JE transmitted among susceptible mosquitoes () due to contact with infected pigs and human populations. Hence, the transmission of JE disease to susceptible mosquitoes should be standardized incidence rate of infected pig and human population and it is in the form of . The interaction between and , and are assumed to be rated B, C respectively. Ultimately, the JE spreads among due to contact with . Hence, the JE transmission to susceptible pigs should be a standard incidence rate of infected mosquitoes and it is in the form of . Also, we consider the proportion of the infected human population become recovered naturally and with medicine control at a rate of and respectively. After that, some of the recovered individuals go to susceptible class at rate . Humans in each compartment are assumed to die with a natural mortality and die with the disease at a rate of . In the SI model, the recruitment rate of mosquitoes and pigs is and respectively. The natural mortality of the mosquito population in each compartment is . Also, the pig populations in each compartment are assumed to die with natural mortality of . Further, we consider that in the whole modeling process the disease transmission is possible only in susceptible individual class, for other classes JE spreads among them after they become susceptible class only. So, we do not consider direct JE transmission from vaccinated humans class to infected human class. The schematic diagram for the dynamics of transmission of JE has appeared in Fig. 1.

The mathematical model is described by following system of ODEs:For our simplicity we can re-write Eqs. (1)–(8) in notations form as-with initial conditions  .

where,

Positivity and boundedness of solutions

Now, we present the positivity and boundedness of solutions of the system of Eqs. (1)–(8). For the system of Eqs. (1)–(8) to be biologically significant and well-proposed, essentially we have to show that all the solutions with non-negative initial conditions will remain positive for all time . In the theoretical study of disease transmission, the boundedness of the system measures that the system is biologically legitimate and well-behaved [16]. Biologically, validation of the model will be appeared by the accompanying theorem.

Theorem 4.1

If are non negative, then are also non-negative for all time . Moreoverwhere ,  ,  ,

Furthermore,and

The proof is excluded. The feasible region for the system of Eqs. (1)–(8) is given bywhere,andis positive invariant.

Dynamical behavior of the system

Equilibrium states and basic reproductive number

Clearly, the system (9)–(16) has following two equilibrium states, namely

(1) Disease or infection-free steady state .

where,The existence of equilibrium is obvious. It demonstrates that if the model has an infected population, then the equilibrium magnitude of susceptible human, mosquito and pig population will reach the value   respectively and vaccinated human population will remain at its equilibrium .

Let us define threshold parameter , asBiologically, is the average number of secondary Culex mosquitoes produced by single Culex mosquitoes [25]. To maintain the Culex mosquitoes in the environment, must be greater than one. To measure the disease transmission potential in such a population, the basic reproduction number can be established by using the concept of the next-generation matrix approach method. Since the system of Eqs. (9)–(16) are consistent, therefore the conditions for new infection and transition terms are given by the matriceswhere,then,For infection-free equilibrium state, the next generation matrix () is calculated as:The potent characteristic root of matrix is the basic reproduction number [26]. It is calculated asBiologically, is the average number of secondary cases generated by one primary JEV infected person, which has been introduced in the susceptible population, in which some individuals have been vaccinated. From Eq. (17), it is found that as the increase with the parameters  and , the value of decreases. Also the parameters B, C, affects the value of significantly. From the expression of , it is also found that vaccination rate affects the value of . It is imperative to take note that, the above parameters assume a significant job in the investigation of global threshold dynamics with regards to . Pronouncedly, if , the disease will stop after some time and disease will stay endemic and irresistible population will arrive at positive continued levels [2, 4, 27, 28]. Clearly, from the perspective of controlling the JE infection, one would normally need to decrease the . In this way, it is advantageous to research how we can successfully diminish by sensible usage of vaccines and different parameters. The variation of w.r.t. some are shown in the Figs. 2, 3, 4 and 5 (Table 3).

Fig. 2

This figure shows the variation of with respect to

Fig. 3

The figure shows the changes in with respect to B

Fig. 4

This figure shows the variation in with respect to

Fig. 5

This figure shows the variation in with respect to

Table 3

The impact of the basic reproduction number () as for parameters and where all the parameter esteem continue as before as in Table 2

Figures	Parameter	R_0	Explanation
2	μ_m>0.05	R_0<1	This signifies that if μ_m>0.05, R_0<1 as R_0 is a decreasing function for μ_m. Thus, if the decay rate of mosquitoes exceeds μ_m=0.05, disease can be eradicated from the population
3	B>1.4	R_0>1	This shows that if B>1.4, R_0>1 as R_0 is a increasing function for B. But if B<1.4, R_0<1 population can be make disease-free
4	μ_h>0.0005	R_0<1	From this, we found that μ_h also help to eradicate the infection from environment when μ_h>0.0005
5	μ_p>0.02	R_0<1	This figure signifies that disease can be control from the population when μ_p>0.02, because if μ_p increases, R_0 decreases

Local stability around

The local stability of infection-free equilibrium state () is determined after the characteristic roots of the Jacobian matrix of system (9)–(16) by following theorem:

Theorem 5.1

The infection-free equilibrium state () is locally asymptotically stable (LAS) if and unstable if  .

See “Appendix A” for proof of this theorem.

Remark 5.1

The local stability of will depend as the sign of and . The infection-free equilibrium is locally asymptotically stable if and . Biologically, implies the absence of Culex mosquito population. In this manner, the infection-free equilibrium must be steady without the Culex mosquito population, in any case, disease despite everything perseveres, regardless of whether the susceptible have no immediate contact with the infected population.

Global stability around

In this section, we have studied the global behavior of infection-free equilibrium point for the system of Eqs. (9)–(16).

Theorem 5.2

Assumed thatIf the the infection-free equilibrium state is globally asymptotically stable (GAS) on D.

See “Appendix B” for proof of this theorem.

(2) Endemic or positive equilibrium .

where,Using the value of Eqs. (19) and (21) in the expression for in the Eq. (20). This shows that the endemic equilibrium of this model satisfy the Eq.where,The root of Eq. (22) correspond to the infection-free equilibrium and the positive roots of the equationcorrespond to endemic equilibria point. From Eq. (23) we haveIf ,  and hence . Thus, the positive equilibrium point is asBut if ,  and hence , this gives imaginary value of . To verify above numerically, we use the set of parametric values from Table 2. Thus, the components of the endemic equilibrium point () are found as

Theorem 5.3

The endemic equilibrium state () is LAS if and only if       and   where,See “Appendix C” for proof of this theorem.

Remark 5.2

To verify Theorem 5.3 numerically, we use the set of parametric values from Table 2. Thus, the conditions are as

Theorem 5.4

The endemic steady state () is GAS in the region D, if the following inequalities are satisfied:See “Appendix D” for proof of this theorem.

Remark 5.3

To verify Theorem 5.4 numerically, we use the set of parametric values from Table 2. The simulations are done in the following seven cases for different initial starts and graphically displayed in Figs. 6 and 7. In this figure, the variation of susceptible humans with susceptible mosquitoes and vaccinated humans with infected humans have been shown respectively. The trajectories start with different initial points reach at the equilibrium point () and (). Therefore, we infer that the system of Eqs. (1)–(8) may be non-linear global asymptotically stable about equilibrium points () and () for the following set of parameters (Table 4).

Fig. 6

The plot represents the non-linear stability of ()

Fig. 7

The plot represents the non-linear stability of ()

Table 4

The components of endemic equilibrium points for non-linear global asymptotically stable

Sh∗(0)	Vh∗(0)	Ih∗(0)	Rh∗(0)	Sm∗(0)	Im∗(0)	Sp∗(0)	Ip∗(0)
275	207	20	30	35	40	17	20
250	207	20	30	25	40	17	20
301	207	20	30	30	40	17	20
301	195	23	30	35	40	17	20
301	200	15	30	35	40	17	20
301	203	17	30	35	40	17	20
301	207	20	30	35	40	17	20

Also for following set of parameter values, the conditions for the global asymptotical stability of the endemic equilibrium are satisfied.

Vaccine-induced reproduction number

From Eq. (17), we define the vaccine-induced reproduction number () asNow,From above it is clear that is a decreasing function of if . This signalize that the effect of vaccination in reducing vaccine-induced reproduction number. Also, when vaccination is zero i.e. then Eq. (25) can be written asFrom Eqs. (25) and (26), we can write the vaccination asIf then , therefore is locally asymptotically stable. Also if , then this is also the case in which locally asymptotically stable. From Eqs. (25) and (26) we can write,If then . From this we find , therefore no amount of vaccination can bring less then one.

Hence, defines the stationary values for the vaccine related reduction rate of infection. Therefore, if the vaccination level , then possibility of disease eradication from population can be maximized (Fig. 8).

Fig. 8

The graph represents the variation of vaccine-induced reproduction number with vaccination rate

Solving Eq. (17) for , we get a critical vaccination level , is given by

Here, is always positive, since . If , as is a decreasing function for . In this manner, in the event that the vaccination coverage surpasses the condition , at that point disease can be eradicated.

From Eq. (19), we also find thatThis critical vaccination level guarantees that the disease annihilation when inoculated people stay at its harmony. Along these lines, if the vaccination rate is with the end goal that , at that point disease annihilation is conceivable in the human population.

Model parameters sensitivity investigation

In this section, we want to examine the relationship between input parameters and the results of the proposed model using the reproduction number . To assist us with discovering the parameters that highly affect , we dissected the to decide if immunization, treatment of infectives and mortality of mosquitoes can prompt the successful disposal or control of the JE in the population. The parameters chosen to perform the sensitivity analysis were B, and  . We evaluate the sensitivity indices for    and all other parameters are same as in Table 2. For sensitivity analysis of the model, we definecalled normalized forward sensitivity index for a variable x and it depends on a parameter s [29]. The sensitivity indices of   w. r. t. the parameters B,    and  are calculated asThe fact is that signify that 1% increase in B, will produce 0.01743% increase in , keeping all other parameters fixed. When the parameters   and  increase by 1%, while keeping all others fixed then the values of decrease by 0.01587%, 0.001987% and 0.00495% respectively. Overall a lower value of is better because it increases the chance of disease eradication in the areas. Therefore, all the above prevention practices should be considered to reduce the parameter B, while increases the parameters   and , instead it should be favored. It is seen that the impact of vaccination on people assumes an overwhelming job in the control of the disease. On the off chance that the contact rate (B) between susceptible people and infected mosquitoes can be managed then Japanese encephalitis can be killed totally from the population. It is additionally seen that the infected human population consistently expanding as for B and it is progressively delicate for the spread of JE disease, (cf. Figs. 9, 10). That’s why we have chosen the parameter B only for the model sensitivity investigation.

Fig. 9

The surface plot shows the variation of with respect to and

Fig. 10

The surface plot shows the variation of with respect to B and

Figure 9 shows the sensitivity plot of the basic reproduction number for and . From this, we found that a higher vaccination rate and a low contact rate will keep . Also, in Fig. 10, the sensitivity plot of the regarding B and has appeared. This implies a higher estimation of and a low contact pace of B will keep .

Simulations and discussions

In this paragraph, we provide some simulations on the dynamics of the system (1)–(8). Moreover, to check how the human, mosquito and pig populations are affected by various parameters, we perform simulations by using the parameter values as shown in Table 2, and numerical results are manifested graphically in Figs. 11, 12, 13, 14, 15, 16, 17 and 18 by using MATLAB(R2014a) software. In Fig. 11, the variation of the human population with vaccination is shown and in Fig. 12, the variation of the human population without vaccination is shown. From these two figures, we can study the effect of vaccination on the human population. In the case of vaccination, the susceptible population is below 50 and the vaccinated population is above 50, whereas in case of without vaccination, the vaccinated level goes to zero and susceptible level increases. We also see that vaccination does not affect the infected and recovered human population. This signifies that vaccination only affects susceptible individuals if individuals infected, vaccination does not affect the infected population. Figures 13 and 14, show the variation of mosquito and pig population at an endemic equilibrium point. These two figures (Figs. 13, 14) show the stability of the endemic equilibrium when is larger than one and the threshold condition of the JE virus larger than one. This reflects that the virus can be dominant when the and are larger than one. We observed in Fig. 15, that due to the control strategy, the number of susceptible human decreases in society. This signifies that the spread of JE infection in society can also be controlled through the proper implementation of vaccines, medicine and some effective prevention. This strategy also shows that the effect of vaccine and medicine control on susceptible humans is almost the same. This reflects that some effective medical treatment and medicine can also control the spread of JE infection. Whereas in Fig. 16, the number of vaccinated human increases in the community. This shows that the community can make infection-free through the effective treatment and vaccination control of a human. From Fig. 17, we observed that due to the control strategy, initially decreases the infected human but after some time it stabilized. This shows that the vaccination strategy has no significant impact on the infected human. Finally, in Fig. 18 we observed that due to controlled strategy, the number of recovered human increases in the community but after some time it stabilized. This shows that infected humans can be immune (but not completely) through some basic supportive medical care and treatment.

Fig. 11

The graph shows the variation of human population with vaccination at endemic equilibrium

Fig. 12

The graph shows the variation of human population without vaccination at endemic equilibrium

Fig. 13

The figure shows the variation of mosquito population at endemic equilibrium

Fig. 14

The figure shows the variation of pig population at endemic equilibrium

Fig. 15

The figure shows the effect of susceptible human with no control, all control, vaccine control and medicine control

Fig. 16

The figure shows the effect of vaccinated human with no control, all control, vaccine control and medicine control

Fig. 17

The plot shows the effect of infected  human  with  no control, all control, vaccine control and medicine control

Fig. 18

The plot shows the effect of recovered human with no control, all control, vaccine control and medicine control

Conclusions

In this article, perspective on human, mosquito and pig population, a nonlinear scientific model for JEV has been proposed and examined to contemplate the impact of vaccination on control of infectious disease, Japanese encephalitis. The model displays two equilibria, to be a specific disease-free and endemic equilibrium. The first equilibrium exists without any conditions but the second equilibrium exists if the threshold parameter . We have discovered another threshold parameter in terms of vaccine, vaccine-induced reproduction number which is, if not as much as unity, the disease ceases to exist, without infected population otherwise the infection is kept up in the population and disease turns out to be progressively endemic. The system of the differential equations has a unique infection-free steady state which is locally asymptotically stable if and unstable if . Furthermore if approaches to 1, becomes unstable but exist. If then exists and it is locally asymptotically stable under certain condition. Also for both the equilibrium point global stability have to show by using Liapunov’s function in the region of attraction under certain conditions. Model analysis shows, if the vaccination program has adequate capacity to mobilize a high degree of population, there is an abundance of mosquitoes in the environment, then the disease can be controlled. Moreover, it is found that vaccine coverage has a substantial effect on the basic reproduction number, . Also, by continuous efforts and effectiveness of vaccine coverage, the disease can be eradicated. Also, found by sensitivity investigation, if 1% increase in the parameters keeping all others parameters fixed, then the value of decreases by 0.01587%, 0.001987% and 0.00495% respectively. Also, if we increase 1% in B then increases by 0.01743%. This signifies that the parameter B is more sensitive to the transmission of JEV. Then again, in the event that the vaccination rate is over a critical level , at that point the JE doesn’t take a type of pestilence gave the density of mosquito population stays at its equilibrium level. During the entire investigation of the JEV model, it is discovered that vaccination is the main long haul control measure against JE in the population. The mosquitoes reduction strategies may incorporate the end of mosquito breeding sites, keeping environmental factors perfect and clean. On the other hand, vaccinating the susceptible individuals in the human population with the end goal that vaccination level is kept up over the critical vaccination level and the vaccine efficacy is sufficiently high so that vaccinated individuals are not infected by infected mosquitoes. Along these lines, our present paper has been investigated various sorts of control procedures and their adequacy in controlling the JE disease by utilizing the steadiness hypothesis of differential conditions and the basic reproduction number. From the clinical perspective, we accept that our model examination is gainful for persons and control of the JE disease.