Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination

  • The current manuscript investigates a model of the spread of polio under the condition of vaccination by using the novel modified Atangana-Baleanu-Caputo (mABC) fractional derivative. This problem has been studied for non-zero solutions under the modified operator. The series-type solution has been obtained through the application of a Laplace transformation, along with the decomposition technique and Adomian polynomial for the nonlinear terms. The qualitative analysis for the solution of the model has been tested by using fixed point theory. The stability of the solution is also crucial for a dynamical system; therefore, it was checked by using the T-Picard method. With the help of the approximate scheme, numerical simulations were conducted for the proposed model by using different fractional orders and transmission parameters. Based on the obtained positivity of the solutions and numerical stability, we have established the analysis of the mABC operator in the field of fractional calculus and other physical sciences.

    Citation: Mati ur Rahman, Mehmet Yavuz, Muhammad Arfan, Adnan Sami. Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination[J]. AIMS Biophysics, 2024, 11(1): 97-120. doi: 10.3934/biophy.2024007

    Related Papers:

    [1] Yasir Nadeem Anjam, Mehmet Yavuz, Mati ur Rahman, Amna Batool . Analysis of a fractional pollution model in a system of three interconnecting lakes. AIMS Biophysics, 2023, 10(2): 220-240. doi: 10.3934/biophy.2023014
    [2] Mohammed Alabedalhadi, Mohammed Shqair, Ibrahim Saleh . Analysis and analytical simulation for a biophysical fractional diffusive cancer model with virotherapy using the Caputo operator. AIMS Biophysics, 2023, 10(4): 503-522. doi: 10.3934/biophy.2023028
    [3] Shaimaa A. M. Abdelmohsen, D. Sh. Mohamed, Haifa A. Alyousef, M. R. Gorji, Amr M. S. Mahdy . Mathematical modeling for solving fractional model cancer bosom malignant growth. AIMS Biophysics, 2023, 10(3): 263-280. doi: 10.3934/biophy.2023018
    [4] Stephen Haller, Rabin Gerrah, Sandra Rugonyi . Towards virtual surgery planning: the modified Blalock-Taussig Shunt. AIMS Biophysics, 2020, 7(3): 169-188. doi: 10.3934/biophy.2020014
    [5] Omar El Deeb, Joseph El Khoury Edde . COVID19 vaccines as boosters or first doses: simulating scenarios to minimize infections and deaths. AIMS Biophysics, 2024, 11(2): 239-254. doi: 10.3934/biophy.2024014
    [6] Etimad Alattar, Eqbal Radwan, Khitam Elwasife . Improvement in growth of plants under the effect of magnetized water. AIMS Biophysics, 2022, 9(4): 346-387. doi: 10.3934/biophy.2022029
    [7] Mary Jane Beilby, Sabah Al Khazaaly . Re-modeling Chara action potential: II. The action potential form under salinity stress. AIMS Biophysics, 2017, 4(2): 298-315. doi: 10.3934/biophy.2017.2.298
    [8] Mehmet Yavuz, Kübra Akyüz, Naime Büşra Bayraktar, Feyza Nur Özdemir . Hepatitis-B disease modelling of fractional order and parameter calibration using real data from the USA. AIMS Biophysics, 2024, 11(3): 378-402. doi: 10.3934/biophy.2024021
    [9] Natarajan Mala, Arumugam Vinodkumar, Jehad Alzabut . Passivity analysis for Markovian jumping neutral type neural networks with leakage and mode-dependent delay. AIMS Biophysics, 2023, 10(2): 184-204. doi: 10.3934/biophy.2023012
    [10] Jean Berthier, Kenneth A. Brakke, David Gosselin, Maxime Huet, Erwin Berthier . Metastable capillary filaments in rectangular cross-section open microchannels. AIMS Biophysics, 2014, 1(1): 31-48. doi: 10.3934/biophy.2014.1.31
  • The current manuscript investigates a model of the spread of polio under the condition of vaccination by using the novel modified Atangana-Baleanu-Caputo (mABC) fractional derivative. This problem has been studied for non-zero solutions under the modified operator. The series-type solution has been obtained through the application of a Laplace transformation, along with the decomposition technique and Adomian polynomial for the nonlinear terms. The qualitative analysis for the solution of the model has been tested by using fixed point theory. The stability of the solution is also crucial for a dynamical system; therefore, it was checked by using the T-Picard method. With the help of the approximate scheme, numerical simulations were conducted for the proposed model by using different fractional orders and transmission parameters. Based on the obtained positivity of the solutions and numerical stability, we have established the analysis of the mABC operator in the field of fractional calculus and other physical sciences.



    In the middle of the 20th century, the disease of Poliomyelitis (polio) was one of the contagious diseases which has been transmitted from the western regions of the Asian continent. Jonas Salk, in the year 1952 developed a vaccination process for this outbreak. Poliovirus became the cause of this outbreak among most societies. It affects the nerves of the human body and may cause whole-body paralysis within a short duration. The associated disease mainly affects children in the first three years of age but it may also affect anyone at any age. The poliovirus is transmitted to the children's bodies through the mouth, eventually reaching the intestines of the digestive system. The poliovirus can be transferred from one person to the other. The starting signs of the associated infection are fever, headache, fatigue, vomiting, and pain produced in the limbs. There is still no effective treatment for the virus, but its spread can be controlled through the process of vaccination and immunizations which commonly save many children's lives.

    Mostly, assumptions regarding dynamic problems of daily life are necessary to manage these problems; therefore, in every mathematical model the infections are only transmitted to the susceptible or healthy class due to the infective population. Several diseases like polio are contagious during the time of incubation. So one can conclude that the interactions between the healthy population classes exposed classes have their importance in the spread of diseases. Thus, scholars have discussed the relationship between the susceptible and the exposed classes. Many scientists have attempted to construct mathematical models of Polio. The authors [1] developed a mathematical system spread of polio after the re-introduction of poliovirus. Garfinkel and Sarewitz et al. [2] briefly studied the conditions that would lead to polio eradication. Furthermore, to prove efficacy of polio vaccines in the control of polio infection, most scholars have constructed various models for the spread of polio [3][5]. The authors of [6] specifically considered the SIR problem with pulse vaccination and have shown that pulse vaccination can lead to eradication under specific conditions related to the magnitude of vaccination proportion. In some of the articles like [7], the authors dealt with a two-dimensional SIS polio system with the addition of vaccination terms for bifurcation in the opposite direction. The authors of [8] developed a polio mathematical system to investigate its dynamical behavior in the presence of vaccinations strategy with nonlinear incidence, which also deals with the optimality of the model that converges the vaccine quantity to the threshold value. Agarwal and Bhadauria [9] developed a polio epidemic model with vaccine quantity to investigate the effect of vaccination once a strategy is applied to the susceptible and exposable populations. In this study, the polio model was adapted from [9], and it has been divided into five subgroups, namely susceptible persons S, exposed persons E, acutely infected persons I, chronically infected persons V and recovered persons R:

    ((1.1))

    with the following initial conditions

    The parameters used in model (1.1) are described as follows: A is the rate of immigration in the society, µ is the rate of death that occurs naturally, β is the probability of disease transmission because of the infectious population, is the rate of chance of infection because of the exposed population, r is the rate of reduction in the transmission of infection by exposed class, ν is the rate of individuals moving from the susceptible to vaccinated class, ν1 is the rate at which exposed individuals are vaccinated, b is the rate at which exposed individuals move to the infective class and α is the rate of death due to the disease.

    The utilization of fractional calculus in the fractional-order model, which combines differentiation and integration, offers a more effective approach to comprehending real-world problems as compared to classical derivatives [10][16]. The concept of fractional derivatives, originally introduced by Riemann-Liouville based on the power law, has recently faced scrutiny from researchers due to the presentation of experimental results questioning the adequacy of a single fractional operator, such as the Caputo, Caputo-Fabrizio and Atangana-Baleanu operators, as a tool to describe complex phenomena in science and engineering [17][21]. These findings suggest that a broader range of fractional operators may be necessary to accurately capture the behavior of diverse systems [22][26]. The authors of [27] developed a model by using a new definition of the constant proportional Caputo operator, which describes the generalized memory effects. A new recursive algorithm has been constructed to solve certain initial value problems involving a fractional differential equation [28]. This new approach Adomian decomposition method (ADM) is based on the application of the ADM, and it involves combining the decomposition with a recurrence formula and utilizing the solutions of the generalized Abel equation. According to the researchers, a new fractional derivative can be defined by utilizing the exponential kernel [29], [30]. The utilization of non-singular kernel fractional derivatives in the modeling of epidemics offers valuable insights into the dynamics of infectious disease outbreaks, particularly when applied to trigonometric and exponential functions [31][33]. Initialization can pose challenges when working with non-singular kernel fractional derivatives, especially in cases involving non-singular kernels. Regardless of their type, equations in the following form have been observed to present similar challenges, as noted by the authors of [34]:

    xaF(x,τ)g(τ)dτ=y(x),axb.

    The condition F(a,a)0 introduces distinct limitations on differential equations involving non-singular kernels, implying that if F(x,τ) and y(x) are continuous and F(a,a)0, then y(a)0. This condition gives rise to peculiar behaviors in these equations, necessitating careful consideration during the development of models in various research fields. As mentioned in [34], this issue remains unsolved. However, the problem associated with the operator introduced in [17] has been solved by various researchers [35], [36]. To overcome the challenges associated with non-singular operators, the authors of [37] proposed a modification to the operator that involves utilizing a Mittag-Leffler kernel. They demonstrated that the resulting fractional differential equations, based on this modified operator, are easier to initialize than those based on non-singular kernels. Additionally, they showcased that the modified operator, known as the modified Atangana-Baleanu-Caputo (mABC) derivative, is capable of solving multiple fractional differential problems, which is not possible for the Atangana-Baleanu-Caputo (ABC) derivative. The mABC derivative possesses an integrable singularity at the origin. In [38], the authors developed a unique numerical method for the mABC derivative, utilizing finite differences; the method streamlines the initialization of the corresponding fractional differential equations. Researchers have also developed modified fractional difference operators by employing the mABC derivative and Mittag-Leffler kernels [39].

    We propose a memory-affected mABC model to model the spread of polio. Although the epidemic model described by Eq (1.1) utilizes classical derivatives and should be considered, our memory-affected mABC model provides additional insights and enhances the accuracy of hepatitis B spread predictions. To convert the ordinary system of equations given by Eq (1.1) into the mABC operator, we make the following transformation:

    ((1.3))

    with the following initial conditions

    The rest of the paper is arranged as follows. In Section 2, we recall some basic definitions and statements from the literature. In Section 3, we give some formulas and show the non-zero solution to the homogeneous fractional initial value problem. With the use of the Laplace Adomian decomposition method, we find the approximate solution for the mABC derivative in Section 4. In the same section, we also establish the stability and uniqueness results for the considered system. The approximate series solution is graphically presented in Section 5. We conclude our work in Section 6.

    Here, we recall basic results from the literature on fractional calculus.

    Definition 2.1. [37] Let f(t)L1(0,T) be a function; then, the mABC derivation is presented as follows:

    ((2.1))

    where Eθ is known as the Mittag-Leffler function for one parameter while Eθ,θ denotes the Mittag-Leffler for two parameters; from the above definition one may prove that $^\mathrm{mABC}\mathbf{D}_t^\vartheta\mathrm{f}=0.

    Definition 2.2. [37] Let f(t)L1(0,T) be a function, then, the integral of mABC operator is given as

    ((2.2))

    Lemma 2.2.1. For fL1(0,) and the order ϑ(0,1), we obtain

    ((2.3))

    The Laplace transform for the mABC is defined as

    ((2.4))

    We state that there is a non-zero solution to the homogeneous fractional initial value problem. To achieve that, we use the following formulas:

    ((3.1))
    ((3.2))

    Lemma 3.0.1. [37] Suppose that the fractional initial value problem is as follows:

    ((3.3))

    where 0<ϑ<1.

    (1) For Ω=M(ϑ)1ϑ, the solution is as follows:

    ((3.4))

    (2). For ΩM(ϑ)1ϑ, the solution is as follows:

    ((3.5))

    where, τϑ=Ω(1ϑ)M(ϑ).

    Proof. (1) Given that

    t0(ty)ϑ1Eϑ,ϑ(µϑ(1y)ϑ)yϑdy=Γ(1ϑ)Eϑ(µtϑ),

    for 0<t, we have

    ((3.7))
    ((3.8))
    ((3.9))
    ((3.10))
    ((3.11))

    which shows that the proof is complete.

    Applying Eq.(3.1) and Eq.(3.2) for t>0, we have

    Thus, the proof is complete.

    Lemma 3.0.2. [37] Assume the following fractional differential equations

    where 0<ϑ<1 and ΩM(ϑ)1ϑ; the solution of the above initial value system is given by

    ((3.12))

    where

    ((3.13))

    and zϑ=M(ϑ+Ω(1ϑ)).

    Proof. Utilizing Eq (3.1) and (3.2) we may verify that

    ((3.14))

    by Eq (2.4), we obtain

    ((3.15))

    Further calculation then yields the following:

    ((3.16))

    By substituting Eq (3.14) into Eq. (3.16), the following results are obtained:

    ((3.17))
    ((3.18))

    This finishes the proof.      �

    Remark 1. Suppose that GiC[0,T], with i=1,2,3,4,; then, we have

    ((3.19))

    subject to the following conditions:

    ((3.20))

    Therefore, ˜S=S0,˜E=E0,˜I=I0 and ˜V=V0; hence, the soultion given by Eq (3.12) presents the condition of continuity. On the basis of the above conditions the solution exists for the proposed model.

    Applying the Laplace transform to the mABC, as given in [37] to the considered model (1.3) gives

    ((4.1))

    Proceeding we get

    ((4.2))

    By application of the initial conditions we get

    ((4.3))

    Decomposing each quantity as yields

    We write the non-linear term as follows:

    where

    Gi=1Γ(i+1)didpi[k=0pkSjVj]|p=0.

    Plugging all of the above values into Eq (4.3) we get

    ((4.4))

    after calculating the above Eq. (4.4) we get

    ((4.5))
    ((4.6))
    ((4.7))
    ((4.8))

    Considering the first three terms and by applying the inverse Laplace transform we get

    ((4.9))
    ((4.10))
    ((4.11))
    ((4.12))

    By putting the values S1, E1, I1, V1 into the second term of each quantity, we get

    ((4.13))
    ((4.14))
    ((4.15))
    ((4.16))

    Lastly, each quantity can be written as follows:

    S(t)=S0+S1(t)+S2(t)\hdots
    E(t)=E0+E1(t)+E2(t)\hdots
    I(t)=I0+I1(t)+I2(t)\hdots
    V(t)=V0+V1(t)+V2(t)\hdots

    Theorem 1. Suppose that a Banach space is denoted by (B,|.|) with a mapping T:BB which satisfies

    TxTyΘTxTy+πxy

    for every x,yB; also, 0Θand0π<1. Then, T is Picard T-stable.

    Theorem 2. Suppose T to be a self map defined as given below:

    ((4.17))

    Thus, the iteration is T-stable in L1(x,y), if we obtained the following results:

    (1a(N1+N4)w1(θ)(η2+µ)w2(θ))<1,(1+a(N1+N4)w1(θ)(α+µ)w3(θ))<1,(1+αw4(θ)(µ+a1+a2)w5(θ))<1,(1+a2w5(θ)(η1+a3+µ)w6(θ))<1,

    Proof. We need to prove that T has a fixed point; for this we use (q,p)N×N

    ((4.19))

    Taking the first equation of Eq (4.19) and calculating the norm on both sides, we have

    ((4.20))

    with the help of the triangular inequality, and by solving Eq (4.20), we obtain

    ((4.21))

    By establishing the relation as follows:

    ((4.22))

    replacing the above relation in Eq (4.21), we get the following relation

    ((4.23))

    Furthermore, the convergent sequence Vq is bounded.

    Additionally, one may obtain different constants N1, N2, N3 and N4 for every t such that

    ((4.24))

    Therefore, considering Eq (4.23) and Eq (4.24), we obtain

    ((4.25))

    here, w1 and w2 are the functions of L1[(1ϑ)(sϑ+μϑ)M(ϑ)sϑL]. Considering the same procedure for the remaining equations, we have

    ((4.26))

    Hence, T has a fixed point. Applying Eq (4.25) and (4.26), we assume that

    ψ=(0,0,0,0),

    Ψ={(1a(N1+N4)w1(θ)(η2+µ)w2(θ)),(1+a(N1+N4)w1(θ)(α+µ)w3(θ)),(1+αw4(θ)(µ+a1+a2)w5(θ)),(1+a2w5(θ)(η1+a3+µ)w6(θ)),

    Therefore, the conditions of Theorem 1 are satisfied; hence, the proof is complete.

    Theorem 3. The aforementioned method gives a unique solution for the considered model.   �

    The graphical representations for the proposed model under the modified operator have been established by using the data taken from [9]. Four studied compartments of the model were tested on different fractional orders by using the obtained numerical scheme. The initial values and corresponding parameter descriptions are given in Table 1.

    Table 1.  Parameters and their numerical values for model (1.3).
    Notation Value Source Notation Value Source
    S0 500 [9] E0 200 [9]
    I0 150 [9] E0 600 [9]
    A 1000 [9] β 0.002, 0.003, 0.0016, 0.05 [9]
    r 0.5 [9] µ 0.5 [9]
    ν 0.6 [9] ν1 0.001 [9]
    b 0.9 [9] α 0.6 [9]

     | Show Table
    DownLoad: CSV
    Figure 1.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.002.

    In Figure 1(a–d), the dynamics of all four compartments are shown for different fractional orders under the modified fractional operator which has extra terms as compared to the ABC operator. The susceptible population gains stability after increasing to some peak value and the stability is best achieved for relatively small fractional orders instead of the higher order. The exposed class declines and then increases to the equilibrium point to achieve stability. During this period, the exposed population is transferred to the infected population, whose density is increasing and then stabilizes, or decreases as a result of individuals recovering from the infection after the vaccination is applied to the population. The vaccination process is terminated once the infection is controlled.

    Figure 2.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.003.

    In Figure 2(a–d), the dynamics of all four compartments are shown for different arbitrary orders under the modified fractional operator with extra informative terms. The susceptible class gains stability after declining to some small value and the stability is achieved faster with small fractional orders. The exposed class here declines quickly and then increases to the equilibrium point to achieve the convergence. During this period, the susceptible population is transferred to the exposed class, and then to the infected class whose density increases and then stabilizes, or decreases as a result of recovering from the infection after the vaccination is injected into the total population. The vaccination process is reduced as the transmission is increased and the infection is controlled after more populations.

    Figure 3.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.0016.

    In Figure 3(a–d), the rate of transmission is decreased which affects the dynamics of all four compartments as shown for different arbitrary orders. The susceptible class increases and gains stability after reaching some high value. The exposed class here declines much more quickly and then converges to its equilibrium point. In this case, after a small increase, the infection is controlled quickly. Here the vaccination rate also increases, i.e. as the vaccination rate increases the infection may be easier to control.

    Figure 4.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.0002.

    In Figure 4(a–d), the rate of transmission is again decreased which affects the dynamics of all the four compartments, as shown for different arbitrary orders. In this case, the dynamics are about the same as observed for the previous dynamics but the dynamics can be controlled more quickly. Such dynamics provide the total density of each compartment in the form of a continuous spectrum lying between 0 and 1.

    This study developed the analysis of the impact of vaccination on the spread of polio in the human population by using a four-compartment system with the generalized novel mABC fractional operator. The dynamical behaviors of all compartments have been successfully examined for different fractional orders to demonstrate the validity of an extra degree of freedom in the selection of the derivative order. The fixed-point theory shows the uniqueness of the solution in the generalized novel modified format. The stability analysis of the model has been conducted by T-Picard-type stability techniques. The numerical solution for the model has been achieved by using the Laplace transform, along with the Adomian technique in the format of the decomposition process. Such a study involved the analysis of the complex geometry in the dynamical system. The study was performed by using various fractional orders and iterations with different transmission rate; also, each curve has been plotted for six different fractional orders and compared with the integer-order result. The study also shows that the impact of vaccination is more significant on the total population. In terms of stability, relatively smaller fractional orders yield high accuracy. In all of the numerical simulations, decreasing the transmission rate led to better control of the polio infection.

    The authors declare that they have not used artificial intelligence tools in the creation of this article.



    Conflict of interest



    The authors declare no conflict of interest.

    [1] Tebbens RJD, Pallansch MA, Kew OM, et al. (2005) A dynamic model of poliomyelitis outbreaks: learning from the past to help inform the future. Am J Epidemiol 162: 358-372. https://doi.org/10.1093/aje/kwi206
    [2] Garfinkel MS, Sarewitz D (2003) Parallel path: poliovirus research in the vaccine era. Sci Eng Ethics 9: 319-338. https://doi.org/10.1007/s11948-003-0028-7
    [3] Anderson RM, May RM (1983) Vaccination against rubella and measles: quantitative investigations of different policies. J Hyg 90: 259-325. https://doi.org/10.1017/S002217240002893X
    [4] Farrington CP (2003) On vaccine efficacy and reproduction numbers. Math Biosci 185: 89-109. https://doi.org/10.1016/s0025-5564(03)00061-0
    [5] Struchiner CJ, Halloran ME, Spielman A (1989) Modeling malaria vaccines I: new uses for old ideas. Math Biosci 94: 87-113. https://doi.org/10.1016/0025-5564(89)90073-4
    [6] Shulgin B, Stone L, Agur Z (1998) Pulse vaccination strategy in the SIR epidemic model. B Math Biol 60: 1123-1148. https://doi.org/10.1016/S0092-8240(98)90005-2
    [7] Kribs-Zaleta CM, Velasco-Hernández JX (2000) A simple vaccination model with multiple endemic states. Math Biosci 164: 183-201. https://doi.org/10.1016/S0025-5564(00)00003-1
    [8] Gumel AB, Moghadas SM (2003) A qualitative study of a vaccination model with non-linear incidence. Appl Math Comput 143: 409-419. https://doi.org/10.1016/S0096-3003(02)00372-7
    [9] Agarwal M, Bhadauria AS (2011) Modeling spread of polio with the role of vaccination. Appl Appl Math 6: 11. https://digitalcommons.pvamu.edu/aam/vol6/iss2/11
    [10] Podlubny I Fractional Differential Equations: an Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Elsevier (1999). https://doi.org/10.1016/s0076-5392(99)x8001-5
    [11] Öztürk BA (2023) Examination of Sturm-Liouville problem with proportional derivative in control theory. Math Model Numer Simul Appl 3: 335-350. http://dx.doi.org/10.53391/mmnsa.1392796
    [12] Abbas M, Aslam S, Abdullah FA, et al. (2023) An efficient spline technique for solving time-fractional integro-differential equations. Heliyon 9: e19307. https://doi.org/10.1016/j.heliyon.2023.e19307
    [13] Awawdeh F, Adawi A, Mustafa Z (2009) Solutions of the SIR models of epidemics using HAM. Chaos Soliton Fract 42: 3047-3052. https://doi.org/10.1016/j.chaos.2009.04.012
    [14] Joshi H, Yavuz M, Stamova I (2023) Analysis of the disturbance effect in intracellular calcium dynamic on fibroblast cells with an exponential kernel law. Bull Biomath 1: 24-39. https://doi.org/10.59292/bulletinbiomath.2023002
    [15] Zhu XH, Xia PF, He QZ, et al. (2023) Coke price prediction approach based on dense GRU and opposition-based learning salp swarm algorithm. Int J Bio-Inspir Comp 21: 106-121. https://doi.org/10.1504/ijbic.2023.130549
    [16] Abdelrazec A Adomian decomposition method: convergence analysis and numerical approximations, thesis, McMaster University (2008).
    [17] Caputo M, Fabrizio M (2015) A new definition of fractional derivative without singular kernel. Prog Fract Differ Appl 1: 73-85. http://dx.doi.org/10.12785/pfda/010201
    [18] Atangana A, Baleanu D (2016) New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model. Therm Sci 20: 763-785. https://doi.org/10.2298/TSCI160111018A
    [19] Evirgen F, Uçar E, Uçar S, et al. (2023) Modelling influenza a disease dynamics under Caputo-Fabrizio fractional derivative with distinct contact rates. Math Model Numer Simul Appl 3: 58-72. http://dx.doi.org/10.53391/mmnsa.1274004
    [20] Anjam YN, Shafqat R, Sarris IE, et al. (2022) A fractional order investigation of smoking model using Caputo-Fabrizio differential operator. Fractal Fract 6: 623. https://doi.org/10.3390/fractalfract6110623
    [21] Mohamed EG, Mohamed N, Adel W (2023) Genocchi collocation method for accurate solution of nonlinear fractional differential equations with error analysis. Math Model Numer Simul Appl 3: 351-375. https://doi.org/10.53391/mmnsa.1373647
    [22] Zguaid K, El Alaoui FZ (2023) On the regional boundary observability of semilinear time-fractional systems with Caputo derivative. Int J Optim Control The Appl 13: 161-170. http://dx.doi.org/10.11121/ijocta.2023.1286
    [23] Li B, Zhang TX, Zhang C (2023) Investigation of financial bubble mathematical model under fractal-fractional Caputo derivative. Fractals 31: 1-13. https://doi.org/10.1142/S0218348X23500500
    [24] Joshi H, Yavuz M, Townley S, et al. (2023) Stability analysis of a non-singular fractional-order covid-19 model with nonlinear incidence and treatment rate. Phys Scripta 98: 045216. https://doi.org/10.1088/1402-4896/acbe7a
    [25] Işık E, Daşbaşı B (2023) A compartmental fractional-order mobbing model and the determination of its parameters. Bull Biomath 1: 153-176. http://dx.doi.org/10.59292/bulletinbiomath.2023008
    [26] Zhang X, Ding Z, Hang J, et al. (2022) How do stock price indices absorb the COVID-19 pandemic shocks?. N Am J Econ Financ 60: 101672. https://doi.org/10.1016/j.najef.2022.101672
    [27] Ur Rehman A, Hua S, Riaz MB, et al. (2023) A fractional study with Newtonian heating effect on heat absorbing MHD radiative flow of rate type fluid with application of novel hybrid fractional derivative operator. Arab J Basic Appl Sci 30: 482-495. https://doi.org/10.1080/25765299.2023.2250063
    [28] Mennouni A, Bougoffa L, Wazwaz AM (2022) A new recursive scheme for solving a fractional differential equation of ray tracing through the crystalline lens. Opt Quant Electron 54: 373. https://doi.org/10.1007/s11082-022-03766-w
    [29] Haidong Q, Ur Rahman M, Arfan M (2023) Fractional model of smoking with relapse and harmonic mean type incidence rate under Caputo operator. J Appl Math Comput 69: 403-420. https://doi.org/10.1007/s12190-022-01747-6
    [30] Khan A, Khan TS, Syam MI, et al. (2019) Analytical solutions of time-fractional wave equation by double Laplace transform method. Eur Phys J Plus 134: 163. https://doi.org/10.1140/epjp/i2019-12499-y
    [31] Li B, Eskandari Z (2023) Dynamical analysis of a discrete-time SIR epidemic model. Journal Franklin I 360: 7989-8007. https://doi.org/10.1016/j.jfranklin.2023.06.006
    [32] Qayyum M, Afzal S, Ahmad E, et al. (2023) Fractional modeling and analysis of unsteady squeezing flow of Casson nanofluid via extended He-Laplace algorithm in Liouville-Caputo sense. Alex Eng J 73: 579-591. https://doi.org/10.1016/j.aej.2023.05.010
    [33] Toufik M, Atangana A (2017) New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models. Eur Phys J Plus 132: 444. https://doi.org/10.1140/epjp/i2017-11717-0
    [34] Atangana A (2018) Non validity of index law in fractional calculus: a fractional differential operator with Markovian and non-Markovian properties. Physica A 505: 688-706. https://doi.org/10.1016/j.physa.2018.03.056
    [35] Caputo M, Fabrizio M (2021) On the singular kernels for fractional derivatives. Some applications to partial differential equations. Prog Fract Differ Appl 7: 1-4. http://dx.doi.org/10.18576/PFDA/070201
    [36] Losada J, Nieto JJ (2021) Fractional integral associated to fractional derivatives with nonsingular kernels. Prog Fract Differ Appl 7: 137-143. http://dx.doi.org/10.18576/pfda/070301
    [37] Al-Refai M, Baleanu D (2022) On an extension of the operator with Mittag-Leffler kernel. Fractals 30: 2240129. http://dx.doi.org/10.1142/S0218348X22401296
    [38] Shiri B, Baleanu D (2023) All linear fractional derivatives with power functions' convolution kernel and interpolation properties. Chaos Soliton Fract 170: 113399. https://doi.org/10.1016/j.chaos.2023.113399
    [39] Mohammed PO, Srivastava HM, Baleanu D, et al. (2022) Modified fractional difference operators defined using Mittag-Leffler kernels. Symmetry 14: 1519. https://doi.org/10.3390/sym14081519
  • This article has been cited by:

    1. Elham Raeisi, Mehmet Yavuz, Mohammadreza Khosravifarsani, Yasin Fadaei, Mathematical modeling of interactions between colon cancer and immune system with a deep learning algorithm, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05111-4
    2. Elhoussine Azroul, Sara Bouda, Mathematical modeling of poliomyelitis virus with vaccination and post‐paralytic syndrome dynamics using Caputo and ABC fractional derivatives, 2024, 0170-4214, 10.1002/mma.10406
    3. Mohamed A. Abdoon, Abdulrahman B. M. Alzahrani, Comparative Analysis of Influenza Modeling Using Novel Fractional Operators with Real Data, 2024, 16, 2073-8994, 1126, 10.3390/sym16091126
    4. Elkhateeb S. Aly, Mohammed A. Almalahi, Khaled A. Aldwoah, Kamal Shah, Criteria of existence and stability of an n-coupled system of generalized Sturm-Liouville equations with a modified ABC fractional derivative and an application to the SEIR influenza epidemic model, 2024, 9, 2473-6988, 14228, 10.3934/math.2024691
    5. Lahcen Boulaasair, Hassane Bouzahir, Mehmet Yavuz, Global mathematical analysis of a patchy epidemic model, 2024, 14, 2146-5703, 365, 10.11121/ijocta.1558
    6. Mehmet Yavuz, Feyza Nur Özdemir, Kübra Akyüz, Naime Büşra Bayraktar, The relationship between colon cancer and immune system: a fractional order modelling approach, 2025, 27, 1301-7985, 126, 10.25092/baunfbed.1498218
    7. Isah Abdullahi, D.G. Yakubu, M.Y. Adamu, Musa Ali, A.M. Kwami, Inclined magnetic fields heat transfer and thermal radiation on fractionalized EMHD Burgers’ fluid flow via bifurcated artery for tumor treatments, 2025, 26668181, 101093, 10.1016/j.padiff.2025.101093
    8. R. Rajaraman, Modeling Immobilized Enzyme Reactions: Nonlinear Kinetics With Fractional‐ and Integer‐Order Analysis, 2025, 0170-4214, 10.1002/mma.10791
    9. Ateq Alsaadi, Advancing water quality management: A synergistic approach using fractional differential equations and neural networks, 2025, 10, 2473-6988, 5332, 10.3934/math.2025246
  • Reader Comments
  • © 2024 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(1277) PDF downloads(78) Cited by(9)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog