Abstract
Objectives: Vector-borne diseases speedily infest the human population. The control techniques must be applied to such ailment and work swiftly. We proposed a compartmental model of dengue disease by incorporating the standard incidence relation between susceptible vectors and infected humans to see the effect of manageable parameters of the model on the basic reproduction number.
Methods: We compute the basic reproduction number by using the next -generation matrix method to study the local and global stability of disease free and endemic equilibrium points along with sensitivity analysis of the model.
Results: Numerical results are explored the global behaviourism of disease-free/endemic state for a choice of arbitrary initial conditions. Also, the biting rate of vector population has more influence on the basic reproduction number as compared the other parameters.
Conclusion: In this paper, shows that controlling the route of transmission of this disease is very important if we plan to restrict the transmission potential.
Introduction
One of the major concerns of the present era is the birth of a new class of diseases, especially vector-borne diseases like Japanese encephalitis, dengue, malaria, yellow fever etc. Earlier, the sources of these diseases were limited. But, new sources have evolved to a similar kind of ailment, which critically hits the human population. Therefore, one of the significant ways to understand the dynamic phenomenon of disease transmission is mathematical modelling. It helps in determining the intervention techniques for disease control and reduce the mortality rate of humans. Nowadays, a common viral disease in the tropical region of the world is dengue. Previously, it is known as dengue hemorrhagic fever. In the 1950s, the Philippines and Thailand countries were plagued, dengue epidemics and recognized for the first time. Currently, Latin, American and Asian countries are the most affected ones causing hospitalization and death of people from almost all age groups. Here, the mosquitoes act as transmitters, and humans act as receiver, describing the life cycle of the dengue virus (DENV). Aedes mosquito is the transmitter of dengue after an infectious bite by anyone of the given dengue virus (DENV-1 to 4). Symptoms of dengue appear in 3–14 days after the infectious bite of mosquitoes. Mostly, adult and young children are affected by viral dengue. Dengue causes grievous fever-like illness due to mosquito-borne viral infection. Sometimes, dengue becomes a potentially lethal convolution. The incidence of dengue has increased 30-fold over the last 50 years. Distinct serotypes (DENV-1 to 4) are the disease origin, belonging to the genus Flavivirus family. If the human recovered from dengue virus, they are never affected again by the virus and receive immunity for a lifetime against that specific virus.
Only partial and momentary protection produced such immunity against the infection due to distinct serotypes of the virus. In comparison to transmitter, mosquitoes have the only selective method to eradicate the dengue virus. According to WHO (World Health Organisation), the estimation of dengue infections worldwide touches the range of 50–100 million, subjecting nearly half of the population all over the world at grave risk. There is no specific vaccine available for dengue. One recent observation in 2012 predicted that 3.9 billion people of around 128 countries are at potential risk of infection with DENV. Another contemplation in 2013 indicates that out of 390 million, 96 million (67–136 million) people infected through dengue per year and manifest clinically (with any severity of disease) (WHO 2020).
To research the mechanisms of dengue transmission, several mathematical models have been proposed and studied from Ghosh, Tiwari, and Chttopadhyay (2019) to Agusto and Khan (2018). In Ghosh, Tiwari, and Chttopadhyay (2019), the authors analyzed a compartmental model and studied dynamic behaviour of the model. Also, discuss ACF (active case finding) for asymptomatic and symptomatic individuals plays a crucial role in reducing disease prevalence. In Zhu and Xu (2019) studied an SIS-SI dengue model associated with a periodic environment to discuss temporal periodicity and spatial heterogeneity. Driessch and Watmough (2002) utilized a system of ODE to formulate a compartmental model for disease transmission and analyzed the stability and instability at R 0 < 1 and R 0 > 1 respectively. In Busenberg and Driessch (1990) and Esteva and Vargas (1999), the authors analyzed a dengue transmission model with varying size and applied competitive system theory to discuss the equilibrium’s global asymptotic stability. Esteva and Vargas (1998, 2009 analyzed the dynamics of dengue disease with the help of interrupted feeding in the role of vertical and mechanical transmission. Also, discussed a transmission model of dengue with varying vector populations and the constant human population. They also established the global stability of endemic equilibrium with the threshold condition’s help to measure the population of vector control parameters. In Linda and Amy (2000), the authors compared the deterministic and stochastic SIS and SIR model in discrete time. Misra and Sharma (2013), evaluate the effect of media awareness by advertising programs on vector-borne diseases. A nonlinear mathematical model is suggested and evaluated in this article. Tewa, Dimi, and Bowong (2009) studied a dengue fever model when only one serotype existed for the dynamics and showed the global asymptotic stable condition of disease-free and endemic equilibrium. In Kar and Jana (2009), the authors proposed a vector-borne disease model with three types of controls for eradicating disease and formulated the optimal control problem. They solved the model by considering control parameters as time-dependent. Syafruddin and Noorani (2013) constructed a model and discussed the system’s local and global stability behaviour. Nur, Rachman, and Abdal (2018) analyzed the ODE’s SIR model system with climate factors for dengue fever transmission characteristics in closed population and explained the equilibrium’s stability using Lyapunov. In Panja, Mondal, and Chattopadhyay (2018), authors identify time-dependent cost-effective control strategies to minimize dengue spread, used the optimal control theory. In this article, the authors proposed a dengue transmission model using a control variable, namely insecticide and vaccination and develop the result based on control strategy (Agusto and Khan 2018). In Bonyah et al. (2017), authors analyze an SEIR Zika epidemic model and investigate the model with basic mathematical results. Treatment is not possible if the number of infectives exceeds a certain threshold, as it is in developing countries such as India, China, Swaziland, and Brazil. In this present paper, susceptible human have categorized by considering two different classes. We have incorporated a modified model, as discussed in Syafruddin and Noorani (2013) and compared it to the single-compartment model of susceptible beings with multiple compartments. Lyapunov function has applied to this model to examine the local and global stability analysis and the basic reproduction number. The application of sensitivity analysis has been analyzing to see the manageable parameters of the model and simulated the variables.
The structure of paper is as follows: the formation and formulation of the model described in Section “Mathematical formulation”. The contained covered in Sections “Positivity and boundedness of the system”, and “Qualitative analysis of model” are analysis of positivity and boundedness of the system and basic reproduction number of the model along with estimation of the disease-free, endemic equilibrium point. Stability analysis of the model and performance of sensitivity analysis are in Sections “Stability analysis”, and “Sensitivity analysis”. The simulation results, as seen in Section “Numerical simulation”. Finally, a detailed discussion and conclusion is presented in Section “Conclusions”.
Mathematical formulation
At any time t > 0, the total population
Variables used in the framework are listed.
Parameters | Significance |
---|---|
N H | Total human population |
|
Susceptible human at high risk |
|
Susceptible human at low risk |
E H | Exposed human |
I H | Infected human |
R H | Recovered human |
N v | Total vector population |
S v | Susceptible vector (mosquito) |
E v | Exposed vector |
I v | Infected vector |

Compartmental model.
On the basis of transmission dynamics mathematical formulation for above model is:
The meaning of parameters used in (1)–(8) are listed in Tables 1 and 2.
Biological significance of parameters.
Parameters | Significance |
---|---|
π v | Mobilization rate of susceptible vector population |
β | Biting rate of vector per individual |
|
Natural death rate of humans, vector respectively |
|
Mobilization rate from infected human to suspectable vector |
γ v | Extrinsic death rate of vector |
|
Intrinsic death rate of humans |
|
Mobilization rate of suspectable human population |
q | Recovery rate of disease |
r | Section of first time admitted person joining the susceptible class at high risk |
α v | Transmission rate from infected mosquitoes to susceptible human |
θ | Relative possibility of low risk susceptible with reference to high risk susceptible |
The above system can also be written as:
with initial conditions,
Positivity and boundedness of the system
System (1)–(8) has positively invariant and bounded solution in the closed set.
Lemma 1
The region of existence for all solutions initiating in the positive region is recommended by set ω
where, Z
1 = max
Proof
System (1)–(8) can be written as:
where,
The vector
Using a standard comparison theorem used in Lakshmikantham, Leela, and Martynyuk (1989) we have,
Assume that Z
2 = max
Assume that Z
4 = max
Therefore, all the feasible solution of system (1)–(8) are entered in the region ω, which means that it is the attracting set for region ω.
Qualitative analysis of model
Basic reproduction number
In this section, the basic reproduction number (R 0) is determined by using the next-generation matrix method as in Driessch and Watmough (2002), and Ullah et al. (2020). It follows that the system will be locally asymptotically stable if R 0 < 1, system will be unstable if R 0 > 1. The process to find the basic reproduction number through Jacobian matrices F and V are:
where,
and
then,
then (F.V −1) is calculated as:
So, the basic reproduction number R 0 of system (9)–(16) is defined by maximum eigenvalue of the matrix (F.V −1) and it is given by
Disease-free equilibrium
For system (9)–(16), the disease-free equilibrium point is
Endemic equilibrium
For Eqs. (1)–(8) the endemic-equilibrium point is
A newer variable Z is defined as:
From the equilibrium equations (23) of the system (1)–(8), we have
the total human population is given by;
Now using Eqs. (19)–(22) in Eq. (18), we get the following quadratic equation in Z:
where,
On applying Descartes’ rule of sign, we discussed the various chances for positive roots of Eq. (23) in Table 3. Here B 0 = 0, therefore it has no effect on evaluation of roots as in Driessche (2017) and Ghosh, Tiwari, and Chttopadhyay (2019).
For R 0 < 1 and R 0 > 1 number of possible positive roots of Eq. (23).
S.N. | B 3 B 2 B 1 B 0 | R 0 | Change of sign | Possibility of positive root |
---|---|---|---|---|
1 | + + + 0 | R 0 < 1 | 0 | 0 |
+ + − 0 | R 0 > 1 | 1 | 1 | |
2 | + − + 0 | R 0 < 1 | 2 | 0, 2 |
+ − − 0 | R 0 > 1 | 1 | 1 | |
3 | − + + 0 | R 0 < 1 | 1 | 1 |
− + − 0 | R 0 > 1 | 2 | 1, 2 | |
4 | − − + 0 | R 0 < 1 | 1 | 1 |
− − − 0 | R 0 > 1 | 0 | 0 |
Stability analysis
In this section, the stability analysis of disease-free and endemic equilibrium is carried out. The basic reproduction number (R 0) is used to discuss both the equilibrium points regarding local stability and global stability. Stability analysis can perform, using the sign of eigenvalues of the Jacobian matrix’s characteristic equation as in Esteva and Vargas (1998) and Lakshmikantham, Leela, and Martynyuk (1989).
Local stability around equilibrium point
Now in this section, we are going to analyze the stability of two equilibrium points W 0 and W 1. The outcome of stability analysis of these equilibria are shown in the following two theorems.
Theorem 1
For R 0 < 1 the disease-free equilibrium (W 0) is locally asymptotically stable and unstable if R 0 > 1.
Proof
The Jacobian matrix W 01 for the model system (9)–(16) corresponding to equilibrium W 0 is given by
The four eigen values of W
01 at disease-free equilibrium are −μ
v
,
where,
v 1 = d 12 + d 14 + d 15 + μ v
v 2 = d 12 d 14 + d 12 d 15 + d 14 d 15 + d 12 μ v + d 14 μ v + d 15 μ v
v 3 = d 12 d 14 d 15 + d 12 d 14 μ v + d 12 d 15 μ v + d 14 d 15 μ v
Here, v
1 > 0, v
1
v
2 − v
3 > 0,
So, if R 0 < 1, all of the above conditions are fulfilled. Hence by Routh–Hurwitz criterion the disease-free equilibrium point W 0 is locally asymptotically stable otherwise, it is unstable.
Theorem 2
The endemic equilibrium W 1 for the system (9)–(16) is locally asymptotically stable under some conditions stated in the proof of this theorem.
Proof
The Jacobian matrix W 10 for the model system (9)–(16) corresponding to equilibrium W 1 is given by
where,
Four eigenvalues of the above matrix are −A
1,
where,
Here, w
1 > 0, w
1
w
2 − w
3 > 0,
For this, if w
4 > 0 and
Global stability around equilibrium point
Here, we are going to analyze the stability of two equilibrium points W 0 and W 1. The outcomes of stability analysis of these equilibrium points are shown in the following two theorems.
Theorem 1
If R 0 ≤ 1, the disease-free equillibrium (W 0) is globally asymptotically stable on ω with assumption
Proof
Consider the Lyapunov function below:
Differentiating with respect to time t, get:
On solving further get:
Using the conditions (26) and (27), also we have
Hence the conditions (26) and (27) ensures that G′(t) ≤ 0. By using Lasalle’s extension to Lyapunov’s method (Bonyah et al. 2017; Ullah, Khan, and Farooq 2018a), the limit set of each solution is contained in the largest invariant set for which
Theorem 2
The endemic equilibrium (W 1) is globally asymptotically stable in ω provided the following inequalities hold
Proof
Consider the following positive definite function, to study the global stability behaviour of the endemic equilibrium (W 1):
Differentiating with respect to time t,
With the help of Eqs. (9)–(16) on rearranging the terms we get:
For system (9)–(16) inside the region of attraction ω, V′(t) can be turned negative definite if conditions (28)–(31) hold as in Baniya and Keval (2020).
Sensitivity analysis
Here, we deal with the sensitivity analysis of the model to see the effect of manageable parameters of the model on the basic reproduction number (R 0) as in Zheng and Nie (2018) and Agusto and Khan (2018). The basic reproduction number behave like an important epidemiological parameter, and sensitivity analysis gives us outputs of the model, divided for quantitative values through variations in different parameter values. Here, a specified space of parameter examines the response model from parameter variation to output variables. The sensitivity index for a variable γ, which depends differentially on a parameter σ, is defined as:
We measure elasticities
where,
From the above expression, we see that

Influence of

Influence of

Influence of β and α v on R 0.

Influence of β and μ v on R 0.

Influence of

Influence of β and q on R 0.

Influence of β and π v on R 0.

Influence of β and γ v on R 0.

PRCC results for significance of parameters in R 0.
Sensitivity indices of R 0 to the parameter values of the model.
Variable | Parameter | Sensitivity index |
---|---|---|
R 0 | β | 0.817806 |
α v | 0.408903 | |
γ v | 0.137216 | |
|
0.189935 | |
μ v | −0.955022 | |
|
−0.412268 | |
π v | 0.408903 | |
q | −0.407538 |
Numerical simulation
This section aims to find the numerical solution of ODEs Eqs. (1)–(8) with the help of a medium order method in a finite interval of time using Matlab (Derouich and Boutayeb 2006; Khan, Ullah, and Farooq 2018; Ullah et al. 2020; Ullah, Khan, and Farooq 2018b), confirming results acquired with analysis and multiple observations give a deeper understanding of the model. For this used parameter values are provided in Table 5. These parameters value satisfy the conditions for the existence of endemic equilibrium, W
1 are obtained as
Epidemiological parameters and their values with unit.
Parameters | Values (Ghosh, Tiwari, and Chttopadhyay 2019) | Unit |
---|---|---|
π v | 350,000 | day−1 |
β | 2 | day−1 |
|
0.75 | day−1 |
μ v | 0.025(Assumed) | day−1 |
γ v | 0.0495(Assumed) | day−1 |
|
273,600 | day−1 |
r | 0.25 | day−1 |
α v | 0.75 | day−1 |
|
0.0154 | day−1 |
|
3.3 | day−1 |
θ | 0.5 | day−1 |
q | 4.5972 | day−1 |
The value of basic reproduction number (R 0) is got as 0.001835 which is less than one as our requirement. Stability conditions are also fully desired for these parameter values. The matrix’s eigenvalues are obtained for W 01 as −0.025, −0.0154 (multiplicity 3), −4.6130, −3.3147, −0.08589 and −0.0140 the polynomial roots satisfied the given stability conditions. It is articulate that all the eigenvalues are negative. Similarly, for the matrix W 10, found the condition for basic reproduction number R 0 > 3.364 and the matrix eigenvalues are −0.0473, −0.0154, −0.1, −0.8211, −4.6380, −3.2928, −0.0665 and −0.0322.
To show the nonlinear stability behaviour of disease prevalence in their respected population. For this, we have taken 110 units for susceptible vector and exposed vector, 800 units for infected vector. We plotted the solution trajectories at different initial starts using the above set of parameter values. These trajectories are shown in Figures 11, 12, and 13 respectively. In Figure 11, the variation of S v w.r.t the time t for the different values of transmission probability α v = 0.55, 0.65, 0.75, 0.85 and 0.95 are shown properly. It suggests that the population of susceptible vectors increase with a decrease in transmission probability. Therefore, to reduce the susceptible vector population’s size, we have to increase the transmission between infected mosquito to susceptible human. Also, in Figure 12, it is interesting to note that for different values of α v , the exposed vector population become zero in 27 unit of time whenever, if we decrease the transmission rate (α v ), it increases up to 40 unit of time. It indicates that if we decrease the rate of mosquito, the exposed vector become increases. Further, the variation of infected vector w.r.t time t for different values of β is given in Figure 13. It is evident that as the value of β increases, the total infected population (I v ) also increases. So, we infer that by multiple biting of mosquitoes, abundant of infected vector rapidly.

The plot of the susceptible vector S v with different value of α v .

The behavior of exposed vectors with α v .

The total number of infected individuals for different values of β.
Conclusions
Vector-borne disease outbreaks have ruined many countries. From which, focused on assessing the dynamics of dengue fever in this article. We proposed a mathematical model, modelling them and incorporating the effect of the variable human with exponential growth. We discussed the positivity and boundedness of the system by Metzler property and found the basic reproduction number (R 0) that control the transmission of disease and the growth rate of the infected human population. The sensitivity of R 0 obtained very sensitive for the model’s parameters corresponding to transmission rate, mobilization rate of vector, biting rate, and human death rate. Numerical simulations with different parameter settings illustrate the succession of epidemics, their amplitudes and justify theoretical findings. So, in this paper, main objective is to understand the dynamics of the dengue virus with variable human and investigate the model stability criteria for the dynamics of the dengue virus. Based on stability results and the basic reproduction number, we can say that our model is very beneficial for us. Also, it discussed the sensitivity and numerical simulation of the system. Thus we can conclude that our model gives us the most reliable results on changing the nonlinear model’s complexity by incorporating more epidemiological parameters, more saturating functions with respect to different virus types of the disease and mode of transmission.
According to the medical point of view, we believe that our model analysis is beneficial for humans also. It would assist in controlling the spreads of dengue disease.
-
Research funding: Not applicable.
-
Author contribution: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
-
Competing interest: Not applicable.
-
Informed consent: Not applicable.
-
Ethical approval: Not applicable.
References
Abate, A., A. Tiwari, and S. Sastry. 2009. “Box Invariance in Biologically-Inspired Dynamical Systems.” Automatica 45 (7): 1601–10. https://doi.org/10.1016/j.automatica.2009.02.028.Search in Google Scholar
Agusto, F., and M. Khan. 2018. “Optimal Control Strategies for Dengue Transmission in Pakistan.” Mathematical Biosciences 305: 102–21. https://doi.org/10.1016/j.mbs.2018.09.007.Search in Google Scholar
Altaf Khana, M., S. Ullah, and M. Farooqb. 2018. “A New Fractional Model for Tuberculosis with Relapse via Atangana–Baleanu Derivative.” Chaos, Solitons & Fractals 116: 227–38. https://doi.org/10.1016/j.chaos.2018.09.039.Search in Google Scholar
Baniya, V., and R. Keval. 2020. “The Influence of Vaccination on the Control of JE with a Standard Incidence Rate of Mosquitoes, Pigs and Humans.” Journal of Applied Mathematics and Computing 64 (1–2): 519–50. https://doi.org/10.1007/s12190-020-01367-y.Search in Google Scholar
Bonyah, E., M. A. Khan, K. O. Okosun, and S. Islam. 2017. “A Theoretical Model for Zika Virus Transmission.” PLoS One 12 (10): e0185540. https://doi.org/10.1371/journal.pone.0185540.Search in Google Scholar
Busenberg, S., and P. V. Driessch. 1990. “Analysis of a Disease Transmission Model in a Population with Varying Size.” Journal of Mathematical Biology 28 (3): 257–70. https://doi.org/10.1007/bf00178776.Search in Google Scholar
Derouich, M., and A. Boutayeb. 2006. “Dengue Fever: Mathematical Modelling and Computer Simulation.” Applied Mathematics and Computation 177 (2): 528–44. https://doi.org/10.1016/j.amc.2005.11.031.Search in Google Scholar
Driessche, P. V. 2017. “Reproduction Number of Infection Disease Model.” Infectious Disease Modelling 2 (3): 288–303. https://doi.org/10.1016/j.idm.2017.06.002.Search in Google Scholar
Driessch, P. V., and J. Watmough. 2002. “Reproduction Number and Subthreshold Endemic Equilibria for Compartmental Models of Disease Transmission.” Mathematical Biosciences 180 (1–2): 29–48. https://doi.org/10.1016/s0025-5564(02)00108-6.Search in Google Scholar
Esteva, L., and C. Vargas. 1998. “Analysis of a Dengue Disease Transnmission Model.” Mathematical Biosciences 150 (2): 131–51. https://doi.org/10.1016/s0025-5564(98)10003-2.Search in Google Scholar
Esteva, L., and C. Vargas. 1999. “A Model for Dengue Disease with Variable Human Population.” Journal of Mathematical Biology 38 (3): 220–40. https://doi.org/10.1007/s002850050147.Search in Google Scholar
Esteva, L., and C. Vargas. 2009. “Influence of Vertical and Mechanical Transmission on the Dynamics of Dengue Disease.” Mathematical Biosciences 167 (1): 51–64.10.1016/S0025-5564(00)00024-9Search in Google Scholar
Ghosh, I., P. K. Tiwari, and J. Chttopadhyay. 2019. “Effect of Active Case Finding on Dengue Control: Implication From a Mathematical Model.” Journal of Theoretical Biology 464: 50–62. https://doi.org/10.1007/978-981-13-7888-1.Search in Google Scholar
Kar, T. K., and S. Jana. 2009. “Application of Three Controls Optimally in a Vector-Borne Disease- A Mathematical Study.” Communications in Nonlinear Science and Numerical Simulation 18 (10): 2868–84. https://doi.org/10.1016/j.cnsns.2013.01.022.Search in Google Scholar
Khan, M. A., and A. Atangana. 2020. “Modeling the Dynamics of Novel Corona-Virus (2019-nCov) with Fractional Derivative.” Alexandria Engineering Journal 59 (4): 2379–89. https://doi.org/10.1016/j.aej.2020.02.033.Search in Google Scholar
Khan, M. A., Z. Hammouch, and D. Baleanu. 2019. “Modeling the Dynamics of Hepatitis E via the Caputo-Fabrizio Derivative.” Mathematical Modelling of Natural Phenomena 14 (3): 311. https://doi.org/10.1051/mmnp/2018074.Search in Google Scholar
Khan, M. K., S. Ullah, and M. Farooq. 2018. “A New Factional Model for Tuberculosis with Relapse via Atangana-Baleanu Derivative.” Chaos, Solitons & Fractals 116: 227–38. https://doi.org/10.1016/j.chaos.2018.09.039.Search in Google Scholar
Lakshmikantham, V., S. Leela, and A. A. Martynyuk. 1989. Stability analysis of Nonlinear Systems. New York: Marcel Dekker.10.1142/1192Search in Google Scholar
Linda, J. A., and M. B. Amy. 2000. “Comparison of Deterministic and Stochastic SIS and SIR Models Discrete-Time.” Mathematical Biosciences 163 (1): 1–33. https://doi.org/10.1016/s0025-5564(99)00047-4.Search in Google Scholar
Misra, A., and A. Sharma. 2013. “A Mathematical Model For Control of Vector-Borne Disease through Media Campaigns.” Discrete & Continuous Dynamical Systems - Series B 18 (7): 1909–27. https://doi.org/10.3934/dcdsb.2013.18.1909.Search in Google Scholar
Nur, W., H. Rachman, and M. Abdal. 2018. “SIR Model Analysis For Transmission of Dengue Fever Disease With Climate Lyapunov Function.” Journal of Physics: Conference Series 1028: 012117. https://doi.org/10.1088/1742-6596/1028/1/012117.Search in Google Scholar
Panja, P., K. Mondal, and J. Chattopadhyay. 2018. “Stability, Bifurcation and Optimal Control Analysis of a Malaria Model in a Periodic Environment.” International Journal of Nonlinear Sciences and Numerical Simulation 19 (6): 627–42. https://doi.org/10.1515/ijnsns-2017-0221.Search in Google Scholar
Syafruddin, S., and M. Noorani. 2013. “Lyapunov Function of SIR and SEIR Model for Transmission of Dengue Fever Disease.” International Journal of Simulation and Process Modelling 8 (2/3): 177. https://doi.org/10.1504/ijspm.2013.057544.Search in Google Scholar
Tewa, J., J. L. Dimi, and S. Bowong. 2009. “Lyapunov Function for a Dengue Disease Transmission Model.” Chaos, Solitons & Fractals 39: 936–41. https://doi.org/10.1016/j.chaos.2007.01.069.Search in Google Scholar
Ullah, S., M. A. Khan, M. Farooq, and E. O. Alzahrani. 2020. “A Fractional Model for the Dynamics of Tuberculosis (TB) using Atangana-Baleanu derivative.” Discrete & Continuous Dynamical Systems - Series B 13 (3): 937–56. https://doi.org/10.3934/dcdss.2020055.Search in Google Scholar
Ullah, S., M. A. Khan, and M. Farooq. 2018a. “Modelling and Analysis the Fractional HBV Model with Atangana- Baleanu Derivative.” The European Physical Journal Plus 133 (8), https://doi.org/10.1140/epjp/i2018-12120-1.10.1140/epjp/i2018-12120-1Search in Google Scholar
Ullah, S., M. A. Khan, and M. Farooq. 2018b. “A New Fractional Model for the Dynamics of Hepatitis B Virus Using Caputo-Fabrizio Derivative.” The European Physical Journal Plus 133 (6), https://doi.org/10.1140/epjp/i2018-12072-4.10.1140/epjp/i2018-12072-4Search in Google Scholar
WHO. 2020. Dengue Control. World Health Organization. https://www.who.int/denguecontrol/human/en.Search in Google Scholar
Zheng, T., and L. Nie. 2018. “Modelling the Transmission Dynamics of Two-Strain Dengue In the Presence Awareness and Vector Control.” Journal of Theoretical Biology 443: 82–91. https://doi.org/10.1016/j.jtbi.2018.01.017.Search in Google Scholar
Zhu, M., and Y. Xu. 2019. “A Time-Periodic Dengue Fever Model in a Heterogeneous Enviroment.” Mathematics and Computers in Simulation 155: 115–29. https://doi.org/10.1016/j.matcom.2017.12.008.Search in Google Scholar
© 2021 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Research Articles
- Determinants of birth-intervals in Algeria: a semi-Markov model analysis
- A simplified approach to bias estimation for correlations
- Gamma frailty model for survival risk estimation: an application to cancer data
- Analysis for transmission of dengue disease with different class of human population
- Quantifying the influence of location of residence on blood pressure in urbanising South India: a path analysis with multiple mediators
- Mixed methods to assess the use of rare illicit psychoactive substances: a case study
- Reliability of fetal–infant mortality rates in perinatal periods of risk (PPOR) analysis
- Sampling from networks: respondent-driven sampling
- Reviewer Acknowledgment
- Reviewer acknowledgment
- Tutorial
- A guide to value of information methods for prioritising research in health impact modelling
Articles in the same Issue
- Research Articles
- Determinants of birth-intervals in Algeria: a semi-Markov model analysis
- A simplified approach to bias estimation for correlations
- Gamma frailty model for survival risk estimation: an application to cancer data
- Analysis for transmission of dengue disease with different class of human population
- Quantifying the influence of location of residence on blood pressure in urbanising South India: a path analysis with multiple mediators
- Mixed methods to assess the use of rare illicit psychoactive substances: a case study
- Reliability of fetal–infant mortality rates in perinatal periods of risk (PPOR) analysis
- Sampling from networks: respondent-driven sampling
- Reviewer Acknowledgment
- Reviewer acknowledgment
- Tutorial
- A guide to value of information methods for prioritising research in health impact modelling