The Solution of Modified Fractional Bergman ’ s Minimal Blood Glucose-Insulin Model

In the present paper, we use analytical techniques to solve fractional nonlinear differential equations systems that arise in Bergman’s minimal model, used to describe blood glucose and insulin metabolism, after intravenous tolerance testing. We also discuss the stability and uniqueness of the solution.


Introduction
Mathematical modeling is a very important branch of applied mathematics.By using this approach, we can convert a real world problem into a mathematical module and then analyze it in a better manner.In the 1990s when researchers were faced with the blood or blood's constituent dynamical transport phenomena, they used classical theory to describe these processes [1][2][3].After that, many researchers conducted research in that direction to remove the complexity of this phenomenon [1][2][3][4][5][6][7].
Regarding research modules , one of the most important modules is Bergman's minimal model [4][5][6].In this model, a body is described as a compartment with a basal concentration of glucose and insulin.Bergman's minimal model has two variations.the first describes glucose kinetics, and the second describes insulin kinetics.The two models have mostly been used to understand the kinetics during IVGTT test (Glucose Tolerance Test) [7,8].
In applied mathematics, one of the most used concepts is the derivative.Derivatives show the rate of change of a function.This is helpful to describe many real phenomena.After this research, mathematicians faced some complex problems of the real world; to solve them, mathematicians introduced the fractional derivative [9][10][11][12][13].The concept of fractional calculus has great importance in many branches and is also important for modeling real world problems [14][15][16][17].
For this reason, many researchers have engaged in a great amount of research work, conferences, and paper publications.Various definitions of fractional derivatives have been given to date.Recently, researchers have described a new fractional derivative operator named the Caputo-Fabrizio fractional derivative [18][19][20][21].In this paper, we use this operator to describe the Bergman's minimal glucoseinsulin model and solve it by the iterative technique.

The Caputo-Fabrizio Fractional Order Derivative
Singularity at the end point of the interval is the main problem that is faced with the definition of the fractional order derivative.To avoid this problem, Caputo and Fabrizio recently proposed a new fractional order derivative that does not have any singularity.The novel fractional derivative given by Caputo and Fabrizio is more suitable to describe the rate of change in concentration of the model because its kernel is non-local and non-singular.The definition is based on the convolution of a first-order derivative and the exponential function, given in the following definition: Then, the new fractional order Caputo derivative is defined as: Here M(α) denotes the normalization function such that M(0) = M(1) = 1; for detail, see [18].If the function does not belong to H 1 (a, b), then the derivative can be written as and As we have defined a new derivative above, then there should be its anti-derivative; the integral of this new fractional derivative is given by Losada and Nieto [19].Definition 2. The fractional integral of order α (0 < α < 1) of the function f is defined below: Remark 2. It is clear from Equation (5) that the fractional integral of order α (0 < α < 1) is an average of function f and its integral of order 1.Hence we get the condition [19]: the above term yields an explicit formula, Due to the above relation, Nieto and Losada [19] anticipated that the new Caputo derivative of order 0 < α < 1 could be written as: Theorem 1.Here f (t) denotes the normalization function such as then, we have For more detail see [18,19].

Bergman's Minimal Model Fractional Module
The minimal model of the glucose insulin kinetics has been proposed to describe the time course of these concentrations.We will use the standard formulation of the minimal model represented by the following system of differential equations: subject to initial conditions, The parameters for the minimal model ( 11) are given in Table 1.

Parameter Unit Description
The effect of active insulin This model can be used to describe the pancreas as the source of insulin.In a healthy individual, a small amount of insulin is always created and cleared [4].This helps to keep the basal concentration I b .The glucose-independent production and clearance of insulin is proportional to the blood insulin concentration.If the insulin level is above basal concentration, clearance increases.On the other hand, if the insulin level is below basal concentration, production increases.When the glucose level gets high, the pancreas reacts by releasing more insulin at a given rate.To explain this mathematically, one has to derive a function describing the reaction of the pancreas.This function was derived by Bergman et al. and adjusted by Gaetano et al. [7,8]

Existence of the Coupled Solutions
By using the Fixed-Point theorem, we define the existence of the solution.First, transform Equation ( 11) into an integral equation as follows: on using the definition defined by Nieto, we get we also have Let us consider the following kernels: Theorem 2. Show that K 1 , K 2 , and K 3 satisfy Lipschiz condition.
Proof.First we prove this condition for K 1 .Let G and G 1 be two functions, then we have on using the Cauchy's inequality, we get or where Additionally, for on using the Cauchy's inequality, we get where Similarly, for K 3 , by Cauchy's inequality or where We consider the following recursive formula and as well as Now the difference between the consecutive terms is . Now take norm on both sides of Equations ( 37)-(39), respectively as well as From Equations ( 40)-(42), by using triangular inequality Since the kernel satisfies the Lipchitz condition, we obtain: Theorem 3. Show that the Bergman's Minimal Model Fractional Module is the minimal model of the glucose insulin kinetics having a solution.
Proof.As we have seen that the above Equations ( 43)-( 45) are bounded, and we have proven that the kernels satisfy Lipschiz condition, therefore following the results obtained in Equations ( 43)-(45) using the recursive technique, we get the following relation and as well as Therefore, the above solutions exist and are continuous.Nonetheless, to show that the above is a solution of Equation ( 11), we get where P n , Q n and R n are remainder terms of series solution.Thus, as well as It follows from the above that: Now, applying the norm on both sides and using the Lipchitz condition, we get similarly, we get and (56) On taking the limit n → ∞ of Equations ( 54)-(56), we get and Equations ( 57)-( 59) is the solution of the system (11); therefore, we can say that a solution exists.

Uniqueness of the Solutions
In this part, we want to show that solutions presented in the above section are unique.
To prove this, we can another solutions for system (11), say G (t) , X (t) , and I (t); then: apply the norm both sides of Equation (60), (61) On using the Lipchitz condition, having the fact in mind that the solution is bounded, we get this is true for any n; hence, Similarly, we get X (t) = X 1 (t) and I (t) = I 1 (t) .
Hence, it shows the uniqueness of the solution of system (11).

Application of Caputo-Fabrizio Derivative to Bergman's Minimal Model
Watugala introduced the Sumudu transform in early 1990s [20].The Sumudu transform is defined over the set of functions: the Sumudu transform is defined by for detail, see [15,16,21].
Theorem 4. Let f (t) be a function for which the Caputo-Fabrizio exists; then, the Sumudu transform of the Caputo-Frabrizio fractional derivative of f (t) is given as:

Conclusions
This paper is an attempt to describe the existence and uniqueness of the Bergman Minimal Model which is extended by Caputo-Fabrizio fractional derivative in the context of glucose and insulin levels in blood.We obtain the approximate solution of the Model and a numerical solution of the system which shows that effect of time on the concentrations G(t), X(t), and I(t).