Mathematical analysis and numerical simulation for fractal-fractional cancer model

: 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 e ff ective 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 di ff erence 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.


Introduction
Cancer is a word used to describe disorders in which aberrant cells divide uncontrollably and can infiltrate neighboring tissues.According to the World Health Organization (2020), cancer is the second leading cause of mortality globally, accounting for approximately one in every six deaths [1].Since the middle of the 1960s, mathematical modeling and nonlinear simulations of the tumor growth process has been researched due to the significant public health issues and the requirement for immediate health measures [2][3][4][5][6][7][8][9].Edward Lorenz, a meteorologist and mathematician, discovered the chaos phenomenon in the unpredictable and irregular behavior of nonlinear dynamical systems in 1963 [10].Chaos can be expressed mathematically via deterministic iterations of nonlinear difference equations or the development of nonlinear ordinary differential equations (ODEs) or partial differential equations (PDEs).The study of chaotic systems has been heralded as one of the most significant scientific accomplishments of the twentieth century.While the field is still in its infancy, there is no doubt that it is becoming increasingly important in various of scientific disciplines.To that end, chaos has been demonstrated to exist in a wide range of systems, including electronics [11], chemistry [12], economics and finance [13,14], biological systems [15][16][17][18] and so on.
It is worth noting that fractional calculus is a vital branch of mathematics.Because of the memory and genetic peculiarity of fractional-order differential equations, several researchers have modeled biological phenomena using fractional calculus derivatives.As a result, it is a very useful tool for describing genuine natural processes.Many papers on fractional-order dynamical models have recently been published [19][20][21].
Atangana [22] presented a new advanced type of fractal fractional derivative in 2017, bridging the gap between fractional and fractal calculus.Fractal-fractional operators contain two components: the fractional order and the fractal dimension (order).Differential equations using the fractal-fractional derivative transform the assumed system's order and dimension into a rational order system.The major goal of defining these derivatives is to examine fractal nonlocal boundary and initial value problems in nature.Certain mathematicians developed various results and designed some fractal-fractional models that exhibit improved simulations for representing mathematical structures in this direction [23][24][25][26].The nonstandard finite deferential numerical methods were first introduced by Mickens in 1994 [27].These methods are well known for maintaining the positivity, boundedness and stability of nonlinear systems' equilibrium points [27,28].
In the paper [15], authors introduced and on studied the following three-dimensional order cancer model where x(t) stands for the number of tumor cells at time t y(t) for the number of healthy host cells at time t, and z(t) for the number of effector immune cells present at time t within the single tumor-site compartment, and x 0 , y 0 and z 0 are the associated initial values of system (1.1).Here, the parameters a, b and c are positive real numbers that indicate the growth rates of populations of x(t), y(t) and z(t), respectively.If α i = α for every i = 1, ..., 3, then the system (1.3) is called commensurate order; otherwise, it is named incommensurate order [29].The fractional version of system (1.1) was considered in the paper [16], and were described by where C 0 D α is the α-order Caputo differential operator.The three-dimensional fractal-fractional-order cancer model is the main topic of this research: where C 0 D α,β is the (α, β) fractal-fractional-order Caputo differential operator.
The rest of this paper is organized as follows.Section 2 provides some fundamental definitions of generalized Banach spaces in the sense of Perov, its properties and fractal fractional operators in the Caputo sense.Section 3 is devoted to the existence and the uniqueness with Perov fixed point theorem.In Section 4, the suggested model's numerical solution was achieved using the Grünwald-Letnikov nonstandard finite difference scheme of Caputo derivative (in short GL-NSFDM) scheme using MATLAB software.Section 5 presents the Lyapunov's stability of the equilibrium points of the proposed system by varying the fractional order and the set of parameter (a, b, c), and by maintaining the fractal dimension (β 1 , β 2 , β 3 ) = (1, 1, 1).Section 6 shows the Ulam-Hyers stability of the Incommensurate fractal-fractional-order cancer model (1.3).Finally, the discussion and the conclusion are given in the last two sections.

Preliminary
We present some basic notation, results of generalized Banach spaces in the sense of Perov, matrices converges to zero and Fractal-Fractional calculus in Caputo sense, which will be essential in the next sections.
We begin with defining on M m×n (R + ) the partial order relation as follow: and we write I n for the identity n × n matrix and O n for the zero n × n matrix.Definition 2.1.Let E be a vector space over The pair (E, ∥ • ∥ G ) is called a generalized normed space.Moreover, (E, ∥ • ∥ G ) is called a generalized Banach space (in short, GBS), if the vector-valued metric space generated by its vector-valued metric δ G (x, y) = ∥x − y∥ G is complete.
Let (E, ∥ • ∥ G ) be a generalized Banach space.In the rest of this article for r = (r 1 , for the open ball centered at ϑ 0 with radius r, and by: for the closed ball centered at ϑ 0 with radius r.If ϑ 0 = 0 we simply denote B r = B(0, r) and Br = B(0, r).Finally, we respectively denote by K and co(K) for the closure and the convex hull of a subset K of E.
The spectral radius of Υ is strictly less than 1.Definition 2.4.Let (E, δ G ) be a generalized metric space and N be an operator from E into itself.N is called Υ-contraction with matrix Υ ∈ M n×n (R + ) that is converges to O n , if for all ϱ, v ∈ E we have In the following, an extension of the Banach contraction principle by Perov is given.Theorem 2.5.[31] Let E be a complete generalized metric space and let N : E −→ E be an M-contraction operator.Then, N has a unique fixed point in E.
Next, we give some important concepts from fractal-fractional calculus in Caputo sense.We refer the reader for the reference [32] for more details.Definition 2.6.Let ϱ be differentiable in opened interval (a, b), if ϱ is fractal differentiable on (a, b) with order β, then the FF-derivative of ϱ of order α in the Caputo sense with power law is given as: where dϱ(t) dt β = lim s→t ϱ(s) − ϱ(t) s β − t β .Lemma 2.7.The Eq (2.1) can be written as follows:

Existence and uniqueness results
Lemma 3.1.(x, y, z) is a solution of the fractal-fractional-order system (1.3), if and only if it is a solution of the following problem Theorem 3.2.Suppose that there is a vector with positive entries fulfills in addition, if the matrix 3) converges to O 3 , where H(α i , β i ) denotes the beta function of α i and β i , then the system (3.1) has a unique solution in the space + where r satisfies the above inequality in (3.2).We recall that the space .
The proof will be broken up into several steps.
Step 1: First, we shall show that the mapping is G-contraction where N defined by the following formula:

And we have for all
by taking the supermum over t we find And we have for all positive numbers ϱ 1 , ϱ 2 and γ ≥ 1 where Step 2: Our objective here is to prove that the operator N maps K into itself.To do so, let X = (x, y, z), ∈ K and for t ∈ [0, T ] we have By the same manner, we find Hence, .
By using Perov fixed point Theorem 2.5, we conclude that the system (3.1) has a unique solution in K. □
Using the GL approximation and the NSFD framework [27], we discretize the first equation in (1.3) as follows: From the first equation in (4.1), we have: Looking that the function F 1 can be written as next: By substituting this latter in (4.2), and using the fact that the nonlinear term nonlocal way, we find that: .
Repeating the same procedure to the second and the third equation of the system (1.3), we conclude that the discretization of system (1.3) using GL-NSFDM can be formulated as follows:

