Qualitative Analysis of Multi-Terms Fractional Order Delay Differential Equations via the Topological Degree Theory

: With the help of the topological degree theory in this manuscript, we develop qualitative theory for a class of multi-terms fractional order differential equations (FODEs) with proportional delay using the Caputo derivative. In the same line, we will also study various forms of Ulam stability results. To clarify our theocratical analysis, we provide three different pertinent examples.


Introduction
Fractional calculus is the fastest growing area for research in the last three decades. Nowadays it has become an important tool due to its wide range of applications in various scientific disciplines such as biology, chemistry, physics, dynamical systems, electrodynamics, etc. (see [1][2][3][4][5][6] and references therein). Its importance can also be explored in other fields like fluid dynamics traffic models, oscillation due to earthquakes, flow in porous media due to seepage, etc. Therefore such problems have been considered from different aspects to check whether the problem of a differential or integral equation that is to be investigated has a solution or not.
To guarantee the answer, the existence theory is used to find the conditions under which the problem under investigation has a solution or not. Therefore, existence and uniqueness are the important aspects of differential equations (DEs), which have been studied very openly by different authors using various approaches (for example, see [7,8], etc.). Classical fixed point theory has been utilized to study existence and uniqueness for certain problems [9,10]. Using these results, one needs to establish strong compact conditions, which shorten the study to some boundary value problems (BVPs). To manage this limitation and to generalize the techniques to greater extent for BVPs, researchers have been looking for a tool of nonlinear analysis. One of the important tools is topological degree theory which needs weak compact conditions instead of strong compact conditions for operation. The suggested method provides very basic criteria for existence results for many problems. Enormous numbers of problems, both linear and nonlinear DEs and FODEs, have been investigated for existence results by researchers. Mawhin [11] applied degree theory to develop appropriate results for the given BVPs to derive existence and uniqueness results: and Isaia [12] has applied the method of degree theory to form some adequate results about existence and uniqueness results to the integral equations given by where under some growth conditions. One of the weighty classes of DEs is known as pantograph equations (PEs), which involve proportional type delay. In the 1960s, the British railways wanted to make the electric locomotive faster. An important construction was the pantograph, which collects current from an overhead wire. Therefore, Ockendon and Tayler studied the motion of the pantograph for an electric locomotive [13]. The above mentioned class of DEs has a large numbers of applications in real-world scientific disciplines such as dynamical systems, quantum mechanics and electrodynamics. Particularly, as mentioned above, the mentioned delay DEs are used to collect current from overhead wire. Therefore several researchers have attempted to develop conditions for existence and uniqueness of solution to the aforesaid DEs. Many authors also have considered delay DEs using analytical and numerical techniques [14][15][16].
Since most nonlinear problems cannot be solved for exact solutions, we need powerful numerical or analytical techniques. For good numerical results one needs stable algorithms and methods. For such needs, the stability theory was founded. This aspect is important in numerical study and optimization procedures. In the literature, there are different type of stability such as exponential, Mittag-Leffler and Lyapunov type. These stabilities were studied for DEs of ordinary order. In the past few years stability results have been generalized for linear and nonlinear FODEs, (for details, see [17,18]). To establish these stabilities for DEs, some of them need a pre-defined Lyapunov function which is sometimes very difficult and also needs much time. On the other hand the exponential and Mittag-Leffler stability involving exponential functions have difficulties during numerical analysis of the problems. To handle these difficulties Ulam [19] in 1940 introduced another kind of stability, known as Ulam-Hyer's (UH) stability which was further studied by Hyer [20] in 1941. For the first time Wang [21] studied the UH stability for the impulsive ordinary DEs in 2012. UH stability for DEs of different orders have been studied by different authors (see [22,23]). Further we state that the stability analysis is one of the basic problems in the fields of systems and signal processing and control.
Since then, the evolution of a physical system in time has been described using initial value problems. However, this is less informative. Therefore, to get more and better information, the initial (local) conditions are replaced nonlocal conditions. In fact, nonlocal conditions give a better effect as compared to local initial conditions and also the measurement due to nonlocal conditions is usually more precise than the one measurement produced by local conditions. Therefore investigation of problems under nonlocal initial or boundary conditions is one of the important areas of research in recent times (for detail we refer to [24,25]). Inspired by the aforesaid work, this research aims to study Equation (5) under generalized nonlocal integral boundary condition as: where 0 < κ ≤ 1, for i = 1, 2, . . . , n, λ i ∈ (0, 1) and g ∈ C[J × R, R], f ∈ C[J × R m , R], for existence and uniqueness of solutions using the mentioned method. Moreover, some adequate results of various UH type stabilities such as UH stability, generalized UH (GUH) stability, UH-Rasaias (UHR) stability and generalized UHR (GUHR) stability are established. Finally the analysis is justified by some examples.

Fundamental Material
Here, we provide some fundamental material about fractional calculus, topological degree theory and UH type stability.

Definition 1 ([10]
). If κ ∈ R + , then integral of fractional order for the function w ∈ L 1 (J, R) is expressed as Definition 2 ([10]). THe derivative of fractional order to a function w on the interval J in Caputo sense is expressed as where r = [κ] + 1 and [κ] represents an integral part of κ.
is given as
In the following Y = C[J, R] will be Banach space with norm w = sup{|w(t)|, t ∈ J} and the family W ⊂ P (Y ) represents all its bounded sets. Below are some notions and results recalled from [12]. Definition 3. The mapping χ : W → R + for Kuratowski measure of non-compactness is defined as where W ∈ W is covered by finite sets with diameter ≤ . Proposition 1. The mapping χ due to non-compactness enjoys the properties given below: Proposition 2. Let Λ, Π : Φ → Y be χ-Lipschitz for constants K and K , respectively, and Λ + Π : Φ → Y also be χ-Lipschitz for constants K + K .
Consequently, Λ has at least one fixed point and the set of the fixed points of Λ lies in W r (0).
The following definitions are recalled from [26]. (5) is UH stable if for every > 0, ∃ C q ∈ R + and w ∈ Y is any solution of

Definition 7. The Equation
there is a unique solutionw ∈ Y of (5), such that |w −w| ≤ C q ζ(t), ∀ t ∈ J. (5) is GUHR stable with respect to ζ ∈ C[J, R + ] if ∃ C q,ζ ∈ R + such that for every > 0 and for any solution w ∈ Y of the inequality (14) there is a unique solutionw ∈ Y of (5), such that |w −w| ≤ C q,ζ ζ(t), ∀ t ∈ J.

Main Results
In this section we study existence results for nonlinear delay FODEs under integral boundary condition, we use J = [0, θ] and Z = C(J).
The above three conditions will be used to show the existence and uniqueness of the solution of (5).
The fixed point of M will insure the existence of the solution of (5).

Proposition 5.
The map A : Z → Z is Lipschitz with constant L. Consequently it is χ-Lipchitz with constant L.

Proof. Consider
and K 1 be such that L < 1.
for every u, w ∈ Z. Thus A is χ-Lipschitz with constant L, proposition (4). Using condition (i), A obeys the result given by: for every w ∈ Z.
Proof. To prove B is a continuous, let {w n } ⊂ Z, w ∈ Z be such that w n − w → 0 as n → ∞.
We must have to show Bw n − Bw → 0 as n → ∞. For > 0, ∃ K ≥ 0 such that w n ≤ K, ∀ n ∈ N, w ≤ K.

Existence Criteria
In this part of our paper we derive results for the existence and uniqueness of the solution to the considered problem.
Proof. The operators A, B, M : Z → Z are bounded and continuous. Moreover, A is χ-Lipschitz with constant L ∈ [0, 1) and B is χ-Lipschitz with zero constant (Propositions 4 and 3). M is a strict χ-contraction with constant L (Proposition 2). Consider the set Θ is bounded in Z. Take w ∈ Θ and µ ∈ [0, 1) such that w = µMw then Thus Θ is bounded in Z for c 1 < 1, c 2 < 1. Therefore, Theorem 1 guarantees that M possesses at least one fixed point, and the set of the fixed points of M is bounded in Z. Hence the considered problem has at least one solution.
Proof. Thanks to the Banach fixed point theorem for u, w ∈ Z, take Hence, problem (19) possesses at most one solution. Consequently our considered problem (5) has at most one solution.

Stability
Now we provide stability results for the problem (5). Here we say that the goal of stability analysis of time delay problems/systems is to find the region in the delay parameter space where the considered problem/system is still stable. In fact, in dynamical problems, we search for a fixed point also called an equilibrium point and its stability. Therefore, investigating UH stability and its different kinds, we do not need an exact equilibrium point (exact solution) but there exists a close exact solution (fixed point) when the system is UH or UHR stable. First we provide a lemma which will help in establishing stability analysis.

Lemma 3. For the perturb problem
the following holds , t ∈ J.
Proof. By Lemma 2 the solution of perturb problem (21) is given by By Remark 1, we get .
Proof. Similarly we can prove Lemma 3.
Proof. Let w ∈ Y be a unique solution of (5) andw be any solution of (21), then Hence the problem (5) is UHR stable. Now using the above inequality with Thus problem (5) is GUHR stable.

Application of Aforesaid Analysis
Here we present some applications for our analysis. ds.
For uniqueness, if f satisfies condition (iv) with L f = 1 25 , then Hence Theorem 3 guarantees that problem (22) has a unique solution.
Example 2. Now, we discuss the same analysis for the problem given below: c D (s + w(s))ds.
Hence Theorem 3 guarantees that problem (23) has a unique solution.