The mathematical oncology has received a lot of interest in recent years since it helps illuminate pathways and provides valuable quantitative predictions, which will shape more effective and focused future therapies. We discuss a new fractal-fractional-order model of the interaction among tumor cells, healthy host cells and immune cells. The subject of this work appears to show the relevance and ramifications of the fractal-fractional order cancer mathematical model. We use fractal-fractional derivatives in the Caputo senses to increase the accuracy of the cancer and give a mathematical analysis of the proposed model. First, we obtain a general requirement for the existence and uniqueness of exact solutions via Perov's fixed point theorem. The numerical approaches used in this paper are based on the Grünwald-Letnikov nonstandard finite difference method due to its usefulness to discretize the derivative of the fractal-fractional order. Then, two types of stabilities, Lyapunov's and Ulam-Hyers' stabilities, are established for the Incommensurate fractional-order and the Incommensurate fractal-fractional, respectively. The numerical results of this study are compatible with the theoretical analysis. Our approaches generalize some published ones because we employ the fractal-fractional derivative in the Caputo sense, which is more suitable for considering biological phenomena due to the significant memory impact of these processes. Aside from that, our findings are new in that we use Perov's fixed point result to demonstrate the existence and uniqueness of the solutions. The way of expressing the Ulam-Hyers' stabilities by utilizing the matrices that converge to zero is also novel in this area.
Keywords: Grünwald-Letnikov nonstandard finite difference method; cancer chaotic mode; existence and uniqueness; fractal-fractional calculus; stability.