images Computer Modeling in
Engineering & Sciences

DOI: 10.32604/cmes.2021.014988


Regarding on the Fractional Mathematical Model of Tumour Invasion and Metastasis

P. Veeresha1, Esin Ilhan2, D. G. Prakasha3, Haci Mehmet Baskonus4 and Wei Gao5,*

1Department of Mathematics, CHRIST (Deemed to be University), Bengaluru, 560029, India
2Faculty of Engineering and Architecture, Kirsehir Ahi Evran University, Kirsehir, 40500, Turkey
3Department of Mathematics, Faculty of Science, Davangere University, Shivagangothri, Davangere, 577007, India
4Department of Mathematics and Science Education, Faculty of Education, Harran University, Sanliurfa, 63050, Turkey
5School of Information Science and Technology, Yunnan Normal University, Kunming, 650500, China
*Corresponding Author: Wei Gao. Email: gaowei@ynnu.edu.cn
Received: 14 November 2020; Accepted: 09 February 2021

Abstract: In this paper, we analyze the behaviour of solution for the system exemplifying model of tumour invasion and metastasis by the help of q-homotopy analysis transform method (q-HATM) with the fractional operator. The analyzed model consists of a system of three nonlinear differential equations elucidating the activation and the migratory response of the degradation of the matrix, tumour cells and production of degradative enzymes by thetumour cells. The considered method is graceful amalgamations of q-homotopy analysis technique with Laplace transform (LT), and Caputo–Fabrizio (CF) fractional operator is hired in the present study. By using the fixed point theory, existence and uniqueness are demonstrated. To validate and present the effectiveness of the considered algorithm, we analyzed the considered system in terms of fractional order with time and space. The error analysis of the considered scheme is illustrated. The variations with small change time with respect to achieved results are effectively captured in plots. The obtained results confirm that the considered method is very efficient and highly methodical to analyze the behaviors of the system of fractional order differential equations.

Keywords: Tumour cell; invasion and metastasis; q-homotopy analysis transform method; Caputo–Fabrizio derivative

1  Introduction

The existence of cancer among the people around the globe has been developed swiftly, and globally it becomes the second foremost cause of demise after cardiovascular diseases [1]. The process of spreading and formation of secondary tumours is known as Metastasis and this nature of cancer cells is the main reason for the death in cancer patients. In addition, the prediction of size, stage, and evolution of a tumour is very critical for the treatment of cancer. Moreover, mathematics plays as an essential tool and aid us to analyze the behaviour of the tumour. The tumour growth has been mathematically modelled by the number of researchers and which are appeared in the literature [24]. The growth of the tumour invasion and metastasis are described by PDEs. Particularly, deterministic diffusion-reaction equations and these equations are employed to model the spatial spread of tumours at the initial development and later invasive moments [5].

The fractional-order derivatives are introduced by Leibnitz soon after the classical order derivatives. As compared to classical calculus, it was soon discovered that fractional calculus (FC) is more appropriate for describing real-world problems [610]. The calculus of arbitrary order turned out one of the most essential tools to describe biological phenomena. The human diseases which are modelled through derivative having fractional-order help us to incorporate the information about its present and past states [1116]. Moreover, these models demonstrate the non-local distributed effects, hereditary properties and system memory. These properties are necessary to describe the biological models. In connection with this, recently many authors established the arbitrary-order model to analyze the diffusion equation and to forecast the effect of the tumour and they applied many powerful methods to find the solution for these models [1732]. The pivotal aim of generalizing the integer to fractional order is to capture consequences related non-locality, long-range memory and time-based properties and also anomalous diffusion aspects. Many real-world problems exemplified with complex and nonlinear models are effectively, systematically and accurately illustrated and investigated by the aid of theory and fundamental concept of FC. Many pioneers nurtured and developed novel and distinct notions of fractional order for both differential and integral operators. Most familiarly hired operators to analyze many models are Riemann, Liouville, Caputo, Fabrizio and others. However, researchers pointed out some limitations while generalizing the system with these notions. The Riemann–Liouville derivative fails to elucidate the essence of initial conditions; the Caputo derivative has overcome this limitation and latter it has been widely applied to the numerous classes of mathematical models exemplifying real-world problems. But this fractional operator is unable to describe the singular kernel of the systems or problems. However, Caputo et al. [33] in 2015 overcome the foregoing obliges and then the number of authors employed CF derivative to investigate and study wide classes of complex and nonlinear problems. It has been proved by many researchers; CF fractional operator as great results compared to other fractional operators.

The study of mathematical models effectively exemplified diverse phenomena. However, as much as important of nurturing the phenomena with the system of equations finding the corresponding solution is also very vital and difficult. In this path, many authors examined diverse phenomena, for instance, the structured predator-prey model with prey refuge [34], COVID-19[3539], Zika virus transmission [40], planar system-masses in an equilateral triangle [41], a harmonic oscillator with position-dependent mass [42], time fractional Burgers equation [43], fractional optimal control problems [44], Emden–Flower type equations [45], and many others [4653]. These investigations help researchers to understand the importance of generalizing the classical concept into fractional operators, and efficiency and difference between diverse schemes.

Many physicists, engineers and mathematicians recently proposed and modified diverse solution procedure with a different approach with respect to increasing in accuracy and methodology, to reduces the complexity, many additional assumptions and consideration, huge time for evaluation and to save computer memory. Moreover, each method is suitable for some specific family of problems and they have their own limitations, including conversion of nonlinear to linear, partial to ordinary differential equations, splitting complex and nonlinear term to simple parts terms. In this connection, with the help of topological concept called homotopy, Liao Shijun who is a Chinese Mathematician proposed algorithm called homotopy analysis method (HAM) and illustrated to confirm it overcomes almost all the limitations raised while we solving nonlinear systems exists in sciences and other disciplines associated to mathematics [54]. The most familiar thing of employing this method by many authors is including it solves nonlinear problems without linearization and perturbation.

As science and technology-enhanced, mankind always expecting new tools or modifications in existing tools to improve the accuracy and reduces the time taken for finding needful. In this regard, some scholars pointed out similar things in HAM and suggested to union with existing and familiar transformation. Authors in [55] modified q-HAM with the help of Laplace transform (LT) and manifest new modified scheme is called q-HATM. This method is perceptible includes all merits which are achieved by HAM and also it attracted many researchers to analyze the diverse class of models and systems. For instance, the model exemplifying three Lakes pollution with the newly proposed fractional operator is investigated by authors in [56], fractional vibration equation is analyzed by authors in [57] with some interesting results, authors in [58] presented the efficiency of the projected scheme while analyzing Swift–Hohenberg equation having arbitrary order, the accuracy of the hired scheme in comparison with existing results is illustrated by authors in [59] with respect to the physical model, the convergence analysis of the considered method for Lienard’s equation is demonstrated in [60], many others analyzed various biological and physical phenomena by the assist of the projected scheme [6162].

In the present study, we find the series solution for a system of nonlinear differential equations describing the model of tumour invasion and metastasis using q-HATM with the help of a novel fractional operator. By using the important results of fixed-point theory for the projected system the existence and uniqueness are demonstrated. The novelty of the projected scheme gives more freedom to choose the initial conditions and the novelty is it offers a simple solution procedure and associated with parameters to provide the swift convergence. Further, it contains the results achieved by other classical methods including ADM, HPM, q-HAM and some other methods [6372]. In the present study, we analyzed the system describing the tumour invasion and metastasis with different time and space for different fractional-order using q-HATM within the frame of the novel fractional operator which can describe the singular kernel. This study can help us to analyze more complex and nonlinear mathematical models describing the deadly virus or diseases.

The rest of the manuscript is organized as follows: The basic and fundamentals are presented in the next section, the hired model is exemplified in Section 3, the basic procedure of the q-HATM is presented in Section 4, and its algorithm is illustrated for the considered model in Section 5. The existence and uniqueness for the archived results and error analysis are respectively presented in Sections 6 and 7. Moreover, with the aid of behaviour captured for the obtained result, the corresponding comments and conclusion are respectively exemplified in Sections 8 and9.

2  Preliminaries

The basic definitions are presented in this segment for the FC and Laplace transform. Specifically, we recall the notions related to Caputo-Fabrizio fractional operator [33,73].

Definition 1.The CF fractional derivative for fH1(a,b) is presented as [33]

Dtα(f(t))=M(α)1αatf(t)exp[αtϑ1α]dϑ,b>a,α[0,1], (1)

where M(α) is a normalization function and admits M(0)=M(1)=1. Further, if fH1(a,b) then we have

Dtα(f(t))=αM(α)1αat(f(t)f(ϑ))exp[αtϑ1α]dϑ. (2)

Definition 2. The CF fractional derivative for fH1(a,b) is presented as [73]

Itα(f(t))=2(1α)(2α)M(α)f(t)+2α(2α)M(α)0tf(ϑ)dϑ,0<α<1,t0. (3)

Note: According to [73], the following must hold

2(1α)(2α)M(α)+2α(2α)M(α)=1,0<α<1, (4)

which gives M(α)=22α. By the assist of above equation researchers in [73] proposed a novel Caputo derivative as follows:

Dtα(f(t))=11α0tf(t)exp[αtϑ1α]dϑ,0<α<1. (5)

Definition 3. The LT for a CF derivative 0CFDtαf(t) is presented as [33] below

L[0CFCDt(α+n)f(t)]=sn+1L[f(t)]snf(0)sn1f(0)f(n)(0)s+(1s)α. (6)

3  Mathematical Model

On the basis of generic solid tumour growth and assuming it is in avascular stage, the mathematical model has been proposed [74]. In this stage, most of the tumours are asymptomatic and further there is a possibility of cells to migrate and escape to the lymph nodes. The considered system exemplifies the interfaces of the surrounding tissue with the tumour and it can be extended to incorporate tumour and the vasculature. Here, the projected system of equations illustrates the interactions of the matrix-degrading enzymes (MDE, signifies by E), extra cellular matrix (ECM, symbolised by C) and tumour cells (denoted by T). With respect to ECM, most of the macromolecules are essential for cell motility, spreading and adhesion. Further, the ECM associated with many macromolecules, for instances collagen, laminin and fibronectin. During the various stages of metastasis, invasion and turn our growth, MDEs play a vital role. The ECM locally degrades by MDEs which are produced by tumour cells. Further, the method wherein they interact with tumour cells, growth factors and inhibitors are highly intricate. The tumour cells in the considered system as haptotaxis and in order to integrate this concept in the model, the hypotactic flux is considered as [75,76]


where χ>0 denotes haptotactic coefficient and which is constant. The random motion is another contribution to tumour cell motility and it helps to study ECM in isolation. Moreover, the flux is defined for the tumour cells with exemplified random motility is


where T is the chemokinetic response, D(C,E) is the function of either the ECM or MDE concentration, or constant.

For the tumour cell density (T), the conservation equation is presented as


and for the cell proliferation absence, the equation describing tumour cell motion is defined as

Tt=(D(C,M)T)χ(TC) (7)

For the notation, the random motility coefficient of tumour cell is considered as D(C,M)=DT and which is constant. Therefore, the degradation process is exemplified by the subsequent equation with positive constant δ

Ct=δEC. (8)

Active MDEs are formed by T, experience some form of decay and diffuse throughout the tissue. The equation modelling the evolution of MDE concentration is presented with MDE diffusion coefficient DE as

Et=DE2E+g(T,E)h(T,E,C), (9)

where g=μT and h=λE, h and g are the functions respectively describing the MDE decay and the production of active MDEs. Moreover, in the surrounding tissues there is a linear relationship between the level of active MDEs and the density of tumour cells.

From the above description, the system is presented as [18,54,57]:



Et=DE2Ediffusion+μTproductionλE.decay (10)

Here, with appropriate initial conditions, Eq. (10) is assumed to satisfy on a region of tissue or domain Ω. Moreover, the model is nurtured so that the MDEs and tumour cells remain inside the domain of tissue within deliberation and hence no-flux boundary conditions are executed on Ω. The terms contained in the above system with ECM density (C0), tumour cell density (T0) and MDE concentration (E0) by setting

T~=TT0,C~=CC0,E~=EE0,x~=xl,t~=tτ, (11)

where l signifies scale length and τ is the time. Then we have a scaled system of equations by dropping the tildes for notational convenience [18,74,77]



Et=dE2E+αTβE, (12)

where dT=DT/D, γ=χC0/D, η=τE0δ,dE=DE/D, α=τμT0/E0 and β=τλ.

The projected model can be protracted to integrate interactions between blood vessels and the tumour cells [74].

Now, we modify the time derivative by the CF derivative in Eq. (12) and given by



0CFDtαE(x,t)=dE2Ex2+αTβE, (13)

where α is fractional order. The associated initial conditions are



E(x,0)=0.5ex2ε. (14)

4  Fundamental Idea of the Considered Scheme

In this section, we hired the differential equation to present the basic procedure of the projected scheme with initial conditions

0CFDtαv(x,t)+Rv(x,t)+Nv(x,t)=f(x,t),n1<αn, (15)


v(x,0)=𝓆(x). (16)

We obtained by applying LT on Eq. (15)

L[v(x,t)]𝓆(x)s+s+(1s)αs{L[Rv(x,t)]+L[Nv(x,t)]L[f(x,t)]}=0. (17)

For φ(x,t;q), N is contracted as follows:


+s+(1s)αs{L[Rφ(x,t;q)]+L[Nφ(x,t;q)]L[f(x,t)]}, (18)

where q[0,1𝓃]. Then, the homotopy is defined by results in [34]

(1𝓃q)L[φ(x,t;q)v0(x,t)]=qN[φ(x,t;q)], (19)

where L is signifying LT. For q=0 and q=1n, the following conditions satisfies

φ(x,t;0)=v0(x,t),φ(x,t;1n)=v(x,t). (20)

By using Taylor theorem we get

φ(x,t;q)=v0(x,t)+m=1vm(x,t)qm, (21)


vm(x,t)=1m!mφ(x,t;q)qm|q=0. (22)

For the proper chaise of v0(x,t),n and the series (13) converges at q=1𝓃. Then

v(x,t)=v0(x,t)+m=1vm(x,t)(1n)m. (23)

After differentiating Eq. (19) m-times with q and multiplying by 1m! and substituting q=0, one can get

L[vm(x,t)kmvm1(x,t)]=Rm(vm1), (24)

where the vectors are defined as

vm={v0(x,t),v1(x,t),,vm(x,t)}. (25)

Eq. (24) reduces after employing inverse LT to

vm(x,t)=kmvm1(x,t)+L1[Rm(vm1)], (26)



+s+(1s)αsL[Rvm1+Hm1], (27)


km={0,m1,n,m>1. (28)

Here, Hm is homotopy polynomial and presented as

Hm=1m![mφ(x,t;q)qm]q=0and φ(x,t;q)=φ0+qφ1+q2φ2+. (29)

By the help of Eqs. (26) and (27), we found


+L1{s+(1s)αsL[Rvm1+Hm1]}. (30)

By the help of q-HATM, the series solution is

v(x,t)=v0(x,t)+m=1vm(x,t)(1n)m. (31)

5  Implementation of the q-Homotopy Analysis Transform Method

Consider the system of equation cited in Eq. (13) describing the tumour invasion and metastasis in CF fractional derivative



0CFDtαE(x,t)dE2Ex2αT+βE=0. (32)

Applying Laplace transform on Eq. (32) and then with the help of Eq. (13), we get



L[E(x,t)]1s(0.5ex2ε)s+(1s)αsL{dE2Ex2+αTβE}=0. (33)

The non-linear operator N defined as

N1[ φ1(x,t;q),φ2(x,t;q),φ3(x,t;q) ]=L[ φ1(x,t;q) ]1s(ex2ε)s+(1s)αsL{ dT2φ1(x,t;q)x2                 γ[ φ1(x,t;q)xφ2(x,t;q)x+2φ2(x,t;q)x2 ] },N2[ φ1(x,t;q),φ2(x,t;q),φ3(x,t;q) ]=L[ φ2(x,t;q) ]1s(10.5ex2ε)                 +s+(1s)αsL{ ηφ3(x,t;q)φ2(x,t;q) },

N3[ φ1(x,t;q),φ2(x,t;q),φ3(x,t;q) ]=L[ φ3(x,t;q) ]1s(0.5ex2ε)s+(1s)αs                 ×L{ dE2φ3(x,t;q)x2+αφ1(x,t;q)βφ3(x,t;q) }. (34)

The m-th order deformation equation by the projected scheme at H(x,t)=1 is given by



L[Em(x,t)kmEm1(x,t)]=L1{R3,m[Tm1,Cm1,Em1]}, (35)


R1,m[Tm1,Cm1,Em1]=L[Tm1(x,t)](1kmn)1s(ex2ε)s+(1s)αsL{dT2Tm1x2γ[i=0m1TixCm1ix+2Cm1x2]},R2,m[Tm1,Cm1,Em1]=L[Cm1(x,t)](1kmn)1s(10.5ex2ε)+s+(1s)αsL{ηi=0m1EiCm1i},R3,m[Tm1,Cm1,Em1]=L[Em1(x,t)](1kmn)1s(0.5ex2ε)s+(1s)αsL{dE2Em1x2+αTm1βEm1}. (36)

On employing inverse LT on Eq. (35), it simplifies to



Em(x,t)=kmEm1(x,t)+L1{R3,m[Tm1,Cm1,Em1]}. (37)

By using T0(x,t),C0(x,t) and E0(x,t) and then solving the forgoing equations, we can obtain the terms of



E(x,t)=E0(x,t)+m=1Em(x,t)(1n)m. (38)

6  Existence and Uniqueness of Solutions

In this section, the existence and uniqueness are illustrated for the considered system with the assist of fixed-point theory. We consider the Eq. (32) as follows:

{0CFDtα[T(x,t)]=G1(x,t,T),0CFDtα[C(x,t)]=G2(x,t,C),0CFDtα[E(x,t)]=G3(x,t,E). (39)

Now, using Eq. (32) and results derived in [53], we obtained

{T(x,t)T(x,0)=0CFItα{dT2Tx2γ(TxCx+2Cx2)},C(x,t)C(x,0)=0CFItα{ηEC},E(x,t)E(x,0)=0CFItα{dE2Ex2+αTβE}. (40)

Then we have from [73] as follows:

{T(x,t)T(x,0)=2(1α)M(α)G1(x,t,T)+2α(2α)M(α)0tG1(x,ζ,T)dζ,C(x,t)C(x,0)=2(1α)M(α)G2(x,t,C)+2α(2α)M(α)0tG2(x,ζ,C)dζ,E(x,t)E(x,0)=2(1α)M(α)G3(x,t,E)+2α(2α)M(α)0tG3(x,ζ,E)dζ. (41)

Theorem 1. The kernel G1 admits the Lipschitz condition and contraction if 0(dTδ2γ(λ1δ+λ2))<1 satisfies.

Proof. Let us consider the two functions u and u1 to prove the theorem, then




(dTδ2γ(λ1δ+λ2))T(x,t)T(x,t1) (42)

where Cxλ2 and 2C2xλ3 be the bounded function. Putting η1=dTδ2γ(λ1δ+λ2) in the above inequality, then we have

G1(x,t,T)G1(x,t,T1)η1T(x,t)T(x,t1). (43)

Eq. (43) provides the Lipschitz condition for G1. Similarly, we can see that if 0dTδ2γ(λ1δ+λ2)<1, then it implies the contraction. Similarly, we can prove

{G2(x,t,C)G2(x,t,C1)η2C(x,t)C(x,t1),G3(x,t,E)G3(x,t,E1)η3E(x,t)E(x,t1). (44)

By the assist of the above equations, Eq. (41) simplifies to

{T(x,t)=T(x,0)+2(1α)(2α)M(α)G1(x,t,T)+2α(2α)M(α)0tG1(x,ζ,T)dζ,C(x,t)=C(x,0)+2(1α)(2α)M(α)G2(x,t,C)+2α(2α)M(α)0tG2(x,ζ,C)dζ,E(x,t)=E(x,0)+2(1α)(2α)M(α)G3(x,t,E)+2α(2α)M(α)0tG3(x,ζ,E)dζ. (45)

Then we get the recursive form as follows:

{Tn(x,t)=2(1α)(2α)M(α)G1(x,t,Tn1)+2α(2α)M(α)0tG1(x,ζ,Tn1)dζ,Cn(x,t)=2(1α)(2α)M(α)G2(x,t,Cn1)+2α(2α)M(α)0tG2(x,ζ,Cn1)dζ,En(x,t)=2(1α)(2α)M(α)G3(x,t,En1)+2α(2α)M(α)0tG3(x,ζ,En1)dζ. (46)

The associated initial conditions are

T(x,0)=T0(x,t),C(x,0)=C0(x,t)andE(x,0)=E0(x,t). (47)

Now, between the terms, the successive difference is defined as






+2α(2α)M(α)0t(G2(x,t,Cn1)G1(x,t,Cn2))dζ, (48)




Notice that

{Tn(x,t)=i=1nϕ1i(x,t),Cn(x,t)=i=1nϕ2i(x,t),En(x,t)=i=1nϕ3i(x,t). (49)

Then we have



+2α(2α)M(α)0t(G1(x,t,Tn1)G1(x,t,Tn2))dζ, (50)

Application of the triangular inequality, Eq. (50) reduces to


+2α(2α)M(α)0t(G1(x,t,Tn1)G1(x,t,Tn2))dζ (51)

The Lipschitz condition satisfied by the kernel T1, so, we have


+2α(2α)M(α)η10tϕ1(n1)(x,t)dζ. (52)

Similarly, we have

ϕ2n(x,t)2(1α)(2α)M(α)η2ϕ2(n1)(x,t)+2α(2α)M(α)η20tϕ2(n1)(x,ζ)dζ,ϕ3n(x,t)2(1α)(2α)M(α)η3ϕ3(n1)(x,t)+2α(2α)M(α)η30tϕ3(n1)(x,ζ)dζ. (53)

By the help of above result, we state the following theorem:

Theorem 2. If we have specific t0, then the solution for Eq. (32) will exist and unique. Further, we have

2(1α)(2α)M(α)ηi+2α(2α)M(α)ηit0<1,for i=1,2 and 3.

Proof. Let T(x,t),C(x,t) and E(x,t) be the bounded functions admitting the Lipschitz condition. Then, we get by Eqs. (52) and (53)

ϕ1i(x,t)Tn(x,0)[2(1α)(2α)M(α)η1+2α(2α)M(α)η1t]n, ϕ2i(x,t)Cn(x,0)[2(1α)(2α)M(α)η2+2α(2α)M(α)η2t]n,

ϕ3i(x,t)En(x,0)[2(1α)(2α)M(α)η3+2α(2α)M(α)η3t]n. (54)

Therefore, for the obtained solutions, continuity and existence are verified. Now, to prove the Eq. (54) is a solution for Eq. (32), we consider

T(x,t)T(x,0)=Tn(x,t)K1n(x,t), C(x,t)C(x,0)=Cn(x,t)K2n(x,t),

E(x,t)E(x,0)=En(x,t)K3n(x,t). (55)

Let us consider





2(1α)(2α)M(α)η1TTn1+2α(2α)M(α)η1TTn1t. (56)

This process gives


Similarly, at t0 we can obtain

K1n(x,t)(2(1α)(2α)M(α)+2α(2α)M(α)t0)n+1η1n+1M. (57)

As n and from Eq. (57), K1n(x,t)0. Similarly, we can verify for K2n(x,t) and K3n(x,t).

Next, for the solution of the projected model, we prove the uniqueness. Suppose T(x,t),C(x,t) and E(x,t) be the set of other solutions, then


+2α(2α)M(α)0t(G1(x,ζ,T)G1(x,ζ,T))dζ. (58)

Now, employing the norm on the above equation we get



2(1α)(2α)M(α)η1T(x,t)T(x,t)+2α(2α)M(α)η1tT(x,t)T(x,t). (59)

On simplification

T(x,t)T(x,t)(12(1α)(2α)M(α)η12α(2α)M(α)η1t)0. (60)

From the above condition, it is clear that T(x,t)=T(x,t), if

(12(1α)(2α)M(α)η12α(2α)M(α)η1t)0. (61)

Hence, Eq. (61) proves our required result.

7  Error Analysis of the q-Homotopy Analysis Transform Method

Theorem 3. Let (B[0,T],) be a Banach space and suppose vn(x,t) and v(x,t) define in the that, then the solution defined in Eq. (31) converges to the solution of Eq. (15), if 0<λ1<1.

Proof: Let {Sn} be a sequence of partial sum of Eq. (31). Then, we need to prove {Sn} is Cauchy sequence in (B[0,T],). Now, consider




For every n,mN(mn), now we have





λ1m+1(1λ1nm1λ1)v0. (62)

But 0<λ1<1, therefore limn,mSnSm=0. Hence, {Sn} is the Cauchy sequence.

Theorem 4. The maximum absolute error for the series solution of the Eq. (15) defined in Eq. (31) is determined as


Proof: By using Eq. (62), we get


But 0<λ1<01λ1nm<1. Hence, we have


This ends the proof.

Table 1: Description of parameters presented in the projected system [74]


8  Results and Discussion

Here, we demonstrate the future scheme is efficient and reliable and evaluate the approximate results for the system of partial differential equations representing a model of tumour invasion and metastasis. In the present study, we find the fourth-order solution to present the nature of the system. In Tab. 1, we present the specific values of the parameters cited in Fig. 1 captures the behaviour of q-HATM solution for tumour cells (T), extra cellular matrix (C), and matrix degrading enzymes (E) in 3D plots by using the Tab. 1 and the combined surface for the three components at the initial stage (i.e., t=0) is cited in Fig. 2. By generalizing the system with a newly nurtured fractional operator, it aids us to capture more interesting consequences associated with singular kernel. In the present work, we demonstrated the nature of q-HATM results for district α both in the change of x and t, and which are presented in Fig. 3. From these curves, we can observe that, as varying in both time and space with fractional order, the obtained results show noticeable vicissitudes in the behaviour. Specifically, extra cellular matrix and matrix degrading enzymes show stimulating behaviour for the change α.


Figure 1: Surfaces of q-HATM solution for (a) tumour cells (T), (b) extra cellular matrix (C),(c)matrix degrading enzymes (E) at n=1,α=1 and =1 and using Tab. 1


Figure 2: Surface of q-HATM solution for Eq. (32) at n=1,α=1,=1 and using Tab. 1


Figure 3: Nature of obtained solution for (a) tumour cells (T), (b) extra cellular matrix (C),(c)matrix degrading enzymes (E) with the change in time (t) for diverse α at n=1,=1 and using Tab. 1


Figure 4: -curves drown for q-HATM solution of (a) T(x,t),(b) C(x,t),(c) E(x,t) for distinct α at t=0.01,x=0.1,n=1 and using Tab. 1

The behaviors have been captured for different fractional Brownian motions and standard motion (α=1) with the change in . In Fig. 4, we drowned the -curves for the obtained solutions for T(x,t),C(x,t) and E(x,t) with the appropriate value of . The -curves aid to adjust and control the convergence province of the achieved results. Fig. 5 presents the 2D plots of an analytic-approximate solution for Eq. (32) at a distinct time. By the plots we can see that, the tumour cells and matrix degrading enzymes are also increases while time increases, but the extra cellular matrix decreases. Moreover, these types of investigation can open the door for analyses the stimulating models exemplifying deadly disease by incorporating diffusion co-efficient.


Figure 5: Response of obtained solution for the considered model with varying in x at (a) initial time (t=0), (b) t=0.1, (c) t=0.5 and (d) t=1 with n=1,α=1,=1 and Tab. 1

9  Conclusion

In the present study, we analyzed and capture the behaviour of the nonlinear fractional model of tumour invasion and metastasis by using the fractional operator and efficient analytical technique. The existences and uniqueness are demonstrated with the assist of a fixed point hypothesis. The plots captured in the present investigation display the stimulating behaviour and these can help scholars for some essential and interesting consequence of the hired system. The present study shows, the phenomena conspicuously be contingent on the time history and the time instant and, these can be proficiently studied using fundamental perceptions of FC and newly proposed fractional operator. The investigations of these types of models can provide new notions to analyze more real-world problems and it opens the door for employing an efficient method to study complex phenomena associated with science and technology.

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.


  1. American Cancer Society (2012). Cancer facts and figures 2012. Atlanta: American Cancer Society.
  2. Chaplain, M. A. J., & Stuart, A. M. (1993). A model mechanism for the chemotactic response of endothelial cells to tumour angiogenesis factor. Mathematical Medicine and Biology: A Journal of the IMA, 10(3), 149-168. [Google Scholar] [CrossRef]
  3. Sherratt, J. A., & Nowak, M. A. (1992). Oncogenes, anti-oncogenes and the immune response to cancer: A mathematical model. Proceedings of the Royal Society of London B, 248(1323), 261-271. [Google Scholar] [CrossRef]
  4. Gatenby, R. A., & Gawlinski, E. T. (1996). A reaction-diffusion model of cancer invasion. Cancer Research, 56(24), 5745-5753. [Google Scholar]
  5. Anderson, A. R. A., & Chaplain, M. A. J. (1998). Continuous and discrete mathematical models of tumor-induced angiogenesis. Bulletin of Mathematical Biology, 60(5), 857-899. [Google Scholar] [CrossRef]
  6. Riemann, G. F. B. (1896). Versucheinerallgemeinen auffassung der integration und differentiation. Gesammelte mathematische werke. Leipzig, Germany: Druck und Verlag.
  7. Caputo, M. (1969). Elasticita e dissipazione. Zanichelli, Bologna, Italy.
  8. Miller, K. S., Ross, B. (1993). An introduction to fractional calculus and fractional differential equations. New York: A Wiley.
  9. Podlubny, I. (1999). Fractional differential equations. New York: Academic Press.
  10. Kilbas, A. A., Srivastava, H. M., Trujillo, J. J. (2006). Theory and applications of fractional differential equations. Amsterdam: Elsevier.
  11. Ionescu, C., Lopes, A., Copot, D., Machado, J. A. T., & Bates, J. H. T. (2017). The role of fractional calculus in modeling biological phenomena. Communications in Nonlinear Science and Numerical Simulation, 51, 141-159. [Google Scholar] [CrossRef]
  12. Veeresha, P., Prakasha, D. G., & Baskonus, H. M. (2019). New numerical surfaces to the mathematical model of cancer chemotherapy effect in caputo fractional derivatives. Chaos, 29, 13119. [Google Scholar] [CrossRef]
  13. Yang, X. J., Baleanu, D., Khan, Y., & Mohyud-Din, S. T. (2013). Local fractional variational iteration method for diffusion and wave equations on cantor set. Romanian Journal of Physics, 59(1–2), 36-48. [Google Scholar]
  14. Veeresha, P., Baskonus, H. M. (2021). A powerful iterative approach for quintic complex Ginzburg–Landau equation within the frame of fractional operator. Fractals. Singapore: World Scientific. DOI 10.1142/S0218348X21400235. [CrossRef]
  15. Merdan, M., Gökdoğan, A., Yıldırım, A., Mohyud-Din, S. T. (2012). Numerical simulation of fractional Fornberg–Whitham equation by differential transformation method. Abstract and Applied Analysis, 1–8. DOI 10.1155/2012/965367. [CrossRef]
  16. Prakasha, D. G., Malagi, N. S., Veeresha, P., & Prasannakumara, B. C. (2021). An efficient computational technique for time-fractional Kaup–Kupershmidt equation. Numerical Methods for Partial Differential Equations, 37(2), 1299-1316. [Google Scholar] [CrossRef]
  17. Jain, R., Arekar, K., & Dubey, R. S. (2017). Study of Bergman’s minimal blood glucose-insulin model by Adomian decomposition method. Journal of Information and Optimization Sciences, 38(1), 133-149. [Google Scholar] [CrossRef]
  18. Meral, G., & Yamanlar, I. C. (2018). Mathematical analysis and numerical simulations for the cancer tissue invasion model. Communications Faculty of Sciences University of Ankara Series A1 Mathematics and Statistics, 68, 371-391. [Google Scholar] [CrossRef]
  19. Usha, S., Abinaya, V., Loghambal, S., & Rajendran, L. (2012). Non-linear mathematical model of the interaction between tumor and on colytic viruses. Applied Mathematics, 3(9), 1089-1096. [Google Scholar] [CrossRef]
  20. Sun, H. G., Zhang, Y., Baleanu, D., Chen, W., & Chen, Y. Q. (2018). A new collection of real world applications of fractional calculus in science and engineering. Communications in Nonlinear Science and Numerical Simulation, 64(48103), 213-231. [Google Scholar] [CrossRef]
  21. Atangana, A., & Baleanu, D. (2016). New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model. Thermal Science, 20(2), 763-769. [Google Scholar] [CrossRef]
  22. Veeresha, P., Prakasha, D. G., & Baskonus, H. M. (2019). Novel simulations to the time-fractional Fisher’s equation. Mathematical Sciences, 13(1), 33-42. [Google Scholar] [CrossRef]
  23. Veeresha, P., & Prakasha, D. G. (2021). Solution for fractional Kuramoto–Sivashinsky equation using novel computational technique. International Journal of Applied and Computational Mathematics, 7(2), 1. [Google Scholar] [CrossRef]
  24. Prakash, A., Veeresha, P., Prakasha, D. G., & Goyal, M. (2019). A homotopy technique for fractional order multi-dimensional telegraph equation via Laplace transform. The European Physical Journal Plus, 134(1), 3698. [Google Scholar] [CrossRef]
  25. Veeresha, P., & Prakasha, D. G. (2020). An efficient technique for two-dimensional fractional order biological population model. International Journal of Modeling, Simulation, and Scientific Computing, 11(1), 2050005. [Google Scholar] [CrossRef]
  26. Atangana, A., & Alkahtani, B. T. (2015). Analysis of the Keller–Segel model with a fractional derivative without singular kernel. Entropy, 17(12), 4439-4453. [Google Scholar] [CrossRef]
  27. Atangana, A., & Alkahtani, B. T. (2016). Analysis of non-homogenous heat model with new trend of derivative with fractional order. Chaos Solitons Fractals, 89(273), 566-571. [Google Scholar] [CrossRef]
  28. Panda, S. K., Abdeljawad, T., & Ravichandran, C. (2020). A complex valued approach to the solutions of Riemann–Liouville integral, Atangana–Baleanu integral operator and non-linear telegraph equation via fixed point method. Chaos Solitons Fractals, 130(2), 109439. [Google Scholar] [CrossRef]
  29. Belmor, S., Ravichandran, C., & Jarad, F. (2020). Nonlinear generalized fractional differential equations with generalized fractional integral conditions. Journal of Taibah University for Science, 14(1), 114-123. [Google Scholar] [CrossRef]
  30. Gao, W., Veeresha, P., Baskonus, H. M., Prakasha, D. G., & Kumar, P. (2020). A new study of unreported cases of 2019-nCOV epidemic outbreaks. Chaos Solitons Fractals, 138(554), 109929. [Google Scholar] [CrossRef]
  31. Valliammal, N., Ravichandran, C., & Nisar, K. S. (2020). Solutions to fractional neutral delay differential nonlocal systems. Chaos Solitons Fractals, 138(1), 109912. [Google Scholar] [CrossRef]
  32. Veeresha, P., Prakasha, D. G., Singh, J., Khan, I., & Kumar, D. (2020). Analytical approach for fractional extended Fisher–Kolmogorov equation with Mittag-Leffler kernel. Advances in Difference Equations, 174, [Google Scholar] [CrossRef]
  33. Caputo, M., & Fabrizio, M. (2015). A new definition of fractional derivative without singular kernel. Progress in Fractional Differentiation and Applications, 1(2), 73-85. [Google Scholar] [CrossRef]
  34. Baishya, C. (2020). Dynamics of fractional stage structured predator prey model with prey refuge. Indian Journal of Ecology, 47(4), 1118-1124. [Google Scholar]
  35. Baba, I. A., & Nasidi, B. A. (2020). Fractional order model for the role of mild cases in the transmission of COVID-19. Chaos Solitons Fractals, 142, 110374. [Google Scholar] [CrossRef]
  36. Ahmed, I., Baba, I. A., Yusuf, A., Kumam, P., & Kumam, W. (2020). Analysis of Caputo fractional-order model for COVID-19 with lockdown. Advances in Difference Equations, 394(1), 119. [Google Scholar] [CrossRef]
  37. Baba, I. A., & Nasidi, B. A. (2021). Fractional order epidemic model for the dynamics of novel COVID-19. Alexandria Engineering Journal, 60(1), 537-548. [Google Scholar] [CrossRef]
  38. Baba, I. A., Baba, B. A., & Esmaili, P. (2020). A mathematical model to study the effectiveness of some of the strategies adopted in curtailing the spread of COVID-19. Computational and Mathematical Methods in Medicine, 2020(1), 1-6. [Google Scholar] [CrossRef]
  39. Baba, I. A., & Baleanu, D. (2020). Awareness as the most effective measure to mitigate the spread of COVID-19 in Nigeria. Computers, Materials & Continua, 65(3), 1945-1957. [Google Scholar] [CrossRef]
  40. Rezapour, S., Mohammadi, H., & Jajarmi, A. (2020). A new mathematical model for Zika virus transmission. Advances in Difference Equations, 589(2020), 479. [Google Scholar] [CrossRef]
  41. Baleanu, D., Ghanbari, B., Asad, J. H., Jajarmi, A., & Pirouz, H. M. (2020). Planar system-masses in an equilateral triangle: Numerical study within fractional calculus. Computer Modeling in Engineering & Sciences, 124(3), 953-968. [Google Scholar] [CrossRef]
  42. Baleanu, D., Jajarmi, A., Sajjadi, S. S., & Asad, J. H. (2020). The fractional features of a harmonic oscillator with position-dependent mass. Communications in Theoretical Physics, 72(5), 055002. [Google Scholar] [CrossRef]
  43. Akram, T. (2020). An efficient numerical technique for solving time fractional Burgers equation. Alexandria Engineering Journal, 59(4), 2201-2220. [Google Scholar] [CrossRef]
  44. Jajarmi, A., & Baleanu, D. (2019). On the fractional optimal control problems with a general derivative operator. Asian Journal of Control, 23(2), 1062-1071. [Google Scholar] [CrossRef]
  45. Iqbal, M. K., Abbas, M., & Wasim, I. (2018). New cubic B-spline approximation for solving third order Emden–Flower type equations. Applied Mathematics and Computation, 331(1), 319-333. [Google Scholar] [CrossRef]
  46. Khalid, N., Abbas, M., Iqbal, M. K., Singh, J., & Ismail, A. I. M. (2020). A computational approach for solving time fractional differential equation via spline functions. Alexandria Engineering Journal, 59(5), 3061-3078. [Google Scholar] [CrossRef]
  47. Sajjadia, S. S., Baleanu, D., Jajarmi, A., & Pirouz, H. M. (2020). A new adaptive synchronization and hyper chaos control of a biological snap oscillator. Chaos Solitons Fractals, 138, 109919. [Google Scholar] [CrossRef]
  48. Gao, W., Veeresha, P., Prakasha, D. G., & Baskonus, H. M. (2020). New numerical simulation for fractional Benney–Lin equation arising in falling film problems using two novel techniques. Numerical Methods for Partial Differential Equations, 37(1), 210-243. [Google Scholar] [CrossRef]
  49. Baishya, C. (2019). A new application of hermite collocation method. International Journal of Mathematical, Engineering and Management Sciences, 4(1), 182-190. [Google Scholar] [CrossRef]
  50. Jajarmi, A., & Baleanu, D. (2020). A new iterative method for the numerical solution of high-order nonlinear fractional boundary value problems. Frontiers in Physics, 8(220), 545. [Google Scholar] [CrossRef]
  51. Khalid, N., Abbas, M., & Iqbal, M. K. (2019). Non-polynomial quintic spline for solving fourth-order fractional boundary value problems involving product terms. Applied Mathematics and Computation, 349(1), 393-407. [Google Scholar] [CrossRef]
  52. Gao, W., Veeresha, P., Prakasha, D. G., Baskonus, H. M., & Yel, G. (2020). New approach for the model describing the deathly disease in pregnant women using Mittag–Leffler function. Chaos Solitons Fractals, 134, 109696. [Google Scholar] [CrossRef]
  53. Owusu-Mensah, I., Akinyemi, L., Oduro, B., & Iyiola, O. S. (2020). A fractional order approach to modeling and simulations of the novel COVID-19. Advances in Difference Equations, 683(1), 211. [Google Scholar] [CrossRef]
  54. Liao, S. J. (1997). Homotopy analysis method and its applications in mathematics. Journal of Basic Science and Engineering, 5(2), 111-125. [Google Scholar]
  55. Singh, J., Kumar, D., & Swroop, R. (2016). Numerical solution of time-and space-fractional coupled Burgers’ equations via homotopy algorithm. Alexandria Engineering Journal, 55(2), 1753-1763. [Google Scholar] [CrossRef]
  56. Prakasha, D. G., & Veeresha, P. (2020). Analysis of Lakes pollution model with Mittag–Leffler kernel. Journal of Ocean Engineering and Science, 5(4), 310-322. [Google Scholar] [CrossRef]
  57. Srivastava, H. M., Kumar, D., & Singh, J. (2017). An efficient analytical technique for fractional model of vibration equation. Applied Mathematical Modelling, 45(3), 192-204. [Google Scholar] [CrossRef]
  58. Veeresha, P., Prakasha, D. G., & Baleanu, D. (2020). Analysis of fractional Swift-Hohenberg equation using a novel computational technique. Mathematical Methods in the Applied Sciences, 43(4), 1970-1987. [Google Scholar] [CrossRef]
  59. Veeresha, P., & Prakasha, D. G. (2020). A reliable analytical technique for fractional Caudrey–Dodd–Gibbon equation with Mittag–Leffler kernel. Nonlinear Engineering, 9(1), 319-328. [Google Scholar] [CrossRef]
  60. Kumar, D., Agarwal, R. P., & Singh, J. (2018). A modified numerical scheme and convergence analysis for fractional model of Lienard’s equation. Journal of Computational and Applied Mathematics, 399(1), 405-413. [Google Scholar] [CrossRef]
  61. Veeresha, P., & Prakasha, D. G. (2019). Solution for fractional Zakharov–Kuznetsov equations by using two reliable techniques. Chinese Journal of Physics, 60(1), 313-330. [Google Scholar] [CrossRef]
  62. Safare, K. M., Betageri, V. S., Prakasha, D. G., Veeresha, P., & Kumar, S. (2021). A mathematical analysis of ongoing outbreak COVID-19 in India through nonsingular derivative. Numerical Methods for Partial Differential Equations, 37(2), 1282-1298. [Google Scholar] [CrossRef]
  63. Veeresha, P., Prakasha, D. G., & Kumar, S. (2020). A fractional model for propagation of classical optical solitons by using non-singular derivative. Mathematical Methodsin the Applied Sciences, 75(4), 125. [Google Scholar] [CrossRef]
  64. Gao, W., Baskonus, H. M., & Shi, L. (2020). New investigation of bats-hosts-reservoir-people coronavirus model and apply to 2019-nCoV system. Advances in Difference Equations, 2020(391), 1-11. [Google Scholar] [CrossRef]
  65. Al-Ghafri, K. S., & Rezazadeh, H. (2019). Solitons and other solutions of (3 + 1)-dimensional space-time fractional modified KdV–Zakharov–Kuznetsov equation. Applied Mathematics and Nonlinear Sciences, 4(2), 289-304. [Google Scholar] [CrossRef]
  66. Durur, H., Ilhan, E., & Bulut, H. (2020). Novel complex wave solutions of the (2 + 1)-dimensional hyperbolic nonlinear schrödinger equation. Fractal and Fractional, 4(3), 41. [Google Scholar] [CrossRef]
  67. Ilhan, E., & Kıymaz, I. O. (2020). A generalization of truncated M-fractional derivative and applications to fractional differential equations. Applied Mathematics and Nonlinear Sciences, 5(1), 171-188. [Google Scholar] [CrossRef]
  68. Gao, W., Veeresha, P., Prakasha, D. G., & Baskonus, H. M. (2020). Novel dynamical structures of 2019-nCoV with nonlocal operator via powerful computational technique. Biology, 9(5), 107. [Google Scholar] [CrossRef]
  69. Yokus, A., & Gulbahar, S. (2019). Numerical solutions with linearization techniques of the fractional harry dym equation. Applied Mathematics and Nonlinear Sciences, 4(1), 35-42. [Google Scholar] [CrossRef]
  70. Gao, W., Veeresha, P., Prakasha, D. G., Baskonus, H. M., & Gulnur, Y. (2020). New numerical results for the time-fractional Phi-four equation using a novel analytical approach. Symmetry, 12(478), 1-16. [Google Scholar] [CrossRef]
  71. Brzeziński, D. W. (2018). Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus. Applied Mathematics and Nonlinear Sciences, 3(2), 487-502. [Google Scholar] [CrossRef]
  72. Gao, W., Gulnur, Y., Baskonus, H. M., & Cattani, C. (2020). Complex solitons in the conformable (2 + 1)-dimensional Ablowitz-Kaup–Newell–Segur equation. AIMS Math, 5(1), 507-521. [Google Scholar] [CrossRef]
  73. Losada, J., & Nieto, J. J. (2015). Properties of the new fractional derivative without singular Kernel. Progress in Fractional Differentiation and Applications, 1, 87-92. [Google Scholar] [CrossRef]
  74. Anderson, A. R. A., Chaplain, M. A. J., Newman, E. L., Steeele, R. J. C., & Thompson, A. M. (2000). Mathematical modelling of tumour invasion and metastasis. Computational and Mathematical Methods in Medicine, 2(2), 129-154. [Google Scholar] [CrossRef]
  75. Stetler-Stevenson, W. G., Hewitt, R., & Corcoran, M. (1996). Matrix metallo-proteinases and tumour invasion, from correlation to causality to the clinic. Cancer Biology, 7(3), 147-154. [Google Scholar] [CrossRef]
  76. Chambers, A. F., & Matrisian, L. M. (1997). Changing views of the role of matrix metalloproteinases in metastasis. Journal of the National Cancer Institute, 89(17), 1260-1270. [Google Scholar] [CrossRef]
  77. Mahiddin, N., Ali, S. A. H. (2014). Approximate analytical solutions for mathematical model of tumour invasion and metastasis using modified Adomian decomposition and homotopy perturbation methods. Journal of Applied Mathematics, 1–13. DOI 10.1155/2014/654978. [CrossRef]
images This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.