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 the tumour 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.
Tumour cellinvasion and metastasisq-homotopy analysis transform methodCaputo–Fabrizio derivativeIntroduction
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 [2–4]. 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 [6–10]. 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 [11–16]. 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 [17–32]. 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[35–39], 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 [46–53]. 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 [61–62].
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 [63–72]. 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 and 9.
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 f∈H1(a,b) is presented as [33]
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]
Jhapto=χT∇C,
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
Jrandom=−D(C,E)∇T,
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
∂T∂t+∇⋅(Jhapto+Jrandom)=0,
and for the cell proliferation absence, the equation describing tumour cell motion is defined as
∂T∂t=∇⋅(D(C,M)∇T)−χ∇⋅(T∇C)
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 δ
∂C∂t=−δEC.
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
∂E∂t=DE∇2E+g(T,E)−h(T,E,C),
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]:
∂T∂t=DT∇2T⏟randommotility−χ∇⋅(T∇C)⏟haptotaxis,
∂C∂t=−δEC⏟,degradation
∂E∂t=DE∇2E⏟diffusion+μT⏟production−λE⏟.decay
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τ,
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]
∂T∂t=dT∇2T−γ∇⋅(T∇C),
∂C∂t=−ηEC,
∂E∂t=dE∇2E+αT−βE,
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αT(x,t)=dT∂2T∂x2−γ[∂T∂x∂C∂x+∂2C∂x2],
0CFDtαC(x,t)=−ηEC,
0CFDtαE(x,t)=dE∂2E∂x2+αT−βE,
where α is fractional order. The associated initial conditions are
T(x,0)=e−x2ε,
C(x,0)=1−0.5e−x2ε,
E(x,0)=0.5e−x2ε.
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
By using T0(x,t),C0(x,t) and E0(x,t) and then solving the forgoing equations, we can obtain the terms of
T(x,t)=T0(x,t)+∑m=1∞Tm(x,t)(1n)m,
C(x,t)=C0(x,t)+∑m=1∞Cm(x,t)(1n)m,
E(x,t)=E0(x,t)+∑m=1∞Em(x,t)(1n)m.
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:
where ‖∂C∂x‖≤λ2 and ‖∂2C∂2x‖≤λ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)‖≤η1‖T(x,t)−T(x,t1)‖.
Eq. (43) provides the Lipschitz condition for G1. Similarly, we can see that if 0≤dTδ2−γ(λ1δ+λ2)<1, then it implies the contraction. Similarly, we can prove
From the above condition, it is clear that T(x,t)=T∗(x,t), if
(1−2(1−α)(2−α)M(α)η1−2α(2−α)M(α)η1t)≥0.
Hence, Eq. (61) proves our required result.
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
‖Sn+1(t)−Sn(t)‖=‖vn+1(x,t)‖
≤λ1‖vn(x,t)‖
≤λ12‖vn−1(x,t)‖≤…≤λ1n+1‖v0(x,t)‖
For every n,m∈N(m≤n), now we have
‖Sn−Sm‖=‖(Sn−Sn−1)+(Sn−1−Sn−2)+…+(Sm+1−Sm)‖
≤‖Sn−Sn−1‖+‖Sn−1−Sn−2‖+…+‖Sm+1−Sm‖
≤(λ1n+λ1n−1+…+λ1m+1)‖v0‖
≤λ1m+1(λ1n−m−1+λ1n−m−2+…+λ1+1)‖v0‖
≤λ1m+1(1−λ1n−m1−λ1)‖v0‖.
But 0<λ1<1, therefore limn,m→∞‖Sn−Sm‖=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
‖v(x,t)−∑n=0Mvn(x,t)‖≤λ1M+11−λ1‖v0(x,t)‖.
Proof: By using Eq. (62), we get
‖v(x,t)−Sn‖=λ1m+1(1−λ1n−m1−λ1)‖v0(x,t)‖.
But 0<λ1<0⇒1−λ1n−m<1. Hence, we have
‖v(x,t)−∑n=0Mvn(x,t)‖≤λ1M+11−λ1‖v0(x,t)‖.
This ends the proof.
Description of parameters presented in the projected system [74]
Parameters
Descriptions
Parameters
Value
DE
MDE diffusion coefficient
ε
0.01
DT
Tumour cell random motility coefficient
dT, dE
0.001
δ
Degradation rate for normal cells
η
10
χ
Haptotaxis coefficient
γ
0.005
μ
Production for MDE
α
0.1
λ
Decay rate for MDE
β
0.5
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 α.
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
Surface of q-HATM solution for Eq. (32) at n=1,α=1,ℏ=−1 and using Tab. 1
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
ℏ-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.
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
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.
ReferencesAmerican Cancer Society (2012). . Atlanta: American Cancer Society.Chaplain, M. A. J., Stuart, A. M. (1993). A model mechanism for the chemotactic response of endothelial cells to tumour angiogenesis factor. ,10(3),149–168. DOI 10.1093/imammb/10.3.149.Sherratt, J. A., Nowak, M. A. (1992). Oncogenes, anti-oncogenes and the immune response to cancer: A mathematical model. ,248(1323),261–271. DOI 10.1098/rspb.1992.0071.Gatenby, R. A., Gawlinski, E. T. (1996). A reaction-diffusion model of cancer invasion. ,56(24),5745–5753.Anderson, A. R. A., Chaplain, M. A. J. (1998). Continuous and discrete mathematical models of tumor-induced angiogenesis. ,60(5),857–899. DOI 10.1006/bulm.1998.0042.Riemann, G. F. B. (1896). Versucheinerallgemeinen auffassung der integration und differentiation. . Leipzig, Germany: Druck und Verlag.Caputo, M. (1969). . Zanichelli, Bologna, Italy.Miller, K. S., Ross, B. (1993). . New York: A Wiley.Podlubny, I. (1999). . New York: Academic Press.Kilbas, A. A., Srivastava, H. M., Trujillo, J. J. (2006). . Amsterdam: Elsevier.Ionescu, C., Lopes, A., Copot, D., Machado, J. A. T., Bates, J. H. T. (2017). The role of fractional calculus in modeling biological phenomena. ,51,141–159. DOI 10.1016/j.cnsns.2017.04.001.Veeresha, P., Prakasha, D. G., Baskonus, H. M. (2019). New numerical surfaces to the mathematical model of cancer chemotherapy effect in caputo fractional derivatives. ,29,13119. DOI 10.1063/1.5074099.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. ,59(1–2),36–48.Veeresha, P., Baskonus, H. M. (2021). A powerful iterative approach for quintic complex Ginzburg–Landau equation within the frame of fractional operator. . Singapore: World Scientific. DOI 10.1142/S0218348X21400235.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. ,1–8. DOI 10.1155/2012/965367.Prakasha, D. G., Malagi, N. S., Veeresha, P., Prasannakumara, B. C. (2021). An efficient computational technique for time-fractional Kaup–Kupershmidt equation. ,37(2),1299–1316. DOI 10.1002/num.22580.Jain, R., Arekar, K., Dubey, R. S. (2017). Study of Bergman’s minimal blood glucose-insulin model by Adomian decomposition method. ,38(1),133–149. DOI 10.1080/02522667.2016.1187919.Meral, G., Yamanlar, I. C. (2018). Mathematical analysis and numerical simulations for the cancer tissue invasion model. ,68,371–391. DOI 10.31801/cfsuasmas.421546.Usha, S., Abinaya, V., Loghambal, S., Rajendran, L. (2012). Non-linear mathematical model of the interaction between tumor and on colytic viruses. ,3(9),1089–1096. DOI 10.4236/am.2012.39160.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. ,64(48103),213–231. DOI 10.1016/j.cnsns.2018.04.019.Atangana, A., Baleanu, D. (2016). New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model. ,20(2),763–769. DOI 10.2298/TSCI160111018A.Veeresha, P., Prakasha, D. G., Baskonus, H. M. (2019). Novel simulations to the time-fractional Fisher’s equation. ,13(1),33–42. DOI 10.1007/s40096-019-0276-6.Veeresha, P., Prakasha, D. G. (2021). Solution for fractional Kuramoto–Sivashinsky equation using novel computational technique. ,7(2),1. DOI 10.1007/s40819-021-00956-0.Prakash, A., Veeresha, P., Prakasha, D. G., Goyal, M. (2019). A homotopy technique for fractional order multi-dimensional telegraph equation via Laplace transform. ,134(1),3698. DOI 10.1140/epjp/i2019-12411-y.Veeresha, P., Prakasha, D. G. (2020). An efficient technique for two-dimensional fractional order biological population model. ,11(1),2050005. DOI 10.1142/S1793962320500051.Atangana, A., Alkahtani, B. T. (2015). Analysis of the Keller–Segel model with a fractional derivative without singular kernel. ,17(12),4439–4453. DOI 10.3390/e17064439.Atangana, A., Alkahtani, B. T. (2016). Analysis of non-homogenous heat model with new trend of derivative with fractional order. ,89(273),566–571. DOI 10.1016/j.chaos.2016.02.012.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. ,130(2),109439. DOI 10.1016/j.chaos.2019.109439.Belmor, S., Ravichandran, C., Jarad, F. (2020). Nonlinear generalized fractional differential equations with generalized fractional integral conditions. ,14(1),114–123. DOI 10.1080/16583655.2019.1709265.Gao, W., Veeresha, P., Baskonus, H. M., Prakasha, D. G., Kumar, P. (2020). A new study of unreported cases of 2019-nCOV epidemic outbreaks. ,138(554),109929. DOI 10.1016/j.chaos.2020.109929.Valliammal, N., Ravichandran, C., Nisar, K. S. (2020). Solutions to fractional neutral delay differential nonlocal systems. ,138(1),109912. DOI 10.1016/j.chaos.2020.109912.Veeresha, P., Prakasha, D. G., Singh, J., Khan, I., Kumar, D. (2020). Analytical approach for fractional extended Fisher–Kolmogorov equation with Mittag-Leffler kernel. ,174. DOI 10.1186/s13662–020-02617-w.Caputo, M., Fabrizio, M. (2015). A new definition of fractional derivative without singular kernel. ,1(2),73–85. DOI 10.12785/pfda/010201.Baishya, C. (2020). Dynamics of fractional stage structured predator prey model with prey refuge. ,47(4),1118–1124.Baba, I. A., Nasidi, B. A. (2020). Fractional order model for the role of mild cases in the transmission of COVID-19. ,142,110374. DOI 10.1016/j.chaos.2020.110374.Ahmed, I., Baba, I. A., Yusuf, A., Kumam, P., Kumam, W. (2020). Analysis of Caputo fractional-order model for COVID-19 with lockdown. ,394(1),119. DOI 10.1186/s13662-020-02853-0.Baba, I. A., Nasidi, B. A. (2021). Fractional order epidemic model for the dynamics of novel COVID-19. ,60(1),537–548. DOI 10.1016/j.aej.2020.09.029.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. ,2020(1),1–6. DOI 10.1155/2020/5248569.Baba, I. A., Baleanu, D. (2020). Awareness as the most effective measure to mitigate the spread of COVID-19 in Nigeria. ,65(3),1945–1957. DOI 10.32604/cmc.2020.011508.Rezapour, S., Mohammadi, H., Jajarmi, A. (2020). A new mathematical model for Zika virus transmission. ,589(2020),479. DOI 10.1186/s13662-020-03044-7.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. ,124(3),953–968. DOI 10.32604/cmes.2020.010236.Baleanu, D., Jajarmi, A., Sajjadi, S. S., Asad, J. H. (2020). The fractional features of a harmonic oscillator with position-dependent mass. ,72(5),055002. DOI 10.1088/1572-9494/ab7700.Akram, T. (2020). An efficient numerical technique for solving time fractional Burgers equation. ,59(4),2201–2220. DOI 10.1016/j.aej.2020.01.048.Jajarmi, A., Baleanu, D. (2019). On the fractional optimal control problems with a general derivative operator. ,23(2),1062–1071. DOI 10.1002/asjc.2282.Iqbal, M. K., Abbas, M., Wasim, I. (2018). New cubic B-spline approximation for solving third order Emden–Flower type equations. ,331(1),319–333. DOI 10.1016/j.amc.2018.03.025.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. ,59(5),3061–3078. DOI 10.1016/j.aej.2020.06.007.Sajjadia, S. S., Baleanu, D., Jajarmi, A., Pirouz, H. M. (2020). A new adaptive synchronization and hyper chaos control of a biological snap oscillator. ,138,109919. DOI 10.1016/j.chaos.2020.109919.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. ,37(1),210–243. DOI 10.1002/num.22526.Baishya, C. (2019). A new application of hermite collocation method. ,4(1),182–190. DOI 10.33889/24557749.Jajarmi, A., Baleanu, D. (2020). A new iterative method for the numerical solution of high-order nonlinear fractional boundary value problems. ,8(220),545. DOI 10.3389/fphy.2020.00220.Khalid, N., Abbas, M., Iqbal, M. K. (2019). Non-polynomial quintic spline for solving fourth-order fractional boundary value problems involving product terms. ,349(1),393–407. DOI 10.1016/j.amc.2018.12.066.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. ,134,109696. DOI 10.1016/j.chaos.2020.109696.Owusu-Mensah, I., Akinyemi, L., Oduro, B., Iyiola, O. S. (2020). A fractional order approach to modeling and simulations of the novel COVID-19. ,683(1),211. DOI 10.1186/s13662-020-03141-7.Liao, S. J. (1997). Homotopy analysis method and its applications in mathematics. ,5(2),111–125.Singh, J., Kumar, D., Swroop, R. (2016). Numerical solution of time-and space-fractional coupled Burgers’ equations via homotopy algorithm. ,55(2),1753–1763. DOI 10.1016/j.aej.2016.03.028.Prakasha, D. G., Veeresha, P. (2020). Analysis of Lakes pollution model with Mittag–Leffler kernel. ,5(4),310–322. DOI 10.1016/j.joes.2020.01.004.Srivastava, H. M., Kumar, D., Singh, J. (2017). An efficient analytical technique for fractional model of vibration equation. ,45(3),192–204. DOI 10.1016/j.apm.2016.12.008.Veeresha, P., Prakasha, D. G., Baleanu, D. (2020). Analysis of fractional Swift-Hohenberg equation using a novel computational technique. ,43(4),1970–1987. DOI 10.1002/mma.6022.Veeresha, P., Prakasha, D. G. (2020). A reliable analytical technique for fractional Caudrey–Dodd–Gibbon equation with Mittag–Leffler kernel. ,9(1),319–328. DOI 10.1515/nleng-2020-0018.Kumar, D., Agarwal, R. P., Singh, J. (2018). A modified numerical scheme and convergence analysis for fractional model of Lienard’s equation. ,399(1),405–413. DOI 10.1016/j.cam.2017.03.011.Veeresha, P., Prakasha, D. G. (2019). Solution for fractional Zakharov–Kuznetsov equations by using two reliable techniques. ,60(1),313–330. DOI 10.1016/j.cjph.2019.05.009.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. ,37(2),1282–1298. DOI 10.1002/num.22579.Veeresha, P., Prakasha, D. G., Kumar, S. (2020). A fractional model for propagation of classical optical solitons by using non-singular derivative. ,75(4),125. DOI 10.1002/mma.6335.Gao, W., Baskonus, H. M., Shi, L. (2020). New investigation of bats-hosts-reservoir-people coronavirus model and apply to 2019-nCoV system. ,2020(391),1–11. DOI 10.1186/s13662-019-2438-0.Al-Ghafri, K. S., Rezazadeh, H. (2019). Solitons and other solutions of (3 + 1)-dimensional space-time fractional modified KdV–Zakharov–Kuznetsov equation. ,4(2),289–304. DOI 10.2478/AMNS.2019.2.00026.Durur, H., Ilhan, E., Bulut, H. (2020). Novel complex wave solutions of the (2 + 1)-dimensional hyperbolic nonlinear schrödinger equation. ,4(3),41. DOI 10.3390/fractalfract4030041.Ilhan, E., Kıymaz, I. O. (2020). A generalization of truncated M-fractional derivative and applications to fractional differential equations. ,5(1),171–188. DOI 10.2478/amns.2020.1.00016.Gao, W., Veeresha, P., Prakasha, D. G., Baskonus, H. M. (2020). Novel dynamical structures of 2019-nCoV with nonlocal operator via powerful computational technique. ,9(5),107. DOI 10.3390/biology9050107.Yokus, A., Gulbahar, S. (2019). Numerical solutions with linearization techniques of the fractional harry dym equation. ,4(1),35–42. DOI 10.2478/AMNS.2019.1.00004.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. ,12(478),1–16. DOI 10.3390/sym12030478.Brzeziński, D. W. (2018). Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus. ,3(2),487–502. DOI 10.2478/AMNS.2018.2.00038.Gao, W., Gulnur, Y., Baskonus, H. M., Cattani, C. (2020). Complex solitons in the conformable (2 + 1)-dimensional Ablowitz-Kaup–Newell–Segur equation. ,5(1),507–521. DOI 10.3934/math.2020034.Losada, J., Nieto, J. J. (2015). Properties of the new fractional derivative without singular Kernel. ,1,87–92. DOI 10.12785/pfda/010202.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. ,2(2),129–154. DOI 10.1080/10273660008833042.Stetler-Stevenson, W. G., Hewitt, R., Corcoran, M. (1996). Matrix metallo-proteinases and tumour invasion, from correlation to causality to the clinic. ,7(3),147–154. DOI 10.1006/scbi.1996.0020.Chambers, A. F., Matrisian, L. M. (1997). Changing views of the role of matrix metalloproteinases in metastasis. ,89(17),1260–1270. DOI 10.1093/jnci/89.17.1260.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. ,1–13. DOI 10.1155/2014/654978.