A singularly perturbed HIV model with treatment and antigenic variation

  • We study the long term dynamics and the multiscale aspects of a within-host HIV model that takes into account both mutation and treatment with enzyme inhibitors. This model generalizes a number of other models that have been extensively used to describe the HIV dynamics. Since the free virus dynamics occur on a much faster time-scale than cell dynamics, the model has two intrinsic time scales and should be viewed as a singularly perturbed system. Using Tikhonov's theorem we prove that the model can be approximated by a lower dimensional nonlinear model. Furthermore, we show that this reduced system is globally asymptotically stable by using Lyapunov's stability theory.

    Citation: Nara Bobko, Jorge P. Zubelli. A singularly perturbed HIV model with treatment and antigenic variation[J]. Mathematical Biosciences and Engineering, 2015, 12(1): 1-21. doi: 10.3934/mbe.2015.12.1

    Related Papers:

    [1] Saima Akter, Zhen Jin . Simulations and fractional modeling of dengue transmission in Bangladesh. Mathematical Biosciences and Engineering, 2023, 20(6): 9891-9922. doi: 10.3934/mbe.2023434
    [2] Vanessa Steindorf, Sergio Oliva, Jianhong Wu . Cross immunity protection and antibody-dependent enhancement in a distributed delay dynamic model. Mathematical Biosciences and Engineering, 2022, 19(3): 2950-2984. doi: 10.3934/mbe.2022136
    [3] Haitao Song, Dan Tian, Chunhua Shan . Modeling the effect of temperature on dengue virus transmission with periodic delay differential equations. Mathematical Biosciences and Engineering, 2020, 17(4): 4147-4164. doi: 10.3934/mbe.2020230
    [4] Biao Tang, Weike Zhou, Yanni Xiao, Jianhong Wu . Implication of sexual transmission of Zika on dengue and Zika outbreaks. Mathematical Biosciences and Engineering, 2019, 16(5): 5092-5113. doi: 10.3934/mbe.2019256
    [5] Zuohuan Zheng, Kaihua Wang, Daiyu Yang, Feifei Yin, Dingwei Sun, Weiyan Yu, Jialun Lin, Ying Liu, Changkuan Chen, Zehui Yang, Biao Wu . Factors affecting the transmission of dengue fever in Haikou city in 2019. Mathematical Biosciences and Engineering, 2023, 20(9): 16045-16059. doi: 10.3934/mbe.2023716
    [6] Bo Zheng, Lihong Chen, Qiwen Sun . Analyzing the control of dengue by releasing Wolbachia-infected male mosquitoes through a delay differential equation model. Mathematical Biosciences and Engineering, 2019, 16(5): 5531-5550. doi: 10.3934/mbe.2019275
    [7] Chen Liang, Hai-Feng Huo, Hong Xiang . Modelling mosquito population suppression based on competition system with strong and weak Allee effect. Mathematical Biosciences and Engineering, 2024, 21(4): 5227-5249. doi: 10.3934/mbe.2024231
    [8] David Cuesta-Frau, Pau Miró-Martínez, Sandra Oltra-Crespo, Antonio Molina-Picó, Pradeepa H. Dakappa, Chakrapani Mahabala, Borja Vargas, Paula González . Classification of fever patterns using a single extracted entropy feature: A feasibility study based on Sample Entropy. Mathematical Biosciences and Engineering, 2020, 17(1): 235-249. doi: 10.3934/mbe.2020013
    [9] Noura Laksaci, Ahmed Boudaoui, Seham Mahyoub Al-Mekhlafi, Abdon Atangana . Mathematical analysis and numerical simulation for fractal-fractional cancer model. Mathematical Biosciences and Engineering, 2023, 20(10): 18083-18103. doi: 10.3934/mbe.2023803
    [10] Salihu Sabiu Musa, Shi Zhao, Hei-Shen Chan, Zhen Jin, Daihai He . A mathematical model to study the 2014–2015 large-scale dengue epidemics in Kaohsiung and Tainan cities in Taiwan, China. Mathematical Biosciences and Engineering, 2019, 16(5): 3841-3863. doi: 10.3934/mbe.2019190
  • We study the long term dynamics and the multiscale aspects of a within-host HIV model that takes into account both mutation and treatment with enzyme inhibitors. This model generalizes a number of other models that have been extensively used to describe the HIV dynamics. Since the free virus dynamics occur on a much faster time-scale than cell dynamics, the model has two intrinsic time scales and should be viewed as a singularly perturbed system. Using Tikhonov's theorem we prove that the model can be approximated by a lower dimensional nonlinear model. Furthermore, we show that this reduced system is globally asymptotically stable by using Lyapunov's stability theory.


    Dengue fever is a viral disease produced by a pathogenic agent and is spread through female Aedes aegypti mosquitoes. Specific temperature areas of the world are affected by dengue infection and are almost spread to 128 countries which brings much more losses in public health and economy [1]. The symptoms include, muscle pain, nausea, red eyes, high fever, joints pain, pain behind eyes, lower back pain, vomiting, severe fatigue and weakness, and mild bleeding etc. The carrier vector becomes infected when they bit an infected host and transmit four different serotypes of dengue to the host through their bites. Rare cases of vertical transmission in dengue infection are reported [2,3]. Due to the increased prominence of dengue viral infection in the past few decades, much effort has been made recently to develop a dengue vaccine. Yet, there is no fully effective vaccine although the partially effective vaccine is accessible in several states of the world [4,5]. Reduction of vector bites and environmental domain for vectors are other adopted control policies for dengue infection.

    Epidemiological models are of great significance to gain insight into the comprehensive dynamics of infectious disease and to provide better control policies. To be more specific, analysis of these models predicts critical factors that play a central part in the prevention and spread of infection [6,7]. A number of the dynamical models have been established and approved to visualize the transmission of dengue infection [8,9,10,11]. The basic model of dengue infection was introduced by Lourdes Esteva and analyzed the stability of the model with variable human population [8,9]. Stability results for the equilibria of dengue infection model were established by in [10,11]. The role of vaccination on the transmission dynamics of dengue is examined in [12,13]. In dengue infection, asymptomatic carriers are frequently reported in [14,15] and are critical in the sense to increase the level of infection in endemic and non-endemic areas of the world.

    It is eminent that vector-borne disease posses information about its previous stages and associative learning mechanism, to be more specific memory has a prominent role in the transmission dynamics of vector-borne disease. The memory in the host population associated to the individual awareness which reduces the contact rate between vector and hosts, while the mosquitoes use their prior experience about the human's location, blood selection, color, and the smell of humans sweat [16,17]. These types of phenomena can easily be captured by a fractional-order system in mathematical modeling of infectious diseases.

    Fractional order system posses information about its previous and present stages to provide more realistic information about the dynamics of viral infections. Therefore, the dynamical behavior of epidemic diseases can be conceptualized and explore more accurately through non-integer order derivatives. A number of fractional-order models are introduced in the literature and successfully showed their positive role in different areas of physics, engineering, mathematics, and biology. For instant, a Maxwell model in fractional derivative is studied in [18]. Dengue dynamics in fractional derivative is considered in [19]. Application of fractional calculus to subdiffusion-reaction process is explored in [20]. An SIR model in fractional derivative is considered in [21]. The wave equation in a delay type fractional derivative is considered [22]. Recently, researchers proved that the fractional-order model can effectively handle the crossover behavior and complexity of communicable diseases, and offers a preferable fit for real data of the infection [23,24]. These properties make fractional-order systems more efficient than integer-order systems. For more related research on fractional order derivative and their applications to science and engineering problems, the authors investigated many piratical problems. For example, the author in [25] studied the fractional derivative and its applications to epidemic models. The corona virus model with lockdown is discussed through a mathematical model in [26]. The CoVID-19 and its impact on HIV and their future observation has been investigated in [27]. The corona virus model with fractional derivative is studied in [28], where the fractional operator is considered the Atangana-Baleanu type. The authors in [29] studded a corona virus model with nonlocal operator and using a power computational method. The M-fractional derivative and its applications to differential equations is discussed in [30]. An engineering problem with a new numerical results is studied in [31]. A fractional Schrodinger-Hirota equation and its optical solution is analyzed in [32]. The computer virus model in fractional derivative is considered in [33]. A review article that describing the harmonic wavelets and its fractal extension is studied in [34]. The solution of heat and diffusion equation with new technology is considered in [35]. A vector-host model with saturated treatment function and its optimal control analysis is discussed in [36]. The solution of Telegraph equation in new fractional operators is considered in [37]. The boundary value problems in fractional derivatives are studied in [38]. A fractional SEIR type model with treatment is studied in [39] while he dynamics of HIV with different transmission rates are analyzed in [40]. Therefore, motivated by these extraordinary properties, we analyze the dynamical behavior of dengue infection in the framework of the fractional-order derivative.

    The article is structured as: In the second section, we put forward a brief summary of the basic concepts of fractional calculus. In the third section, we construct a fractional-order model for dengue infection with new assumptions, in addition, we proved basic results for our proposed fractional-order model. We established stability results for the steady-state in section four and carried out sensitivity analysis through the PRCC technique to point out the importance of input parameter on the output of $ R_0 $. We have shown the influence of different input factors on the basic reproduction number $ R_0 $, numerically. The fractional-order system is inspected numerically to observe the influence of $ \vartheta $ on the dynamics of dengue in Section 5. In the last section of the article, ending remarks and conclusion of the overall analysis is presented.

    Here, we will present a brief summary of the Caputo fractional derivative for further analysis of the dengue model. Caputo fractional derivative has the benefit of dealing with initial value problems and involves integration which posses information about the previous state of the system.

    Definition 2.1. Let $ \mathrm{g}:\mathbb{R}^+\rightarrow \mathbb{R} $, then the fractional integral of $ \mathrm{g} $ is given by

    $ Iϑt(g(t))=1Γ(ϑ)t0(ty)ϑ1g(y)dy.
    $
    (2.1)

    where $ \vartheta $ indicates the order of fractional integral and Gamma function is indicated by $ \Gamma $.

    Definition 2.2. The fractional derivative of order $ \vartheta $ for a given function $ g $ in the Caputo form is given by

    $ CaDϑt(g(t))=InϑDng(t)=1Γ(nϑ)ta(ty)nϑ1gn(y)dy,
    $
    (2.2)

    where $ n-1 < \vartheta < n, n\in N $.

    Let $ \eta_1, \eta_2 $ be two positive numbers, then the Mittag-Leffler function is given by

    $ Eη1,η2(s)=k=k=0skΓ(η1k+η2).
    $
    (2.3)

    Let $ ^{C}_0 D^{\vartheta}_t $ be the Caputo fractional derivative (CFD) of order $ \vartheta $, then its laplace transform is define by

    $ L[C0Dϑtg(t)]=uϑG(u)n1k=0gk(0)uϑk1,
    $
    (2.4)

    further, for the function $ t^{\eta_2-1}E_{\eta_1, \eta_2}(\pm \omega t^{\vartheta}) $, the laplace transform is given by

    $ L[tη21Eη1,η2(±ωtϑ)]=uη1η2uη1±ω.
    $
    (2.5)

    The following equation hold true by Mittage-Leffler function given in [41]:

    $ E_{\eta_1,\eta_2}(s) = s . E_{\eta_1,\eta_1+\eta_2}(s)+\frac{1}{\Gamma(\eta_2)}. $

    Lemma 2.1. [42]. Consider a function $ g:R^+ \times R^5 \rightarrow R^5 $, which fulfills the conditions:

    The function $ g(t, Y(t)) $ is Lebesgue measurable in $ R^+ $ with respect to t;

    The function $ g(t, Y(t)) $ is continuous on $ R^5 $ with respect to $ Y(t) $;

    The function $ \frac{\partial g(t, Y(t))}{\partial Y} $ is continuous on $ R^5 $ with respect to $ Y(t) $;

    and $ \| g(t, Y(t)) \| \leq \aleph +\hbar \| Y \|, \forall t \in R^+, Y \in R^5, $ where $ \aleph, \hbar $ are positive constants.

    Then the fractional order system

    $ CaDϑtY(t)=g(t,Y(t)),Y(0)=Y0,
    $
    (2.6)

    where $ 0 < \vartheta\leq 1 $, has a unique solution.

    Lemma 2.2. [43]. Let a function $ g(t) \in C(r, s] $ such that $ ^{C}_r D^{\vartheta}_t g(t) \in C[r, s] $ for $ \vartheta \in (0, 1] $, then we have

    $ g(t) = g(r)+ \frac{1}{\Gamma(\vartheta)} ^{C}_r D^{\vartheta}_t g(\epsilon)(t-r)^{\vartheta}, \ \ r \lt \epsilon \lt t, \forall t \in (r,s]. $

    Remark 1. Let a function $ g \in C(r, s] $, such that $ ^{C}_r D^{\vartheta}_t g(t) \in C[r, s] $ for $ \vartheta \in (0, 1] $. Then by Lemma (2.2) if $ ^{C}_r D^{\vartheta}_t g(t) \leq 0 $, $ \forall \; t \in (r, s) $, then $ g(t) $ is non-increasing function for all $ t \in [r, s] $, and if $ ^{C}_r D^{\vartheta}_t g(t) \geq 0 $, $ \forall \; t \in (r, s) $, then $ g(t) $ is non-decreasing function for all $ t \in [r, s] $.

    In construction of the model, we represent the mosquitoes and humans population by $ N_v $ and $ N_h $. The total mosquitoes size is divided into susceptible ($ S_v $) and infected ($ I_v $) compartments, while the total humans size is divided into susceptible ($ S_h $), infected ($ I_h $), carrier ($ I_{hA} $), partially immune ($ P $), and recovered ($ R_h $) compartments.

    We indicate the natural birth (or death) rate of mosquito and human by $ \mu_v $ and $ \mu_h $, respectively, and assumed that the disease-induced death rate to be negligible. Furthermore, we assumed that a fraction $ \upsilon $ of susceptible hosts goes to partially immune class after implementation of vaccination and $ \tau $ indicates treatment of infected host. The symptomatic fraction of infected class from susceptible and partially immune classes are indicated by $ \psi $ and $ \omega $, respectively. the parameter $ \rho $ denotes the fraction of treated humans which join the partially immune class $ P $ while remanning enter to recovered class $ R_h $. Parameter $ \gamma_h $ is used for recovery of symptomatic and asymptomatic humans while $ \phi $ denotes the fraction of recovered humans which join the partially immune class $ P $. The terms ($ \frac{b \beta_1}{N_h} I_v $) and ($ \frac{b \beta_2 }{N_h} I_v $) indicate the infection rate per susceptible host and partially immune host, respectively, while ($ \frac{b \beta_3 }{N_h} I_h $) denotes the infection rate of per susceptible vector. In these terms $ b $ accounts the bitting rate, $ \beta_1 $ and $ \beta_2 $ are the rates at which infected vectors transmit disease to susceptible and partially immune humans respectively, whereas $ \beta_3 $ represents the disease transmission rate from infected human to susceptible vector. The dynamics of dengue infection is given by

    $ {dShdt=μhNhβ1bNhShIvυShμhSh,dIhdt=ψβ1bNhShIv+ωβ2bNhPIv(μh+τ+γh)Ih,dIhAdt=(1ψ)β1bNhShIv+(1ω)β2bNhPIv(μh+γh)IhA,dPdt=υSh+ρτIh+ϕγh(Ih+IhA)β2bNhPIvμhP,dRhdt=(1ρ)τIh+(1ϕ)γh(Ih+IhA)μhRh,dSvdt=μvNvβ3bNhSv(Ih+IhA)μvSv,dIvdt=β3bNhSv(Ih+IhA)μvIv,
    $
    (3.1)

    with the following initial condition

    $ S_h(0)\geq 0, I_h(0)\geq 0, I_{hA}(0)\geq 0, P(0)\geq 0, R_h(0)\geq 0, S_v(0)\geq 0, I_v(0) \geq 0. $

    We denote the total human population size by $ N_h $, and vector by $ N_v $ so that, $ N_h = S_h+I_h+I_{hA}+P+R_h $ and $ N_v = S_v+I_v $. As the system (3.1) is independent of the state-variable $ R_h $, thus, the proposed model (3.1) of dengue infection can be re-written in the following form without $ R_h $:

    $ {dShdt=μhNhβ1bNhShIvυShμhSh,dIhdt=ψβ1bNhShIv+ωβ2bNhPIv(μh+τ+γh)Ih,dIhAdt=(1ψ)β1bNhShIv+(1ω)β2bNhPIv(μh+γh)IhA,dPdt=υSh+ρτIh+ϕγh(Ih+IhA)β2bNhPIvμhP,dSvdt=μvNvβ3bNhSv(Ih+IhA)μvSv,dIvdt=β3bNhSv(Ih+IhA)μvIv,
    $
    (3.2)

    further, assume that

    $ x_1 = \frac{S_h}{N_h}, x_2 = \frac{I_h}{N_h}, x_3 = \frac{I_{hA}}{N_h}, x_4 = \frac{P}{N_h}, x_5 = \frac{I_v}{N_v} \ and \ x_6 = \frac{S_v}{N_v}, $

    then, we have $ x_1+x_2+x_3+x_4 = 1 $ and $ x_5+x_6 = 1 $. After replacing $ x_6 = 1-x_5 $ in the last equation of normalized system we obtained the following five dimensional reduced dengue model:

    $ {dx1dt=μhbβ1NvNhx1x5υx1μhx1,dx2dt=ψbβ1NvNhx1x5+ωbβ2NvNhx4x5(μh+τ+γh)x2,dx3dt=(1ψ)bβ1NvNhx1x5+(1ω)bβ2NvNhx4x5(μh+γh)x3,dx4dt=υx1+ρτx2+ϕγh(x2+x3)bβ2NvNhx4x5μhx4,dx5dt=bβ3(1x5)(x2+x3)μvx5.
    $
    (3.3)

    As fractional-order models describe the non-local behavior of biological systems and posses hereditary property, moreover, it provides information about its past and present state for the future, therefore, we represent the dynamical system (3.3) of dengue infection in the framework of fractional order Caputo's derivative to conceptualize the transmission of dengue fever in a more accurate way. Thus, the system consist of fractional derivatives is presented by

    $ {C0Dϑtx1=μϑhbϑax1x5υϑx1μϑhx1,C0Dϑtx2=ψbϑax1x5+ωbϑcx4x5(μϑh+τϑ+γϑh)x2,C0Dϑtx3=(1ψ)bϑax1x5+(1ω)bϑcx4x5(μϑh+γϑh)x3,C0Dϑtx4=υϑx1+ρτϑx2+ϕγϑh(x2+x3)bϑcx4x5μϑhx4,C0Dϑtx5=bϑe(1x5)(x2+x3)μϑvx5,
    $
    (3.4)

    where $ ^{C}_0 D_t ^{\vartheta} $ indicates Caputo's fractional derivative of order $ \vartheta $, the order $ \vartheta $ indicates the index of memory in the system. In addition, the values

    $ a = \frac{{ \beta_1 N_v}}{N_h} , c = \frac{{ \beta_2 N_v}}{N_h}, \ \ and \ \ e = \beta_3 . $

    In Caputo's sense, the derivative of constant is equal to zero which is another advantage to make the system more reliable and flexible for analysis. Next, we will analyze the biologically feasible region of the fractional-order dengue model (3.4).

    Theorem 3.1. The proposed fractional-order dengue infection system of equations (3.4) has a unique solution.

    Proof. For the required result, we first prove that the fractional system (3.4) has a unique solution for all initial conditions in $ R^5 $. Clearly, the first three conditions of Lemma (2.1) are hold by the vector function $ g $ of the system (3.4). Next, to prove the last condition of Lemma (2.1), we rewrite system (3.4) as

    $ ^{C}_0 D^{\vartheta}_t x(t) = \aleph+ A_1 x(t)+A_2 x_1(t) x(t)+A_3 x_2(t) x(t)+A_4 x_3(t) x(t)+A_5 x_4(t) x(t)+A_6 x_5(t) x(t), $

    where

    $ \aleph = \left[{μϑh0000
    }\right], A_1 = \left[{(μϑh+υϑ)00000(μϑh+τϑ+γϑh)00000(μϑh+γϑh)00υϑρτϑ+ϕγϑhϕγϑhμϑh00bϑebϑe0μϑv
    }\right] , $
    $ A_2 = \left[{000000000ψbϑa0000(1ψ)bϑa0000000000
    }\right] , A_3 = \left[{000000000000000000000000bϑe
    }\right], A_4 = \left[{000000000000000000000000bϑe
    }\right] , $
    $ A_5 = \left[{000000000ωbϑc0000(1ω)bϑc0000000000
    }\right], A_6 = \left[{bϑa00000000000000000bϑc000000
    }\right]. $

    Symbolize the above

    $ g(t,x(t)) = \aleph+ A_1 x(t)+A_2 x_1(t) x(t)+A_3 x_2(t) x(t)+A_4 x_3(t) x(t)+A_5 x_4(t) x(t)+A_6 x_5(t) x(t), $

    and taking norm, we get

    $ g(t,x(t))=+A1x(t)+A2x1(t)x(t)+A3x2(t)x(t)+A4x3(t)x(t)+A5x4(t)x(t)+A6x5(t)x(t)+A1x(t)+A2x1(t)x(t)+A3x2(t)x(t)+A4x3(t)x(t)+A5x4(t)x(t)+A6x5(t)x(t)+A1x(t)+A2x(t)+A3x(t)+A4x(t)+A5x(t)+A6x(t)=+(A1+A2+A3+A4+A5+A6)x(t)=+x(t).
    $

    Hence, the requirements of Lemma (2.1) are fulfilled, therefore the system of equations (3.4) has a unique solution.

    Theorem 3.2. The closed set $ \Omega = \{ (x_1, x_2, x_3, x_4, x_5) \in R^5_+: 0 \leq x_1+x_2+x_3+x_4 \leq M_1, 0 \leq x_5\leq M_2 \} $ is a positive invariant set for the proposed fractional order system (3.4).

    Proof. To prove that the system of equations (3.4) has a non-negative solution, the system of equations (3.4) implies

    $ {C0Dϑtx1x1=0=μϑh>0,C0Dϑtx2x2=0=ψbϑax1x5+ωbϑcx4x50,C0Dϑtx3x3=0=(1ψ)bϑax1x5+(1ω)bϑcx4x50,C0Dϑtx4x4=0=υϑx1+ρτϑx2+ϕγϑh(x2+x3)0,C0Dϑtx5x5=0=bϑe(x2+x3)0.
    $
    (3.5)

    Thus, the fractional system (3.4) has non-negative solutions. In the end, from the first four equations of the fractional system (3.4), we obtain

    $ C0Dϑt(x1+x2+x3+x4)μϑhμϑhx1((1ρ)τϑ+(1ϕ)γϑh+μϑh)x2((1ϕ)γϑh+μϑh)x3μϑhx4,μϑhW(x1+x2+x3+x4)
    $
    (3.6)

    where $ \mathcal{W} = min(\mu_h^\vartheta, ((1-\rho)\tau^\vartheta+(1-\phi)\gamma_h^\vartheta+\mu_h^\vartheta), ((1-\phi)\gamma_h^\vartheta+\mu_h^\vartheta), \mu_h^\vartheta) $. Solving the above inequality, we obtain

    $ \bigg(x_1(t)+x_2(t)+x_3(t)+x_4(t) \bigg) \leq \bigg(x_1(0)+x_2(0)+x_3(0)+x_4(0)-\frac{\mu_h^\vartheta}{\mathcal{W}} \bigg) E_{\vartheta}(-\mathcal{W} t^{\vartheta})+\frac{\mu_h^\vartheta}{\mathcal{W}}, $

    so by the asymptotic behavior of Mittag-Leffler function [44], we obtain

    $ \bigg(x_1(t)+x_2(t)+x_3(t)+x_4(t) \bigg) \leq \frac{\mu_h^\vartheta}{\mathcal{W}}\cong M_1, $

    taking the same steps for the last equation of system (3.4), we get $ x_5(t) \leq M_2 $, where $ M_2 = \frac{b^\vartheta e M_1}{b^\vartheta e M_1+\mu_v^\vartheta} $. Hence, the closed set $ \Omega $ is a positive invariant region for the fractional-order dengue model (3.4).

    The biological meaningful equilibria of fractional system (3.4) are DFE and EE, depending on infected classes in both the populations. To obtain the infection-free equilibrium, we set the fractional derivative $ ^{C}_0 D^{\vartheta}_t x_1, ^{C}_0 D^{\vartheta}_t x_2, ^{C}_0 D^{\vartheta}_t x_3, ^{C}_0 D^{\vartheta}_t x_4, $ and $ ^{C}_0 D^{\vartheta}_t x_5 $ to zero of the fractional system (3.4) without infection, and get

    $ E_0(x_1^0,x_2^0,x_3^0,x_4^0,x^0_5) = \left(\frac{\mu_h^\vartheta}{\mu_h^\vartheta+\upsilon^\vartheta},0,0,\frac{\upsilon^\vartheta}{\mu_h^\vartheta+\upsilon^\vartheta},0\right). $

    The detailed concept of the basic reproduction number is presented in [45], and is normally denoted by $ R_0 $, and is calculated as

    $ \mathcal{F} = \left[{ψbϑax1x5+ωbϑcx4x5(1ψ)bϑax1x5+(1ω)bϑcx4x5bϑe(1x5)(x2+x3)
    }\right] and \ \ \mathcal{V} = \left[{(μϑh+τϑ+γϑh)x2(μϑh+γϑh)x3μϑvx5
    }\right] , $

    this is because of three infected comportment in the system, i.e., $ m = 3 $, which provides

    $ F = \left[{00ψbϑax01+ωbϑcx0400(1ψ)bϑax01+(1ω)bϑcx04bϑebϑe0
    }\right] and \ \ V = \left[{(μϑh+τϑ+γϑh)000(μϑh+γϑh)000μϑv
    }\right] , $

    which gives

    $ FV^{-1} = \left[{00ψbϑax01+ωbϑcx04μϑv00(1ψ)bϑax01+(1ω)bϑcx04μϑvbϑe(μϑh+τϑ+γϑh)bϑe(μϑh+γϑh)0
    }\right] . $

    The $ R_0 $ of the fractional-order dengue model (3.4) is obtained by using the next generation matrix as

    $ \rho(FV^{-1}) = \sqrt{ \frac{b^\vartheta e}{\mu_v^\vartheta(\mu_h^\vartheta+\upsilon^\vartheta)} \bigg[ \frac{(1-\psi) b^\vartheta a \mu_h^\vartheta + (1-\omega) b^\vartheta c \upsilon^\vartheta}{\mu_h^\vartheta+\gamma_h^\vartheta}+ \frac{\psi b^\vartheta a \mu_h^\vartheta+ \omega b^\vartheta c \upsilon^\vartheta}{\mu_h^\vartheta+\gamma_h^\vartheta+\tau^\vartheta}\bigg]}, $
    $ \Rightarrow R_0 = \sqrt{ \frac{b^\vartheta e}{\mu_v^\vartheta(\mu_h^\vartheta+\upsilon^\vartheta)} \bigg[ \frac{(1-\psi) b^\vartheta a \mu_h^\vartheta + (1-\omega) b^\vartheta c \upsilon^\vartheta}{\mu_h^\vartheta+\gamma_h^\vartheta}+ \frac{\psi b^\vartheta a \mu_h^\vartheta+ \omega b^\vartheta c \upsilon^\vartheta}{\mu_h^\vartheta+\gamma_h^\vartheta+\tau^\vartheta}\bigg]}. $

    Assume a fractional-order linear homogenous system of the following form

    $ C0Dϑty(t)=By(t),y(0)=y0,
    $
    (4.1)

    where $ B\in M_{m\times m}(R) $ and $ 0 < \vartheta \leq 1 $. The following theorems are on the stability of linear homogenous system (4.1).

    Theorem 4.1. [46]. The origin of the fractional dynamical system (4.1) is asymptotically stable $ \Leftrightarrow \mid arg(\lambda_i) \mid > \frac{\vartheta \pi}{2} $ is fulfilled for all eigenvalues $ \lambda_i $ of matrix B.

    Theorem 4.2. [47]. The steady state of fractional system (2.6) is locally asymptotically stable (LAS) if $ \mid arg(\lambda_i) \mid > \frac{\vartheta \pi}{2} $ for all eigenvalues $ \lambda_i $'s of $ g(y) $ at steady state, otherwise unstable.

    In order to prove the $ LAS $ of infection-free steady state $ E_0 $ it is sufficient to confirm that all of the eigenvalues of the matrix of dynamics given by $ J(E_0) $ lie outside the closed angular sector as shown in [46]. We state the following theorem for the desired result.

    Theorem 4.3. The $ DFE $ of the fractional-order system (3.4) is LAS if $ R_0 < 1 $ and $ \mid arg(\lambda_i) \mid > \frac{\vartheta \pi}{2} $ for all eigenvalues $ \lambda_i $'s of the Jacobian matrix $ J(E_0) $ at steady state, otherwise it is unstable.

    Proof. To obtain the demanded result for the fractional system (3.4), we take the Jacobian matrix of the system at infection-free equilibrium as

    $ J(E_0) = \left[{(υϑ+μϑh)000bϑax010(μϑh+γϑh+τϑ)00ψbϑax01+ωbϑcx0400(μϑh+γϑh)0(1ψ)bϑax01+(1ω)bϑcx04υϑρτϑ+ϕγϑhϕγϑhμϑhbϑcx040bϑebϑe0μϑv
    }\right] . $

    The characteristic equation in term of $ \Lambda $ of the above matrix $ J_{\mathcal{E}_0} $ which is given below:

    $ (λ+μϑh)(λ+(μϑh+υϑ))(λ3+a1λ2+a2λ+a3)=0,
    $
    (4.2)

    where the associated coefficients are as follow:

    $ a1=(μϑh+γϑh+τϑ)+(μϑh+γϑh)+μϑv,a2=(μϑh+γϑh+τϑ)(μϑh+γϑh)+(μϑh+γϑh+τϑ)μϑv+(μϑh+γϑh)μϑvbϑeBbϑeA,a3=(μϑh+γϑh+τϑ)(μϑh+γϑh)μϑv[1R20],
    $

    , where $ A = \psi b^\vartheta a x_1^0 + \omega b^\vartheta c x_4^0 $ and $ B = (1-\psi) b^\vartheta a x_1^0 + (1-\omega)b^\vartheta c x_4^0 $. From (4.2) it is clear the argument of the first two eigenvalues $ \lambda = -\mu_h^\vartheta $ and $ \lambda = -(\mu_h^\vartheta+\upsilon^\vartheta) $ satisfy the necessary condition given (i.e., $ \mid arg(\lambda_i) \mid > \frac{\vartheta \pi}{2} $ for all $ \vartheta\in(0, 1) $). Furthermore, if $ R_0 < 1, $ then all $ a_i > 0 $, for $ i = 1, 2, 3, $ and it is easy to show that $ a_1 a_2-a_3 > 0 $ Thus, following [46], the arguments of all eigenvalues satisfy the necessary Matignon condition. Consequently, the infection-free steady-state is $ LAS $.

    Theorem 4.4. The infection-free steady state of the fractional-order system (3.4) is GAS without vaccination, if $ R_0 < 1 $.

    Proof. Let $ \upsilon = 0 $, and $ (x_1, x_2, x_3, x_4, x_5) $ be the solution of system (3.4) with suitable initial conditions $ (x_1(0), x_2(0), x_3(0), x_4(0), x_5(0)) $ in $ \Omega $. Here, it is clear that the fractional system has only one equilibrium $ E_0 $ on the boundary of $ \Omega $. Therefore, to achieve the target, it is enough to prove that the solution $ (x_1, x_2, x_3, x_4, x_5) $ tends to the infection-free equilibrium as time tends to infinity. Since $ x_i \leq x_i^0 $, for $ i = 1, 4 $, then the system implies that

    $ C0Dϑtx2ψbϑax01x5+ωbϑcx04x5(μϑh+τϑ+γϑh)x2,C0Dϑtx3(1ψ)bϑax01x5+(1ω)bϑcx04x5(μϑh+γϑh)x3,C0Dϑtx5bϑe(x2+x3)μϑvx5.
    $

    Taking the auxiliary system

    $ C0Dϑty1=ψbϑax01y3+ωbϑcx04y3(μϑh+τϑ+γϑh)y1,C0Dϑty2=(1ψ)bϑax01y3+(1ω)bϑcx04y3(μϑh+γϑh)y2,C0Dϑty3=bϑe(y1+y2)μϑvy3,
    $

    this further implies

    $ C0DϑtY=(FV)Y
    $

    the coefficient matrix of the above fractional system is $ F-V $, and if $ R_0 = \rho(FV^{-1}) < 1 $, then the eigenvalues of $ F-V $ lies in the left half plane. Consequently, each positive solution of the fractional system (4.3) fulfills $ \lim \limits_{t\rightarrow\infty}y_1 = 0 $, $ \lim \limits_{t\rightarrow\infty}y_2 = 0 $ and $ \lim \limits_{t\rightarrow\infty}y_3 = 0 $ by Theorem 4.1. By the comparison theory of fractional differential equations [48], we have $ \lim \limits_{t\rightarrow\infty}x_2 = 0 $, $ \lim \limits_{t\rightarrow\infty}x_3 = 0 $ and $ \lim \limits_{t\rightarrow\infty}x_5 = 0 $. Then from system (3.4), we have

    $ C0Dϑtx1=μϑhμϑhx1,C0Dϑtx4=0μϑhx4,
    $

    this can be further converted into

    $ ^{C}_0 D^{\vartheta}_t X = \Lambda - B X, $

    having the solution

    $ X(t)=tϑEϑ,ϑ+1(Btϑ)Λ+Eϑ,1(Btϑ)X0,
    $

    taking the asymptotic behavior of Mittag-Leffler function [44] as

    $ E_{\vartheta,1}(-\omega t^{\vartheta})\approx _{t\rightarrow\infty} \frac{t^\vartheta}{\omega \Gamma(1-\vartheta)}; \ \ \ 0 \lt \vartheta \lt 1 \ \ and \ \ \omega \gt 0, $

    also the real parts of eigenvalues of $ -B $ are negative, so it can be observed that $ X(t) \rightarrow X^0 $ as t tends to infinity. Hence, the infection-free steady-state of the system (3.4) is GAS without vaccination.

    Next, we will demonstrate the persistence of infection in the fractional-order system. It describes the level of endemicity of infection in the system. Biologically speaking, the infection persists in the system if the level of infected fraction stays at a higher level for t large enough.

    Theorem 4.5. In the fractional-order system (3.4), the dengue infection is uniformly persistent at $ \vartheta = 1 $ if $ R_0 > 1 $.

    Proof. The threshold parameter $ R_0 $ is of high significance, which measures the level of infection in the system. If $ R_0 $ less than one then the asymptomatic stability of infection-free equilibrium$ (E_0) $ prevent the system from the persistent of the infection. If $ R_0 $ greater than one then the sufficiently close solution of fractional system (3.4) to $ E_0 $ moves away from it due to its instability in the interior of $ \Omega $ by Theorem (4.3). For uniform persistence of our system (3.4) at $ \vartheta = 1 $, we will prove that it fulfills the criteria of Theorem 4.3 in [49] for $ R_0 > 1 $. Select $ X = R^5 $ and $ E = \Omega $. Here, the maximal invariant set on the boundary of $ \Omega $ is the set $ \{E_0\} $ and is isolated. So, the criteria of Theorem 4.3 in [49] satisfies for the system. Thus, the system (3.4) is uniformly persistent if $ E_0 $ is unstable for $ R_0 > 1 $.

    In the section, we present the global sensitivity analysis in order to determine those model parameters which are most influential on the disease dynamics. The main goal of this analysis is to detect and measure the influence of input parameters on the output of the system. To be more specific, it is used to know how the input parameters and initial values contribute to the output of a system. Mostly, when there is a little uncertainty in initial conditions and input parameters partial derivative of output functions are computed with respect to the input factors around the base values. This method is named as the local sensitivity analysis and relies on the variations of parameters close to the base values. This technique is not most suitable for epidemiological models due to the uncertainty in the input of the system. Therefore, global sensitivity analysis is preferred to perform this analysis and to provide more accurate results.

    Here, we used the PRCC method [50] for sensitivity analysis to point out the input parameters that highly influence the results of $ R_0 $. It is an effective method and can successfully measure the monotonic, nonlinear relationship between input and output values of the system. PRCC analysis provides PRCC and p values for each factor involve therein, with which we can measure the contribution of each factor. More specifically, the input factors with sizeable PRCC and negligible p-values are considered to be highly effective factors in the system. In our analysis, we investigate all the parameter presented in Table 1 to know their contribution to the outcomes of $ R_0 $ and listed all the associated PRCC and p-values provided by the PRCC significance test. Figure 1 and Table 1 demonstrate that the parameter $ b $ is highly influential with PRCC value 0.7022. After that, the parameter $ \beta_1 $, and $ \vartheta $ are highly influential with PRCC values 0.6427, and -0.5112, respectively. Sensitivity analysis illustrated that $ b, \beta_1 $, and $ \vartheta $ are critical parameters in $ R_0 $, which can remarkably influence the infection. Therefore by controlling these factors can greatly decrease and prevent the level of new dengue cases in the community.

    Table 1.  Sensitivity results of $ R_0 $ with PRCC and corresponding p values.
    Parameter Interpretation PRCC values p values
    $ \vartheta $ index of memory or fractional order -0.5112 0.0000
    $ \mu_h $ recruitment rate of humans host 0.1095 0.0023
    $ \upsilon $ vaccinated fraction of susceptible host + 0.3133 0.0000
    $ \psi $ fraction of incidence rate from $ S_h $ to $ I_h $ + 0.1804 0.0010
    $ \beta_{1} $ transmission probability from vectors to susceptible host + 0.6427 0.0000
    $ \beta_{2} $ transmission probability from vectors to partially immune host + 0.5094 0.0000
    $ \tau $ treatment fraction of infected host individuals -0.2110 0.0000
    $ \gamma_h $ recovery rate of host individuals -0.3361 0.0000
    $ \beta_3 $ transmission probability from infected hosts to susceptible vectors +0.2320 0.0000
    $ \mu_v $ recruitment rate of mosquitoes vector -0.3197 0.0000
    $ \omega $ fraction of incidence rate from $ P $ to $ I_h $ + 0.3596 0.0000
    $ b $ bitting rate of vectors +0.7022 0.0000

     | Show Table
    DownLoad: CSV
    Figure 1.  Sensitivity test for input factors of $ R_0 $ and PRCC outcomes for $ R_0 $.

    Next, we illustrate the influence of vaccination and treatment in on the threshold $ R_0 $ of the proposed system 2 (a). We noted that vaccination and treatment decrease the threshold parameter, moreover, we observed that vaccination reduces the susceptibility of host individuals in the system while the treatment reduces the infectivity of host individuals. In Figure 2 (b), we showed the influence of $ \beta_3 $ and $ \gamma_h $ on the basic reproduction number $ R_0 $ of the proposed system of dengue fever.

    Figure 2.  Illustration of basic reproduction number $ R_0 $ (a) with the variation of vaccination $ \upsilon $ and treatment $ \tau $, (b) and with the variation of $ \beta_3 $ and $ \gamma_h $.

    We present in this section the solution of the model (3.4) by using the algorithm presented in [51]. To have a numerical scheme, we write the model (3.4) in the following form:

    $ {CDϑtg(t)=G(t,g(t)),g(0)=g0,0<T<,
    $
    (5.1)

    where $ g = (x_1, x_2, x_3, x_4, x_5)\in \mathbb{R_+}^5 $, $ \mathcal{G}\Big(t, g(t)\Big) $ is used for a continuous real valued vector function, which additionally satisfies the Lipschitz condition and $ g_0 $ stands for initial state vector. Taking Caputo integral on both sides of (5.1) we get

    $ g(t)=g0+1Γ(ϑ)t0(tλ)ϑ1G(λ,g(λ))dλ.
    $
    (5.2)

    In order to formulate an iterative scheme, we consider a uniform grid on $ [0, \mathcal{T}] $ with $ h = \frac{\mathcal{T}-0}{m} $ is the step size and $ m\in \mathbb{N} $. Thus, the equation (5.2) gets the structure as follows after make use of the Euler method [52]

    $ {gn+1=g0+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)G(tj,g(tj)),n=0,1,2,,m.
    $
    (5.3)

    Thus, utilizing the above scheme (5.3), we deduced the following iterative formulae for the corresponding classes of the model (3.4)

    $ x1n+1=x10+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)(μϑhbϑax1jx5jυϑx1jμϑhx1j),x2n+1=x20+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)(ψbϑax1jx5j+ωbϑcx4jx5j(μϑh+τϑ+γϑh)x2j),x3n+1=x30+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)((1ψ)bϑax1jx5j+(1ω)bϑcx4jx5j(μϑh+γϑh)x3j),x4n+1=x40+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)(υϑx1j+ρτϑx2j+ϕγϑh(x2j+x3j)bϑcx4jx5jμϑhx4j),x5n+1=x50+hϑΓ(ϑ+1)nj=0((nj+1)ϑ(nj)ϑ)(bϑe(1x5j)(x2j+x3j)μϑvx5j).
    $
    (5.4)
    Table 2.  Interpretation of input parameters with values used in numerical simulations.
    Parameter Interpretation Values Reference
    $ \mu_v $ natural mortality and recruitment rate of vector 0.032300 & 0.029410 [53,54]
    $ \beta_{1} $ transmission probability from vectors to susceptible host 0.75 [10]
    $ \beta_{2} $ transmission probability from vectors to partially immune host 0.375 Assumed
    $ \vartheta $ index of memory or fractional order variable Assumed
    $ b $ bitting rate of vectors 0.50 [10]
    $ \tau $ treatment fraction of infected host 0.40 Assumed
    $ \upsilon $ vaccinated fraction of susceptible host 0.20 Assumed
    $ \beta_{3} $ transmission probability from infected humans to mosquitoes 0.75 [10]
    $ \psi $ fraction of incidence rate from $ S_h $ to $ I_h $ 0.48 Assumed
    $ \gamma_h $ recovery rate of host individuals 0.3288330 [53,54]
    $ \omega $ fraction of incidence rate from $ P $ to $ I_h $ 0.54 Assumed
    $ \mu_h $ natural mortality and recruitment rate of humans 0.004500 & 0.000046 [53,54]
    $ \rho $ fraction of treated humans enter to $ P $ class 0.6 Assumed

     | Show Table
    DownLoad: CSV

    We used the above approximation for the solution of our fractional system. In Figure (3), we demonstrate the dynamics of infected host ($ I_h $), asymptomatic host ($ I_{hA} $), partially immune $ (P) $ host and infected vector ($ I_v $) with the variation of fractional-order $ \vartheta $. We noticed that the variation of fractional-order $ \vartheta $ has a great influence on the infection level of dengue in both the population. In other words, it can highly reduce the level of dengue fever in the community. In Figure (4) and Figure (5), we represent the dynamics of dengue fever with the variation of vaccination and treatment. In these simulations, we observed that vaccination and treatment decrease susceptibility and infectivity of the host population, and have a slight influence on the level of infection. We demonstrated the effect of biting rate $ b $ of the mosquitoes on the dynamics of dengue and observed that the peak of infection can be greatly decreased by decreasing the biting rate $ b $ in Figure (6). The mosquitoes biting and its generation further can be decreased by spraying or wasting the standing water around home or inside the home, which has great influence on population of mosquitoes that end up biting humans. Using bed-nets, avoid to visit areas prone to mosquitoes, using mosquito repellent, covering legs and arms by wearing long-sleeves and long pants are useful to prevent himself from the biting of mosquitoes. These scenarios predict that the infection can be controlled and prevented by decreasing the index of memory and biting rate of vectors in the community.

    Figure 3.  Illustration of the path moment of (a) Infected host $ I_h $, (b) Carrier host(Asymptomatic) $ I_{hA} $, (c) Partially immune $ P $ and (d) Infected host $ I_v $ individuals with different values of fractional-order $ \vartheta $, i.e., $ \vartheta = 0.6, 0.8, 1.0 $.
    Figure 4.  Illustration of the path moment of (a) Infected host $ I_h $, (b) Carrier host(Asymptomatic) $ I_{hA} $, (c) Partially immune $ P $ and (d) Infected vector $ I_v $ individuals with different vaccination rates $ \upsilon $, i.e., $ \upsilon = 0.2, 0.4, 0.6 $.
    Figure 5.  Illustration of the path moment of (a) Infected host $ I_h $, (b) Carrier host(Asymptomatic) $ I_{hA} $, (c) Partially immune $ P $ and (d) Infected vector $ I_v $ individuals with different treatment rates $ \tau $, i.e., $ \tau = 0.4, 0.6, 0.8 $.
    Figure 6.  Illustration of the path moment of (a) Infected host $ I_h $, (b) Carrier host(Asymptomatic) $ I_{hA} $, (c) Partially immune $ P $ and (d) Infected vector $ I_v $ individuals with variation of biting rate of mosquitoes $ b $, i.e., $ b = 0.4, 0.6, 0.8 $.

    In this article, we constructed a new epidemic model for the transmission of dengue infection with non-integer derivative. Basic results of the suggested fractional system are investigated through analytic skills. The basic reproduction number of the system is calculated through the next-generation method, indicated by $ R_0 $. We established stability results of infection-free steady-state of the system. Sensitivity analysis of $ R_0 $ is carried out in order to know the contribution of input factors in the results of $ R_0 $ and observed that $ b, \mu_v $, and $ \beta_3 $ are the most critical parameters that highly contribute in the control and subsequent spread of dengue infection. We showed that the dengue infection is uniformly persistent in the system for $ R_0 > 1 $. We obtained feasible results for the dynamics of dengue infection with the variation of memory index $ \vartheta $ and suggested that the index of memory has a dominant influence on the system. We concluded that the fractional-order model can explore more accurately the dengue epidemic disease transmission model rather than the integer-order derivative models. The numerical results for the suggested parameter can be considered useful for the possible eliminations of dengue infection in the community. We suggest that fractional-order(index of memory) $ \vartheta $ and biting rate $ b $ can remarkably control and greatly decrease the level of disease in the society.

    No conflict of interest exists regarding the publications of this work.

    The authors are thankful to the anonymous reviewers for the careful checking and suggestions that improved the presentation of the paper greatly. Fatmawati and Muhammad Altaf Khan are supported by Universitas Airlangga through Hibah Riset Mandat (No. 342/UN3.14/LT/2019). The funding of this paper was made from the project Universitas Airlangga through Hibah Riset Mandat (No. 342/UN3.14/LT/2019). Yasir Khan extends their appreciation to the Deanship of Scientific Research, University of Hafr Al Batin for funding this work through the research group project no.(G-108-2020).

    [1] Nature, 333 (1988), 514-519.
    [2] Proceedings of the Royal Society of London. Series B: Biological Sciences, 270 (2003), 1651-1657.
    [3] Mathematical Biosciences and Engineering, 10 (2013), 499-521.
    [4] (2010).
    [5] Proceedings of the National Academy of Sciences, 94 (1997), 6971-6976.
    [6] Proceedings of the National Academy of Sciences, 83 (1986), 9159-9163.
    [7] Ministério da Saúde, 2013.
    [8] Journal of Differential Equations, 31 (1979), 53-98.
    [9] JAIDS Journal of Acquired Immune Deficiency Syndromes, 7 (1994), 236-244.
    [10] PLoS computational biology, 9 (2013), e1002971, 12 pp.
    [11] Current HIV research, 2 (2004), 23-37.
    [12] Nature, 373 (1995), 123-126.
    [13] Cambridge University Press, 1998.
    [14] vol. 114, Springer New York, 1996.
    [15] AMS notices, 43 (1996), 191-202.
    [16] Bulletin of Mathematical Biology, 66 (2004), 879-883.
    [17] Proceedings of the Royal Society of London. Series B: Biological Sciences, 268 (2001), 847-854.
    [18] Journal of Clinical Investigation, 105 (2000), R1-R8.
    [19] Immunology and cell biology, 85 (2006), 6-15.
    [20] Mathematical Biosciences, 179 (2002), 73-94.
    [21] Oxford University Press, 2000.
    [22] Science, 272 (1996), 74-79.
    [23] (2012).
    [24] Mac Graw Hill, 1978.
    [25] Ph.D. thesis, IMPA, 2005.
    [26] Mathematical biosciences, 114 (1993), 81-125.
    [27] SIAM review, 41 (1999), 3-44.
    [28] The Hopkins HIV Report, 8 (1996), no. 3.
    [29] Nature Reviews Genetics, 5 (2004), 52-61.
    [30] Journal of molecular evolution, 40 (1995), 249-259.
    [31] no. Ed. 4, WB Saunders, 1995.
    [32] Ph.D. thesis, AIMS, 2012.
    [33] Nature Reviews Microbiology, 1 (2003), 181-190.
    [34] SIAM Journal on Applied Mathematics, 63 (2003), 1313-1327.
    [35] Science, 242 (1988), 1554-1557.
    [36] Journal of mathematical biology, 68 (2014), 1269-1293.
    [37] Bull. Math. Biol., 73 (2011), 609-625.
    [38] Journal of Theoretical Biology, 203 (2000), 285-301.
    [39] Journal of differential equations, 92 (1991), 252-281.
    [40] Springer-Verlag Berlin, 1984.
    [41] Moskva: Vysshaya Shkola, 1990 (Russian).
    [42] Advances in Complex Systems, 10 (2007), 495-503.
    [43] Dover Publications, Inc., New York, 1987.
    [44] (2013).
  • This article has been cited by:

    1. Zahir Shah, Rashid Jan, Poom Kumam, Wejdan Deebani, Meshal Shutaywi, Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative, 2021, 26, 1420-3049, 1806, 10.3390/molecules26061806
    2. Rashid Jan, Salah Boulaaras, Sultan Alyobi, Muhammad Jawad, Transmission dynamics of Hand–Foot–Mouth Disease with partial immunity through non-integer derivative, 2023, 16, 1793-5245, 10.1142/S1793524522501157
    3. Olumuyiwa James Peter, Festus Abiodun Oguntolu, Mayowa M Ojo, Abdulmumin Olayinka Oyeniyi, Rashid Jan, Ilyas Khan, Fractional order mathematical model of monkeypox transmission dynamics, 2022, 97, 0031-8949, 084005, 10.1088/1402-4896/ac7ebc
    4. Tahajuddin Sk, Santosh Biswas, Tridip Sardar, The impact of a power law-induced memory effect on the SARS-CoV-2 transmission, 2022, 165, 09600779, 112790, 10.1016/j.chaos.2022.112790
    5. Ajeet M. Wagle, Smita R. Hegde, Srinivasan Sanjay, Kah-Guan Au Eong, Ophthalmic manifestations in seropositive dengue fever patients during epidemics caused by predominantly different dengue serotypes, 2022, 2, 26673762, 100049, 10.1016/j.aopr.2022.100049
    6. Din Prathumwan, Inthira Chaiya, Kamonchat Trachoo, Study of Transmission Dynamics of Streptococcus suis Infection Mathematical Model between Pig and Human under ABC Fractional Order Derivative, 2022, 14, 2073-8994, 2112, 10.3390/sym14102112
    7. Suthep Suantai, Zulqurnain Sabir, Muhammad Umar, Watcharaporn Cholamjiak, Scaled Conjugate Gradient for the Numerical Simulations of the Mathematical Model-Based Monkeypox Transmission, 2023, 7, 2504-3110, 63, 10.3390/fractalfract7010063
    8. Salah Boulaaras, Ziad Ur Rehman, Farah Aini Abdullah, Rashid Jan, Mohamed Abdalla, Asif Jan, Coronavirus dynamics, infections and preventive interventions using fractional-calculus analysis, 2023, 8, 2473-6988, 8680, 10.3934/math.2023436
    9. H. M. Srivastava, Rashid Jan, Asif Jan, Wejdan Deebani, Meshal Shutaywi, Fractional-calculus analysis of the transmission dynamics of the dengue infection, 2021, 31, 1054-1500, 053130, 10.1063/5.0050452
    10. Muhammad Sinan, Hijaz Ahmad, Zubair Ahmad, Jamel Baili, Saqib Murtaza, M.A. Aiyashi, Thongchai Botmart, Fractional mathematical modeling of malaria disease with treatment & insecticides, 2022, 34, 22113797, 105220, 10.1016/j.rinp.2022.105220
    11. Rashid Jan, Salah Boulaaras, Syed Azhar Ali Shah, Fractional-calculus analysis of human immunodeficiency virus and CD4+ T-cells with control interventions, 2022, 74, 0253-6102, 105001, 10.1088/1572-9494/ac7e2b
    12. Tao-Qian Tang, Zahir Shah, Rashid Jan, Wejdan Deebani, Meshal Shutaywi, A robust study to conceptualize the interactions of CD4+ T-cells and human immunodeficiency virus via fractional-calculus, 2021, 96, 0031-8949, 125231, 10.1088/1402-4896/ac2d7b
    13. Joshua Kiddy K. Asamoah, Ernest Yankson, Eric Okyere, Gui-Quan Sun, Zhen Jin, Rashid Jan, , Optimal control and cost-effectiveness analysis for dengue fever model with asymptomatic and partial immune individuals, 2021, 31, 22113797, 104919, 10.1016/j.rinp.2021.104919
    14. Saima Rashid, Fahd Jarad, Sobhy A. A. El-Marouf, Sayed K. Elagan, Global dynamics of deterministic-stochastic dengue infection model including multi specific receptors via crossover effects, 2023, 8, 2473-6988, 6466, 10.3934/math.2023327
    15. Shah Hussain, Elissa Nadia Madi, Naveed Iqbal, Thongchai Botmart, Yeliz Karaca, Wael W. Mohammed, Fractional Dynamics of Vector-Borne Infection with Sexual Transmission Rate and Vaccination, 2021, 9, 2227-7390, 3118, 10.3390/math9233118
    16. Salah Boulaaras, Rashid Jan, Amin Khan, Muhammad Ahsan, Dynamical analysis of the transmission of dengue fever via Caputo-Fabrizio fractional derivative, 2022, 8, 25900544, 100072, 10.1016/j.csfx.2022.100072
    17. Zahir Shah, Ebenezer Bonyah, Ebraheem Alzahrani, Rashid Jan, Nasser Aedh Alreshidi, Mariya Gubareva, Chaotic Phenomena and Oscillations in Dynamical Behaviour of Financial System via Fractional Calculus, 2022, 2022, 1099-0526, 1, 10.1155/2022/8113760
    18. Rashid Jan, Zahir Shah, Wejdan Deebani, Ebraheem Alzahrani, Analysis and dynamical behavior of a novel dengue model via fractional calculus, 2022, 15, 1793-5245, 10.1142/S179352452250036X
    19. Tao-Qian Tang, Zahir Shah, Rashid Jan, Ebraheem Alzahrani, Modeling the dynamics of tumor–immune cells interactions via fractional calculus, 2022, 137, 2190-5444, 10.1140/epjp/s13360-022-02591-0
    20. Abilasha Balakumar, Sumathi Muthukumar, Veeramani Chinnadurai, Analyzing nonlinear dynamics of dengue epidemics: A fractional order model with alert state and optimal control strategies, 2024, 45, 0143-2087, 1404, 10.1002/oca.3103
    21. Afeez Abidemi, Olumuyiwa James Peter, An optimal control model for dengue dynamics with asymptomatic, isolation, and vigilant compartments, 2024, 10, 27726622, 100413, 10.1016/j.dajour.2024.100413
    22. Khaled A. Aldwoah, Mohammed A. Almalahi, Kamal Shah, Muath Awadalla, Ria H. Egami, Dynamics analysis of dengue fever model with harmonic mean type under fractal-fractional derivative, 2024, 9, 2473-6988, 13894, 10.3934/math.2024676
    23. M. Prakash Raj, A. Venkatesh, K. Arun Kumar, M. Manivel, Mathematical Modeling of the Co‐Infection Dynamics of Dengue and Malaria Using Delay Differential Equations, 2024, 2513-0390, 10.1002/adts.202400609
    24. Rashid Jan, Salah Boulaaras, Muhammad Jawad, Karthikeyan Rajagopal, Effect of virotherapy treatment on the dynamics of tumor growth through fractional calculus, 2023, 45, 0142-3312, 2981, 10.1177/01423312231164451
    25. Jéssica V.L. Macêdo, Alberto G.S. Júnior, Maria D.L. Oliveira, César A.S. Andrade, Systematic review and meta-analysis: assessing the accuracy of rapid immunochromatographic tests in dengue diagnosis, 2024, 109, 07328893, 116227, 10.1016/j.diagmicrobio.2024.116227
    26. Jiraporn Lamwong, Puntani Pongsumpun, A fractional derivative model of the dynamic of dengue transmission based on seasonal factors in Thailand, 2025, 457, 03770427, 116256, 10.1016/j.cam.2024.116256
    27. Azhar Iqbal Kashif Butt, Muhammad Imran, Brett A. McKinney, Saira Batool, Hassan Aftab, Mathematical and Stability Analysis of Dengue–Malaria Co-Infection with Disease Control Strategies, 2023, 11, 2227-7390, 4600, 10.3390/math11224600
    28. Muhammad Younas Khan, Saif Ullah, Muhammad Farooq, Muhammad Bilal Riaz, Modeling the dynamics of dengue fever with double susceptibility and optimal control strategies, 2024, 2363-6203, 10.1007/s40808-024-02152-7
    29. Ashfaq Ahmad, Rashid Ali, Ijaz Ahmad, Muhammad Ibrahim, Fractional view analysis of the transmission dynamics of norovirus with contaminated food and water, 2024, 17, 1793-5245, 10.1142/S1793524523500729
    30. A. Refaie Ali, Hammad Alotaibi, Rashid Jan, M. Abu-Shady, Nesreen A. Yaseen, An Exact Solution of the Fractional Transient Electromagnetic Field Inside an Atmospheric Duct, 2024, 31, 1776-0852, 10.1007/s44198-024-00235-8
    31. Yuyan Qin, Lixin Yang, Jia Li, Analyzing the dynamics and optimal control of a vector-borne model with dual vertical transmission and multiple serotypes, 2024, 99, 0031-8949, 105285, 10.1088/1402-4896/ad7bf9
    32. Muhammad Farman, Ali Hasan, Changjin Xu, Kottakkaran Sooppy Nisar, Evren Hincal, Computational techniques to monitoring fractional order type-1 diabetes mellitus model for feedback design of artificial pancreas, 2024, 257, 01692607, 108420, 10.1016/j.cmpb.2024.108420
    33. R. Prem Kumar, G.S. Mahapatra, Sanjoy Basu, P.K. Santra, Global stability and sensitivity analysis of dengue transmission using four host and three vector classes along with control strategies, 2024, 0020-7160, 1, 10.1080/00207160.2024.2431683
    34. Syed Khayyam Shah, Muhammad Sarwar, Manel Hleili, Mohammad Esmael Samei, Thabet Abdeljawad, Generalized (ψ,φ)-contraction to investigate Volterra integral inclusions and fractal fractional PDEs in super-metric space with numerical experiments, 2025, 14, 2192-8029, 10.1515/nleng-2024-0032
    35. Zeshan Faiz, Iftikhar Ahmed, Dumitru Baleanu, Shumaila Javeed, A Novel Fractional Dengue Transmission Model in the Presence of Wolbachia Using Stochastic Based Artificial Neural Network, 2024, 139, 1526-1506, 1217, 10.32604/cmes.2023.029879
    36. Patience Pokuaa Gambrah, Abdul-Samad Abdul-Rahaman, Gratien Twagirumukiza, Herman Matondo Mananga, Louis Ndinyun Tawam, Arpan Hazra, A Fractional‐Order Peer Influence Mathematical Model, 2025, 2025, 1110-757X, 10.1155/jama/9649567
    37. G. M. Vijayalakshmi, M. Ariyanatchi, Fractional order modelling of Wolbachia-carrying mosquito population dynamics for dengue control, 2025, 11, 2363-6203, 10.1007/s40808-025-02360-9
    38. Ahmed S. Rashed, Mahy M. Mahdy, Samah M. Mabrouk, Rasha Saleh, Fractional Order Mathematical Model for Predicting and Controlling Dengue Fever Spread Based on Awareness Dynamics, 2025, 13, 2079-3197, 122, 10.3390/computation13050122
    39. Jiraporn Lamwong, Puntani Pongsumpun, Fractional-order modeling of dengue dynamics: exploring reinfection mechanisms with the Atangana–Baleanu derivative, 2025, 11, 2363-6203, 10.1007/s40808-025-02441-9
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3235) PDF downloads(650) Cited by(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog