Investigation of Fractional Order Tumor Cell Concentration Equation Using Finite Difference Method


  • Kalyanrao Takale
  • Uttam Kharde Department of Mathematics, R.N.C. Arts, J.D.B. Commerce and N.S.C. Science College, Nashik, S.P. Pune University, India.
  • Shrikisan Gaikwad Department of Mathematics, New Arts, Commerce and Science College, Ahmednagar, S.P. Pune University, India.



Caputo Derivative, Finite Difference Scheme, Fractional Differential Equation, Python, Tumor cell.



The fractional differential equation provides a powerful mathematical framework for modeling and understanding a wide range of complex, non-local and memory-dependent phenomena in various scientific, engineering and real-world applications. Liver metastasis is a secondary cancerous tumor that develops in the liver due to the spread of cancer cells from a primary cancer originating in another part of the human body. The primary focus of this study is to understand tumor growth in the human liver, both in the presence and absence of medication therapy. To achieve this, a temporal fractional-order parabolic partial differential equation is utilized, and its analysis is carried out using numerical methods. The Caputo derivative is employed to explore the impact of medication therapy on tumor growth. To numerically solve the mathematical model, the Crank-Nicolson finite difference method has been developed. This method is selected for its remarkable attributes, including unconditional stability and second-order accuracy in both spatial and temporal dimensions. The outcomes and insights derived from this study are effectively communicated through various graphical representations. These visual aids serve as invaluable tools for comprehending the profound impact of medication therapy on the growth of liver tumors. Through the medium of these graphical depictions, one can glean a clearer and more intuitive understanding of the complex dynamics at play. The numeric solution to this intricate problem is achieved through the implementation of algorithms meticulously crafted using versatile and powerful Python programming language. Python’s flexibility, extensive libraries, and robust numerical computing capabilities make it a good choice for handling the complexities of this study.


Mubarak S, Khanday MA, Lone AUH. Mathematical analysis based on eigenvalue approach to study liver metastasis disease with applied drug therapy. Netw Model Anal Health Inform Bioinform. 2020; 9(1). DOI:

Khanday MA, Nazir K. Mathematical and numerical analysis of thermal distribution in cancerous tissues under the local heat therapy. Int J Biomath. 2017; 10(7): 1–10.

Swanson KR, Bridge C, Murray JD, Alvord EC. Virtual and real brain tumors: Using mathematical modeling to quantify glioma growth and invasion. J Neurol Sci. 2003; 216(1): 1–10.

Filipovic N, Djukic T, Saveljic I, Milenkovic P, Jovicic G, Djuric M. Modeling of liver metastatic disease with applied drug therapy. Comput Methods Programs Biomed. 2014; 115(3): 162–70. Available from:

Ghode K, Takale K, Gaikwad S. Traveling Wave Solutions of Fractional Differential Equations Arising in Warm Plasma. Baghdad Sci J. 2023; 20(1(SI)): 0318–0318.

Sonawane J, Sontakke B, Takale K. Approximate Solution of Sub diffusion Bio heat Transfer Equation. Baghdad Sci J.; 20(1(SI)): 0394.

Joshi H, Jha BK. Fractional-order mathematical model for calcium distribution in nerve cells. Comput Appl Math. 2020; 39(2): 1–22.

Khan AA, Amin R, Ullah S, Sumelka W, Altanji M. Numerical simulation of a Caputo fractional epidemic model for the novel coronavirus with the impact of environmental transmission. Alexandria Eng J. 2022; 61(7): 5083–95.

Alzubaidi AM, Othman HA, Ullah S, Ahmad N, Alam MM. Analysis of Monkeypox viral infection with human to animal transmission via a fractional and Fractal-fractional operators with power law kernel. Math Biosci Eng. 2023; 20(4): 6666–90.

Bolton L, Cloot AHJJ, Schoombie SW, Slabbert JP. A proposed fractional-order Gompertz model and its application to tumour growth data. Math Med Biol. 2015; 32(2): 187–207.

Iyiola OS, Zaman FD. A fractional diffusion equation model for cancer tumor. AIP Adv. 2014; 4(10): 107121.

Rihan FA, Velmurugan G. Dynamics of fractional-order delay differential model for tumor-immune system. Chaos Solitons Fractals. 2020; 132: 109592.

Abaid Ur Rehman M, Ahmad J, Hassan A, Awrejcewicz J, Pawlowski W, Karamti H, et al. The Dynamics of a Fractional-Order Mathematical Model of Cancer Tumor Disease. Symmetry (Basel). 2022; 14(8): 1694.

Owolabi KM, Atangana A. Numerical Methods for Fractional Differentiation. Springer Series in Computational Mathematics 54. Singapore. 1st Ed. 2019. 328 p.

Kadalbajoo MK, Awasthi A. A numerical method based on Crank-Nicolson scheme for Burgers’ equation. Appl Math Comput. 2006; 182(2): 1430–42.

Ali AH, Jaber AS, Yaseen MT, Rasheed M, Bazighifan O, Nofal TA. A Comparison of Finite Difference and Finite Volume Methods with Numerical Simulations: Burgers Equation Model. Complexity. 2022: 1–9.

Podlubny I. An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. Math Sci Eng. 1999; 198: 340.

Kharde U, Takale K, Gaikwad S. Crank-Nicolson Method For Time Fractional Drug Concentration Equation in Central Nervous System. Adv Appl Math Sci. 2022; 22(2): 407–433.





How to Cite

Investigation of Fractional Order Tumor Cell Concentration Equation Using Finite Difference Method. Baghdad Sci.J [Internet]. [cited 2024 May 28];21(9). Available from: