Numerical simulation for fractional-order differential system of a Glioblastoma Multiforme and Immune system

In this paper, we present a numerical simulation to study a fractional-order differential system of a glioblastoma multiforme and immune system. This numerical simulation is based on spectral collocation method for tackling the fractional-order differential system of a glioblastoma multiforme and immune system. We introduce new shifted fractional-order Legendre orthogonal functions outputted by Legendre polynomials. Also, we state and derive some corollaries and theorems related to the new shifted fractional order Legendre orthogonal functions. The shifted fractional-order Legendre–Gauss–Radau collocation method is developed to approximate the fractional-order differential system of a glioblastoma multiforme and immune system. The basis of the shifted fractional-order Legendre orthogonal functions is adapted for temporal discretization. The solution of such an equation is approximated as a truncated series of shifted fractional-order Legendre orthogonal functions for temporal variable, and then we evaluate the residuals of the mentioned problem at the shifted fractionalorder Legendre–Gauss–Radau quadrature points. The accuracy of the novel method is demonstrated with several test problems.

A large amount of work on modeling biological systems has been restricted to fractionalorder ordinary differential equations [9][10][11][12]. In the last decades, fractional calculus theory [13] has been developed rapidly. It has been applied in several scientific areas such as engineering, viscoelasticity, physics, diffusion processes, rheology, etc. [14][15][16][17]. Actually, fractional calculus theory is considered as a substitution model to classical calculus theory. Fractional differential equations including those of fixed, variable, distributed, tempered order have been used as perfect models for real world problems including earthquake analysis, viscoelasticity, bio-chemical, controller design, electric circuits, diffusion processes, signal processing, etc. [18][19][20]. Therefore, the urgent necessity to find the exact solutions or merely the approximate ones to these problems has emerged. Since the finding of the exact solutions is not possible for most fractional differential equations [21]; the numerical methods have been developed to obtain the approximate solutions to them.
Some local numerical techniques have been introduced for solving fractional-order biological systems [22,23], and this method may become computationally heavy due to the nonlocal property of fractional differential operator. Moreover, the local methods listed the approximate solution at specific points, while the global methods give the approximate solution in the whole mentioned interval. Hence, the global behavior of the solution can be naturally taken into account. In many areas of science such as engineering, biology, economics, physics, and others, several high-order numerical methods have been developed to deal with the related problems. Among these algoriths, spectral methods [24,25] have been rapidly developed in the last four decades. Spectral methods [26,27] are widely applied as powerful techniques in the construction of numerical solution for differential and fractional differential equations.
The spectral collocation method, a global numerical technique [28][29][30], is a particular kind of famous spectral methods that is widely applicable for almost all types of differential equations. The convergence speed is one of the major advantages of the spectral method. Spectral methods have exponential convergence rates as well as high accuracy level. Thus, the spectral collocation method is a more reliable, suitable, and accurate technique for treating such problems. In many areas of science, such as engineering, biology, economics, physics, and others, several high-order numerical methods have been developed to deal with the related problems. Recently, spectral methods have been known as efficient and highly accurate schemes. The spectral method is classified into four kinds, namely collocation [28][29][30], Tau [31], Galerkin [32], and Petrov-Galerkin [33] methods.
Here, a shifted fractional-order Legendre-Gauss-Radau collocation (SFL-GR-C) method is developed to approximate the fractional-order differential system of a glioblastoma multiforme and immune system (GBM-IS). In this paper, we want to numerically solve the fractional-order differential system of GBM-IS [34][35][36] Here, the temporal derivative is taken in the Caputo sense with fractional order μ, 0 < μ < 1. Where the parameters β 1 , β 2 , β 3 , τ 1 , τ 2 , γ , δ, , p, d 1 , d 2 , d 3 , d 4 , k 1 , k 2 , k 3 , r 1 , r 2 , and r 3 are positive constants. S(t), R(t), N(t), and Z(t) are used for sensitive tumor cell, resistant tumor cell, activated macrophages, and macrophages, respectively, where p equals division rate of the sensitive tumor cell, k 1 and k 2 are the carrying capacities of the sensitive (including necrotic part) and resistant cell, respectively. Parameters β 1 and β 2 are necessary to construct logistic differential equations. r 1 and r 2 are the growth rate of the sensitive cell and the resistant cell, respectively. γ is the converting rate of the sensitive tumor cell to the resistant tumor cell. Parameters d 1 and d 2 are their dead rate caused from drugs, respectively. δ is the conversion rate of macrophages to active macrophages. d 3 is the natural death of active macrophages and d 4 is the natural death of macrophages. r 3 is the growth rate of macrophages and the parameter β 3 is selected in view of logistic differential equations. k 3 gives the carrying capacity of macrophages. is the conversion rate of the activated macrophages to macrophages. τ 1 and τ 2 are the destroying rate caused by the activated macrophages.
The rest of this paper is organized as follows. In the next section, we introduce several analytical and numerical prerequisites. In Sect. 3, we construct a spectral numerical scheme. In Sect. 4, three numerical examples are presented to support our theoretical results and to demonstrate the significant gain in accuracy. Finally, we give a brief conclusion in Sect. 5.
Next, denote by χ and (χ, v) the norm and the inner product of space L 2 [-1, 1]. The set of k (t) is a complete orthogonal system in L 2 [-1, 1] where h i = 2 2i+1 and δ jk is the Dirac function. Thus, for any Let S N [-1, 1] be the set of all polynomials of degree at most N (N ≥ 0). Thus, for any where The norm and the discrete inner product are defined as follows:

Numerical simulations
The aim of this paper is devoted to the development of a new scheme based on spectral collocation method for solving the fractional-order differential system of GBM-IS (1). The accuracy of the proposed methods is proved by convergence analysis. According to the rigorous mathematical proof provided by convergence analysis of fractional Legendre spectral collocation method [37][38][39][40] for differential equations, we can expect the spectral convergence rate for the fractional-order differential system of GBM-IS. Therefore, we can achieve very high accuracy by merely using a small number of fractional Legendre collocation nodes. We will discuss error analysis of our proposed methods to show that spectral convergence is obtained for sufficiently smooth solution functions. In this paper, shifted fractional-order Legendre orthogonal functions (SFOLOF) in the interval [0, L] are outputted of the classical Legendre polynomial. Now, we define the SFOLOF outputted Legendre polynomial. The SFOLOF is offered by Also, we list and derive some facts related to the SFOLOF like orthogonality, corresponding nodes, and Christoffel numbers of the shifted fractional-order Legendre- where h (ε) Proof Based on the orthogonality property of Legendre polynomials, we have Consequently, this yields The corresponding nodes and corresponding Christoffel numbers of the shifted fractional Legendre-Gauss (Gauss-Radau or Gauss-Lobatto) interpolation in the interval [0, L] can be given by where x K,s , and K,s , 0 ≤ s ≤ M, are the nodes and the Christoffel numbers of the standard Legendre-Gauss (Gauss-Radau or Gauss-Lobatto) interpolation in the interval [-1, 1].
Proof Using Legendre-Gauss quadrature property, we have The proof is completed.
The solution of Eq. (1) is approximated as follows: Moreover, the Caputo fractional derivative of any order μ of the shifted fractional-order Legendre orthogonal function is derived. The Caputo fractional derivative of order μ is defined as follows: where m is the ceiling function of μ. The analytical form of (ε) L,j (t) is given by By means of Eq. (15), we find The Caputo fractional derivative of order μ of the approximate solution is then estimated as follows: Spectral collocation technique is addressed to solve the previous problem subject to the initial conditions. Firstly, a collocation method based on fractional-order Legendre basis functions is introduced for the fractional-order differential system of GBM-IS with initial conditions. Consequently, we obtain a system of algebraic equations. Based on the above details, we obtain Using the initial conditions, we obtain In the proposed method, the residual of (17) is set to zero at the M points Taking s = 1, 2, . . . , M, we get a system of (4M + 4) algebraic equations which can be solved for the unknown coefficients. The rigor of the novel method is demonstrated with several test problems. We can indicate that our numerical method can also accommodate other methods. For example, the spectral Tau approach for nonsmooth temporal solution may be fallen apart. Employing fractional-order Legendre functions instead of the classical Legendre functions stopped this deterioration. Suppose that D kε Z(t) ∈ C[0, L] for k = 0, 1, . . . , M. If Z M,ε (t) is the best approximation to Z(t) from M , then the error bound is presented as follows: . (20) Proof. Since Z M,ε (t) is the best approximation to Z(t) from M , then by the definition of the best approximation, we have Based on the generalized Taylors formula [41], we obtain V M (t) = M k=0 t kε (kε+1) D kε Z(0 + ). Thus, an upper bound of the absolute error is where E ε = MAX{D kε Z(ξ ), ξ ∈ [0, L]}. Then we conclude that the following is obtained for the approximate solutions. The convergence of the proposed method depends basically on the above error bound.

Numerical results
Using the algorithm presented in the previous section, we give in this section some numerical results. We discuss the fractional-order differential system of a GBM-IS interaction.  Fig. 1 the numerical solution curves of the fractional-order differential system of GBM-IS. Moreover, we obtain the numerical solutions of the fractional-order differential system of GBM-IS as follows:     Fig. 2 the numerical solution curves of the fractional-order differential system of GBM-IS. Moreover, we obtain the numerical solutions of the fractional-order differential system of GBM-IS as follows:   given parameter values, we plot in Fig. 3 the numerical solution curves of the fractionalorder differential system of GBM-IS. Moreover, we obtain the numerical solutions of the Figure 3 Curves of numerical solutions of the fractional-order differential system of GBM-IS fractional-order differential system of GBM-IS as follows:

Conclusions
This paper adopted a fully fractional-order Legendre collocation method to study a fractional-order version of the differential system of GBM-IS. The powerful numerical scheme gave way to a number of impressive numerical results that prove high efficiency of the algorithm. The results of the algorithm pave way to conduct further additional research in this field to display additional results in the future. Also, we provided numerical simulations exhibiting dynamical behavior and stability around equilibria of the system, and on the basis of simulations and comparison with the experimental data, it is concluded that the fractional-order version of the differential system of GBM-IS is a better representative of the system than its integer-order form. In fact, our approach has three advantages: (i) lower computational cost, (ii) ease of implementation, and (iii) exponential accuracy.