Research article Special Issues

A new model of dengue fever in terms of fractional derivative

  • It is eminent that the epidemiological patterns of dengue are threatening for both the global economy and human health. The experts in the field are always in search to have better mathematician models in order to understand the transmission dynamics of epidemics models and to suggest possible control or the minimization of the infection from the community. In this research, we construct a new fractional-order system for dengue infection with carrier and partially immune classes to visualize the intricate dynamics of dengue. By using the basics of fractional theory, we determine the fundamental results of the proposed fractional-order dengue model. We obtain the basic reproduction number R0 by next generation method and present the results based on it. The stability results are established for the infection-free state of the system. Moreover, sensitivity of R0 is analyzed through partial rank correlation coefficient(PRCC) method to show the importance of different parameters in R0. The influence of different input factors is shown on the output of basic reproduction number R0 numerically. Our result showed that the threshold parameter R0 can be decreased by increasing vaccination and treatment in the system. Finally, we illustrate the solution of the suggested dengue system through a numerical scheme to notice the influence of the fractional-order ϑ on the system. We observed that the fractional-order dynamics can explain the complex system of dengue infection more precisely and accurately rather than the integer-order dynamics. In addition, we noticed that the index of memory and biting rate of vector can play a significant part in the prevention of the infection.

    Citation: Fatmawati, Rashid Jan, Muhammad Altaf Khan, Yasir Khan, Saif ullah. A new model of dengue fever in terms of fractional derivative[J]. Mathematical Biosciences and Engineering, 2020, 17(5): 5267-5287. doi: 10.3934/mbe.2020285

    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
  • It is eminent that the epidemiological patterns of dengue are threatening for both the global economy and human health. The experts in the field are always in search to have better mathematician models in order to understand the transmission dynamics of epidemics models and to suggest possible control or the minimization of the infection from the community. In this research, we construct a new fractional-order system for dengue infection with carrier and partially immune classes to visualize the intricate dynamics of dengue. By using the basics of fractional theory, we determine the fundamental results of the proposed fractional-order dengue model. We obtain the basic reproduction number R0 by next generation method and present the results based on it. The stability results are established for the infection-free state of the system. Moreover, sensitivity of R0 is analyzed through partial rank correlation coefficient(PRCC) method to show the importance of different parameters in R0. The influence of different input factors is shown on the output of basic reproduction number R0 numerically. Our result showed that the threshold parameter R0 can be decreased by increasing vaccination and treatment in the system. Finally, we illustrate the solution of the suggested dengue system through a numerical scheme to notice the influence of the fractional-order ϑ on the system. We observed that the fractional-order dynamics can explain the complex system of dengue infection more precisely and accurately rather than the integer-order dynamics. In addition, we noticed that the index of memory and biting rate of vector can play a significant part in the prevention of the infection.


    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 R0. We have shown the influence of different input factors on the basic reproduction number R0, numerically. The fractional-order system is inspected numerically to observe the influence of ϑ 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 g:R+R, then the fractional integral of g is given by

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

    where ϑ indicates the order of fractional integral and Gamma function is indicated by Γ.

    Definition 2.2. The fractional derivative of order ϑ 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 n1<ϑ<n,nN.

    Let η1,η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 C0Dϑt be the Caputo fractional derivative (CFD) of order ϑ, 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η21Eη1,η2(±ωtϑ), 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η1,η2(s)=s.Eη1,η1+η2(s)+1Γ(η2).

    Lemma 2.1. [42]. Consider a function g:R+×R5R5, 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 R5 with respect to Y(t);

    The function g(t,Y(t))Y is continuous on R5 with respect to Y(t);

    and g(t,Y(t))+Y,tR+,YR5, where , are positive constants.

    Then the fractional order system

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

    where 0<ϑ1, has a unique solution.

    Lemma 2.2. [43]. Let a function g(t)C(r,s] such that CrDϑtg(t)C[r,s] for ϑ(0,1], then we have

    g(t)=g(r)+1Γ(ϑ)CrDϑtg(ϵ)(tr)ϑ,  r<ϵ<t,t(r,s].

    Remark 1. Let a function gC(r,s], such that CrDϑtg(t)C[r,s] for ϑ(0,1]. Then by Lemma (2.2) if CrDϑtg(t)0, t(r,s), then g(t) is non-increasing function for all t[r,s], and if CrDϑtg(t)0, t(r,s), then g(t) is non-decreasing function for all t[r,s].

    In construction of the model, we represent the mosquitoes and humans population by Nv and Nh. The total mosquitoes size is divided into susceptible (Sv) and infected (Iv) compartments, while the total humans size is divided into susceptible (Sh), infected (Ih), carrier (IhA), partially immune (P), and recovered (Rh) compartments.

    We indicate the natural birth (or death) rate of mosquito and human by μv and μh, respectively, and assumed that the disease-induced death rate to be negligible. Furthermore, we assumed that a fraction υ of susceptible hosts goes to partially immune class after implementation of vaccination and τ indicates treatment of infected host. The symptomatic fraction of infected class from susceptible and partially immune classes are indicated by ψ and ω, respectively. the parameter ρ denotes the fraction of treated humans which join the partially immune class P while remanning enter to recovered class Rh. Parameter γh is used for recovery of symptomatic and asymptomatic humans while ϕ denotes the fraction of recovered humans which join the partially immune class P. The terms (bβ1NhIv) and (bβ2NhIv) indicate the infection rate per susceptible host and partially immune host, respectively, while (bβ3NhIh) denotes the infection rate of per susceptible vector. In these terms b accounts the bitting rate, β1 and β2 are the rates at which infected vectors transmit disease to susceptible and partially immune humans respectively, whereas β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

    Sh(0)0,Ih(0)0,IhA(0)0,P(0)0,Rh(0)0,Sv(0)0,Iv(0)0.

    We denote the total human population size by Nh, and vector by Nv so that, Nh=Sh+Ih+IhA+P+Rh and Nv=Sv+Iv. As the system (3.1) is independent of the state-variable Rh, thus, the proposed model (3.1) of dengue infection can be re-written in the following form without Rh:

    {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

    x1=ShNh,x2=IhNh,x3=IhANh,x4=PNh,x5=IvNv and x6=SvNv,

    then, we have x1+x2+x3+x4=1 and x5+x6=1. After replacing x6=1x5 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 C0Dϑt indicates Caputo's fractional derivative of order ϑ, the order ϑ indicates the index of memory in the system. In addition, the values

    a=β1NvNh,c=β2NvNh,  and  e=β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 R5. 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

    C0Dϑtx(t)=+A1x(t)+A2x1(t)x(t)+A3x2(t)x(t)+A4x3(t)x(t)+A5x4(t)x(t)+A6x5(t)x(t),

    where

    =[μϑh0000],A1=[(μϑh+υϑ)00000(μϑh+τϑ+γϑh)00000(μϑh+γϑh)00υϑρτϑ+ϕγϑhϕγϑhμϑh00bϑebϑe0μϑv],
    A2=[000000000ψbϑa0000(1ψ)bϑa0000000000],A3=[000000000000000000000000bϑe],A4=[000000000000000000000000bϑe],
    A5=[000000000ωbϑc0000(1ω)bϑc0000000000],A6=[bϑa00000000000000000bϑc000000].

    Symbolize the above

    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),

    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 Ω={(x1,x2,x3,x4,x5)R5+:0x1+x2+x3+x4M1,0x5M2} 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 W=min(μϑh,((1ρ)τϑ+(1ϕ)γϑh+μϑh),((1ϕ)γϑh+μϑh),μϑh). Solving the above inequality, we obtain

    (x1(t)+x2(t)+x3(t)+x4(t))(x1(0)+x2(0)+x3(0)+x4(0)μϑhW)Eϑ(Wtϑ)+μϑhW,

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

    (x1(t)+x2(t)+x3(t)+x4(t))μϑhWM1,

    taking the same steps for the last equation of system (3.4), we get x5(t)M2, where M2=bϑeM1bϑeM1+μϑv. Hence, the closed set Ω 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 C0Dϑtx1,C0Dϑtx2,C0Dϑtx3,C0Dϑtx4, and C0Dϑtx5 to zero of the fractional system (3.4) without infection, and get

    E0(x01,x02,x03,x04,x05)=(μϑhμϑh+υϑ,0,0,υϑμϑh+υϑ,0).

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

    F=[ψbϑax1x5+ωbϑcx4x5(1ψ)bϑax1x5+(1ω)bϑcx4x5bϑe(1x5)(x2+x3)]and  V=[(μϑh+τϑ+γϑh)x2(μϑh+γϑh)x3μϑvx5],

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

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

    which gives

    FV1=[00ψbϑax01+ωbϑcx04μϑv00(1ψ)bϑax01+(1ω)bϑcx04μϑvbϑe(μϑh+τϑ+γϑh)bϑe(μϑh+γϑh)0].

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

    ρ(FV1)=bϑeμϑv(μϑh+υϑ)[(1ψ)bϑaμϑh+(1ω)bϑcυϑμϑh+γϑh+ψbϑaμϑh+ωbϑcυϑμϑh+γϑh+τϑ],
    R0=bϑeμϑv(μϑh+υϑ)[(1ψ)bϑaμϑh+(1ω)bϑcυϑμϑh+γϑh+ψbϑaμϑh+ωbϑcυϑμϑh+γϑh+τϑ].

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

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

    where BMm×m(R) and 0<ϑ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 ⇔∣arg(λi)∣>ϑπ2 is fulfilled for all eigenvalues λi of matrix B.

    Theorem 4.2. [47]. The steady state of fractional system (2.6) is locally asymptotically stable (LAS) if arg(λi)∣>ϑπ2 for all eigenvalues λi's of g(y) at steady state, otherwise unstable.

    In order to prove the LAS of infection-free steady state E0 it is sufficient to confirm that all of the eigenvalues of the matrix of dynamics given by J(E0) 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 R0<1 and arg(λi)∣>ϑπ2 for all eigenvalues λi's of the Jacobian matrix J(E0) 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(E0)=[(υϑ+μϑ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].

    The characteristic equation in term of Λ of the above matrix JE0 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=ψbϑax01+ωbϑcx04 and B=(1ψ)bϑax01+(1ω)bϑcx04. From (4.2) it is clear the argument of the first two eigenvalues λ=μϑh and λ=(μϑh+υϑ) satisfy the necessary condition given (i.e., arg(λi)∣>ϑπ2 for all ϑ(0,1)). Furthermore, if R0<1, then all ai>0, for i=1,2,3, and it is easy to show that a1a2a3>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 R0<1.

    Proof. Let υ=0, and (x1,x2,x3,x4,x5) be the solution of system (3.4) with suitable initial conditions (x1(0),x2(0),x3(0),x4(0),x5(0)) in Ω. Here, it is clear that the fractional system has only one equilibrium E0 on the boundary of Ω. Therefore, to achieve the target, it is enough to prove that the solution (x1,x2,x3,x4,x5) tends to the infection-free equilibrium as time tends to infinity. Since xix0i, 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 FV, and if R0=ρ(FV1)<1, then the eigenvalues of FV lies in the left half plane. Consequently, each positive solution of the fractional system (4.3) fulfills limty1=0, limty2=0 and limty3=0 by Theorem 4.1. By the comparison theory of fractional differential equations [48], we have limtx2=0, limtx3=0 and limtx5=0. Then from system (3.4), we have

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

    this can be further converted into

    C0DϑtX=ΛBX,

    having the solution

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

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

    Eϑ,1(ωtϑ)ttϑωΓ(1ϑ);   0<ϑ<1  and  ω>0,

    also the real parts of eigenvalues of B are negative, so it can be observed that X(t)X0 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 ϑ=1 if R0>1.

    Proof. The threshold parameter R0 is of high significance, which measures the level of infection in the system. If R0 less than one then the asymptomatic stability of infection-free equilibrium(E0) prevent the system from the persistent of the infection. If R0 greater than one then the sufficiently close solution of fractional system (3.4) to E0 moves away from it due to its instability in the interior of Ω by Theorem (4.3). For uniform persistence of our system (3.4) at ϑ=1, we will prove that it fulfills the criteria of Theorem 4.3 in [49] for R0>1. Select X=R5 and E=Ω. Here, the maximal invariant set on the boundary of Ω is the set {E0} 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 E0 is unstable for R0>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 R0. 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 R0 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 β1, and ϑ are highly influential with PRCC values 0.6427, and -0.5112, respectively. Sensitivity analysis illustrated that b,β1, and ϑ are critical parameters in R0, 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 R0 with PRCC and corresponding p values.
    Parameter Interpretation PRCC values p values
    ϑ index of memory or fractional order -0.5112 0.0000
    μh recruitment rate of humans host 0.1095 0.0023
    υ vaccinated fraction of susceptible host + 0.3133 0.0000
    ψ fraction of incidence rate from Sh to Ih + 0.1804 0.0010
    β1 transmission probability from vectors to susceptible host + 0.6427 0.0000
    β2 transmission probability from vectors to partially immune host + 0.5094 0.0000
    τ treatment fraction of infected host individuals -0.2110 0.0000
    γh recovery rate of host individuals -0.3361 0.0000
    β3 transmission probability from infected hosts to susceptible vectors +0.2320 0.0000
    μv recruitment rate of mosquitoes vector -0.3197 0.0000
    ω fraction of incidence rate from P to Ih + 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 R0 and PRCC outcomes for R0.

    Next, we illustrate the influence of vaccination and treatment in on the threshold R0 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 β3 and γh on the basic reproduction number R0 of the proposed system of dengue fever.

    Figure 2.  Illustration of basic reproduction number R0 (a) with the variation of vaccination υ and treatment τ, (b) and with the variation of β3 and γ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=(x1,x2,x3,x4,x5)R+5, G(t,g(t)) is used for a continuous real valued vector function, which additionally satisfies the Lipschitz condition and g0 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,T] with h=T0m is the step size and mN. 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
    μv natural mortality and recruitment rate of vector 0.032300 & 0.029410 [53,54]
    β1 transmission probability from vectors to susceptible host 0.75 [10]
    β2 transmission probability from vectors to partially immune host 0.375 Assumed
    ϑ index of memory or fractional order variable Assumed
    b bitting rate of vectors 0.50 [10]
    τ treatment fraction of infected host 0.40 Assumed
    υ vaccinated fraction of susceptible host 0.20 Assumed
    β3 transmission probability from infected humans to mosquitoes 0.75 [10]
    ψ fraction of incidence rate from Sh to Ih 0.48 Assumed
    γh recovery rate of host individuals 0.3288330 [53,54]
    ω fraction of incidence rate from P to Ih 0.54 Assumed
    μh natural mortality and recruitment rate of humans 0.004500 & 0.000046 [53,54]
    ρ 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 (Ih), asymptomatic host (IhA), partially immune (P) host and infected vector (Iv) with the variation of fractional-order ϑ. We noticed that the variation of fractional-order ϑ 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 Ih, (b) Carrier host(Asymptomatic) IhA, (c) Partially immune P and (d) Infected host Iv individuals with different values of fractional-order ϑ, i.e., ϑ=0.6,0.8,1.0.
    Figure 4.  Illustration of the path moment of (a) Infected host Ih, (b) Carrier host(Asymptomatic) IhA, (c) Partially immune P and (d) Infected vector Iv individuals with different vaccination rates υ, i.e., υ=0.2,0.4,0.6.
    Figure 5.  Illustration of the path moment of (a) Infected host Ih, (b) Carrier host(Asymptomatic) IhA, (c) Partially immune P and (d) Infected vector Iv individuals with different treatment rates τ, i.e., τ=0.4,0.6,0.8.
    Figure 6.  Illustration of the path moment of (a) Infected host Ih, (b) Carrier host(Asymptomatic) IhA, (c) Partially immune P and (d) Infected vector Iv 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 R0. We established stability results of infection-free steady-state of the system. Sensitivity analysis of R0 is carried out in order to know the contribution of input factors in the results of R0 and observed that b,μv, and β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 R0>1. We obtained feasible results for the dynamics of dengue infection with the variation of memory index ϑ 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) ϑ 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] Center for Disease Control and Prevention. Dengue, Accessed January 14, 2013. URL http://www.cdc.gov/dengue/.
    [2] D. J. Gubler, Dengue and dengue hemorrhagic fever, Clin. Microbiol. Rev., 11(1998), 480-496.
    [3] M. G. Guzman, G. Kouri, J. Bravo, L. Valdes, S. Vasquez, S. B. Halstead, Effect of age on outcome of secondary dengue 2 infections, Int. J. Infect. Dis., 6(2002), 118-124.
    [4] East, Susie (6 April 2016). World's first dengue fever vaccine launched in the Philippines. CNN. Archived from the original on 18 October 2016. Retrieved 17 October 2016.
    [5] Epidemiology, https://www.cdc.gov/dengue/epidemiology/index.html. Retrieved 7 March 2018.
    [6] Y. Xiao, T. Zhao, S. Tang, Dynamics of an infectious disease with media/ psychology induced non-smooth incidence, Math. Biosci. Eng., 10(2013), 445-461.
    [7] A. Wang, Y. Xiao, A Filippov system describing media effects on the spread of infectious diseases, Nonlinear Analysis, Hybrid Syst., 11(2014), 84-97.
    [8] L. Esteva, C. Vargas, Analysis of a dengue disease transmission model, Math. Biosci., 150(1998), 131-151.
    [9] L. Esteva, C. Vargas, A model for dengue disease with variable human population, J. Math. Biol., 38(1999), 220-240.
    [10] M. Derouich, A. Boutayeb, E. H. Twizell, A model of Dengue fever, BioMed. Eng. Online, 2(2003).
    [11] J. J. Tewa, J. L. Dimi, S. Bowang, Lypaunov functions for a dengue disease transmission model, Chaos Solit. Fract., 39(2009), 936-941.
    [12] H. S. Rodrigues, M. T. T. Monteiro, D. F. M. Torres, Vaccination models and optimal control strategies to Dengue, Math. Biosci., 5(2014), 1-12.
    [13] R. Jan, Y. Xiao, Effect of partial immunity on transmission dynamics of dengue disease with optimal control, Math. Method Appl. Sci., 42(2019), 1967-1983.
    [14] D. S. Burke, A. Nisalak, D. E. Johnson, R. M. Scott, A prospective study of dengue infections in Bangkok, Am. J. Trop. Med. Hyg., 38(1988), 172-180.
    [15] T. P. Endy, S. Chunsuttiwat, A. Nisalak, D. H. Libraty, S. Green, A. L. Rothman, et al., Epidemiology of inapparent and symptomatic acute dengue virus infection: A prospective study of primary school children in Kamphaeng Phet, Thailand Am. J. Epidemiol., 156(2002), 40-51.
    [16] L. F. Chaves, L. C. Harrington, C. L. Keogh, A. M. Nguyen, U. D. Kitron, Blood feeding patterns of mosquitoes: Random or structured?, Front. Zool., 7(2010), 1-11.
    [17] C. Vinauger, L. Buratti, C. R. Lazzari, Learning the way to blood: First evidence of dual olfactory conditioning in a blood sucking insect, Rhodnius prolixus. I. Appetitive learning, J. Exp. Biol., 214(2011), 3032-3038.
    [18] Y. Carrera, G. A. Rosa, E. J. Vernon-Carter, A. Ramirez, A fractional-order Maxwell model for non-Newtonian fluids, Phys. A, 482(2017), 276-285.
    [19] R. Jan, M. A. Khan, P. Kumam, P. Thounthong, Modeling the transmission of dengue infection through fractioal derivatives, Chaos Solit. Fract., 127(2019), 189-216.
    [20] T. Kosztolowicz, K. D. Lewandowska, Application of fractional differential equations in modelling the subdiffusion-reaction process, Math. Model Nat. Pheno., 8(2013), 44-54.
    [21] C. N. Angstmann, B. I. Henry, A. V. McGann, A fractional-order infectivity SIR model, Phys. A, 452(2016), 86-93.
    [22] M. Alquran, I. Jaradat, Delay-asymptotic solutions for the time fractional delay-type wave equation, Phys. A, 527(2019), 121-275.
    [23] M. A. Khan, S. Ullah, M. Farooq, A new fractional model for tuberculosis with relapse via Atangana-Baleanu derivative, Chaos Solit. Fract., 116(2018), 227-238.
    [24] S. Qureshi, A. Atangana, Mathematical analysis of dengue fever outbreak by novel fractional operators with field data, Phys. A, 526(2019), 121-127.
    [25] A. Atangana, Application of fractional calculus to epidemiology, Fractional Dynamic, (eds. C. Cattani, H. M. Srivastava, X. Yang), De Gruyter Open, (2015), 174-190.
    [26] A. Atangana, Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination?. Chaos Solit. Fract., 136(2020), 109860.
    [27] E. F. Doungmo Goufo, Y. Khan, Q. A. Chaudhry, HIV and shifting epicenters for COVID-19, an alert for some countries,Chaos Solit. Fract., 139(2020), 110030.
    [28] M. A. Khan, A. Atangana, Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative, Alex. Eng. J., (2020).
    [29] W. Gao, P. Veeresha, D. G. Prakasha, H. M. Baskonus, Novel dynamic structures of 2019-nCoV with nonlocal operator via powerful computational technique, Biology, 9(2020), 107.
    [30] E. Ilhan, I. O. Kiymaz, A generalization of truncated M-fractional derivative and applications to fractional differential equations, Appl. Math. Non. Sci., 5(2020), 171-188.
    [31] W. Gao, P. Veeresha, D. G. Prakasha, H. M. Baskonus, G. Yel, New numerical results for the time-fractional Phi-four equation using a novel analytical approach, Symmetry, 12(2020), 478.
    [32] T. A. Sulaiman, H. Bulut, S. S. Atas, Optical solitons to the fractional Schrodinger-Hirota equation, Appl. Math. Non. Sci., 4(2019),535-542.
    [33] J. Singh, D. Kumar, Z. Hammouch, A. Atangana, A fractional epidemiological model for computer viruses pertaining to a new fractional derivative, Appl. Math. Comput., 316(2018), 504-515.
    [34] C. Cattani, A review on Harmonic wavelets and their fractional extension, J. Adv. Eng. Comput., 2(2018), 224-238.
    [35] X. J. Yang, F. Gao, A new technology for solving diffusion and heat equations, Therm. Sci., 21(2017), 133-140.
    [36] M. A. Khan, N. Iqbal, Y. Khan, E. Alzahrani, A biologoical mathematical model of vector-host disease with saturated treatment function and optimal control strategies, Math. Biosci. Eng., 17 (2020), 3972.
    [37] S. K. Panda, T. Abdeljawad, C. Ravichandran, A complex valued approach to the solutions of Riemann-Liouville integral, Atangana-Baleanu integral operator and non-linear Telegraph equation via fixed point method, Chaos Solit. Fract., 130(2020), 109439.
    [38] X. J. Yang, D. Baleanu, M. P. Lazaveric, M. S. Cajic, Fractal boundary value problems for integral and differential equations with local fractional operators, Therm. Sci., 19(2015), 959-966.
    [39] R. Almeida, Analysis of a fractional SEIR model with treatment, Appl. Math. Lett., 84(2018), 56-62.
    [40] R. M. A. Carvalho, C. M. A. Pinto, Immune response in HIV epidemics for distinct transmission rates and for saturated CTL response, Math. Model. Nat. Phenom., 14(2019), 13.
    [41] K. Diethelm, The analysis of fractional differential equations: An application-oriented exposition using operators of Caputo type, Springer,2004.
    [42] S. Mao, R. Xu, Y. Li, A fractional order SIRS model with standard incidence rate, J. Beihua Univ., 12(2012), 379-382.
    [43] Z. M. Odibat, N. T. Shawagteh, Genralized Taylor's formula, Appl. Math. Comput., 186(2007), 286-293.
    [44] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and application of fractional differential equations, Elsevier, The Netherlands, 2006.
    [45] C. Castillo-Chavez, Z. Feng, W. Huang, On the computation of R0 and its role on global stability, Mathematical Approach for Emerging and Reemerging Infectious Diseases: An Introduction, Springer, p. 229, 2002.
    [46] D. Matignon, Stability results for fractional differential equations with applications to control processing, Comput. Eng. Syst. Appl., 2(1996), 963-968.
    [47] C. Li, Y. Ma. Fractional dynamical system and its linearization theorem, Nonlin. Dyn., 71(2013), 621-633.
    [48] Z. Wang, D. Yang, T. Ma, N. Sun, Stability analysis for nonlinear fractional-order system based on comparison principle, Nonlin. Dyn., 75(2014), 387-402.
    [49] H. I. Freedman, M. X. Tang, S. G. Ruan, Uniform persistence and flows near a closed positively invariant set, J. Dyn. Differ. Equ., 6(1994), 583-600.
    [50] S. Marino, I. B. Hogue, C. J. Ray, D Kirschner, A methodolojgy for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254(2008), 178-196.
    [51] A. Jajarmi, D. Baleanu, A new fractional analysis on the interaction of HIV with CD4+ T-cells, Chaos Solit. Fract., 113(2018), 221-229.
    [52] C. Li, F. Zeng, Numerical methods for fractional calculus, Chapman and Hall/CRC, (2015).
    [53] KKM. 2007, Health Facts 2007.
    [54] KKRI, Dengue fever is still high in south sulawesi. Tribun Timur Makassar, Newspaper fact, 2009.
  • 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
  • Reader Comments
  • © 2020 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(6031) PDF downloads(389) Cited by(37)

Figures and Tables

Figures(6)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog