Qualitative analysis of nonlinear coupled pantograph differential equations of fractional order with integral boundary conditions

In this research article, we develop a qualitative analysis to a class of nonlinear coupled system of fractional delay differential equations (FDDEs). Under the integral boundary conditions, existence and uniqueness for the solution of this system are carried out. With the help of Leray–Schauder and Banach fixed point theorem, we establish indispensable results. Also, some results affiliated to Ulam–Hyers (UH) stability for the system under investigation are presented. To validate the results, illustrative examples are given at the end of the manuscript.


Introduction
Theoretical and applied aspects of fractional calculus can be found comprehensively in the literature. This area of mathematics has numerous applications of modeling various phenomena and processes in nonlinear oscillations of earthquakes, nanotechnology, engineering discipline, and other different scientific fields (see [1][2][3][4][5]). These fields played a vital role for the birth of some new subareas of research like fractional dynamics, fractional fuzzy calculus, etc. Different aspects related to qualitative and numerical analysis corresponding to boundary and initial value problems of fractional order differential equations (FODEs) have been widely explored. The literature related to this area contains significant contributions (we refer to [6][7][8][9]).
There is a rapidly growing body of literature on a delay type of differential equations including discrete, continuous, and proportional delay. The area of delay FODEs has many applications in mathematical modeling of real world problems and processes. Among the delay problems, there are those containing a proportional delay, are called pantograph differential equations(PDEs). The existing work about PDEs indicates that such problems have many real world applications. Initially, pantograph was introduced to govern the drawing and scaling aspects of a device. Later on this device was more refined and now scientists use it in electric trains [10,11], material modeling [12], and in modeling lasers, especially quantum dot lasers [13]. Basically, PDEs form a subclass of differential equations called delay, which gives change in terms of dependent variable at previous time [14]. Not long ago, in this area, some beneficial work has been carried out (see [15][16][17]).
FODEs have been extensively studied with respect to different angles. Among these, stability analysis in the UH sense is an important aspect that gained proper attention from researchers [18,19]. Based on the fundamental definition of UH stability of a system, the notion was later modified to more general types, and their results were successfully applied to various problems (we refer to [20,21] and the references therein). The UH-type stability has been developed very well in the sense of differential equations in the last twenty years [22][23][24].
To the best of our knowledge, qualitative study to a coupled system of FDDEs under integral boundary conditions has not been investigated properly. Therefore, using results from fixed point theory, we study the qualitative aspects of a coupled system of FDDEs under integral conditions given as follows: where t ∈ I 1 = [0, 1], β, γ ∈ (1, 2], and λ, p, q ∈ (0, 1) = I 2 . Also the functions f 1 , f 2 : I 1 × R 3 → R are nonlinear continuous functions, φ : I 2 → [0, ∞) is a bounded function, and D is a Riemann-Liouville derivative. Further, the Banach space is defined by (W, , ) with the norm w = max t∈I 1 |w|. Consequently, U = W × W is a Banach space with norms denoted by (w 1 , w 2 ) = w 1 + w 2 or (w 1 , w 2 ) = max{ w 1 , w 2 }.

Preliminaries
Related results and definitions for the current work, which may be found in [25][26][27][28], are discussed here.

Definition 1
Integral with fractional order γ ∈ R + for a function x ∈ L[0, 1] is defined as follows: provided that the integral exists on the right-hand side on (0, ∞).

Definition 2
Riemann-Liouville-type fractional order γ ∈ R + derivative for a function x ∈ L[0, 1] is recalled as follows: where m > γ and the integral on R + is pointwise defined.

Lemma 2
For the FODE, the following result holds: with d k ∈ R and m = [γ ] + 1.

Theorem 1 ([29])
For the two operators S 1 , and if the following matrix tends to zero, subsequently fixed points of (2) are UH-type stable.
Some needed assumptions for our work are given as follows: (M 4 ) For positive real number A and ∀t ∈ I 1 ,

Results about the existence of solutions
Theorem 3 Let r(t) ∈ C(I 1 , R) and u(t) ∈ L(I 1 , R), the solution for the problem is as follows: where G 1 (t, z) is given below: Proof Taking integral to both sides of (7) and using Lemma (2) implies With the help of boundary conditions in (7), we obtain Hence, (10) implies In the same way, for s(t) ∈ C(I 1 , R) and u(t) ∈ L(I 1 , R), the solution for the problem can be obtained in the form where G 2 (t, z) is given below: Corollary 1 In light of Theorem 3, the solution of system (1) is given as follows: where G(t, z) = (G 1 (t, z), G 2 (t, z)) is called Green's function of (1).
Theorem 4 Let f 1 , f 2 , f 3 be continuous functions, subsequently (r, s) ∈ U is a solution of (1) if and only if (r, s) is the solution of the integral equations given in (15).
Proof Assuming that (r, s) is the solution of (15), then differentiating both sides of (15), we get (1). Anyway, if (r, s) is a solution of (1), then (r, s) is a solution of (15).
Let F 1 , F 2 : U → U such that

φ(t)s(t) dt
and F(r, s) = F 1 (r,s) F 2 (r,s) . Thus, solutions of (15) are fixed points of F. Proof Let and We define B a closed subset of E : For arbitrary (r, s) ∈ B, we have

r(λz), s(z) , D p s(z) dz
Obviously, the right-hand side in inequality (19) tends to zero on t 1 → t 2 . Also F 1 is bounded and continuous. Thus it is uniformly bounded. Hence In a similar way, one can also show that Consequently, F : E → E is equicontinuous and so by Arzelá-Ascoli theorem, F is completely continuous.
Proof Taking r,r, s,s ∈ W, and for every t ∈ I 1 , let where Similarly, where Hence, from (20) and (21), one has where d = max t∈I 1 {d 1 , d 2 }. Hence, F is a contraction, thus F has a unique fixed point, which implies that the concerned system (1) has a unique solution.

Results regarding stability
where .
So, from (24) and (25), we get F 1 (r, s) -F 1 (r,s) ≤ 1 r -r + 2 s -s , From (26), we can get the matrix M as follows: Because it is given that the matrix tends to zero, the solution of (1) is UH-type stable.

Conclusion
We have performed a qualitative analysis of the solutions of (1). The problem has been studied with integral boundary conditions having proportional delay called pantograph. Through fixed point theory, the problem in hand has been investigated. Some results about the stability of UH have been discussed via nonlinear analysis. Then, some relevant examples were presented, the results were also discussed. The considered problem is random and includes many applied problems of fluid mechanics and dynamics as special cases [31].