Lyapunov's stability of the incommensurate fractional-order cancer model
In this section we analyze the dynamics of the incommensurate by taking the initial conditions (x 0 , y 0 , z 0 ) = (0.4,0.5, 0.5), (β 1 , β 2 , β 3 ) = (1, 1, 1) and by selecting different values of the fractional-orders α 1 , α 2 , α 3 and varying the set of parameter (a, b, c).Definition 5.1.[36] The equilibrium point E is called a saddle point of index one (two) if the Jacobian matrix evaluated at point E has exactly one (two) eigenvalue with non-negative real part.Scrolls are generally created only around the saddle points of index two.
The Following Lemma gives the sufficient condition to exhibit the equilibrium point E a stability nature.Lemma 5.2.[37] The equilibrium point E of the fractional-order system is locally asymptotically stable in the Lyapunov sense if the following condition is satisfied: where ∆(λ) = det(J − diag(λ δα 1 , ..., λ δα k )) and J = ( ∂F i ∂x j ) i, j=1,...,k is the Jacobian matrix evaluated at E.
The parameter δ is the least common multiple of the denominators q i s of α i s, where α i = p i q i , (p i , q i ) = 1, If the condition (5.1) does not satisfy, we are in the following state.
The Jacobian matrix of system (1.3) is .
In [39], the authors established that the system (1.3) has five real equilibrium points, where four of them are obtained analytically and can be described as follows: 1) E 0 = (0, 0, 0), , Because they have negative coordinates, the equilibrium points E 1 , E 2 , and E 3 are irrelevant to the ensuing dynamics (negative populations are not defined and, consequently, the dynamics must take place in the positive octant).The equilibrium point E 0 relates to a situation in which there is no cell at all.The fifth equilibrium point changes according to the set of parameters (a, b, c).The following Table 1 gives the index of saddle points (ISP), and the IMFOSs of three sets selected parameters and different fractional-orders.
Table 1.The IMFOSs and the index of saddle points (ISP) of the incommensurate fractalfractional-order cancer model (1.3) for different fractional-orders and system parameters, and According to Table 1 the IMFOSs for the equilibrium points E * and E * * are negative numbers, which implies that E * and E * * are stables.Therefore, for the given derivative orders, the systems in case 1 and case 2 do not have the necessary condition to exhibit chaos.Numerical simulation results in Figures 1  and 2, respectively, confirm this conclusion.In the third case, Table 1 shows that the IMFOS is non-negative number, and the equilibrium point E * * * is a saddle point of index 2.This implies that the system (1.3) in case 3 in The Table 1 satisfies the necessary condition for exhibiting a 1-scroll attractor.As shown in Figure 3, numerical simulation results confirm this conclusion.

Ulam-Hyers stability of the incommensurate fractal-fractional-order cancer model
Here, we are going to demonstrate the stability of Ulam-Hyers sense of the proposed model.We adopt the following definitions from [40].Definition 6.1.Let (X, d G ) be a generalized metric space and F : X → X be an operator.Then, the fixed point equation is said to be generalized Ulam-Hyers stable if there exists an increasing function ψ : R m + → R m + , continuous in 0 R m with ψ(0) = 0, such that, for any ε := (ε 1 , . . ., ε m ) with ε i > 0 for i ∈ {1, . . ., m} and any solution there exists a solution X * of (6.1) such that Lemma 6.2.The solution of the perturbed model fulfills the relation given below Proof.The solution of (6.3) is given by Then, we have sup Repeating the same procedure to the second and the third equations of the system (6.3),we have sup Hence, the proof is completed.□ Theorem 6.3.If the matrix Θ (3.3) converges to O 3 , then (1.3) is generalized Ulam-Hyers stable.
Proof.Let X = (x, y, z) be any solution of the inequality (6.4), and let X * = (x * , y * , z * ) be the unique solution of (1.3), then By the same manner, we find Since the matrix Θ converges to zero, then we have Hence, the solution of the proposed problem is generalized Ulam-Hyers stable.□ Example 6.4.Consider the following fractal-fractional-order cancer model Note that for β = (1, 1, 1) the system (6.6) was stated in the third case in Table 1 and exhibit a chaotic behavior, as it have shown in Figure 3.By taking Υ = (1, 1, 1) we find Furthermore, Then, the eigenvalues of matrix Θ are as follows: hence, the spectral radius of Θ is ρ = 1.6069 × 10 −2 < 1.As a result, the system (6.6) has a unique solution that is generalized Ulam-Hyers stable according Theorem 3.2 and Theorem 6.3.Figure 4 illustrates the conclusion.

Discussion
This section is devoted to numerical simulations of the proposed model under investigation in the present paper.As described in Table 1 and Example 6.4, the approximate solutions of the fractal-fractional system (1.3) are given in Figures 1-4 with varying values of fractional-order parameters (α 1 , α 2 , α 3 ), the fractal dimension (β 1 , β 2 , β 3 ) and the set parameters (a, b, c).We briefly presented the simulation of this model using the GL-NSFDM numerical method as given in (4.3).The time interval is [0, 3000] and ϕ(t) = exp(t) − 1. MATLAB computer language was used to accomplish all computations in this work.
By maintaining the fractal dimension (β 1 , β 2 , β 3 ) = (1, 1, 1), the phase plots in Figures 1 and 2 presented above indicated that the system exhibits a non-chaotic behavior, counter to Figure 3, where the phase portrait shows a chaotic behavior.By varying the fractal dimension of the previous case to β = (0.88, 0.8, 0.9), the fractal-fractional cancer model widths a stable behavior in Figure 4.This demonstrates how the fractal dimension -which is absent in both the fractional and classical models-can turn the behavior of the solutions from chaotic into a stable state and vice versa.By returning to the proposed model, we observe that the parameters (a, b, c) are constants.However, from the biological point of view, these parameters may show a randomness behavior; that is a limitation of our study.In future research, we are focusing on replacing these parameters with Ornstein-Uhlenbeck process to make the model more realistic.

Conclusions
A mathematical study of the growth of tumor has been discussed in this paper.The contribution is based on describing the cancer process by a novel fractal-fractional order model.This is inspired by population dynamics and contains terms that refer to tumor cells, effector immune cells and healthy tissue cells.The study in [16] is a special case from the present paper.Perov's fixed point theorem showed the existence and the uniqueness result.Besides, the numerical simulations were received with the Grünwald-Letnikovv nonstandard finite difference scheme.The dynamics of the proposed Incommensurate fractal-fractional cancer model were analyzed by varying the value of the fractional order, the fractal dimension and the values of the system parameters.The obtained results are also compatible with theoretical analysis.The proposed model could describe a wide range of biologically observed tumor states, including stable and chaotic states.

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

Figure 1 .
Figure 1.Numerical simulation for the system in (1.3) of the case 2 as stated in the Table 1 using the GL-NSFDM scheme.(a): Time behaviors of the three state variables: x(t), y(t) and z(t).(b)-(d): The corresponding projection in xy; yz and xz planes, respectively.(e):Behavior of the model in xyz-plane.

Figure 2 .
Figure 2. Results of the numerical simulation for the system in (1.3) of the case 2 as stated in the Table 1 using the GL-NSFDM scheme.(a): Time behaviors of the three state variables: x(t), y(t) and z(t).(b)-(d): The corresponding projection in xy; yz and xz planes, respectively.(e):Behavior of the model in xyz-plane.

Figure 3 .
Figure 3. Results of the numerical simulation for the system in (1.3) of the case 3 as stated in the Table 1 using the GL-NSFDM scheme.(a): Time behaviors of the three state variables: x(t), y(t) and z(t).(b)-(d): The corresponding projection in xy; yz and xz planes, respectively.(e):Behavior of the model (1.3) in xyz-plane.

Figure 4 .
Figure 4. Numerical simulation for the system in (1.3) of the case 3 as stated in the Table 1 using the GL-NSFDM scheme.(a): Time behaviors of the three state variables: x(t), y(t) and z(t).(b)-(d): The corresponding projection in xy; yz and xz planes, respectively.(e):Behavior of the model in xyz-plane.