Research article Special Issues

Investigating of an immune system-cancer mathematical model with Mittag-Leffler kernel

  • Cancer that is difficult to treat, is a very common disease today and there are many types of cancer such as lung, colon, stomach. When cancer settles in the body, the immune system tries to resist it. In this study, the mathematical model of the interaction between immune system components and cancer is discussed and is modified by using Atangana-Baleanu derivative. After investigating the existence and uniqueness of the solution of the fractional immune system-cancer model, numerical simulations are given via predictor-corrector scheme.

    Citation: Necati Özdemir, Esmehan Uçar. Investigating of an immune system-cancer mathematical model with Mittag-Leffler kernel[J]. AIMS Mathematics, 2020, 5(2): 1519-1531. doi: 10.3934/math.2020104

    Related Papers:

    [1] Ridha Dida, Hamid Boulares, Bahaaeldin Abdalla, Manar A. Alqudah, Thabet Abdeljawad . On positive solutions of fractional pantograph equations within function-dependent kernel Caputo derivatives. AIMS Mathematics, 2023, 8(10): 23032-23045. doi: 10.3934/math.20231172
    [2] Saeed M. Ali, Mohammed S. Abdo, Bhausaheb Sontakke, Kamal Shah, Thabet Abdeljawad . New results on a coupled system for second-order pantograph equations with ABC fractional derivatives. AIMS Mathematics, 2022, 7(10): 19520-19538. doi: 10.3934/math.20221071
    [3] Abdelkader Moumen, Hamid Boulares, Tariq Alraqad, Hicham Saber, Ekram E. Ali . Newly existence of solutions for pantograph a semipositone in Ψ-Caputo sense. AIMS Mathematics, 2023, 8(6): 12830-12840. doi: 10.3934/math.2023646
    [4] Ayub Samadi, Chaiyod Kamthorncharoen, Sotiris K. Ntouyas, Jessada Tariboon . Mixed Erdélyi-Kober and Caputo fractional differential equations with nonlocal non-separated boundary conditions. AIMS Mathematics, 2024, 9(11): 32904-32920. doi: 10.3934/math.20241574
    [5] Hui Huang, Kaihong Zhao, Xiuduo Liu . On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Mathematics, 2022, 7(10): 19221-19236. doi: 10.3934/math.20221055
    [6] Cuiying Li, Rui Wu, Ranzhuo Ma . Existence of solutions for Caputo fractional iterative equations under several boundary value conditions. AIMS Mathematics, 2023, 8(1): 317-339. doi: 10.3934/math.2023015
    [7] Mohamed Houas, Kirti Kaushik, Anoop Kumar, Aziz Khan, Thabet Abdeljawad . Existence and stability results of pantograph equation with three sequential fractional derivatives. AIMS Mathematics, 2023, 8(3): 5216-5232. doi: 10.3934/math.2023262
    [8] Ahmed M. A. El-Sayed, Wagdy G. El-Sayed, Kheria M. O. Msaik, Hanaa R. Ebead . Riemann-Liouville fractional-order pantograph differential equation constrained by nonlocal and weighted pantograph integral equations. AIMS Mathematics, 2025, 10(3): 4970-4991. doi: 10.3934/math.2025228
    [9] Isra Al-Shbeil, Abdelkader Benali, Houari Bouzid, Najla Aloraini . Existence of solutions for multi-point nonlinear differential system equations of fractional orders with integral boundary conditions. AIMS Mathematics, 2022, 7(10): 18142-18157. doi: 10.3934/math.2022998
    [10] Yujun Cui, Chunyu Liang, Yumei Zou . Existence and uniqueness of solutions for a class of fractional differential equation with lower-order derivative dependence. AIMS Mathematics, 2025, 10(2): 3797-3818. doi: 10.3934/math.2025176
  • Cancer that is difficult to treat, is a very common disease today and there are many types of cancer such as lung, colon, stomach. When cancer settles in the body, the immune system tries to resist it. In this study, the mathematical model of the interaction between immune system components and cancer is discussed and is modified by using Atangana-Baleanu derivative. After investigating the existence and uniqueness of the solution of the fractional immune system-cancer model, numerical simulations are given via predictor-corrector scheme.


    The key to solving the general quadratic congruence equation is to solve the equation of the form x2amodp, where a and p are integers, p>0 and p is not divisible by a. For relatively large p, it is impractical to use the Euler criterion to distinguish whether the integer a with (a,p)=1 is quadratic residue of modulo p. In order to study this issue, Legendre has proposed a new tool-Legendre's symbol.

    Let p be an odd prime, the quadratic character modulo p is called the Legendre's symbol, which is defined as follows:

    (ap)={1, if a is a quadratic residue modulo p;1, if a is a quadratic non-residue modulo p;0, if pa.

    The Legendre's symbol makes it easy for us to calculate the level of quadratic residues. The basic properties of Legendre's symbol can be found in any book on elementary number theory, such as [1,2,3].

    The properties of Legendre's symbol and quadratic residues play an important role in number theory. Many scholars have studied them and achieved some important results. For examples, see the [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21].

    One of the most representative properties of the Legendre's symbol is the quadratic reciprocal law:

    Let p and q be two distinct odd primes. Then, (see Theorem 9.8 in [1] or Theorems 4–6 in [3])

    (pq)(qp)=(1)(p1)(q1)4.

    For any odd prime p with p1mod4 there exist two non-zero integers α(p) and β(p) such that

    p=α2(p)+β2(p). (1)

    In fact, the integers α(p) and β(p) in the (1) can be expressed in terms of Legendre's symbol modulo p (see Theorems 4–11 in [3])

    α(p)=12p1a=1(a3+ap)andβ(p)=12p1a=1(a3+rap),

    where r is any integer, and (r,p)=1, (rp)=1, (p)=χ2 denote the Legendre's symbol modulo p.

    Noting that Legendre's symbol is a special kind of character. For research on character, Han [7] studied the sum of a special character χ(ma+ˉa), for any integer m with (m,p)=1, then

    |p1a=1χ(ma+ˉa)|2=2p+(mp)p1a=1χ(a)p1b=1(b(b1)(a2b1)p),

    which is a special case of a general polynomial character sums N+Ma=N+1χ(f(a)), where M and N are any positive integers, and f(x) is a polynomial.

    In [8], Du and Li introduced a special character sums C(χ,m,n,c;p) in the following form:

    C(χ,m,n,c;p)=p1a=0p1b=0χ(a2+nab2nb+c)e(mb2ma2p),

    and studied the asymptotic properties of it. They obtained

    p1c=1|C(χ,m,n,c;p)|2k={p2k+1+k23k22p2k+O(p2k1),ifχ is the Legendre symbol modulo p;p2k+1+k23k22p2k+O(p2k1/2),ifχ is a complex character modulo p.

    Recently, Yuan and Zhang [12] researched the question about the estimation of the mean value of high-powers for a special character sum modulo a prime, let p be an odd prime with p1mod6, then for any integer k0, they have the identity

    Sk(p)=13[dk+(d+9b2)k+(d9b2)k],

    where

    Sk(p)=1p1p1r=1Ak(r),
    A(r)=1+p1a=1(a2+rˉap),

    and for any integer r with (r,p)=1.

    More relevant research on special character sums will not be repeated. Inspired by these papers, we have the question: If we replace the special character sums with Legendre's symbol, can we get good results on p1mod4?

    We will convert β(p) to another form based on the properties of complete residues

    β(p)=12p1a=1(a+nˉap),

    where ˉa is the inverse of a modulo p. That is, ˉa satisfy the equation xa1modp for any integer a with (a,p)=1.

    For any integer k0, G(n) and Kk(p) are defined as follows:

    G(n)=1+p1a=1(a2+nˉa2p)andKk(p)=1p1p1n=1Gk(n).

    In this paper, we will use the analytic methods and properties of the classical Gauss sums and Dirichlet character sums to study the computational problem of Kk(p) for any positive integer k, and give a linear recurrence formulas for Kk(p). That is, we will prove the following result.

    Theorem 1. Let p be an odd prime with p1mod4, then we have

    Kk(p)=(4p+2)Kk2(p)8(2α2p)Kk3(p)+(16α416pα2+4p1)Kk4(p),

    for all integer k4 with

    K0(p)=1,K1(p)=0,K2(p)=2p+1,K3(p)=3(4α22p),

    where

    α=α(p)=p12a=1(a+ˉap).

    Applying the properties of the linear recurrence sequence, we may immediately deduce the following corollaries.

    Corollary 1. Let p be an odd prime with p1mod4. Then we have

    1p1p1n=111+p1a=1(a2+nˉa2p)=16α2p28α28p2+14p16α416α2p+4p1.

    Corollary 2. Let p be an odd prime with p1mod4. Then we have

    1p1p1n=1p1m=0(1+p1a=1(a2+nˉa2p))e(nm2p)=p.

    Corollary 3. Let p be an odd prime with p1mod4. Then we have

    1p1p1n=1p1m=0[1+p1a=1(a2+nˉa2p)]2e(nm2p)=(4α22p)p.

    Corollary 4. Let p be an odd prime with p1mod8. Then we have

    p1n=1(1+p1a=1(a2+nˉa2p))p1m=0e(nm4p)=p(1+B(1))p,

    where

    B(1)=p1m=0e(m4p).

    If we consider such a sequence Fk(p) as follows: Let p be a prime with p1mod8, χ4 be any fourth-order character modulo p. For any integer k0, we define the Fk(p) as

    Fk(p)=p1n=11Gk(n),

    we have

    Fk(p)=116α416α2p+4p1Fk4(p)(4p+2)16α416α2p+4p1Fk2(p)+4(4α22p)16α416α2p+4p1Fk1(p).

    Lemma 1. Let p be an odd prime with p1mod4. Then for any fourth-order character χ4modp, we have the identity

    τ2(χ4)+τ2(¯χ4)=2pα,

    where

    τ(χ4)=p1a=1χ4(a)e(ap)

    denotes the classical Gauss sums, e(y)=e2πiy,i2=1, and α is the same as in the Theorem 1.

    Proof. See Lemma 2.2 in [9].

    Lemma 2. Let p be an odd prime. Then for any non-principal character ψ modulo p, we have the identity

    τ(ψ2)=ψ2(2)τ(χ2)τ(ψ)τ(ψχ2),

    where χ2=(p) denotes the Legendre's symbol modulo p.

    Proof. See Lemma 2 in [12].

    Lemma 3. Let p be a prime with p1mod4, then for any integer n with (n,p)=1 and fourth-order character χ4modp, we have the identity

    p1a=1(a2+nˉa2p)=1χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)).

    Proof. For any integer a with (a,p)=1, we have the identity

    1+χ4(a)+χ2(a)+¯χ4(a)=4,

    if a satisfies ab4modp for some integer b with (b,p)=1 and

    1+χ4(a)+χ2(a)+¯χ4(a)=0,

    otherwise. So from these and the properties of Gauss sums we have

    p1a=1(a2+nˉa2p)=p1a=1(a2p)(a4+np)=p1a=1χ2(a4)χ2(a4+n)=p1a=1(1+χ4(a)+χ2(a)+¯χ4(a))χ2(a)χ2(a+n)=p1a=1(1+χ4(na)+χ2(na)+¯χ4(na))χ2(na)χ2(na+n)=p1a=1χ2(a)χ2(a+1)+p1a=1χ4(na)χ2(a)χ2(a+1) (2)
    +p1a=1χ2(na)χ2(a)χ2(a+1)+p1a=1¯χ4(na)χ2(a)χ2(a+1)=p1a=1χ2(1+ˉa)+p1a=1χ4(na)χ2(a)χ2(a+1)+p1a=1χ2(n)χ2(a+1)+p1a=1¯χ4(na)χ2(a)χ2(a+1).

    Noting that for any non-principal character χ,

    p1a=1χ(a)=0

    and

    p1a=1χ(a)χ(a+1)=1τ(ˉχ)p1b=1p1a=1ˉχ(b)χ(a)e(b(a+1)p).

    Then we have

    p1a=1χ2(1+ˉa)=1,p1a=1χ2(a+1)=1,
    p1a=1χ4(a)χ2(a)χ2(a+1)=1τ(χ2)p1b=1p1a=1χ2(b)χ4(a)χ2(a)e(b(a+1)p)=1τ(χ2)p1b=1¯χ4(b)e(bp)p1a=1χ4(ab)χ2(ab)e(abp) (3)
    =1τ(χ2)τ(¯χ4)τ(χ4χ2).

    For any non-principal character ψ, from Lemma 2 we have

    τ(ψ2)=ψ2(2)τ(χ2)τ(ψ)τ(ψχ2). (4)

    Taking ψ=χ4, note that

    τ(χ2)=p,  τ(χ4)τ(¯χ4)=χ4(1)p,

    from (3) and (4), we have

    p1a=1χ4(a)χ2(a)χ2(a+1)=¯χ42(2)τ(χ24)τ(χ2)τ(¯χ4)τ(χ2)τ(χ4)=χ2(2)τ(χ2)τ2(¯χ4)τ(χ4)τ(¯χ4)=χ2(2)pτ2(¯χ4)χ4(1)p (5)
    =χ2(2)τ2(¯χ4)χ4(1)p.

    Similarly, we also have

    p1a=1¯χ4(a)χ2(a)χ2(a+1)=χ2(2)τ2(χ4)χ4(1)p. (6)

    Consider the quadratic character modulo p, we have

    (2p)=χ2(2)={1,if p±1mod8;1,if p±3mod8. (7)

    And when p1mod8, we have χ4(1)=1; when p5mod8, we have χ4(1)=1. Combining (2) and (5)–(7) we can deduce that

    p1a=1(a2+nˉa2p)=1χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)).

    This prove Lemma 3.

    Lemma 4. Let p be an odd prime with p1mod4. Then for any integer k4 and n with (n,p)=1, we have the fourth-order linear recurrence formula

    Gk(n)=(4p+2)Gk2(n)+8(p2α2)Gk3(n)+[(4α22p)2(2p1)2]Gk4(n),

    where

    α=α(p)=12p1a=1(a3+ap)=p12a=1(a+ˉap),

    (p)=χ2 denotes the Legendre's symbol.

    Proof. For p1mod4, any integer n with (n,p)=1, and fourth-order character χ4 modulo p, we have the identity

    χ44(n)=¯χ44(n)=χ0(n),  χ24(n)=χ2(n),

    where χ0 denotes the principal character modulo p.

    According to Lemma 3,

    p1a=1(a2+nˉa2p)=1χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)),
    G(n)=1+p1a=1(a2+nˉa2p).  

    We have

    G(n)=χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)), (8)
    G2(n)=[χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4))]2=12χ2(n)1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4))+1p(χ2(n)τ4(¯χ4)+χ2(n)τ4(χ4)+2p2)=12χ2(n)1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4))+1p(χ2(n)(τ4(¯χ4)+τ4(χ4))+2p2).

    According to Lemma 1, we have

    (τ2(χ4)+τ2(¯χ4))2=τ4(¯χ4)+τ4(χ4)+2p2=4pα2.

    Therefore, we may immediately deduce

    G2(n)=12(χ2(n)(G(n)+χ2(n))+1p(χ2(n)(τ4(¯χ4)+τ4(χ4))+2p2)=12χ2(n)(G(n)+χ2(n)) (9)
    +1p[χ2(n)((τ2(¯χ4)+τ2(χ4))22p2)+2p2]=2p12χ2(n)G(n)+(4α22p)χ2(n),
    G3(n)=[χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4))]3=(2p12χ2(n)G(n)+(4α22p)χ2(n))G(n) (10)
    =(4α22p)χ2(n)G(n)+(2p+3)G(n)(4p2)χ2(n)2(4α22p)

    and

    [G2(n)(2p1)]2=[χ2(n)(4α22p)2χ2(n)G(n)]2,

    which implies that

    G4(n)=(4p+2)G2(n)+8(p2α2)G(n)+[(4α22p)2(2p1)2]. (11)

    So for any integer k4, from (8)–(11), we have the fourth-order linear recurrence formula

    Gk(n)=Gk4(n)G4(n)=(4p+2)Gk2(n)+8(p2α2)Gk3(n)+[(4α22p)2(2p1)2]Gk4(n).

    This proves Lemma 4.

    In this section, we will complete the proof of our theorem.

    Let p be any prime with p1mod4, then we have

    K0(p)=1p1p1n=1G0(n)=p1p1=1. (12)
    K1(p)=1p1p1n=1G1(n)=1p1p1n=1(χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)))=0, (13)
    K2(p)=1p1p1n=1G2(n)=1p1p1n=1(χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)))2=2p+1, (14)
    K3(p)=1p1p1n=1G3(n)=1p1p1n=1(χ2(n)+1p(χ4(n)τ2(¯χ4)+¯χ4(n)τ2(χ4)))3=3(4α22p). (15)

    It is clear that from Lemma 4, if k4, we have

    Kk(p)=1p1p1n=1Gk(n)=(4p+2)Kk2(p)8(2α2p)Kk3(p)+(16α416pα2+4p1)Kk4(p). (16)

    Now Theorem 1 follows (12)–(16). Obviously, using Theorem 1 to all negative integers, and that lead to Corollary 1.

    This completes the proofs of our all results.

    Some notes:

    Note 1: In our theorem, know n is an integer, and (n,p)=1. According to the properties of quadratic residual, χ2(n)=±1, χ4(n)=±1.

    Note 2: In our theorem, we only discussed the case p1mod8. If p3mod4, then the result is trivial. In fact, in this case, for any integer n with (n,p)=1, we have the identity

    G(n)=1+p1a=1(a2+nˉa2p)=1+p1a=1(a4p)(a4+np)=1+p1a=1(ap)(a+np)=1+p1a=1(a2+nap)=1+p1a=1(1+nˉap)=p1a=0(1+nap)=0.

    Thus, for all prime p with p3mod4 and k1, we have Kk(p)=0.

    The main result of this paper is Theorem 1. It gives an interesting computational formula for Kk(p) with p1mod4. That is, for any integer k, we have the identity

    Kk(p)=(4p+2)Kk2(p)8(2α2p)Kk3(p)+(16α416pα2+4p1)Kk4(p).

    Thus, the problems of calculating a linear recurrence formula of one kind special character sums modulo a prime are given.

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

    The authors are grateful to the anonymous referee for very helpful and detailed comments.

    This work is supported by the N.S.F. (11971381, 12371007) of China and Shaanxi Fundamental Science Research Project for Mathematics and Physics (22JSY007).

    The authors declare no conflicts of interest.



    [1] L. Marsha, K. R. Conroy, J. L. Davis, et al. Atlas Pathophysiology, Lippincott Williams & Wilkins, 2010.
    [2] V. Kumar, A. Abbas, J. Aster, Robbins and cotran pathologic basis of disease, Canada: Elsevier, 2014.
    [3] M. Zanetti, Tapping CD4 T cells for cancer immunotherapy: The choice of personalized genomics, J. Immunol., 194 (2015), 2049-2056. doi: 10.4049/jimmunol.1402669
    [4] D. Cassell, J. Forman, Linked recognition of helper and cytotoxic antigenic determinants for he generation of cytotoxic T lymphocytes, Ann. N. Y. Acad. Sci., 532 (1998), 51-60.
    [5] H. Choudhry, N. Helmi, W. H. Abdulaal, et al. Prospects of IL-2 in cancer immunotherapy, BioMed Res. Int., 2018 (2018), 9056173.
    [6] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, 2006.
    [7] D. Baleanu, K. Diethelm, E. Scalas, et al. Fractional calculus models and numerical methods, World Scientific, 2012
    [8] N. Özdemir, D. Karadeniz, B. B. Iskender, Fractional optimal control problem of a distributed system in cylindrical coordinates, Phys. Lett. A, 373 (2009), 221-226. doi: 10.1016/j.physleta.2008.11.019
    [9] F. Evirgen, N. Özdemir, Multistage adomian decomposition method for solving NLP problems over a nonlinear fractional dynamical system, J. Comput. Nonlinear Dyn., 6 (2011), 21003.
    [10] F. Evirgen, Analyze the optimal solutions of optimization problems by means of fractional gradient based system using VIM, An International Journal of Optimization and Control: Theories & Applications (IJOCTA), 6 (2016), 75-83.
    [11] Z. Hammouch, T. Mekkaoui, Circuit design and simulation for the fractional-order chaotic behavior in a new dynamical system, Complex Intell. Syst., 4 (2018), 251-260. doi: 10.1007/s40747-018-0070-3
    [12] E. Bonyah, A. Atangana, M. A. Khan, Modeling the spread of computer virus via Caputo fractional derivative and the beta derivative, Asia Pacific Journal on Computational Engineering, 4 (2017), 1-15. doi: 10.1186/s40540-016-0019-1
    [13] N. Özdemir, M. Yavuz, Numerical solution of fractional Black-Scholes equation by using the multivariate pade approximation, Acta Phys. Pol. A., 132 (2017), 1050-1053. doi: 10.12693/APhysPolA.132.1050
    [14] E. Uçar, N. Özdemir, E. Altun, Fractional order model of immune cells influenced by cancer cells, Math. Model. Nat. Phenom., 14 (2019), 308.
    [15] A. Atangana, D. Baleanu, New fractional derivatives with non-local and non-singular kernel: theory and applications to heat transfer model, Therm. Sci., 20 (2016), 763-769. doi: 10.2298/TSCI160111018A
    [16] M. Yavuz, N. Özdemir, H. M. Baskonus, Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel, Eur. Phys. J. Plus, 133 (2018), 215.
    [17] V. F. Morales-Delgadoa, J. F. Gomez-Aguilar, M. A. Taneco-Hernandez, et al. Mathematical modeling of the smoking dynamics using fractional differential equations with local and nonlocal kernel, J. Nonlinear Sci. Appl., 11 (2018), 994-1014. doi: 10.22436/jnsa.011.08.06
    [18] N. A. Asif, Z. Hammouch, M. B. Riaz, et al. Analytical solution of a Maxwell fluid with slip effects in view of the Caputo-Fabrizio derivative, Eur. Phys. J. Plus, 133 (2018), 272.
    [19] I. Koca, Analysis of rubella disease model with non-local and non-singular fractional derivatives, Int. J. Optim. Control Theor. Appl. IJOCTA, 8 (2018), 17-25.
    [20] D. Avcı A. Yetim, Analytical solutions to the advection-diffusion equation with the AtanganaBaleanu derivative over a finite domain, J. BAUN Inst. Sci. Technol., 20 (2018), 382-395.
    [21] S. Uçar, E. Uçar, N. Özdemir, et al. Mathematical analysis and numerical simulation for a smoking model with Atangana-Baleanu derivative, Chaos, Solitons & Fractals, 118 (2019), 300-306.
    [22] D. Baleanu, A. Fernandez, On some new properties of fractional derivatives with Mittag-Leffler kernel, Commun. Nonlinear Sci. Numer. Simulat., 59 (2018), 444-462. doi: 10.1016/j.cnsns.2017.12.003
    [23] A. Fernandez, D. Baleanu, H. M. Srivastava, Series representations for fractional-calculus operators involving generalised Mittag-Leffler functions, Commun. Nonlinear Sci. Numer. Simulat., 67 (2019), 517-527. doi: 10.1016/j.cnsns.2018.07.035
    [24] S. Uçar, Existence and uniqueness results for a smoking model with determination and education in the frame of non-singular derivatives, Discrete Continuous Dyn. Syst. Ser. S, in press.
    [25] F. Evirgen, S. Uçar, N. Özdemir, et al. System response of an alcoholism model under the effect of immigration via non-singular kernel derivative, Discrete Continuous Dyn. Syst. Ser. S, in press.
    [26] J. E. Solis-Perez, J. F. Gomez-Aguilar, A. Atangana, A factional mathematical model of breast cancer competition model, Chaos, Solitons and Fractals, 127 (2019), 38-54. doi: 10.1016/j.chaos.2019.06.027
    [27] V. F. Morales-Delgado, J. F. Gomez-Aguilar, K. Saad, et al. Application of the Caputo-Fabrizio and Atangana-Baleanu fractional derivatives to mathematical model of cancer chemotherapy effect, Math. Methods Appl. Sci., 42 (2019), 1167-1193. doi: 10.1002/mma.5421
    [28] P. Vereesha, D. G. Prakasha, H. M. Baskonus, New numerical surfaces to the mathematical model of cancer chemotherapy effect in Caputo fractional derivatives, CHAOS, 29 (2019), 1-13.
    [29] A. Minelli, F. Topputo, F. Bernelli, Controlled drug delivery in cancer immunotherapy: Stability, optimization and monte carlo analysis, SIAM J. Appl. Math., 71 (2011), 2229-2245. doi: 10.1137/100815190
    [30] L. G. De Pillis, A. Radunskaya, A mathematical tumour model with immune resistance and drug therapy: An optimal control approach, Journal of Theoretical Medicine, 3 (2001), 79-100. doi: 10.1080/10273660108833067
    [31] F. Castiglione, B. Piccoli, Cancer immunotheraphy, mathematical modeling and optimal control, J. Theor. Biol., 247 (2007), 723-732. doi: 10.1016/j.jtbi.2007.04.003
    [32] D. Baleanu, A. Jajarmi, M. Hajipour, On the nonlinear dynamical systems within the generalized fractional derivatives with Mittag-Leffler kernel, Nonlinear Dyn., 94 (2018), 397-414. doi: 10.1007/s11071-018-4367-y
  • This article has been cited by:

    1. Mounia Mouy, Hamid Boulares, Saleh Alshammari, Mohammad Alshammari, Yamina Laskri, Wael W. Mohammed, On Averaging Principle for Caputo–Hadamard Fractional Stochastic Differential Pantograph Equation, 2022, 7, 2504-3110, 31, 10.3390/fractalfract7010031
    2. Sabbavarapu Nageswara Rao, Manoj Singh, Ahmed Hussein Msmali, Abdullah Ali H. Ahmadini, Existence of Positive Solutions for a Coupled System of p-Laplacian Semipositone Hadmard Fractional BVP, 2023, 7, 2504-3110, 499, 10.3390/fractalfract7070499
  • 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(4999) PDF downloads(593) Cited by(35)

Figures and Tables

Figures(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog