## 1 Introduction

Fuzzy set theory was introduced by Zadeh [34]. Kandel and Byatt [23] studied fuzzy differential equations, where other studies about fuzzy differential equations and the concept of the strongly generalized derivative of higher order fuzzy differential equation can be shown in [6, 7, 21, 24]. Many terms we use randomly in daily life often have a fuzzy structure. Verbal and numerical expressions that we use when describing something, explaining an event, commanding and in many other cases contain fuzziness. As for the fuzzy set, it has two basic features. The first is the modeling of systems whose mathematical model is uncertain or whose behavior can be estimated approximately. The second can decide when there is incomplete and uncertain information. Because of these features, new mathematical concepts have emerged, new research issues have emerged and engineering practices have been designed. Especially interesting applications in the field of medicine and artificial intelligence began to emerge. Nowadays, diseases can now be expressed and solved mathematically. This leads to better ways to treat diseases. Fuzzy numbers can give closer results than classical mathematics. For this reason, in our study, we have more accurate results by using these numbers and to be more effective in the diagnosis and treatment of the disease. One of these diseases is cancer that is one of the greatest killers in the world. The control of tumor growth requires special attention [17] and interdisciplinary studies, like biology, medicine and mathematics, many fields are attracted by modeling the spread of this disease. The typical approach for treating GBM(Glioblastoma Multiforme) involves surgical resection, that is followed by radiation treatment and chemotherapy [20]. Works about modeling of multi subpopulations can be shown in [7, 9, 15, 27, 29, 31]. Differential equations have high importance in biological modeling. In the last years and through the use of different type of models, F. Bozkurt has used a system of differential equations as a model for the brain tumor, and its interaction with the immune system(IS) of the human body [10, 11]. Also, another interesting model is the GBM model constructed by F.Bozkurt [12], which explains the interaction between cancer and the human body. Furthermore, in [2], Ö. Akın and Ö. Oruç have studied a predator-prey model with fuzzy initial values. Then, by considering a second-order initial values Ö. Akın, T. Khaniyev, Ö. Oruç and I. B. Turksen [3] have generalized the model. Works about the prey-predator model of fuzzy numbers can be seen in [1, 26, 28]. Lastly, Benli and Keskin [8] considered a model which has a predator-prey structure between the monoclonal tumor and the macrophages. In this study, a new model upon the work of F. Bozkurt [12] is modeled, where fuzzy initial values are used. Furthermore, Allee threshold functions are embedded in the system to investigate the threshold effect on the system and to explain the extinction conditions. A novel multi-resolution method proposed by Seyedi [35] for solving partial differential equations. He tested this method for the solution of well-known viscous Burger’s equation and the obtained results showed superior accuracy in comparison to the finite difference and boundary element methods. Owing to the importance of studying some of equations, several distinct techniques have been proposed that some of them are: 3D-block-pulse functions method [36], homotopy method [37, 38], Hirota bilinear method [39], and Dynamic buckling in lamina plates and dynamic analysis [40, 41].

The paper is constructed as follows: The model and the preliminary definitions are in Section 2., in Section 3 nonlinear fuzzy differential equations are given, a numerical study of the model are given in Section 4. The last Section we have the conclusion part to summarize the study in this paper.

## 2 The Model and Preliminary Definitions

The model is constructed as follows:

where t≥0 denotes the time and E is an Allee constant, the parameters *α*_{1}, *β*_{1}, *τ*_{1}, *τ*_{2}, *p*, *d*_{1}, *d*_{2}, *K*_{1}, *K*_{2}, *r*_{1} and *r*_{2} (see Table 1) are positive numbers. x(t) is used for the GB which represents the monoclonal brain tumor, y(t) denotes the activated macrophages in the system [18]. Differential equations are indispensable for modeling real- world phenomena, unfortunately, every time uncertainty can intervene with real worlds problems the uncertainty can arise from deficient data, measurement errors or when determining initial conditions. Fuzzy set theory is a powerful tool to overcome these problems. Followings are some definitions that are needed to have a fundamental idea of the work.

Values of the parameter of system (1).

E | Allee constant | 1.23 |

r_{1} | The growth rate of the macrophage | 0.55 |

r_{2} | The growth rate of the macrophage | 0.5 |

p | Division rate of the sensitive cells | 0.148 |

K_{1} | Carriying capacity of the tumor cells | 4.704 |

K_{2} | Carriying capacity of the macrophages | 1.132 |

α_{1} | Logistic population rate of tumor cell population | α_{1}ϵ [0.5; 0.95] |

τ_{i} | Destroying rate caused from the interaction | 0.01 |

d_{1} | Causes of drug treatment to the tumor cells | 0.55 |

d_{2} | Causes of drug treatment to the macrophages | 0.06 |

β_{1} | Logistic population rate of macrophages | β_{1}ϵ [0.05; 0.25] |

A fuzzy set A in a universe set X is a mapping *A* : *X* → [0, 1]. We think of A as assigning to each element *x* ∈ *X* a degree of membership, 0 ≤ *A*(*x*) ≤ 1. Let us denote by 𝓕 the class of fuzzy subsets of the real numbers, *A* : *X* → [0, 1] satisfying the following properties:

- A is a convex fuzzy set, i.e.
*A*(*rλ*+ (1 –*λ*)*s*) ≥ min [*A*(*r*),*A*(*s*)],*λ*∈ [0, 1] and*r*,*s*∈*X*; - A is normal, i.e. ∃
*x*_{0}∈*X*with*A*(*x*_{0}) = 1; - A is upper semicontinuous, i.e.
$\begin{array}{}A\left({x}_{0}\right)\ge {lim}_{x\to {x}_{0}^{\mp}}A(x);\end{array}$ - [
*A*]^{0}= =*sup**p*(*A*){ is compact, where*x*∈*R*|*A*(*x*) ≥ 0} denotes the closure of A.*A*

Then 𝓕 is called the space of fuzzy numbers.

If A is a fuzzy set, we define [*A*]^{α} = {*x* ∈ *X*|*μ _{A}* ≥

*α*}, the

*α*–level (cut) of A, with 0 <

*α*≤ 1. For

*u*,

*v*

*ϵ*𝓕 and

*λ*∈

*R*the sum

*u*⊕

*v*and the product

*λ*⊙

*u*are defined by [

*u*⊕

*v*]

^{α}= [

*u*]

^{α}+ [

*v*]

^{α},

*λ*⊙

*u*=

*λ*[

*u*]

^{α}, ∀

*α*∈ [0, 1]. Additionally,

*u*⊕

*v*=

*v*⊕

*u*,

*λ*⊙

*u*=

*u*⊙

*λ*. Also if

*u*∈ 𝓕 the

*α*–cut of

*u*, denoted by [

*u*]

^{α}= [

__,__

*u*^{α}*u*

^{α}

*α*∈ [0, 1].

Let *D* : 𝓕*x*𝓕 → *R*_{+} ∪ {0}, *D*(*u*, *v*) = *sup*_{α∈[0,1]} *max*{|__ u^{α}__,

__|, |__

*v*^{α}*u*

^{α}

*v*

^{α}

*u*]

^{α}= [

__,__

*u*^{α}*u*

^{α}

*v*]

^{α}= [

__,__

*v*^{α}*v*

^{α}

The following properties are well- known [18, 33].

and (𝓕, *D*) is a complete metric space.

(H- Difference) Let ∀ *u*, *v* ∈ 𝓕. If there exists *ω* ∈ 𝓕 such that *u* = *v* ⊕ *ω*, then *ω* is called the H- difference of *u* and *v* and is denoted by *u* ⊖ *v*.

(Hukuhara Derivative) [29] Consider a fuzzy mapping *F* : (*a*, *b*) → 𝓕 and *t*_{0} ∈ (*a*, *b*). We say that *F* is H-differentiable at *t*_{0} ∈ (*a*, *b*) if there exists an element *F*′(*t*_{0}) ∈ 𝓕 such that for all *h* > 0 sufficiently small *F*(*t*_{0} + *h*) ⊖ *F*(*t*_{0}), *F*(*t*_{0}) ⊖ *F*(*t*_{0} – *h*), and the limits (in the metric *D*)

exist and are equal to *F*′(*t*_{0}).

Note that the definition of the derivative above is restrictive; for instance in [6, 7] it is shown that, if *F*(*t*) = *c.g*(*t*) where c is a fuzzy number and *g* : [*a*, *b*] → *R*^{+} is a function with *g*′(*t*) < 0, then F is not differentiable. To avoid this difficulty, a more general definition of the derivative for fuzzy mappings are given in [6, 7] which is as following:

(Generalized Fuzzy Derivative) [6, 7] Let *F* : (*a*, *b*) → 𝓕 and *t*_{0} ∈ (*a*, *b*). We say that *F* is strongly generalized differentiable at *t*_{0} if there exists an element *F*′(*t*_{0}) ∈ 𝓕 such that

- for
*h*> 0 sufficiently small ∃*F*(*t*_{0}+*h*) ⊖*F*(*t*_{0}),*F*(*t*_{0}) ⊖*F*(*t*_{0}–*h*) and the limits satisfy$$\begin{array}{}{\displaystyle \underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}+h)\ominus F({t}_{0})}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0})\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0}-h)}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$ - for
*h*> 0 sufficiently small ∃*F*(*t*_{0}) ⊖*F*(*t*_{0}+*h*),*F*(*t*_{0}–*h*) ⊖*F*(*t*_{0}) and the limits satisfyor$$\begin{array}{}{\displaystyle \underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0})\ominus F({t}_{0}+h)}{\left(-h\right)}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}-h)\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0})}{\left(-h\right)}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$ - for
*h*> 0 sufficiently small ∃*F*(*t*_{0}+*h*) ⊖*F*(*t*_{0}),*F*(*t*_{0}–*h*) ⊖*F*(*t*_{0}), and the limits satisfyor$$\begin{array}{}{\displaystyle \underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}+h)\ominus F({t}_{0})}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}-h)\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0})}{\left(-h\right)}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$ - for
*h*> 0 sufficiently small ∃*F*(*t*_{0}) ⊖*F*(*t*_{0}+*h*),*F*(*t*_{0}) ⊖*F*(*t*_{0}–*h*) and the limits satisfy$$\begin{array}{}{\displaystyle \underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0})\ominus F({t}_{0}+h)}{\left(-h\right)}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}-h)\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0})}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$

Definition 5 is equivalent to the Definition 6 that we will use in the paper.

Let *F* : (*a*, *b*) → 𝓕 and *t*_{0} ∈ (*a*, *b*).

- for
*h*> 0 sufficiently small, ∃*F*(*t*_{0}+*h*) ⊖*F*(*t*_{0}),*F*(*t*_{0}) ⊖*F*(*t*_{0}–*h*) andor$$\begin{array}{}{\displaystyle \underset{h\to {0}^{+}}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}+h)\ominus F({t}_{0})}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to {0}^{+}}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0})\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0}-h)}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$ - for
*h*> 0 sufficiently small ∃*F*(*t*_{0}+*h*) ⊖*F*(*t*_{0}),*F*(*t*_{0}) ⊖*F*(*t*_{0}–*h*) and$$\begin{array}{}{\displaystyle \underset{h\to {0}^{-}}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0}+h)\ominus F({t}_{0})}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}\underset{h\to {0}^{-}}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{F({t}_{0})\phantom{\rule{thinmathspace}{0ex}}\ominus \phantom{\rule{thinmathspace}{0ex}}F({t}_{0}-h)}{h}\phantom{\rule{thinmathspace}{0ex}}=\phantom{\rule{thinmathspace}{0ex}}{F}^{\prime}({t}_{0}),\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}}\end{array}$$

The following theorem is very important for us to solve fuzzy differential equations.

*[14*, *22] Let* *F* : *T* → 𝓕 *be a function and set* [*F*(*t*)]^{α} = [*f*_{α}(*t*), *g*_{α}(*t*)] *for each* *α* ∈ [0, 1]. *Then*

*If**F**is differentiable following the form (1) in Definition 6 then**f*_{α}(*t*)*and**g*_{α}(*t*)*are differentiable functions and*[*F*^{′}(*t*)]^{α}= [*f*^{′}(*t*),*g*^{′}(*t*)].*If**F**is differentiable following the form (2) in Definition 6 then**f*_{α}(*t*)*and**g*_{α}(*t*)*are differentiable functions and*[*F*^{′}(*t*)]^{α}= [*g*^{′}_{α}(*t*),*f*^{′}_{α}(*t*)].

## 3 Solving Fuzzy Differential Equations with Fuzzy Initial Values

Consider the following equation with fuzzy initial values

where *F* : [0, *α*] *x* 𝓕 → 𝓕 and *x*_{0} is a fuzzy number, [*x*(*t*)]^{α} = [*u*_{α}(*t*), *v*_{α}(*t*)], *F*(*t*, *x*(*t*))]^{α} = [*f*_{α}(*t*, *u*_{α}(*t*), *v*_{α}(*t*)), *g*_{α}(*t*, *u*_{α}(*t*), *v*_{α}(*t*))].

Then, we get the following alternatives for solving the initial value problem (2):

- If we consider
*x*^{′}(*t*) by using the derivative in the first form (1), then from Theorem 1 [*x*^{′}(*t*)]^{α}= So we have the following equalities:$\begin{array}{}{\displaystyle [{u}_{\alpha}^{{}^{\mathrm{\prime}}}(t),{v}_{\alpha}^{{}^{\mathrm{\prime}}}(t)].}\end{array}$ By solving the above system for$$\begin{array}{}{\displaystyle {u}_{\alpha}^{{}^{\mathrm{\prime}}}(t)={f}_{\alpha}(t,{u}_{\alpha}(t),{v}_{\alpha}(t)),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{\alpha}(0)={u}_{\alpha}^{0}}\\ {\displaystyle {v}_{\alpha}^{{}^{\mathrm{\prime}}}(t)={g}_{\alpha}(t,{u}_{\alpha}(t),{v}_{\alpha}(t)),{v}_{\alpha}(0)={v}_{\alpha}^{0}.}\end{array}$$ *u*_{α}and*v*_{α}, we get the fuzzy solution [*x*(*t*)]^{α}= [*u*_{α}(*t*),*v*_{α}(*t*)]. Finally we ensure that [*x*(*t*)]^{α}= [*u*_{α}(*t*),*v*_{α}(*t*)] and are valid level sets.$\begin{array}{}{\displaystyle [{x}^{{}^{\mathrm{\prime}}}(t){]}^{\alpha}=[{u}_{\alpha}^{{}^{\mathrm{\prime}}}(t),{v}_{\alpha}^{{}^{\mathrm{\prime}}}(t)]}\end{array}$ - If we consider
*x*^{′}(*t*) by using the derivative in the second form (2), then from Theorem 1 [*x*^{′}(*t*)]^{α}= So we get following:Solving the above system for$$\begin{array}{}{\displaystyle {u}_{\alpha}^{{}^{\mathrm{\prime}}}(t)={g}_{\alpha}(t,{u}_{\alpha}(t),{v}_{\alpha}(t)),{u}_{\alpha}(0)={u}_{\alpha}^{0}}\\ {\displaystyle {v}_{\alpha}^{{}^{\mathrm{\prime}}}(t)={f}_{\alpha}(t,{u}_{\alpha}(t),{v}_{\alpha}(t)),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{v}_{\alpha}(0)={v}_{\alpha}^{0}.}\end{array}$$ *u*_{α}and*v*_{α}, we get the fuzzy solution [*x*(*t*)]^{α}= [*u*_{α}(*t*),*v*_{α}(*t*)]. Finally we ensure that [*x*(*t*)]^{α}= [*u*_{α}(*t*),*v*_{α}(*t*)] and are valid level sets.

## 4 Numerical Study with Allee Threshold For A GB-IS Interaction with Fuzzy Initial Values

Important research for population models by the Allee [4], who demonstrated that Allee effect occurs when the population growth rate is reduced at low population size. It is well known that the logistic model assumes that per-capita growth rate declines monotonically when the density increase; however, it is shown that for population subjected to an Allee effect, per-capita growth rate gives a humped curve increasing at low density, up to a maximum intermediate density and then declines again [32]. Many theoretical and laboratory studies have demonstrated the importance of the Allee effect in the dynamics of small populations, see for example [5, 16, 19, 25, 30]. This section presents the numerical study for the following Allee threshold tumor model, which considers fuzzy initial values. Since the critical point x = 0 is unstable, we need a threshold to explain the extinction of the tumor population. Biological facts led to assume the Allee function as follows:

- If
*N*= 0, then*α*(*N*) = 0; that is, there is no reproduction without partners,*α*is the Allee function and N is the populations density at time t, *α*^{′}(*N*)>0 for*N*∈(0,∞); that is, Allee effect decreases as density increases,- lim
_{N→∞}*α*(*N*) = 1; that is, Allee effect vanishes at high density [13].

Considering system (3), we can see that we embed to the monoclonal tumor population on Allee threshold model.

where x(t) and y(t) are the density of the tumor and macrophages at time t, respectively. Table 1 shows the values of the parameters in Equation (1), which gives us (4).

ordering and simplifying system (3), we have

The crisp solutions for the problem (4) are shown in Figure 1.

Let the initial values be fuzzy, that is, *x*(0) =
*y*(0) =
*α*-level sets be as follows;

Let the *α*–level sets of *x*(*t*, *α*) be [*x*(*t*, *α*)]^{α} = [*u*(*t*, *α*), *v*(*t*, *α*)] and for simplicity denote them as [*u*, *v*], similarly [*y*(*t*, *α*)]^{α} = [*r*(*t*, *α*), *s*(*t*, *α*)] = [*r*, *s*]

If *x*(*t*, *α*) and y(*t*, *α*) are (1) differentiable according to Definition 6, system (4) becomes

Hence for *α* = 0 the following initial value problem derives from (4):

*u*(0), *v*(0), *r*(0), *s*(0) are fuzzy initial conditions for the system. Now if *x*(*t*, *α*) and y(*t*, *α*) are (2) differentiable according to Definition 6, system (4) becomes:

In Figure 2, we can see the graphical solution of all cases for *α* = 0.

In Figure 2; according to Definition 6, (a) means that *x*(*t*, *α*) and *y*(*t*, *α*) are (1) differentiable, (b) means that *x*(*t*, *α*) is (1) differentiable and *y*(*t*, *α*) is (2) differentiable, (c) means that *x*(*t*, *α*) is (2) differentiable and *y*(*t*, *α*) is (1) differentiable and (d) means that *x*(*t*, *α*) and *y*(*t*, *α*) are (2) differentiable. Now, if Figure 2 is analyzed, we see that when *x*(*t*, *α*) and *y*(*t*, *α*) are (2) differentiable graphical solution (Fig.(d)) is biologically meaningful. In addition, the graphical solution is compatible with a crisp solution. In contrast, when *x*(*t*, *α*) and *y*(*t*, *α*) are differentiable as in (a), (b) and (c), the graphical solutions are incompatible with biological facts. *x*(*t*) and *y*(*t*) are the density of tumors and macrophages at time t respectively. Being compatible with biological facts means that the tumor is increasing or decreasing and the macrophages are increasing or decreasing. It does not give us information about tumors and macrophages when it is differentiated as in (a), (b) and (c), so these conditions are biologically meaningless.

Now, we will focus on the situation when *x*(*t*, *α*) and *y*(*t*, *α*) are (2) differentiable. When the crisp graphical solution and the fuzzy graphical solution *x*(*t*, *α*) and *y*(*t*, *α*) are (2) differentiable, we will plot their graphs on the same graph for *α* = 0 and *α*∈[0, 1]. The fuzzy solution for *α* = 0 and the crisp solution are given in Table 2 and Figure 3.

Numerical values of Figure 3.

Time(t) | v | x | u | s | y | r |
---|---|---|---|---|---|---|

0.0 | 2.0000 | 1.6000 | 1.2000 | 0.3000 | 0.2000 | 0.10000 |

0.1 | 1.9800 | 1.6800 | 1.3762 | 0.3037 | 0.20950 | 0.11499 |

0.2 | 2.0005 | 1.7767 | 1.5461 | 0.30812 | 0.2194 | 0.13010 |

0.3 | 2.0589 | 1.8931 | 1.7179 | 0.31326 | 0.22972 | 0.14535 |

0.4 | 2.1543 | 2.0329 | 1.8995 | 0.31909 | 0.24046 | 0.16078 |

0.5 | 2.2876 | 2.2001 | 2.0980 | 0.3256 | 0.25164 | 0.17639 |

0.6 | 2.4605 | 2.3992 | 2.3203 | 0.33279 | 0.26324 | 0.19222 |

0.7 | 2.6760 | 2.6350 | 2.5731 | 0.34064 | 0.27528 | 0.20826 |

0.8 | 2.9375 | 2.9123 | 2.8626 | 0.34914 | 0.28775 | 0.22454 |

0.9 | 3.2487 | 3.2361 | 3.1948 | 0.35826 | 0.30065 | 0.24105 |

1.0 | 3.6133 | 3.6107 | 3.5750 | 0.36799 | 0.31397 | 0.25780 |

## 5 Conclusion

In this paper, we consider a model which has a predator-prey structure between the monoclonal tumor and the macrophages. Building upon the work of F. Bozkurt [12], we include fuzzy initial values to study about the interaction of a monoclonal brain tumor and the macrophages to see the extinction conditions for the tumor population. And also by including on Allee threshold function to the model we have seen that by using fuzzy initial values, the uniqueness of the solution is lost. Furthermore, using the strongly generalized derivative biologically we obtain more realistic behavior that explains the interaction phenomena. In Figure 2, the graphical solution is compatible with a crisp solution. Being compatible with biological facts means that the tumor is increasing or decreasing and the macrophages are increasing or decreasing. It does not give us information about tumors and macrophages when it is differentiated as in Figure 2 (a), (b) and (c), so these conditions are biologically meaningless.

**Funding**: This work is supported by the Scientific Research Center at the Erciyes University with the project code FBA-2014-4808.

## References

- [2]↑
Ö. Akın Ö. Oruç A prey predator model with fuzzy initial values Hacettepe Journal of Mathematics and Statistics 41 (2012) 387-395

- [3]↑
Ö. Akın et al. An Algorithm for the solution of second order Fuzzy Initial Value Problem Expert system and Applications 40 (2013) 953-957.

- [4]↑
W. C. Allee Animal Aggregations:A Study in General Sociology University of Chicago Press Chicago (1931).

- [5]↑
M. A. Asmussen Ö. Akın et al. Density-Dependent Selection II. The Allee effect Am. Nat. 114 (1979) 796-809.

- [6]↑
B. Bede S. Gal Generalizations of differentiability of fuzzy number valued functions with applications to fuzzy differential equation Fuzzy Sets and Systems 151 (2005) 581-599.

- [7]↑
B. Bede et al. First order linear fuzzy differential equations under generalized differentiability Inform. Sci. 177 (2007) 1648-1662.

- [8]↑
F. B. Benli Ö. Keskin A Tumor-Macrophage Interaction Model with Fuzzy Initial Values IAIP Proceedings 1648 (2015) Article ID 370007.

- [9]↑
B.G.Birkhuadetal A mathematical model of the development of drug resistance to cancer chemotherapy Europ. J. of Cancer and Clinical Oncology 23 (1987) 14211427.

- [10]↑
F. Bozkurt Modeling a Tumor Growth with Piecewise Constant Arguments Discrete Dynamics in Nature and Society 841764 (2013) 1-8.

- [11]↑
F. Bozkurt Stability analysis of a fractional-order differential equation system of a GBM-IS interaction depending on the density Applied Mathematics and Information Science 8(3)(2014) 1-8.

- [12]↑
F. Bozkurt A Mathematical Model of the brain tumor GBM and IS Interaction International Journal of Mathematics and Computation 22(1) (2014) 58-65.

- [13]↑
C. Çelik H. Merdan O. Duman and Ö. Akın Allee effects on population dynamics with delay Chaos Solitons and Fractals 37 (2008) 65-74.

- [14]↑
Chalco-Cano Y. and Roman-Flores H. On the new solution of fuzzy differantial equations Chaos Solitons and Fractals 38 (2006) 112-119.

- [15]↑
A. J. Coldman J. H. Goldiec A mathematical model for relating the drug sensitivity of tumors to their spontaneous mutation rate Cancer Treat. Rep. 63 (1979) 1727-1731.

- [16]↑
B. Dennis Allee-effect: population growth critical density and chance of extinction Nat. Resource Model 3 (1989) 481-538.

- [17]↑
A. El-Gohary Chaos and optimal control of cancer self-remission and tumor system steady states Chaos Solutions and Fractals 37 (2008) 1305-1316.

- [18]↑
Gal S. G. Approximation theory in fuzzy setting in: G. A. Anastassiou(Ed.) Handbook of Analytic-Computional Methods in Applted Mathematics A Chapman and Hall/CRC Press Boca Raton FL (2000) 617-666.

- [19]↑
Hanski I. Single-species metapopulation dynamics: concepts models and observations Biol. J. Linnean Soc. 42 (1991) 17-38.

- [20]↑
E.C. Holland Glioblastoma multiforme: The Terminator Proc. of the Nat. Acad. of Sci. 97 (2000) 6242-6244.

- [21]↑
E. Hüllermeier An approach to modeling and simulation of uncertain dynamical systems International Journal of Uncertainty Fuzziness and Knowledge-Based Systems 5(2) (1997) 117-137.

- [23]↑
A. Kandel W. J. Byatt Fuzzy Differential Equations Proceedings of the International Conference on Cybernetias and Society (1978) 1213-1216.

- [24]↑
A. Khastan F. Bahrami K. Ivaz New results on multiple solutions for nth-order fuzzy differential equations under generalized differentiability Boundary Value Problems (2009) Article ID 395714.

- [25]↑
R. Lande Extinction thresholds in demographic models of territorial populations American Naturalist 130 (1987) 624-635.

- [26]↑
X. Liu Y. Lou Global Dynamics of a Predator-Prey Model Journal of Mathematical Analysis and Applications 371(1) 323-340 2010.

- [27]↑
J.C. Patetta A mathematical model of drug resistance: Heterogeneous tumors Math. Biosci. 147 (1998) 41-61.

- [28]↑
M. S. Peixoto L. C. Barros R. C. Baasanezi Predator-Prey Model Ecological Modeling 214 39-44 2008.

- [29]↑
J. E. Schmitz A.R. Kansal and S. Torquato A cellular Automaton of Brain Tumor Treatment and Resistance Journal of Theoretical Medicine 4(4) (2002) 223-239.

- [30]↑
Stephan T. Wissel C. Stochastic extinction models discrete in time Ecol. Model 75-76 (1994) 183-192.

- [31]↑
De Vladar H. p. Gonzalez J.A. Dynamic reponse of cancer under the influence of immunological activity and therapy Journal of Theoretical Biology 227 (2004) 335-348.

- [32]↑
Gang Wang Xue-Gong Liang and Fei-Zhi Wang The competitive dynamics of populations subject to an Allee effect Ecological Modeling 124 (1999) 183-192.

- [33]↑
C. Wu Z. Gong On Henstock integral of fuzzy-number-valued functions (I) Fuzzy Sets and Systems 120 (2001) 523-532.

- [35]↑
S. H. Seyyedi Multiresolution Solution of Burgers Equation with B-spline Wavelet Basis arXiv:1812.10117 2018.

- [36]↑
J. Manafian P. Bolghar Numerical solutions of nonlinear three-dimensional Volterra integral-differential equations with 3D-block-pulse functions Mathematical Methods in the Applied Sciences 41(12) (2018) 4867-4876.

- [37]↑
C. T. Sindi J. Manafian An optimal homotopy asymptotic method applied to the nonlinear thin film flow problems International Journal of Numerical Methods for Heat and Fluid Flow (2018)

- [38]↑
M. Dehgha J. Manafian A. Saadatmandi Study of the wave-breaking's qualitative behavior of the Fornberg-Whitham equation via quasi-numeric approaches International Journal of Numerical Methods for Heat and Fluid Flow Vol. 22 No. 5 2012 pp. 537-553.

- [39]↑
M. R. Foroutan J. Manafian A. Ranjbaran Lump solution and its interaction to (3+1)-D potential-YTSF equation Nonlinear Dynamics 92 (2018) 2077-2092.

- [40]↑
M. Cinefra Dynamic buckling in lamina plates reinforced by CNTs in each layer using exact solution International Journal of Hydromechatronics 2018 1(4): 415-426.

- [41]↑
S. K. Panda R. Kolahchi Dynamic analysis in three layered conical shells utilising numerical methods International Journal of Hydromechatronics 2018 1(4): 427-446.