Research article Special Issues

On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies

  • In this work, an epidemic model of a susceptible, exposed, infected and recovered SEIR-type is established for the distinctive dynamic compartments and epidemic characteristics of COVID-19 as it spreads across a population with a heterogeneous rate. The proposed model is investigated using a novel approach of fractional calculus known as piecewise derivatives. The existence theory is demonstrated through the establishment of sufficient conditions. In addition, result related to Hyers-Ulam stability is also derived for the considered model. A numerical method based on modified Euler procedure is also constructed to simulate the approximate solutions of the proposed model by employing various values of fractional orders. We testified the numerical results by using real available data of Japan. In addition, some results for the SEIR-type model are also presented graphically using the stochastic process, and the obtained results are discussed.

    Citation: Mdi Begum Jeelani, Kamal Shah, Hussam Alrabaiah, Abeer S. Alnahdi. On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies[J]. AIMS Mathematics, 2023, 8(11): 27268-27290. doi: 10.3934/math.20231395

    Related Papers:

    [1] Fawaz K. Alalhareth, Seham M. Al-Mekhlafi, Ahmed Boudaoui, Noura Laksaci, Mohammed H. Alharbi . Numerical treatment for a novel crossover mathematical model of the COVID-19 epidemic. AIMS Mathematics, 2024, 9(3): 5376-5393. doi: 10.3934/math.2024259
    [2] Abdon Atangana, Saima Rashid . Analysis of a deterministic-stochastic oncolytic M1 model involving immune response via crossover behaviour: ergodic stationary distribution and extinction. AIMS Mathematics, 2023, 8(2): 3236-3268. doi: 10.3934/math.2023167
    [3] Nadiyah Hussain Alharthi, Mdi Begum Jeelani . Analyzing a SEIR-Type mathematical model of SARS-COVID-19 using piecewise fractional order operators. AIMS Mathematics, 2023, 8(11): 27009-27032. doi: 10.3934/math.20231382
    [4] Muhammad Farman, Ali Akgül, Kottakkaran Sooppy Nisar, Dilshad Ahmad, Aqeel Ahmad, Sarfaraz Kamangar, C Ahamed Saleel . Epidemiological analysis of fractional order COVID-19 model with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(1): 756-783. doi: 10.3934/math.2022046
    [5] Hasib Khan, Jehad Alzabut, J.F. Gómez-Aguilar, Praveen Agarwal . Piecewise mABC fractional derivative with an application. AIMS Mathematics, 2023, 8(10): 24345-24366. doi: 10.3934/math.20231241
    [6] Tao Xie, Wenqing Zheng . Robustness analysis of Cohen-Grossberg neural network with piecewise constant argument and stochastic disturbances. AIMS Mathematics, 2024, 9(2): 3097-3125. doi: 10.3934/math.2024151
    [7] 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. AIMS Mathematics, 2023, 8(3): 6466-6503. doi: 10.3934/math.2023327
    [8] Ishtiaq Ali, Sami Ullah Khan . Dynamics and simulations of stochastic COVID-19 epidemic model using Legendre spectral collocation method. AIMS Mathematics, 2023, 8(2): 4220-4236. doi: 10.3934/math.2023210
    [9] Ramsha Shafqat, Ateq Alsaadi . Artificial neural networks for stability analysis and simulation of delayed rabies spread models. AIMS Mathematics, 2024, 9(12): 33495-33531. doi: 10.3934/math.20241599
    [10] Jinji Du, Chuangliang Qin, Yuanxian Hui . Optimal control and analysis of a stochastic SEIR epidemic model with nonlinear incidence and treatment. AIMS Mathematics, 2024, 9(12): 33532-33550. doi: 10.3934/math.20241600
  • In this work, an epidemic model of a susceptible, exposed, infected and recovered SEIR-type is established for the distinctive dynamic compartments and epidemic characteristics of COVID-19 as it spreads across a population with a heterogeneous rate. The proposed model is investigated using a novel approach of fractional calculus known as piecewise derivatives. The existence theory is demonstrated through the establishment of sufficient conditions. In addition, result related to Hyers-Ulam stability is also derived for the considered model. A numerical method based on modified Euler procedure is also constructed to simulate the approximate solutions of the proposed model by employing various values of fractional orders. We testified the numerical results by using real available data of Japan. In addition, some results for the SEIR-type model are also presented graphically using the stochastic process, and the obtained results are discussed.



    The deadly COVID-19 disease broke out at the tail end of 2019 in Wuhan, a prominent Chinese metropolis. The aforementioned disease has rapidly spread around the world. By the end of March 2020, the World Health Organization (WHO) proclaimed it a global pandemic [30,68,72]. Globally, more than six million deaths linked to COVID-19 have been recorded. Infected individuals in the same queue number from 600 to 700. People's health, the way of life and the situation of the economy have all been significantly disturbed. To find the best cure for the aforementioned illness, experts and researchers are working around the clock [69]. The COVID-19 is a persistent major health problem and concern with a global focus. Breakout poses a new type of worldwide threat, regardless of the fact that there are major gaps in our knowledge of COVID-19 epidemiology, transmission dynamics, research techniques and management. For instance, several US states and other countries around the world adopted lockdown and reopening procedures. Hence, various management strategies have been used to control the epidemic's development.

    Technology has advanced epidemiology to the point where several infectious diseases are now examined for treatment, control, curing and other outcomes [52]. It should be mentioned that the study of many diseases also heavily incorporates mathematical biology. As a result, over the past several decades, there has been a tremendous advancement in the mathematical modeling of infectious diseases [58,63] and mathematical modeling has become more common in research. In order to successfully control a variety of diseases, such as those listed in citations [35,36], secure public health approaches are developed with the help of mathematical models. The dynamic behaviour of infections can be studied using these mathematical models, as well as spatiotemporal patterns, academics have studied COVID-19 from a variety of angles during the past three years, with the help of mathematical models. To create effective methods for managing this illness, researchers in this field are taking a range of approaches (several recent studies are included as [9,21,28,70]). A mathematical model was recently used to examine the effects of vaccinations in nursing homes [37]. Researchers [46] examined mathematical modeling and practical COVID-19 epidemic intervention approaches. It is worth mentioning that the area devoted to model COVID-19 using various analysis and tools to investigate the transmission dynamics of the said disease has attracted very well in the last three years. For instance, global stability and cost-effectiveness analysis of COVID-19 with impact of the environment by using data from Ghana has been studied [5]. In the same way, authors [6] have studied optimal control and comprehensive cost-effectiveness analysis for COVID-19. Proceeding with the same process, sensitivity and optimal control analysis have been performed [7,8] for COVID-19 mathematical models. Different numerical results have been presented there. Moreover, authors [61] have studied the role of asymptomatic infections in the COVID-19 epidemic via complex networks and stability analysis. Also, authors [31] have used stochastic concepts along with environmental white noise to investigate a COVID-19 mathematical model. Researchers [64] have used mathematical models for forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan. In all the mentioned studies researchers have used classical or stochastic type models for their study.

    As far as we are aware, the concept of the classical derivative has been intensively investigated in the subject of epidemiology. As we know, the classical differential operators have local nature and a number of inherited, long-term and short-term memory processes. To better comprehend the aforementioned process, fractional calculus has therefore attracted a lot more interest in recent years. It has grown in popularity as a result of the dynamic properties that have shown a variety of uses in practical contexts, including biological and physical processes [64]. Like normal calculus, fractional calculus has a lengthy history [47]. Several authors have explored the topic from various perspectives [48,49,62] and it has been found that it has numerous applications in science and technology [50,55]. The fractional order derivative has a higher degree of freedom [24] due to its non-locality. Therefore, in the mathematical modeling of infectious diseases, the aforementioned derivative might be preferred to the standard order derivative. In the past, a number of authors have produced insightful work, such as the existence theory of solutions to fractional differential equations [4] and qualitative findings [51,59], respectively. In order to investigate fractional order differential and integral equations (FODIEs) for approximate or analytical results, a variety of tools and methods have been developed (such as the fractional visco-elasticity model in [2], the fractional calculus in mechanical system modeling [1] and the fractional model for creep/recovery testing of asphalt mixtures [22].

    It is authenticating that the area devoted to fractional order epidemiological models has been extensively considered in the last two decades. Researchers have investigated numerous diseases including dengue infection, malaria disease, TB, HBV and HCV, AIDS, cancer and so on. Researchers have extended classical susceptible, infected and recovered (SIR), susceptible, exposed, infected and recovered (SEIR), susceptible, exposed, infected, asymptomatic and recovered (SEIAR) type models to fractional orders [23]. For recent work in fractional orders models of infectious diseases, we refer to [3,25]. In addition, problems devoted to describe some physical phenomenons have also been well investigated using the concept of fractional order derivatives [57,60].

    The majority of real-world issues are unpredictable to some extent, which traditional mathematical models are unable to account for. The idea of stochastic mathematical differential equations has been put forth and used extensively in recent decades, with remarkable outcomes. Other issues, however, demonstrate non-locality tendencies instead of randomness, including long-range dependence, fractal processes, power law processes and crossover behaviors, suggesting that physical occurrences exhibit a variety of behaviors. A class of fractional derivatives was proposed to overcome these problems. However, these operators still do a poor job of describing crossover behavior. The idea of short memory with a real or complex order derivative was developed to characterize the aforementioned behavior for the first time. The piecewise idea has been proved to be more powerful than the stated, despite the fact that fractional derivatives have extended memory capability [66]. Various concepts including fractal-fractional derivative, fractional order derivative with singular and non-singular kernels and various additional derivative operator have been defined recently to look at the crossover qualities and important work on nonlocal operators and their applications [13,29,40] refers to fractional dynamics of cellulose degradation, [10] refers to existence and uniqueness with applications to epidemiology and [12] refers to local and nonlocal operators. The crossover dynamical behavior has not been addressed, despite the fact that stochastic equation considerations lead to more realistic conclusions [11]. This behavior is seen in a lot of real-world process models, including heat flow, fluid flow, and many complex advection problems [71]. In fractional calculus, the timing of crossovers cannot be determined using the exponential or Mittag-Leffler mappings because the typical fractional order derivative fails to accurately describe crossover behavior that is common in real-world problems. The phenomena such as earthquakes, pendulum motion, the current economic instability in less developed nations, etc are few examples have crossover behavior. This crossover behavior is exemplified by piecewise equations with fractional order derivatives. Through the analysis of numerous models, certain crucial elements in this regard have recently been discovered [16]. Along with several applications, the authors created classical and global piecewise derivatives. As shown in papers [19,38,39,43,45,67], a number of infectious illness models have lately been studied using non-singular and power-law type operators.

    With the aforementioned significance in mind, we planned to concentrate on these fundamental issues in our work by using a model motivated from [56] that has been specially tailored to reflect the dynamics of said disease as well as the curbs of human response to it. A usual SEIR model that allows for lengthy incubation, we first recreated the epidemic dynamics within a single community with a particular social pattern. Here, using piecewise derivative, we formulated our model by including birth rate and natural death rate in the model investigated as

    PCC0Drt(S)(t)=aβS(I+qE)NdSPCC0Drt(E)(t)=βS(I+qE)NEδdEPCC0Drt(I)(t)=EδIγdIPCC0Drt(R)(t)=IγdR, (1.1)

    where PCC0Drt stands for piecewise Caputo derivative, which can be described for any function say y as

    PCC0Drt(y(t))={C0Dt(y(t))=dydt,  0<tt1,C0Drt(y(t))=1Γ(1r)tt1(tη)ry(η)dη,  t1<tT, (1.2)

    where 0Drt represents the usual Caputo fractional order derivative. The flow chart of our model is given in Figure 1 and the parameters are described in Table 1.

    Table 1.  Nomenclature and their discerption.
    Nomenclature Description
    S Compartment of uninfected population
    E Compartment of exposed population
    I Compartment of infected population
    R Compartment of recovered population
    N0 Total initial papulation
    N Total population at time t
    a Recruitment (birth) rate
    β is the typical number of contacts made by each person each time
    γ is the time required for a person to pass away or recover after entering the infectious stage
    δ is the typical time for a person who has been exposed to an infection
    q is contact rate
    d natural death rate

     | Show Table
    DownLoad: CSV
    Figure 1.  Flow chart of our established model (1.1).

    The birth rate is a, and the carrier can either be infected or exposed. The fractions SIN and SEN indicate the probability of a random contact between a susceptible person and an infected person, and the probability of a contact between a susceptible person and an exposed person respectively. The scaling factor q in the model is used for the rate of interaction with an exposed person, which results in a different likelihood of transmission than contact with an infected person. Following the transition state, individuals may contract an infectious class at a rate determined by 1δ. The same way, 1γ represents the rate at which people pass away or recover from an infection.

    It is notable that such dynamical concerns are given some unique perspectives by existence theory with piecewise derivatives of fractional orders. The hypothesis suggests that such physical problems can be resolved. In recent years, a variety of numerical techniques have proven to be particularly successful for classical fractional order systems [32]. For instance, the Range-Kutta method was used to address a number of fractional order problems [18], and the authors [44] also took use of a groundbreaking parameter estimate method, according to researchers. Problems of fractional order were solved by a novel numerical method in [33], and in [27] studied a fractional order system by using the finite difference approach. For a different set of non-integer order problems, improved finite-difference techniques were also used [17,42,65]. We used actual data of Japan for infected, recovered cases whose sources have listed as [73,74,75] and a numerical approach to elaborate the numerical analysis of the model under consideration at various fractional orders. In fact, it was stated that the memory rate was highlighted by the order of the time derivative, and the memory function was the kernel of the fractional-order derivative [14]. The boundedness were established for the existence of solution and feasible region was determined. The existence and uniqueness of approximation solutions were then investigated for the aforementioned model using the Banach and Schauder fixed point theorems. In addition, the results devoted to Hyers-Ulam stability were also derived. The concerned stability was investigated about the best approximate or exact solution. For the existence theory, we used research from Arzelá-Ascoli [41] and the Schauder fixed point theorem [15].

    It is remarkable that in the actual world, everything is susceptible to chance, including people and animal movement. Stochastic models are used to mathematically interpret this effect. Mathematical models have been investigated by using the afore said concepts. It is worth mentioning that there are significant applications in modeling real world process by using stochastic calculus. Hence, due to this fact, researchers have recently investigated COVID-19 mathematical models by using the aforementioned area. The said differential equations have been increasingly used to model other infectious disease as well. Here we referenced to articles [20,34,54]. Motivated from the said importance of stochastic calculus, we attempted to simulate the proposed model (1.1) under the stochastic white noise as

    dr(S)(t)=[aβS(I+qE)NdS]dtr+ρ1Sdb1(t)dr(E)(t)=[βS(I+qE)NEδdE]dtr+ρ2Edb2(t)dr(I)(t)=[EδIγdI]dtr+ρ3Idb3(t)dr(R)(t)=[IγdR]dtr+ρ4Rdb4(t), (1.3)

    where bi denotes Brownian motion with bi(0)=0. Further, ρi>0, for i=1,2,3,4 and denotes the intensity of white noise. Here, it should be kept in mind that we simulated our model (1.2) using the numerical scheme developed in [53].

    This manuscript is arranged as follows: A substantial introduction to our work is contained in Part one, and crucial results that we need for this work are given in Section two. We established the existence theory for a rough solution to the proposed model using the fixed point theory in Section three. Section four discusses the numerical approach for a rough solution to the suggested model. Section five is devoted to illustrating our conclusions. Finally, Section six offers a brief summary and explanation of the numerical outcomes.

    Recollecting some basic results as follows:

    Definition 2.1. [16] If ϖ is a differentiable function with r>0, then the piecewise integral is described by considering I=[0,T], I1=[0,t1], I2=(t1,T] as

    PC0Irtϖ(t)={t10ϖ()d,  tI1,1Γ(r)tt1(t)r1ϖ()d(), tI2,

    where PC0It stands for classical integration in I1 and represents Riemann-Liouville integration in I2.

    Definition 2.2. [16] Let 0<r1, and if ϖC(I) is differentiable, then the classical and fractional order piecewise derivative is defined as

    PCC0Drtϖ(t)={ϖ(t),  tI1,C0Drtϖ(t),  tI2.

    Lemma 2.3. [16] Let ϖL(I)C(I) and ψL(I), then the solution of the given problem

    PCC0Drtϖ(t)=ψ(t),  0<r1

    is derived as

    ϖ(t)={ϖ0+t0ψ()d, tI1,ϖ(t1)+1Γ(r)tt1(t)r1ψ()d(),  tI2.

    For (1.1), we derived some axillary results. The feasible region is given in Remark 2.4.

    Remark 2.4. Let N be the total papulation at any time t, we have

    N=S+E+I+R. (2.1)

    Applying the piecewise derivative defined in (2.1) w.r.t 't' from model (1.1), we get

    PCC0DrtN(t)=adN, (2.2)

    which on taking Laplace transform and using t yields

    N(t)ad.

    Hence, the feasible region is described as

    Φ={(S,E,I,R)R4+:Nad}.

    Putting left hand sides of model (1.1) equal to zero and solving the equations, the disease free equilibrium is obtained as follows [32]. From the first two equations of (1.1), one has

    aβS(I+qE)NdS=0βS(I+qE)NEδdE=0EδIγdI=0IγdR=0,

    and on solving at disease free equilibrium state when I0=0, we have S0=ad, E0=0, R0=0. Hence, local equilibrium is given by

    E0=(S0,0,0,0)=(ad,0,0,0).

    In the same line, we can also compute the endemic equilibria as

    S=aδγd(1dδ)(1+γd)IdδγE(t)=1+γdγII=γ[βaδ(γ+q(1+γd))d(1δd)(1+γd)]β(1δd)(1+γd)(γ+q(1+γd))R=Idγ.

    In addition, the fundamental reproductive number can be computed as R0=βqδ1+δd. Clearly, if R0<1, the local equilibrium point will be locally asymptotically stable. Further, the 3D profile of the said number is given in the following Figure 2.

    Figure 2.  3D profile of R0 for (1.1).

    Here, we developed sufficient results for the qualitative theory of existence and uniqueness using the fixed point approach. Assume that F:I×RR is a nonlinear continuous function. In view of Lemma 2.3, the solution of

    PCC0Drtϖ(t)=F(t,ϖ),  0<r1,ϖ(0)=ϖ0 (3.1)

    is given as

    ϖ(t)={ϖ0+t0F(,ϖ())d, tI1,ϖ(t1)+1Γ(r)tt1(t)r1F(,ϖ())d(),  tI2, (3.2)

    where

    ϖ(t)={S(t)E(t)I(t)R(t),   ϖ0={S0E0I0R0,   ϖ(t1)={S(t1)E(t1)I(t1)R(t1),   F(t,ϖ(t))={F1(t,ϖ(t))={F1(t,ϖ(t)), tI1,F1(t,ϖ(t)), tI2,F2(t,ϖ(t))={F2(t,ϖ(t)), tI1,F2(t,ϖ(t)), tI2,F3(t,ϖ(t))={F3(t,ϖ(t)), tI1,F3(t,ϖ(t)), tI2,F4(t,ϖ(t))={F4(t,ϖ(t)),tI1,F4(t,ϖ(t)),tI2. (3.3)

    Let >t2t>t1>0 with Banach space defined as H=C(I)×C(I)×C(I)×C(I) under the norm

    ϖ=maxtI|ϖ(t)|.

    These hypotheses hold for further results.

    (C1) If LF>0 is real number and ϖ, ˉϖH, then

    |F(t,ϖ)F(t,ˉϖ)|LF|ϖˉϖ|.

    (C2) For real values CF>0 and MF>0, we have

    |F(t,ϖ(t))|CF|ϖ|+MF.

    Theorem 3.1. Under the hypotheses (C1) and (C2), if there exists a closed bounded subset B={ϖH: ϖR1,2, R1,2>0} of H, where

    R1,2max{|ϖ0|+t1MF1t1CF, tI1,|ϖ(t1)|Γ(r+1)+TrMF(Γ(r+1)TrCF, tI2,

    then (3.1) has at least one solution. Consequently, (1.1) has at least one solution.

    Proof. Considering BH is closed and bounded, such that

    B={ϖH: ϖR1,2, R1,2>0}.

    If Q:BB, then the operator is defined as

    Q(ϖ)={ϖ0+t0F(,ϖ())d, tI1,ϖ(t1)+1Γrtt1(t)σ1F(,ϖ())d(),  tI2. (3.4)

    For ϖB we have

    |Q(ϖ)(t)|{|ϖ0|+t10|F(,ϖ())|d,|ϖ(t1)|+1Γ(r)tt1(t)r1|F(ϖ())|d(),{|ϖ0|+t10[CG|ϖ|+MF]d,|ϖ(t1)|+1Γ(r)tt1(t)r1[CF|ϖ|+MF]d(),{|ϖ0|+t1[CFR1,2+MF]R1,2, tI1,|ϖ(t1)|+TrΓ(r+1)[CFR1,2+MF]R1,2, tI2,

    where for tI2, we put |(t1)r(t2)r|Tr. Hence, we have that Q(ϖ)R1,2 which yields that Q(B)B. Thus, Q maps bounded set to bounded and Q is bounded operator. Since F is continuous function, Q is also continuous operator. Next, for complete continuity, consider tm<tnI1, then

    |Q(ϖ)(tn)Q(ϖ)(tm)|=|tn0F(,ϖ())dtm0F(,ϖ())d|tntm|F(,ϖ())|dtntm[CF|ϖ|+MF]d(CFR1,2+MF)[tntm]. (3.5)

    From (3.5), we see that tmtn, then

    |Q(ϖ)(tn)Q(ϖ)(tm)|0, as tmtn.

    Also Q is bounded operator. So,

    Q(ϖ)(tn)Q(ϖ)(tm)0, as tmtn.

    Hence, Q is equi-continuous in this case. In addition, take tm<tnI2 and consider

    |Q(ϖ)(tn)Q(ϖ)(tm)|=|1Γ(r)tn0(tn)r1F(,ϖ())d1Γ(r)tm0(tm)r1F(,ϖ())d|1Γ(r)tm0[(tm)r1(tn)r1]|G(,ϖ())|d+1Γ(r)tntm(tn)r1|F(,ϖ())|d1Γ(r)[tm0[(tm)r1(tn)r1]d+tntm(tn)r1d](CF|ϖ|+MF)(CFR1,2+MF)Γ(r+1)[trntrm+2(tntm)r]. (3.6)

    Further from (3.6), we see that

    |Q(ϖ)(tn)Q(ϖ)(tm)|0, as tmtn.

    Additionally, Q is bounded over I2 so it is uniformly continuous. Hence,

    Q(ϖ)(tn)Q(ϖ)(tm)0, as tmtn.

    Therefore, Q is equi-continuous in the I2 interval and Q is equi-continuous mapping over I1I2. Thus Q is a relatively compact operator. By using the Arzelá-Ascoli theorem stated in [41], operator Q is completely continuous. In view of the Schauder's fixed point theorem [15], (3.1) has at least one solution, meaning that (1.1) has at least one solution.

    Theorem 3.2. Inview of Hypothesis (C1), and if max{TLF,TrΓ(r+1)LF}<1 holds, then (1.1) has a unique solution.

    Proof. Let Q:HH be the mapping defined as

    Q(ϖ)={ϖ0+t0F(,ϖ())d, tI1,ϖ(t1)+1Γrtt1(t)σ1F(,ϖ())d(),  tI2.

    Let ϖ, ˉϖH, then over I1, one has

    Q(ϖ)Q(ˉϖ)=maxtI1|t0F(,ϖ())dt10F(,ˉϖ())d|TLFϖˉϖ. (3.7)

    From (3.7), we have

    Q(ϖ)Q(ˉϖ)TLFϖˉϖ. (3.8)

    By the same fashion for tI2, we have

    Q(ϖ)Q(ˉϖ)=maxtI2|1Γ(r)tt1(t)r1F(,ϖ())d1Γ(r)tt1(t)r1F(,ˉϖ())d|TrΓ(r+1)LFϖˉϖ. (3.9)

    From (3.9), we have

    Q(ϖ)Q(ˉϖ)TrΓ(r+1)LFϖˉϖ. (3.10)

    Hence, from (3.8) and (3.10), we see that Q is a contraction operator and (3.1) has a unique solution. Consequently, (1.1) has a unique solution.

    Remark 3.3. Let there exist a non-decreasing function ςC[0,T], such that

    (i) |ς(t)|ε, tI.

    In addition, the solution of the problem

    PCC0Drtϖ(t)=F(t,ϖ)+ς(t),  0<r1,ϖ(0)=ϖ0 (3.11)

    is given as

    ϖ(t)={ϖ0+t0F(,ϖ())d+t0ς()d, tI1,ϖ(t1)+1Γ(r)tt1(t)r1F(,ϖ())d()+1Γ(r)tt1(t)r1ς()d(),  tI2, (3.12)

    In view of (i), we have from (3.12) that

    |ϖ(t)(ϖ0+t0F(,ϖ())d)|εt1, if tI1, (3.13)

    and

    |ϖ(t)(ϖ(t1)+1Γ(r)tt1(t)r1F(,ϖ())d())|+εt1Γ(r+1), if tI2. (3.14)

    Theorem 3.4. Under the hypothesis (C1), and Remark 3.3, the solution of problem (3.1) is Hyers-Ulam stable if max{LFt1,εt1LFΓ(r+1)}<1.

    Proof. Let ϖ be any solution of (3.1), for which we have a unique solution ˆϖ. Then for tI1 we have

    |ϖ(t)ˆϖ(t)|=|ϖ(t)(ϖ0+t0F(,ˆϖ())d)||ϖ(t)(ϖ0+t0F(,ϖ())d)|+|t0F(,ϖ())dt0F(,ˆϖ())d|εt1+LFt1ϖˆϖ,

    which further yields that

    ϖˆϖεt11LFt1. (3.15)

    In the same way, if tI2, repeating the same process, one has

    ϖˆϖεt1Γ(r+1)1εt1LFΓ(r+1). (3.16)

    Hence, from (3.15) and (3.16), we have that the solution is Hyers-Ulam stable. Consequently, the solution of (1.1) is Hyers-Ulam stable.

    For numerical presentation of (1.1), we constructed a numerical method for the two sub-intervals of I. The numerical scheme for the piecewise problem is like an integer order numerical scheme as established in [14] by using ϖ=(S,E,I,R) as

    S(tn+1)={Sn1(tn1)+h2F1[tn1+h2,ϖn1(tn1)+K12],  tI1,Sn(tn)+hrΓ(r+1)F1(tn,ϖn(tn))+hr2Γ(r+1)[K2+K3],  tI2, (4.1)

    where h=tn+1tn, and

    K1=F1(tn1,yn1(tn1)),  K2=F1(tn,ϖn(tn)),K3=F1(tn+2hrΓ(r+1)Γ(2r+1),ϖ(tn)+2hrΓ(r+1)Γ(2r+1)F1(tn,yn(tn))). (4.2)

    In the same way, we established for other compartments as given by

    E(tn+1)={Sn1(tn1)+h2F1[tn1+h2,ϖn1(tn1)+K12],  tI1,En(tn)+hrγ(r+1)F1(tn,ϖn(tn))+hr2Γ(r+1)[K2+K3],  tI2, (4.3)
    I(tn+1)={Sn1(tn1)+h2F1[tn1+h2,ϖn1(tn1)+K12],  tI1,In(tn)+hrγ(r+1)F1(tn,ϖn(tn))+hr2Γ(r+1)[K2+K3],  tI2, (4.4)

    and for the last class as

    R(tn+1)={Sn1(tn1)+h2F1[tn1+h2,ϖn1(tn1)+K12],  tI1,Rn(tn)+hrγ(r+1)F1(tn,ϖn(tn))+hr2Γ(r+1)[K2+K3],  tI2. (4.5)

    In this section, we present the numerical simulation in Figures 311 using the obtained scheme under classical and piecewise derivative concepts. We divided the whole interval into two sub-intervals and checked the first interval for integer order derivative while the second interval was tested on different fractional orders in the sense of the Caputo derivative by using the data given in Table 2.

    Table 2.  Description and specification of real values for the variables used in (1.1).
    Parameters Numerical value
    S 76.013223 in millions[73]
    E 8.206109 in million [73]
    I 21.150371 in million[73]
    R 20.330297 in million [73]
    a 0.09 per day assumed
    β 0.08 [56]
    γ 10 days [56]
    δ 6 days [56]
    q 0.001 [56]
    d 0.0012 [75]

     | Show Table
    DownLoad: CSV
    Figure 3.  Graphical presentations of approximate solutions of S and E corresponding to the given fractional orders for the proposed model.
    Figure 4.  Graphical presentations of approximate solutions of I and R corresponding to the given fractional orders for the proposed model.
    Figure 5.  Graphical presentations of approximate solutions of S and E corresponding to the given fractional orders for the proposed model.
    Figure 6.  Graphical presentations of approximate solutions of I and R corresponding to the given fractional orders for the proposed model.
    Figure 7.  Graphical presentations of approximate solutions of S and E corresponding to the given fractional orders for the proposed model.
    Figure 8.  Graphical presentations of approximate solutions of I and R corresponding to the given fractional orders for the proposed model.
    Figure 9.  Presentation of real data plots for confirmed and recovered cases for the given months.
    Figure 10.  Graphical presentations of approximate solutions of S and E corresponding to the given values of ρi for the proposed model in stochastic concepts.
    Figure 11.  Graphical presentations of approximate solutions of I and R corresponding to the given values of ρi for proposed model in stochastic concepts.

    Figure 3 expressed the susceptible and exposed populations which declined, and then became stable as the remaining classes increased on both intervals as presented in Figure 4. The single curve was for the first interval and it showed integer order classical behavior from [0,t1]=[0,50], while the four different curves showed the global order derivative behavior on (t1,t2]=(500,400]. Here, we have considered the first set of fractional orders values in (0,0.75].

    Here, we considered another set of fractional orders values in (0.75,0.95]. Figure 5 contains the graphical presentations of susceptible, and exposed classes while in Figure 6 presents the dynamical behaviors of infected and recovered classes respectively corresponding to given fractional orders.

    Furthermore, Figures 7 and 8 show the dynamical behaviors of different classes of the proposed model by considering another set of fractional orders values in (0.90,1]. The crossover behaviors in dynamics of different classes can be observed clearly at t1=50.

    From Figures 38 we have observed the crossover behavior at the given points. The numerical simulation of all four compartments has been shown for three sets of different fractional orders. The bending effects are also shown on t1 describing the piecewise derivative dynamics. Since the population of the infected class is decreasing, as a result the recovered class population is rising. Here, we have presented the real data fit for the recorded values of confirmed infected and recovered cases per day in Japan from January 2022 to August 2022 in Figure 9.

    Here, we considered the same numerical values given in Table 2 to simulated the aforesaid considered model under the stochastic concepts described in (1.2). We simulate (1.2) corresponding to two different values of white noise. Figures 10 and 11 show, the dynamics of S,E and I,R of the proposed model using stochastic concepts. Here, we have taken the fractional order as one and simulated the results against two distinct values of white noise, ρi=0,0.02 for i=1,2,3,4. In addition to keeping in mind the importance of stochastic calculus, we have also attempted to simulate the proposed model under the mentioned concept in Figures 10 and 11, respectively. The concerned behavior due to stochastic is more realistic than the usual or traditional fractional derivatives. The crossover properties of both stochastic and fractional orders can be handled by our suggested investigation.

    In the presented research work, we have established a model under the concept of piecewise equations with the fractional order Caputo derivative. By using Schauder and Banach's fixed point theory, we have established the existence theory of solution to the proposed model and a numerical scheme has been established by using the modified Euler method. The results have been presented by using real data graphically for various fractional orders. Additionally, concerned results have been demonstrated and compared with real data in the case of reported infected individuals. Since many real-world problems are subject to abrupt changes in their state of rest or uniform motion which is also called crossover behavior traditional derivatives of either classical or fractional form cannot demonstrate this effect properly. The mentioned effect can be excellently explained by using piecewise equations with fractional order derivatives. Further, fractional order derivatives are more flexible and keep a greater degree of freedom. In the future, such a type of analysis can be extended to more complex dynamical problems involving Mittag-Leffler and fractal-fractional type derivatives. In addition, the aforesaid model will be investigated under the stochastic fractional order differential equations using non-singular differential operators.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-RG23124).

    The authors declare that they have no conflict of interest.



    [1] Z. Ali, F. Rabiei, K. Shah, T. Khodadadi, Fractal-fractional order dynamical behavior of an HIV/AIDS epidemic mathematical model, Eur. Phys. J. Plus, 136 (2021), 36. https://doi.org/10.1140/epjp/s13360-020-00994-5 doi: 10.1140/epjp/s13360-020-00994-5
    [2] M. M. Amirian, Y. Jamali, The concepts and applications of fractional order differential calculus in modeling of viscoelastic systems: a primer, Crit. Rev. Biomed. Eng., 47 (2019), 249–276. https://doi.org/10.1615/CritRevBiomedEng.2018028368 doi: 10.1615/CritRevBiomedEng.2018028368
    [3] A. J. Arenas, G. Gonzalez-Parra, B. M. Chen-Charpentier, Construction of nonstandard finite difference schemes for the SI and SIR epidemic models of fractional order, Math. Comput. Simulat., 121 (2016), 48–63. https://doi.org/10.1016/j.matcom.2015.09.001 doi: 10.1016/j.matcom.2015.09.001
    [4] M. Arfan, K. Shah, A. Ullah, Fractal-fractional mathematical model of four species comprising of prey-predation, Phys. Scripta, 96 (2021), 124053. DOI 10.1088/1402-4896/ac2f37 doi: 10.1088/1402-4896/ac2f37
    [5] J. K. K. Asamoah, M. A. Owusu, Z. Jin, F. T. Oduro, A. Abidemi, E. O. Gyasi, Global stability and cost-effectiveness analysis of COVID-19 considering the impact of the environment: using data from Ghana, Chaos Soliton. Fract., 140 (2020), 110103. https://doi.org/10.1016/j.chaos.2020.110103 doi: 10.1016/j.chaos.2020.110103
    [6] J. K. K. Asamoah, E. Okyere, A. Abidemi, S. E. Moore, G. Q. Sun, Z. Jin, et al., Optimal control and comprehensive cost-effectiveness analysis for COVID-19, Results Phys., 33 (2022), 105177. https://doi.org/10.1016/j.rinp.2022.105177 doi: 10.1016/j.rinp.2022.105177
    [7] J. K. K. Asamoah, Z. Jin, G. Q. Sun, B. Seidu, E. Yankson, A. Abidemi, et al., Sensitivity assessment and optimal economic evaluation of a new COVID-19 compartmental epidemic model with control interventions, Chaos Soliton. Fract., 146 (2021), 110885. https://doi.org/10.1016/j.chaos.2021.110885 doi: 10.1016/j.chaos.2021.110885
    [8] J. K. K. Asamoah, C. S. Bornaa, B. Seidu, Z. Jin, Mathematical analysis of the effects of controls on transmission dynamics of SARS-CoV-2, Alex. Eng. J., 59 (2020), 5069–5078. https://doi.org/10.1016/j.aej.2020.09.033 doi: 10.1016/j.aej.2020.09.033
    [9] J. K. K. Asamoah, Fatmawati, A fractional mathematical model of heartwater transmission dynamics considering nymph and adult amblyomma ticks, Chaos Soliton. Fract., 174 (2023), 113905. https://doi.org/10.1016/j.chaos.2023.113905 doi: 10.1016/j.chaos.2023.113905
    [10] J. K. K. Asamoah, E. Okyere, E. Yankson, A. A. Opoku, A. Adom-Konadu, E. Acheampong, et al., Non-fractional and fractional mathematical analysis and simulations for Q fever, Chaos Soliton. Fract., 156 (2022), 111821. https://doi.org/10.1016/j.chaos.2022.111821 doi: 10.1016/j.chaos.2022.111821
    [11] J. K. K. Asamoah, Fractal-fractional model and numerical scheme based on Newton polynomial for Q fever disease under Atangana-Baleanu derivative, Results Phys., 34 (2022), 105189. https://doi.org/10.1016/j.rinp.2022.105189 doi: 10.1016/j.rinp.2022.105189
    [12] A. Atangana, S. I. Araz, Nonlinear equations with global differential and integral operators:existence, uniqueness with application to epidemiology, Results Phys., 20 (2021), 103593. https://doi.org/10.1016/j.rinp.2020.103593 doi: 10.1016/j.rinp.2020.103593
    [13] A. Atangana, D. Baleanu, New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.48550/arXiv.1602.03408 doi: 10.48550/arXiv.1602.03408
    [14] M. A. Abdulwasaa, M. S. Abdo, K. Shah, T. A. Nofal, S. K. Panchal, S. V. Kawale, et al., Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India, Results Phys., 20 (2021), 103702. https://doi.org/10.1016/j.rinp.2020.103702 doi: 10.1016/j.rinp.2020.103702
    [15] R. P. Agarwal, S. Arshad, D. Regan, V. Lupulescu, A Schauder fixed point theorem in semilinear spaces and applications, Fixed Point Theory Appl., 2013 (2013), 306. https://doi.org/10.1186/1687-1812-2013-306 doi: 10.1186/1687-1812-2013-306
    [16] A. Atangana, S. I. Araz, New concept in calculus: piecewise differential and integral operators, Chaos Soliton. Fract., 145 (2021), 110638. https://doi.org/10.1016/j.chaos.2020.110638 doi: 10.1016/j.chaos.2020.110638
    [17] A. J. Arenas, G. González-Parra, B. M. Chen-Charpentier, Construction of nonstandard finite difference schemes for the SI and SIR epidemic models of fractional order, Math. Comput. Simulat., 121 (2016), 48–63. https://doi.org/10.1016/j.matcom.2015.09.001 doi: 10.1016/j.matcom.2015.09.001
    [18] M. S. Arshad, D. Baleanu, M. B. Riaz, M. Abbas, A novel 2-stage fractional Runge-Kutta method for a time fractional logistic growth model, Discrete Dyn. Nat. Soc., 2020 (2020), 1020472. https://doi.org/10.1155/2020/1020472 doi: 10.1155/2020/1020472
    [19] S. Boccaletti, W. Ditto, G. Mindlin, A. Atangana, Modeling and forecasting of epidemic spreading: The case of Covid-19 and beyond, Chaos Soliton. Fract., 135 (2020), 109794. https://doi.org/10.1016/j.chaos.2020.109794 doi: 10.1016/j.chaos.2020.109794
    [20] S. Banihashemi, H. Jafari, A. Babaei, A stable collocation approach to solve a neutral delay stochastic differential equation of fractional order, J. Comput. Appl. Math., 403 (2022), 113845. https://doi.org/10.1016/j.cam.2021.113845 doi: 10.1016/j.cam.2021.113845
    [21] I. I. Bogoch, A. Watts, A. Thomas-Bachli, C. Huber, M. U. Kraemer, K. Khan, Pneumonia of unknown aetiology in Wuhan, China: potential for international spread via commercial air travel, J. Travel Med., 27 (2020), taaa008. https://doi.org/10.1093/jtm/taaa008 doi: 10.1093/jtm/taaa008
    [22] C. Celauro, C. Fecarotti, A. Pirrotta, A. C. Collop, Experimental validation of a fractional model for creep/recovery testing of asphalt mixtures, Constr. Build. Mater., 36 (2012), 458–466. https://doi.org/10.1016/j.conbuildmat.2012.04.028 doi: 10.1016/j.conbuildmat.2012.04.028
    [23] Y. Chen, F. Liu, Q. Yu, T. Li, Review of fractional epidemic models, Appl. Math. Model., 97 (2021), 281–307. https://doi.org/10.1016/j.apm.2021.03.044 doi: 10.1016/j.apm.2021.03.044
    [24] A. Carpinteri, F. Mainardi, Fractals and fractional calculus in continuum mechanics, Vienna: Springer, 1997. https://doi.org/10.1007/978-3-7091-2664-6
    [25] C.T. Deressa, G. F. Duressa, Analysis of Atangana-Baleanu fractional-order SEAIR epidemic model with optimal control, Adv. Differ. Equ., 2021 (2021), 174. https://doi.org/10.1186/s13662-021-03334-8 doi: 10.1186/s13662-021-03334-8
    [26] Z. Dai, Y. Peng, H. A. Mansy, R. H. Sandler, T. J. Royston, A model of lung parenchyma stress relaxation using fractional viscoelasticity, Med. Eng. Phys., 37 (2015), 752–758. https://doi.org/10.1016/j.medengphy.2015.05.003 doi: 10.1016/j.medengphy.2015.05.003
    [27] Z. J. Fu, Z. C. Tang, H. T. Zhao, P. W. Li, T. Rabczuk, Numerical solutions of the coupled unsteady nonlinear convection-diffusion equations based on generalized finite difference method, Eur. Phys. J. Plus, 134 (2019), 272. https://doi.org/10.1140/epjp/i2019-12786-7 doi: 10.1140/epjp/i2019-12786-7
    [28] A. B. Gumel, S. Ruan, T. Day, J. Watmough, F. Brauer, P. van den Driessche, et al., Modelling strategies for controlling SARS out breaks, Proc. R. Soc. Lond. B., 271 (2004), 2223–2232. https://doi.org/10.1098/rspb.2004.2800 doi: 10.1098/rspb.2004.2800
    [29] E. F. D. Goufo, Application of the Caputo-Fabrizio fractional derivative without singular kernel to Korteweg-de Vries-Burgers equation, Math. Model. Anal., 21 (2016), 188–198. https://doi.org/10.3846/13926292.2016.1145607 doi: 10.3846/13926292.2016.1145607
    [30] D. S. Hui, E. I. Azhar, T. A. Madani, F. Ntoumi, R. Kock, O. Dar, et al., The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health–The latest 2019 novel coronavirus outbreak in Wuhan, China, B. Math. Biol., 91 (2020), 264–66. https://doi.org/10.1016/j.ijid.2020.01.009 doi: 10.1016/j.ijid.2020.01.009
    [31] S. Hussain, E. N. Madi, H. Khan, H. Gulzar, S. Etemad, S. Rezapour, et al., On the stochastic modeling of COVID-19 under the environmental white noise, J. Funct. Space., 2022 (2022), 4320865. doilinkhttps://doi.org/10.1155/2022/4320865 doi: 10.1155/2022/4320865
    [32] A. A. Hamou, E. Azroul, Z. Hammouch, A. L. Alaoui, On dynamics of fractional incommensurate model of Covid-19 with nonlinear saturated incidence rate, MedRxiv, 2021 (2021), 07, https://doi.org/10.1101/2021.07.18.21260711 doi: 10.1101/2021.07.18.21260711
    [33] M. T. Hoang, O. F. Egbelowo, Dynamics of a fractional-order hepatitis B epidemic model and its solutions by nonstandard numerical schemes, In: Mathematical modelling and analysis of infectious diseases, Cham: Springer, 2020,127–153. https://doi.org/10.1007/978-3-030-49896-2_5
    [34] G. Jumarie, Stochastic differential equations with fractional Brownian motion input, Int. J. Syst. Sci., 24 (1993), 1113–1131. https://doi.org/10.1080/00207729308949547 doi: 10.1080/00207729308949547
    [35] S. Kumar, A. Kumar, B. Samet, H. Dutta, A study on fractional host-parasitoid population dynamical model to describe insect species, Numer. Meth. Part. D. E., 37 (2021), 1673–1692. https://doi.org/10.1002/num.22603 doi: 10.1002/num.22603
    [36] M. M. Khalsaraei, An improvement on the positivity results for 2-stage explicit Runge-Kutta methods, J. Comput. Appl. Math., 235 (2010), 137–143. https://doi.org/10.1016/j.cam.2010.05.020 doi: 10.1016/j.cam.2010.05.020
    [37] R. Kahn, I. Holmdahl, S. Reddy, J. Jernigan, M. J. Mina, R. B. Slayton, Mathematical modeling to inform vaccination strategies and testing approaches for voronavirus disease 2019 (COVID-19) in nursing homes, Clin. Infect. Dis., 74 (2022), 597–603. https://doi.org/10.1093/cid/ciab517 doi: 10.1093/cid/ciab517
    [38] M. A. Khan, A. Atangana, Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative, Alex. Eng. J., 59 (2020), 2379–2389. https://doi.org/10.1016/j.aej.2020.02.033 doi: 10.1016/j.aej.2020.02.033
    [39] M. A. Khan, A. Atangana, E. Alzahrani, E. Fatmawati, The dynamics of COVID-19 with quarantined and isolation, Adv. Differ. Equ., 2020 (2020), 425. https://doi.org/10.1186/s13662-020-02882-9 doi: 10.1186/s13662-020-02882-9
    [40] A. M. Lopes, J. T. Machado, Fractional order models of leaves, J. Vib. Control, 20 (2014), 998–1008. https://doi.org/10.1177/1077546312473323 doi: 10.1177/1077546312473323
    [41] R. Li, S. Zhong, C. Swartz, An improvement of the Arzela-Ascoli theorem, Topol. Appl., 159 (2012), 2058–2061. http://doi.org/10.1016/j.topol.2012.01.014 doi: 10.1016/j.topol.2012.01.014
    [42] R. Lewandowski, Z. Pawlak, Dynamic analysis of frames with viscoelastic dampers modelled by rheological models with fractional derivatives, J. Sound Vib., 330 (2011), 923–936. https://doi.org/10.1016/j.jsv.2010.09.017 doi: 10.1016/j.jsv.2010.09.017
    [43] B. Li, H. Liang, L. Shi, Q. He, Complex dynamics of Kopel model with nonsymmetric response between oligopolists, Chaos Soliton. Fract., 156 (2022), 111860. https://doi.org/10.1016/j.chaos.2022.111860 doi: 10.1016/j.chaos.2022.111860
    [44] F. Liu, K. Burrage, Novel techniques in parameter estimation for fractional dynamical models arising from biological systems, Comput. Math. Appl., 62 (2011), 822–833. https://doi.org/10.1016/j.camwa.2011.03.002 doi: 10.1016/j.camwa.2011.03.002
    [45] B. Li, H. Liang, Q. He, Multiple and generic bifurcation analysis of a discrete Hindmarsh-Rose model, Chaos Soliton. Fract., 146 (2021), 110856. https://doi.org/10.1016/j.chaos.2021.110856 doi: 10.1016/j.chaos.2021.110856
    [46] J. Mondal, S. Khajanchi, Mathematical modeling and optimal intervention strategies of the COVID-19 outbreak, Nonlinear Dyn., 109 (2022), 177–202. https://doi.org/10.1007/s11071-022-07235-7 doi: 10.1007/s11071-022-07235-7
    [47] J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
    [48] F. C. Meral, T. J. Royston, R. Magin, Fractional calculus in viscoelasticity: an experimental study, Commun. Nonlinear Sci., 15 (2010), 939–945. https://doi.org/10.1016/j.cnsns.2009.05.004 doi: 10.1016/j.cnsns.2009.05.004
    [49] R. L. Magin, Fractional calculus in bioengineering, Crit. Rev. Biomed. Eng., 32 (2004), 1–104. https://doi.org/10.1615/critrevbiomedeng.v32.i1.10 doi: 10.1615/critrevbiomedeng.v32.i1.10
    [50] R. L. Magin, Fractional Calculus in bioengineering, Redding: Begell House, 2006.
    [51] F. Mainardi, An historical perspective on fractional calculus in linear viscoelasticity, Fract. Calc. Appl. Anal., 15 (2012), 712–717. https://doi.org/10.2478/s13540-012-0048-6 doi: 10.2478/s13540-012-0048-6
    [52] I. Nesteruk, Statistics based predictions of coronavirus 2019-nCoV spreading in mainland China, MedRxiv, 4 (2020), 1988–1989. https://doi.org/10.1101/2020.02.12.20021931 doi: 10.1101/2020.02.12.20021931
    [53] O. A. Omar, R. A. Elbarkouky, H. M. Ahmed, Fractional stochastic modelling of COVID-19 under wide spread of vaccinations: Egyptian case study, Alex. Eng. J., 61 (2022), 8595–8609. https://doi.org/10.1016/j.aej.2022.02.002. doi: 10.1016/j.aej.2022.02.002
    [54] J. C. Pedjeu, G. S. Ladde, Stochastic fractional differential equations: Modeling, method and analysis, Chaos Soliton. Fract., 45 (2012), 279–293. https://doi.org/10.1016/j.chaos.2011.12.009 doi: 10.1016/j.chaos.2011.12.009
    [55] A. Y. Rossikhin, M. V. Shitikova, Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids, Appl. Mech. Rev., 50 (1997), 15–67. https://doi.org/10.1115/1.3101682 doi: 10.1115/1.3101682
    [56] A. Radulescu, C. Williams, K. Cavanagh, Management strategies in a SEIR-type model of COVID 19 community spread, Sci. Rep., 10 (2020), 21256. https://doi.org/10.1038/s41598-020-77628-4 doi: 10.1038/s41598-020-77628-4
    [57] Y. B. Sang, Critical Kirchhoff-Choquard system involving the fractional p-Laplacian operator and singular nonlinearities, Topol. Method. Nonl. An., 58 (2021), 233–274. https://doi.org/10.12775/TMNA.2020.070 doi: 10.12775/TMNA.2020.070
    [58] K. Shah, R. Din, W. Deebani, P. Kumam, Z. Shah, On nonlinear classical and fractional order dynamical system addressing COVID-19, Results Phys., 24 (2021), 104069. https://doi.org/10.1016/j.rinp.2021.104069 doi: 10.1016/j.rinp.2021.104069
    [59] M. Shimizu, W. Zhang, Fractional calculus approach to dynamic problems of viscoelastic materials, JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, 42 (1999), 825–837. https://doi.org/10.1299/jsmec.42.825 doi: 10.1299/jsmec.42.825
    [60] Y. B. Sang, S. H. Liang, Fractional Kirchhoff-Choquard equation involving Schrodinger term and upper critical exponent, J. Geom. Anal., 32 (2022), 5. https://doi.org/10.1007/s12220-021-00747-5 doi: 10.1007/s12220-021-00747-5
    [61] L. Stella, A. P. Martínez, D. Bauso, P. Colaneri, The role of asymptomatic infections in the COVID-19 epidemic via complex networks and stability analysis, SIAM J. Control Optim., 60 (2022), S119–S144. https://doi.org/10.1137/20M1373335 doi: 10.1137/20M1373335
    [62] D. Valério, J. T. Machado, V. Kiryakova, Some pioneers of the applications of fractional calculus, Fract. Calc. Appl. Anal., 17 (2014), 552–578. https://doi.org/10.2478/s13540-014-0185-1 doi: 10.2478/s13540-014-0185-1
    [63] Y. Wang, Z. Wei, J. Cao, Epidemic dynamics of influenza-like diseases spreading in complex networks, Nonlinear Dyn., 101 (2020), 1801–1820. https://doi.org/10.1007/s11071-020-05867-1 doi: 10.1007/s11071-020-05867-1
    [64] J. T. Wu, K. Leung, G. M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study, Lancet, 395 (2020), 689–697. https://doi.org/10.1016/S0140-6736(20)30260-9 doi: 10.1016/S0140-6736(20)30260-9
    [65] B. Wang, L. Li, Y. Wang, An efficient nonstandard finite difference scheme for chaotic fractional-order Chen system, IEEE Access, 8 (2020), 98410–98421. https://doi.org/10.1109/ACCESS.2020.2996271 doi: 10.1109/ACCESS.2020.2996271
    [66] G. C. Wu, M. Luo, L. L. Huang, S. Banerjee, Short memory fractional differential equations for new memristor and neural network design, Nonlinear Dyn., 100 (2020), 3611–3623. https://doi.org/10.1007/s11071-020-05572-z doi: 10.1007/s11071-020-05572-z
    [67] A. Zeb, A. Atangana, Z. A. Khan, S. Djillali, A robust study of a piecewise fractional order COVID-19 mathematical model, Alex. Eng. J., 61 (2022), 5649–5665. https://doi.org/10.1016/j.aej.2021.11.039 doi: 10.1016/j.aej.2021.11.039
    [68] S. Zhao, Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, et al., Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak, Int. J. Infect. Dis., 92 (2020), 214–217. https://doi.org/10.1016/j.ijid.2020.01.050 doi: 10.1016/j.ijid.2020.01.050
    [69] S. Zhao, S. S. Musa, Q. Lin, J. Ran, G. Yang, W. Wang, et al., Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020: A data-driven modelling analysis of the early outbreak, J. Clin. Med., 9 (2020), 388. http://doi.org/10.3390/jcm9020388 doi: 10.3390/jcm9020388
    [70] P. Zhou, X. L. Yang, X. G. Wang, B. Hu, L. Zhang, W. Zhang, et al., A pneumonia outbreak associated with a new coronavirus of probable bat origin, Nature, 579 (2020), 270–273. https://doi.org/10.1038/s41586-020-2012-7 doi: 10.1038/s41586-020-2012-7
    [71] Y. Zhang, H. Sun, H. H. Stowell, M. Zayernouri, S. E. Hansen, A review of applications of fractional calculus in Earth system dynamics, Chaos Soliton. Fract., 102 (2017), 29–46. https://doi.org/10.1016/j.chaos.2017.03.051 doi: 10.1016/j.chaos.2017.03.051
    [72] Naming the coronavirus disease (COVID-19) and the virus that causes it, Available from: World Health Organization (WHO), 2019. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it.
    [73] World Health Organization, Japan situation, 2020. Available from: https://covid19.who.int/region/wpro/country/jp.
    [74] Japan COVID coronavirus statistics, 2023. Available from: https://www.worldometers.info/coronavirus/country/japan/.
    [75] Japan COVID cases, 2021. Available from: https://www.nytimes.com/interactive/2021/world/japan-covid-cases.html.
  • This article has been cited by:

    1. Muhammad Umer Saleem, Muhammad Farman, Kottakkaran Sooppy Nisar, Aqeel Ahmad, Zainab Munir, Evren Hincal, Muhammad Aqeel, Investigation and application of a classical piecewise hybrid with a fractional derivative for the epidemic model: Dynamical transmission and modeling, 2024, 19, 1932-6203, e0307732, 10.1371/journal.pone.0307732
    2. Xiaojie Mu, Daqing Jiang, Dynamics caused by the mean-reverting Ornstein–Uhlenbeck process in a stochastic predator–prey model with stage structure, 2024, 179, 09600779, 114445, 10.1016/j.chaos.2023.114445
  • Reader Comments
  • © 2023 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(1707) PDF downloads(58) Cited by(2)

Figures and Tables

Figures(11)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog