On a system of fractional coupled hybrid Hadamard differential equations with terminal conditions

In this manuscript, we study the existence of solutions for a coupled system of nonlinear hybrid differential equations of fractional order involving Hadamard derivative with nonlocal boundary conditions. By using suitable fixed point theorems we establish sufficient conditions for the existence result. An example is provided to illustrate our main result.


Introduction
Differential Equations of Fractional Order (DEFO) gain huge attention among scientists due to the applications which are not possible with integer order ordinary or partial differential equations. In recent days, researchers in various fields pay more attention on the existence of solutions of the aforesaid equations. In this regard, many literature surveys, some notable developments (refer [13,14]), and significant studies are available in the theory of coupled systems with coupled boundary conditions in FDEs. For further references, one can refer to [1, 3, 6-9, 12, 15, 16, 19]. This is because the aforesaid equations include the perturbations of original differential equations in nonidentical ways.
The nonlinear hybrid Hadamard differential equations with initial conditions were discussed by Ahmad and Ntouyas in [4]. In 2015, Wang [22] proved the existence results for the impulsive systems with Hadamard derivative. For more details, the interested readers can refer to Wang et al. [21], where the authors used some classical methods to analyze the Ulam-Hyers stability and proved existence results for Hadamard fractional differential equations. Further, Samko et al. [20] proved the existence and uniqueness of the solution of Cauchy problems for FDE involving Hadamard derivatives in a nonsequential setting. In [17], the author studied the Hadamard-type fractional calculus.
We believe that only few literature studies are available which investigate coupled systems of FDEs of Hadamard type; see [23,24] and the references therein. This paper will definitely fill up the gap between the theories and applications of Hadamard FDEs.
In [4], Ahmad et al. studied In 2015, the authors of [5] investigated the nonlocal boundary value problems for hybrid FDEs and inclusions of Hadamard type as given below: Ahmad et al. [2] in 2017 focused on the issues related to initial and boundary value problems for fractional differential equations of Hadamard type and studied mixed initial value problems involving Hadamard derivatives. The above mentioned works motivated us to study the following coupled Hadamard FDEs with the boundary conditions: where D σ 1 , D δ 1 are the Hadamard fractional derivatives of order σ 1 and δ 1 , respectively; The rest of the paper is organized as follows. In Sect. 2, some basic definitions and lemmas related to Hadamard fractional derivatives are presented. In Sect. 3, we establish sufficient conditions for the existence of the coupled solutions to (3), and in Sect. 4, we present a numerical example to illustrate our main result.

Auxiliary results
In this section we list some basic definitions and lemmas of fractional calculus.
We define the nonhomogeneous boundary value problems of coupled hybrid fractional differential equations of Hadamard type of the function space X = C([1, e], R) of continuous real valued functions f 1 : which is again a Banach algebra with respect to the multiplication "·" defined by The product space E = X × X is a Banach space under the norm Then the normed linear space (E, (·, ·) ) is a Banach space which further becomes a Banach algebra with respect to the multiplication "·" defined by for all t ∈ [1, e], where (ϑ, ν), (ū,v) ∈ X × X = E. The following result concerning this fact of algebraic structure of the product space E = X × X is proved by Dhage in [11]. But first we need where n = [β 1 ] + 1, [β 1 ] denotes the integer part of the real number β 1 and ln(·) = log e (·).

Definition 2.2 ([18])
The Hadamard fractional integral of order β 1 for a function h 1 : provided the integral exists. [11]) Let E be a Banach space. An operator T : E → E is called Lipschitz if there exists a constant L T > 0 such that
Proof Let (ϑ, ν) and (ū,v) be any two elements of E then, by definition of the norm (·, ·) and the multiplication "·" in E, we obtain and Now from (8) and (9), it follows that This shows that (E, (·, ·) ) is a Banach algebra and the proof of lemma is complete.
We employ the following fixed point theorem as a basic tool for proving our main existence result for a coupled system of the proposed problem (3).
Theorem 2.5 (Dhage [10]) Let S be a nonempty subset of a Banach algebra, which is bounded, convex and closed. Assume also that the operators A and B are such that A : X → X and B : S → X, and satisfy the following conditions: Then the operator equation ϑ = AϑBϑ has a solution in S.

Existence result
To prove the main existence result for a mild coupled solution of problem (3), we need to prove the following lemma.
is given by Proof The solution of equation (10) can be written as where the arbitrary constants c 1 , c 2 ∈ R and the boundary conditions given in (10) give Substituting these values of c 1 and c 2 into (12), we get If these conditions hold, then the coupled system ( and Choose Clearly, S is a closed, convex and bounded subset of the Banach space E = X × X. Define the operators A = (A 1 , A 2 ) : E → E and B = (B 1 , Now, the system of equations in (13) and (14) becomes which further, in view of the multiplication (7) of two elements in E, yields This further implies that Now we prove that the operators A and B satisfy the conditions of Theorem 2.5 in a series of steps.
Step I. First we show that A = (A 1 , A 2 ) is Lipschitz on E with constant (L f 1 + L f 2 ). Let (ϑ, ν), (θ,ν) ∈ E be arbitrary. Then, using (H1), we have Taking the supremum over t, we obtain Therefore, by definition of the operator A, we obtain Step II. In this step we prove that B = (B 1 , B 2 ) from S into E is a compact and continuous operator. For continuity of B, take any sequence of points (ϑ n , ν n ) in S converging to a point (ϑ, ν) ∈ S. Then, by Lebesgue dominated convergence theorem, we have for all t ∈ [1, e]. Similarly, we prove for all t ∈ [1, e]. Hence B(ϑ n , ν n ) = (B 1 (ϑ n , ν n ); B 2 (ϑ n , ν n )) converges to B(ϑ, ν) pointwise on [1, e]. This shows that B is continuous. Next, we prove B is a compact operator on S. Let (ϑ, ν) ∈ S be any point. Then, using (H4), we have Taking the supremum for the above inequality, we get Hence B 1 is a uniformly bounded operator on S. In a similar way, we can prove B 2 is also uniformly bounded operator on S. Hence B is a uniformly bounded operator on S. Next to prove that {B(ϑ n , ν n )} is an equicontinuous sequence of functions in E, assume τ 1 , τ 2 ∈ [1, e], provided τ 1 < τ 2 , then uniformly for all (ϑ, ν) ∈ E. Similarly we have uniformly for all (ϑ, ν) ∈ S. This shows that B(S) ⊂ E is uniformly bounded, equicontinuous and compact, in view of Arzela-Ascoli theorem. Consequently, B is a compact and continuous operator on E.
Step III. Now we prove that the third condition (a 3 ) of Theorem 2.5 holds. Let (ϑ, ν) and Then, we have s ds s ds s ds Thus, by taking the supremum of the left-hand side, we get Similarly, proceeding as before, we obtain Adding inequalities (20) and (21), we obtain As (ϑ, ν) = ϑ + ν , we have that (ϑ, ν) ≤ ρ and so condition (a 3 ) of Theorem 2.5 holds.
Step IV. We now prove that M B < 1. Consider From the above estimate, we conclude that condition (a 4 ) of Theorem 2.5 also holds. Thus, the operators A and B satisfy all the conditions of Theorem 2.5 and so, for the considered values of h 1 and h 2 , the coupled system of equations (3) has at least one solution on [1, e].

An example
Hence, all the conditions of Theorem 3.3 are satisfied; consequently, the problem (22) and (23) has a solution on [1, e].

Conclusions
We have presented the existence result for boundary value problems of coupled hybrid differential systems involving the Hadamard fractional derivative. The proof of the existence result is based on the fixed point theorem due to Dhage [10]. The present work can be extended to coupled systems of nonlinear fractional differential equations involving Caputo-Hadamard fractional derivative.