Next Article in Journal
Harris Hawk Optimization-Based Deep Neural Networks Architecture for Optimal Bidding in the Electricity Market
Next Article in Special Issue
New Hermite–Hadamard and Ostrowski-Type Inequalities for Newly Introduced Co-Ordinated Convexity with Respect to a Pair of Functions
Previous Article in Journal
Sufficient Conditions for the Existence and Uniqueness of the Solution of the Dynamics of Biogeochemical Cycles in Coastal Systems Problem
Previous Article in Special Issue
A Survey on the Oscillation of Solutions for Fractional Difference Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequential Fractional Hybrid Inclusions: A Theoretical Study via Dhage’s Technique and Special Contractions

1
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz 53751-71379, Iran
2
Department of Mathematics, University of Ioannina, 451 10 Ioannina, Greece
3
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
4
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 404332, Taiwan
5
Intelligent and Nonlinear Dynamic Innovations Research Center, Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(12), 2090; https://doi.org/10.3390/math10122090
Submission received: 11 May 2022 / Revised: 13 June 2022 / Accepted: 14 June 2022 / Published: 16 June 2022

Abstract

:
The most important objective of the present research is to establish some theoretical existence results on a novel combined configuration of a Caputo sequential inclusion problem and the hybrid integro-differential one in which the boundary conditions are also formulated as the hybrid multi-order integro-differential conditions. In this respect, firstly, some inequalities are proven in relation to the corresponding integral equation. Then, we employ some newly defined theoretical techniques with the help of the product operators on a Banach algebra and also with the aid of some special functions including α - ψ -contractions and α -admissible mappings to extract the existence criteria corresponding to the given mixed sequential hybrid BVPs. Some important useful properties such as the approximate endpoint property, ( C α ) -property, and the compactness play a key role in this regard. The final part of the manuscript is devoted to formulating and computing two applicable examples to guarantee the correctness of the obtained results.

1. Introduction

Human beings need to recognize different interesting phenomena more than ever before. An appropriate approach to meet this demand is to employ the methods and techniques that are available in fractional calculus and, specifically, the fractional operators in modeling of events and processes. Many operators of this fractional type have appeared in recent years, and their consistency and flexibility are becoming known to mathematicians everyday. In such a way, it is convenient that we design complicated and general abstract mathematical models of processes in the format of applicable fractional BVPs. Several examples of the usability of such operators can be seen in branches of science including bio-mathematics, medical science, engineering, and so on; see [1,2,3,4]. All of these points imply that a wide range of scientists are attracted to work on various angles of applicability of such fractional BVPs along with some dynamical behaviors of solutions of these fractional systems. In the mentioned context, most mathematicians have turned toward studying advanced fractional modelings and relevant theoretical findings and graphical behaviors of solutions for these kinds of BVPs; see [5,6,7,8,9,10,11].
Miller with the help of Ross [12] introduced sequential structures of derivatives, which are illustrated by a product of the given derivatives. Later, new findings about different forms of these operators resulted in the publishing of several articles on sequential BVPs in fractional settings. Recently, Alsaedi et al. [13] addressed the primary version of the sequential fractional boundary problem of the Caputo type given by
( C D 0 + η * + μ * C D 0 + η * 1 ) ( z ) = g ˘ * ( z , ( z ) ) , ( 0 ) = 0 , ( 0 ) = 0 , ( σ ) = γ R L I 0 + θ * ( δ )
so that η * ( 2 , 3 ] , δ ( 0 , 1 ) , σ ( δ , 1 ) with θ * > 0 , z [ 0 , 1 ] , γ , μ * R + , and g ˘ * C [ 0 , 1 ] × R , R . The authors in relation to the existence results used some usual notions of the nonlinear analysis theory [13]. For more details, check [14].
During the past few years, a new framework of problems called hybrid differential ones supplemented with abstract types of boundary conditions have been in the spotlight of most mathematicians; see, for example, [15,16,17,18,19,20]. In 2010, the study of this category of equations began with a brilliant paper by Lakshmikantham with the aid of Dhage [21]. Both of them regarded a newly designed hybrid differential equation and obtained its extremal possible solutions by proving several inequalities [21]. After this work, Zhao et al. [22] presented a fractional general extension of the aforementioned BVP in [21] and considered an FBVP consisting of fractional hybrid differential equations. Gradually, researchers have been attracted to this category of modern differential equations. For example, Ahmad et al. [23] carried out a pure mathematical analysis to find the necessary conditions for the existence results of the following Caputo nonlocal hybrid inclusion problem illustrated by
C D 0 + η * ϖ ( z ) j = 1 k R L I 0 + θ j * h j * ( z , ϖ ( z ) ) g * ( z , ϖ ( z ) ) O ( z , ϖ ( z ) ) ,
equipped with conditions ϖ ( 0 ) = β ( z ) and ϖ ( 1 ) = α R , where η * ( 1 , 2 ] and R L I 0 + θ * stands for the RL integral of order θ * > 0 with θ * { θ 1 * , θ 2 * , , θ k * } .
In 2020, Baleanu, Etemad, and Rezapour [24] developed a novel hybrid version of the thermostat fractional model, which describes its performance. Indeed, this model checks an amount of ambient heat and controls it according to the temperature assessed by active sensors. We can review the mentioned thermostat hybrid model as
C D 0 + η * μ ( s ) w * ( s , μ ( s ) ) + Φ ˜ ( s , μ ( s ) ) = 0 , ( η * ( 1 , 2 ] , s [ 0 , 1 ] )
subject to the fractional hybrid boundary conditions
C D 0 + 1 μ ( s ) w * ( s , μ ( s ) ) | s = 0 = 0 , λ * C D 0 + η * 1 μ ( s ) w * ( s , μ ( s ) ) | s = 1 + μ ( s ) w * ( s , μ ( s ) ) | s = δ = 0 ,
where λ * is a positive parameter, δ [ 0 , 1 ] , and η * 1 ( 0 , 1 ] . Furthermore, C D 0 + 1 = d d z , C D 0 + γ * illustrates the Caputo derivative of order γ * { η * , η * 1 } and Φ ˜ , w * C ( [ 0 , 1 ] × R , R ) with w * 0 . The analytical techniques utilized by the authors rely on some fixed point concepts on self-maps and multifunctions.
In 2021, Abbas and Ragusa [25] dealt with a novel category of hybrid FDEs via proportional derivatives depending to a certain function as
ϱ D a + δ , ϑ μ ( s ) w * ( s , μ ( s ) ) = Φ ˜ ( s , μ ( s ) ) , ϱ I a + 1 δ , ϑ μ ( s ) w * ( s , μ ( s ) ) | s = a = b R ,
where ϱ , δ ( 0 , 1 ] and ϱ D a + δ , ϑ displays the δ th-proportional derivative depending on the increasing mapping ϑ with ϑ ( s ) > 0 . They discussed the continuity of solutions of this BVP in terms of some inputs.
By considering and mixing the ideas of the above problems, we formulate a new combined structure of a Caputo sequential hybrid integro-differential inclusion as follows
μ 1 * C D 0 + η * + μ 2 * C D 0 + η * 1 ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) O ˜ ( z , ϖ ( z ) )
furnished with three-point hybrid multi-order integro-differential conditions:
ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , C D 0 + 1 ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , R L I 0 + θ 1 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 1 + R L I 0 + θ 2 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = σ = 0 ,
where z T : = [ 0 , 1 ] , η * [ 2 , 3 ) , η * 1 [ 1 , 2 ) , σ ( 0 , 1 ) , μ 1 * , μ 2 * , θ 1 * , θ 2 * > 0 with θ 1 * 1 > 0 and θ 2 * 1 > 0 . Furthermore, for k N , δ 1 * , δ 2 * , , δ k * > 0 . Here, C D 0 + ( · ) and R L I 0 + ( · ) display the Caputo derivative and the RL integral. Further, notice that C D 0 + 1 = d d z . The continuous real-valued function S ˜ * 0 is defined on T × R k + 1 , and the set-valued map O ˜ is assumed to be considered on T × R having some properties that are pointed out in the sequel. In the next stage, we turn to a special nonhybrid form of the above FBVP as follows:
μ 1 * C D 0 + η * + μ 2 * C D 0 + η * 1 ϖ ( z ) O ˜ ( z , ϖ ( z ) ) , ϖ ( 0 ) = 0 , ϖ ( 0 ) = 0 , R L I 0 + θ 1 * ϖ ( 1 ) + R L I 0 + θ 2 * ϖ ( σ ) = 0 ,
where η * [ 2 , 3 ) , η * 1 [ 1 , 2 ) , z T : = [ 0 , 1 ] , 0 < σ < 1 , μ 1 * , μ 2 * R + , and R L I 0 + ( · ) illustrates the RL integral of both orders θ 1 * , θ 2 * > 0 with θ 1 * 1 > 0 and θ 2 * 1 > 0 . As you see, it is notable that we derive the Caputo sequential inclusion problem (3) if we take S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) = 1 in the given problem (1) and (2). The readers should pay attention to this point that the combined construction of a sequential problem and a hybrid one is a novel and unique fractional modeling, and this kind of mixed sequential hybrid integro-differential inclusion problem has not been discussed in any literature so far. We organize the rest of the research as follows. In Section 2, we review some primitive notions briefly. In Section 3, we develop some newly defined theoretical methods to derive the required existence criteria corresponding to the given mixed BVPs by (1)–(3). Some existing important properties on operators and the existing space including the approximate endpoint property, the ( C α ) -property and the compactness play a key role in this regard. The last section of this manuscript is devoted to formulating and computing two simulation examples to validate the correctness of the results.

2. Preliminaries

This part of the paper is devoted to recalling some basic notions and auxiliary concepts briefly. First of all, we assume that η * > 0 . Then, one can recall the definition of the Riemann–Liouville integral for ϖ : [ 0 , + ) R in the format of R L I 0 + η * ϖ ( z ) = 0 z ( z q ) η * 1 Γ ( η * ) ϖ ( q ) d q such that there exists the R.H.S. value finitely [26,27]. Hereafter, we regard η * ( n 1 , n ) . For ϖ AC R ( n ) ( [ 0 , + ) ) , the derivative in the Caputo sense can be defined by the following framework:
C D 0 + η * ϖ ( z ) = 0 z ( z q ) n η * 1 Γ ( n η * ) ϖ ( n ) ( q ) d q
such that the R.H.S integral exists [26,27]. In addition, for an arbitrary sufficiently smooth map ϖ : [ 0 , + ) R , one can define the sequential fractional derivative as follows:
D 0 + η * ϖ ( z ) = ( D 0 + η 1 * D 0 + η 2 * D 0 + η n * ) ϖ ( z ) ,
where η * = ( η 1 * , η 2 * , , η n * ) is a multi-index [12]. It is notable that the sequential derivative operator D 0 + η * can be considered one of the Caputo, RL, Hadamard, Caputo–Hadamard, etc. In the present research, we shall utilize the sequential derivative in the Caputo setting, which is illustrated by the following. For n = 1 + [ η * ] , the Caputo sequential derivative of the given map ϖ : [ 0 , + ) R is arranged as
C D 0 + η * ϖ ( z ) = D 0 + ( n η * ) d d z n ϖ ( z ) ,
where D 0 + ( n η * ) ϖ ( z ) = R L I 0 + ( n η * ) ϖ ( z ) stands for the RL integral of order n η * [26]. In a monograph presented by Miller and Ross [12], it was demonstrated that the general solution of C D 0 + η * ϖ ( z ) = 0 is illustrated by ϖ ( z ) = c 0 + c 1 z + c 2 z 2 + + c n 1 z n 1 , and we obtain
R L I 0 + η * ( C D 0 + η * ϖ ( z ) ) = ϖ ( z ) + k = 0 n 1 c k z k = ϖ ( z ) + c 0 + c 1 z + c 2 z 2 + + c n 1 z n 1 ,
in which c i are real numbers with n = [ η * ] + 1 . Assuming the normed space ( W , · W ) , we introduce the notations P ( W ) , P c l s ( W ) , P b n d ( W ) , P c m p ( W ) , and P c v x ( W ) to illustrate the set of all nonempty, closed, bounded, compact, and convex subsets of W , respectively. A metric PH d W : P ( W ) × P ( W ) R * is named as the Pompeiu–Hausdorff metric if
PH d W ( E 1 , E 2 ) = max { sup e 1 E 1 d W ( e 1 , E 2 ) , sup e 2 E 2 d W ( E 1 , e 2 ) }
so that d W ( E 1 , e 2 ) = inf e 1 E 1 d W ( e 1 , e 2 ) and d W ( e 1 , E 2 ) = inf e 2 E 2 d W ( e 1 , e 2 ) [28]. We represent all selections of O ˜ at point ϖ C R ( T ) by
( SEL ) O ˜ , ϖ : = { κ ^ L R 1 ( T ) : κ ^ ( z ) O ˜ ( z , ϖ ( z ) ) }
for a.e. z T : = [ 0 , 1 ] [28,29]. Further, ( SEL ) O ˜ , ϖ if dim ( W ) < [28].
A novel category of nonnegative nondecreasing mappings as ψ : [ 0 , ) [ 0 , ) was constructed by Samet et al. [30] in which n = 1 ψ n ( z ) < . This category of functions is illustrated by the notation Ψ . The most important property of these functions is that ψ ( z ) < z for each z > 0 [30]. After this work, Mohammadi et al. introduced a generalized framework for set-valued maps as follows [31]. A multifunction O ˜ : W P b n d , c l s ( W ) is named an α - ψ -contraction whenever for each ϖ , ϖ W , we have
α ( ϖ , ϖ ) PH d W ( O ˜ ϖ , O ˜ ϖ ) ψ ( d W ( ϖ , ϖ ) ) .
We also say that the normed space W possesses the ( C α ) -property if for each convergent sequence { ϖ n } W with ϖ n ϖ and α ( ϖ n , ϖ n + 1 ) 1 for any n N , there exists a subsequence { ϖ n j } of { ϖ n } provided that α ( ϖ n j , ϖ ) 1 for each j N . In the same direction, we say that O ˜ is α -admissible if for every ϖ W and ϖ O ˜ ( ϖ ) with α ( ϖ , ϖ ) 1 , we have α ( ϖ , ϖ ) 1 for all ϖ O ˜ ( ϖ ) [31]. Eventually, an element ϖ W is said to be the endpoint of O ˜ : W P ( W ) if the equality O ˜ ( ϖ ) = { ϖ } holds [32]. In addition, O ˜ is a set-valued map having an approximate endpoint property (APPX-endpoint property) if inf ϖ W sup ϱ O ˜ ϖ d W ( ϖ , ϱ ) = 0 [32]. We need the following lemmas to establish theoretical results about the existence criteria of solutions in this research.
Lemma 1
([33]). The Banach space W is assumed to be separable. Suppose that the multi-function O ˜ : [ 0 , 1 ] × W P c m p , c v x ( W ) is L 1 -Carathéodory and the linear mapping P ˜ * : L W 1 ( [ 0 , 1 ] ) C W ( [ 0 , 1 ] ) is continuous. Then, P ˜ * ( SEL ) O ˜ : C W ( [ 0 , 1 ] ) P c m p , c v x ( C W ( [ 0 , 1 ] ) ) defines an operator in C W ( [ 0 , 1 ] ) × C W ( [ 0 , 1 ] ) via action
ϖ ( P ˜ * ( SEL ) O ˜ ) ( ϖ ) = P ˜ * ( ( SEL ) O ˜ , ϖ )
which has a closed graph.
Lemma 2
([34]). Let W be a Banach algebra and some items be valid for A 1 * : W W and A 2 * : W P c m p , c v x ( W ) including:
(1S)
A 1 * is Lipschitz with l * > 0 ;
(2S)
A 2 * has the upper semi-continuity and the compactness properties;
(3S)
2 l * O ^ < 1 , provided that O ^ = A 2 * ( W ) .
In that case, either ( a ) Σ * = { v * W | α 0 v * ( A 1 * v * ) ( A 2 * v * ) , α 0 > 1 } is unbounded or ( b ) a member exists in W with ϖ ( A 1 * ϖ ) ( A 2 * ϖ ) .
Lemma 3
([31]). Let the metric space ( W , d W ) be complete and α be a nonnegative mapping defined on W 2 , ψ Ψ be a map that increases strictly, and O ˜ : W P c l s , b n d ( W ) be α-admissible and an α-ψ-contraction via α ( ϖ , ϖ ) 1 for some ϖ W and ϖ O ˜ ( ϖ ) . In this case, O ˜ includes a fixed point if W contains the ( C α ) -property.
Lemma 4
([32]). Let the metric space ( W , d W ) be complete and ψ : [ 0 , ) [ 0 , ) involve the upper semi-continuity specification via ψ ( z ) < z and lim inf z ( z ψ ( z ) ) > 0 for all z > 0 . Besides, we assume that O ˜ : W P c l s , b n d ( W ) is such that PH d W ( O ˜ ϖ , O ˜ ϖ ) ψ ( d W ( ϖ , ϖ ) ) for each ϖ , ϖ W . Then, O ˜ involves an endpoint uniquely iff O ˜ possesses the APPX-endpoint property.

3. Main Results

In two previous sections, we assembled some auxiliary and useful notions to achieve our main goals. Now, in the following, we establish other required lemmas to derive the main existence items. To do this, we first regard the sup norm given by ϖ W = sup z T | ϖ ( z ) | on the space W = { ϖ ( z ) : ϖ ( z ) C R ( T ) } . In this case, the Banach space ( W , · W ) along with the multiplication action defined as ( ϖ · ϖ ) ( z ) = ϖ ( z ) ϖ ( z ) is a Banach algebra for all ϖ , ϖ W . In addition, we specify the following constant for the sake of simplicity in computation:
Ω ^ * : = μ 2 * Γ ( θ 1 * + 2 ) 1 Γ ( θ 1 * + 1 ) + μ 2 * σ θ 2 * + 1 Γ ( θ 2 * + 2 ) σ θ 2 * Γ ( θ 2 * + 1 ) + 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) e μ 2 * q d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) e μ 2 * q d q 0 .
On this basis, we explore an integral framework for the possible solution of the proposed Caputo sequential hybrid boundary problem of the inclusion version (1) and (2).
Lemma 5.
Let h ˜ * W , η * [ 2 , 3 ) , η * 1 [ 1 , 2 ) , σ ( 0 , 1 ) , μ 1 * , μ 2 * , θ 1 * , θ 2 * > 0 with θ 1 * 1 > 0 and θ 2 * 1 > 0 . Furthermore, for k N , δ 1 * , δ 2 * , , δ k * > 0 . Then, a solution for the sequential hybrid FDE:
μ 1 * C D 0 + η * + μ 2 * C D 0 + η * 1 ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) = h ˜ * ( z )
furnished with hybrid multi-order integro-differential conditions:
ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , C D 0 + 1 ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , R L I 0 + θ 1 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 1 + R L I 0 + θ 2 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = σ = 0 ,
is displayed as ϖ ( z ) iff
ϖ ( z ) = S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q ] ) ,
where Ω ^ * is given in (4).
Proof. 
We suppose that the function ϖ ˜ 0 * satisfies the given Caputo sequential hybrid Equation (5). Then,
μ 1 * C D 0 + η * 1 + μ 2 * C D 0 + 1 ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) = h ˜ * ( z ) .
In the next step, by taking the η * -th order integral of the Riemann–Liouville type on the above equality, the following non-homogeneous integro-differential equation results:
μ 1 * 1 + μ 2 * C D 0 + 1 ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) = R L I 0 + η * h ˜ * ( z ) + c 0 + c 1 z + c 2 z 2 .
We try to find these unknowns c 0 , c 1 , c 2 R provided
μ 1 * [ ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) ] = R L I 0 + η * h ˜ * ( z ) + c 0 + c 1 z + c 2 z 2 μ 1 * μ 2 * 0 z ϖ ˜ 0 * ( q ) S ˜ * q , ϖ ˜ 0 * ( q ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( q ) , , R L I 0 + δ k * ϖ ˜ 0 * ( q ) d q .
The latter equality implies that μ 1 * ϖ ˜ 0 * ( 0 ) S ˜ * 0 , ϖ ˜ 0 * ( 0 ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( 0 ) , , R L I 0 + δ k * ϖ ˜ 0 * ( 0 ) = c 0 . In addition, by taking the Caputo derivative of Equation (8) of the first order with respect to z, we obtain
μ 1 * C D 0 + 1 ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) = R L I 0 + η * 1 h ˜ * ( z ) + c 1 + 2 c 2 z μ 1 * μ 2 * ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) .
By multiplying both sides of the above equality by e μ 2 * z , one can write
μ 1 * C D 0 + 1 ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) e μ 2 * z = e μ 2 * z R L I 0 + η * 1 h ˜ * ( z ) + c 1 e μ 2 * z + 2 c 2 z e μ 2 * z μ 1 * μ 2 * ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) e μ 2 * z .
After performing some direct computations and necessary simplifications, we obtain
ϖ ˜ 0 * ( z ) S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d q + c 0 μ 1 * e μ 2 * z + c 1 μ 1 * μ 2 * ( 1 e μ 2 * z ) + 2 c 2 μ 1 * μ 2 * 2 ( μ 2 * z 1 + e μ 2 * z ) .
Here, in light of the hybrid multi-order integro-differential conditions (6), we obtain c 0 = c 1 = 0 and
c 2 = μ 2 * 2 2 Ω ^ * 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q μ 2 * 2 2 Ω ^ * 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q ,
where we have mentioned before the nonzero constant Ω ^ * in (4). Eventually, we insert all three values c 0 , c 1 , and c 2 into (9) and (9) becomes
ϖ ˜ 0 * ( z ) = S ˜ * z , ϖ ˜ 0 * ( z ) , R L I 0 + δ 1 * ϖ ˜ 0 * ( z ) , , R L I 0 + δ k * ϖ ˜ 0 * ( z ) ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q ] ) .
This indicates that the function ϖ ˜ 0 * can be regarded as a solution for (7), and the argument is finished. The converse is evident. □
Next, we present the following inequalities, which are useful in the sequel.
Lemma 6.
Let h ˜ * : T R be continuous via h ˜ * = sup z T | h ˜ * ( z ) | . Then, the following inequalities are valid:
(E1)
| 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d q | 1 e μ 2 * μ 2 * Γ ( η * ) h ˜ * ;
(E2)
| 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q | μ 2 * + e μ 2 * 1 μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) h ˜ * ;
(E3)
| 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q | σ θ 2 * + η * 1 ( μ 2 * σ + e μ 2 * σ 1 ) μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) h ˜ * .
Proof. 
( E 1 ) In the first stage, a simple computation yields
0 q ( q r ) η * 2 Γ ( η * 1 ) d r = ( q r ) η * 1 Γ ( η * ) | 0 q = q η * 1 Γ ( η * ) .
In the next step, one can write
0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) d r d q = 0 z e μ 2 * ( z q ) q η * 1 Γ ( η * ) d q z η * 1 Γ ( η * ) 0 z e μ 2 * ( z q ) d q 1 Γ ( η * ) 1 μ 2 * e μ 2 * ( z q ) | 0 z = 1 μ 2 * Γ ( η * ) 1 e μ 2 * z .
Thus, we obtain
| 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d q | 1 e μ 2 * μ 2 * Γ ( η * ) h ˜ * .
( E 2 ) By some easy calculations, we have
0 m ( m r ) η * 2 Γ ( η * 1 ) d r = ( m r ) η * 1 Γ ( η * ) | 0 m = m η * 1 Γ ( η * ) .
In addition, one may write
0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) d r d m = 0 q e μ 2 * ( q m ) m η * 1 Γ ( η * ) d m q η * 1 Γ ( η * ) 0 q e μ 2 * ( q m ) d m = q η * 1 Γ ( η * ) 1 μ 2 * e μ 2 * ( q m ) | 0 q = q η * 1 μ 2 * Γ ( η * ) 1 e μ 2 * q .
Then, by mixing the above results, we obtain
| 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q | h ˜ * 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) q η * 1 ( 1 e μ 2 * q ) μ 2 * Γ ( η * ) d q h ˜ * 1 μ 2 * Γ ( θ 1 * ) Γ ( η * ) 0 1 ( 1 e μ 2 * q ) d q = h ˜ * 1 μ 2 * Γ ( θ 1 * ) Γ ( η * ) 1 + 1 μ 2 * e μ 2 * 1 μ 2 * = h ˜ * 1 μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) μ 2 * + e μ 2 * 1 .
( E 3 ) By similar computations and according to (11) and (12), we have
0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) d r d m q η * 1 μ 2 * Γ ( η * ) 1 e μ 2 * q .
Thus,
| 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) h ˜ * ( r ) d r d m d q | h ˜ * 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) q η * 1 ( 1 e μ 2 * q ) μ 2 * Γ ( η * ) d q h ˜ * 1 μ 2 * Γ ( θ 2 * ) Γ ( η * ) 0 σ σ θ 2 * 1 σ η * 1 ( 1 e μ 2 * q ) d q h ˜ * σ θ 2 * + η * 1 μ 2 * Γ ( θ 2 * ) Γ ( η * ) 0 σ ( 1 e μ 2 * q ) d q = h ˜ * σ θ 2 * + η * 1 μ 2 * Γ ( θ 2 * ) Γ ( η * ) σ + 1 μ 2 * e μ 2 * σ 1 μ 2 * = h ˜ * σ θ 2 * + η * 1 μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) μ 2 * σ + e μ 2 * σ 1 .
Definition 1.
A function ϖ AC R ( T ) is a solution of the suggested sequential multi-order hybrid inclusion boundary problem (1) and (2) if some function κ ^ L R 1 ( T ) exists with κ ^ ( z ) O ˜ z , ϖ ( z ) for almost all z T that satisfies the following hybrid multi-order integro-differential boundary conditions:
ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , C D 0 + 1 ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 0 = 0 , R L I 0 + θ 1 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = 1 + R L I 0 + θ 2 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | z = σ = 0 ,
and
ϖ ( z ) = S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ] )
for any z T .
For now, by taking into account Lemmas 5 and 6, we provide some existence theorems for two Caputo sequential multi-order BVPs (1)–(3).
Theorem 1.
Assume that a continuous real mapping S ˜ * 0 is defined on product space T × W k + 1 and O ˜ : T × W P c v x , c m p ( W ) . Along with these, the following statements are valid:
(HPS1)
There exists a bounded positive-valued mapping υ ˜ : T R + so that for each member ϖ 1 , , ϖ k + 1 , ϖ 1 , , ϖ k + 1 W and z T , we have
| S ˜ * z , ϖ 1 ( z ) , , ϖ k + 1 ( z ) S ˜ * z , ϖ 1 ( z ) , , ϖ k + 1 ( z ) | υ ˜ ( z ) i = 1 k + 1 | ϖ i ( z ) ϖ i ( z ) | ;
(HPS2)
O ˜ : T × W P c v x , c m p ( W ) is supposed to be L 1 -Caratheodory;
(HPS3)
There is a positive function x * L R + 1 ( T ) such that for any ϖ W and for almost all z T , we have
O ˜ ( z , ϖ ) = sup { | κ ^ | : κ ^ O ˜ ( z , ϖ ( z ) ) } x * ( z ) ;
(HPS4)
There is a number r * R + provided that
r * > S ˜ * M * x * 1 υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) M * x * ,
where x * = sup z T { | x ( z ) | } , S ˜ * = sup z T | S ˜ * ( z , 0 , , 0 ) | , υ ˜ * = sup z T | υ ˜ ( z ) | and
M * = | 1 e μ 2 * | μ 1 * μ 2 * Γ ( η * ) + | 1 e μ 2 * | + μ 2 * μ 1 * | Ω ^ * | | μ 2 * + e μ 2 * 1 | μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) + σ θ 2 * + η * 1 | μ 2 * σ + e μ 2 * σ 1 | μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) .
Then, the Caputo sequential hybrid multi-order integro-differential inclusion BVP (1) and (2) has a solution if an inequality υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) M * x * < 1 2 holds strictly.
Proof. 
First of all, for every ϖ W , we construct the following collection of selections for O ˜ as follows:
( SEL ) O ˜ , ϖ = κ ^ L 1 ( T ) : κ ^ ( z ) O ˜ z , ϖ ( z )
for almost all z T . In addition, we consider a set-valued map G * : W P ( W ) by G * ( ϖ ) = { g * W : g * ( z ) = ρ * ( z ) for z T } , where
ρ * ( z ) = S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ] )
for some κ ^ ( SEL ) O ˜ , ϖ . It is evident that g 0 * W is considered as a solution for the Caputo sequential hybrid multi-order inclusion BVP (1) and (2) iff g 0 * is a fixed point of G * . To begin the main proof, with due attention to Lemma 5, we formulate two different structural mappings A 1 * : W W by ( A 1 * ϖ ) ( z ) = S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) and A 2 * : W P ( W ) by
( A 2 * ϖ ) ( z ) = { v W : v ( z ) = φ ( z ) for z T } ,
where
φ ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ]
for some κ ^ ( SEL ) O ˜ , ϖ . Then, one can represent the product operator equation:
G * ( ϖ ) = A 1 * ϖ A 2 * ϖ .
The basic aim at present is to verify that A 1 * and A 2 * settle all hypotheses of Lemma 2. We in the first place intend to confirm that A 1 * is Lipschitz on space W . Let ϖ 1 , ϖ 2 W be arbitrary. By proceeding under the hypothesis ( HPS 1 ), it is deduced that
| ( A 1 * ϖ 1 ) ( z ) ( A 1 * ϖ 2 ) ( z ) | = | S ˜ * z , ϖ 1 ( z ) , R L I 0 + δ 1 * ϖ 1 ( z ) , , R L I 0 + δ k * ϖ 1 ( z ) S ˜ * z , ϖ 2 ( z ) , R L I 0 + δ 1 * ϖ 2 ( z ) , , R L I 0 + δ k * ϖ 2 ( z ) | υ ˜ ( z ) | ϖ 1 ( z ) ϖ 2 ( z ) | + 1 Γ ( δ 1 * + 1 ) | ϖ 1 ( z ) ϖ 2 ( z ) | + + 1 Γ ( δ k * + 1 ) | ϖ 1 ( z ) ϖ 2 ( z ) | = υ ˜ ( z ) ( 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) ) | ϖ 1 ( z ) ϖ 2 ( z ) |
for all z T : = [ 0 , 1 ] . Hence, we obtain
A 1 * ϖ 1 A 1 * ϖ 2 W υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) ϖ 1 ϖ 2 W
for all ϖ 1 , ϖ 2 W . It follows that the single-valued mapping A 1 * is Lipschitz via constant υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) . In the subsequent stage, we proceed to determine that A 2 * involves convex values. To emphasize the correctness of such a claim, let ϖ 1 , ϖ 2 A 2 * ϖ be arbitrary. Choose functions κ 1 ^ , κ 2 ^ ( SEL ) O ˜ , ϖ such that
ϖ l ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ l ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ l ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ l ( r ) d r d m d q ]
for all z T (a.e.) and for l = 1 , 2 . Further, let λ ( 0 , 1 ) . In this case, one may write
λ ϖ 1 ( z ) + ( 1 λ ) ϖ 2 ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) λ κ ^ 1 ( r ) + ( 1 λ ) κ ^ 2 ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) λ κ ^ 1 ( r ) + ( 1 λ ) κ ^ 2 ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) λ κ ^ 1 ( r ) + ( 1 λ ) κ ^ 2 ( r ) d r d m d q ]
for all z T (a.e.). By assumption, we know that O ˜ is convex-valued; thus, we immediately realize that ( SEL ) O ˜ , ϖ is a convex set. Therefore, we obtain that λ κ ^ 1 ( z ) + ( 1 λ ) κ ^ 2 ( z ) ( SEL ) O ˜ , ϖ for each z T , and so, A 2 * w belongs to collection P c v x ( W ) for each ϖ W .
By continuing the proof process, we are going to investigate the complete continuity of A 2 * on W . To reach this aim, we have to prove two notions for the set A 2 * ( W ) including equi-continuity and uniform boundedness. We first check that A 2 * relatesevery bounded set to a bounded one in W . For an arbitrary number r ˜ * R + , we build a bounded ball V r ˜ * = { ϖ W : ϖ W r ˜ * } . In that phase, for every ϖ V r ˜ * and v A 2 * ϖ , there is a function κ ^ ( SEL ) O ˜ , ϖ provided that
v ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ]
for each z T . Then, the estimate of v is implemented as follows:
| v ( z ) | 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q ] 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d m d q ] 1 e μ 2 * μ 1 * μ 2 * Γ ( η * ) + | 1 e μ 2 * | + μ 2 * μ 1 * | Ω ^ * | μ 2 * + e μ 2 * 1 μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) + σ θ 2 * + η * 1 ( μ 2 * σ + e μ 2 * σ 1 ) μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) x * = M * x * ,
so that M * is illustrated by (14). Therefore, v M * x * , which confirms the uniform boundedness of A 2 * ( W ) . In the subsequent step, we verify that A 2 * maps every bounded set to an equi-continuous subset. Let ϖ V r ˜ * and v A 2 * ϖ . We select a function κ ^ ( SEL ) O ˜ , ϖ so that
v ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ]
for all z T . Let us suppose that z 1 , z 2 T along with z 1 < z 2 . Then, we can write
| v ( z 2 ) v ( z 1 ) | | 1 μ 1 * 0 z 2 e μ 2 * ( z 2 q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) r ˜ * d r d q 1 μ 1 * 0 z 1 e μ 2 * ( z 1 q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) r ˜ * d r d q | + ( e μ 2 * z 1 e μ 2 * z 2 ) + μ 2 * ( z 2 z 1 ) μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) r ˜ * d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) r ˜ * d r d m d q ] .
In light of the above arguments, we realize that the limit value of the right-hand expressions equals zero without depending on ϖ V r ˜ * whenever z 1 z 2 . Consequently, we are in a position that one can refer to as the Arzela–Ascoli theorem to deduce that the operator A 2 * : C R ( T ) P ( C R ( T ) is completely continuous. In the sequel, we intend to prove that A 2 * is an operator having a closed graph, which implies its upper semi-continuity property. Let ϖ n V r ˜ * and v n A 2 * ϖ n so that ϖ n ϖ * and v n v * . We claim that the inclusion v * A 2 * ϖ * is valid. For each n 1 and v n A 2 * ϖ n , choose κ ^ n ( SEL ) O ˜ , ϖ n so that
v n ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d m d q ]
for any z T . We have to check the existence of a function κ ^ * ( SEL ) O ˜ , ϖ * provided that
v * ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d m d q ]
for each z T . To achieve this aim, we construct the continuous linear map P ˜ * : L R 1 ( T ) W = C R ( T ) as follows:
P ˜ * ( κ ^ ) ( z ) = ϖ ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ]
for any z T . One can directly confirm
v n ( z ) v * ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) κ ^ * ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) κ ^ * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) κ ^ * ( r ) d r d m d q ] 0 .
Hence, by the assumptions of Lemma 2, it is realized that P ˜ * ( SEL ) O ˜ has a closed graph. For the sake of holding the inclusion v n P ˜ * ( ( SEL ) O ˜ , ϖ n ) and also ϖ n ϖ * , there exists a function κ ^ * ( SEL ) O ˜ , ϖ * provided that
v * ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ * ( r ) d r d m d q ]
for all z T . Hence, v * A 2 * ϖ * and so A 2 * is an operator having a closed graph. From this point, we find that A 2 * is upper semi-continuous. Meanwhile, notice that by assumption, we know that A 2 * is compact-valued. In the following, by taking into account the hypothesis ( HPS 3 ) and continuing a similar argument, we have
O ^ = A 2 * ( W ) = sup z T { | A 2 * ϖ | : ϖ W } 1 e μ 2 * μ 1 * μ 2 * Γ ( η * ) + | 1 e μ 2 * | + μ 2 * μ 1 * | Ω ^ * | μ 2 * + e μ 2 * 1 μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) + σ θ 2 * + η * 1 ( μ 2 * σ + e μ 2 * σ 1 ) μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) x * = M * x * .
As we can observe, an inequality O ^ M * x * holds. Therefore, one can write
υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) ) O ^ υ ˜ * ( 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) M * x * < 1 2 .
By setting l * = υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) , clearly l * O ^ < 1 2 . Until now, it is seen that Lemma 2 is settled on A 1 * and A 2 * . In this position, it is sufficient to check that one of (a) or (b) holds. Our claim is that (b) is not valid. To confirm this claim, with due attention to Lemma 2 and hypothesis ( HPS 4 ), one can suppose that ϖ is an arbitrary member belonging to Σ * with ϖ = r * . Obviously, α 0 ϖ ( z ) A 1 * ϖ ( z ) A 2 * ϖ ( z ) for any α 0 > 1 . By selecting a suitable function κ ^ ( SEL ) O ˜ , ϖ , one can write
ϖ ( z ) = 1 α 0 S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ] )
for each α 0 > 1 and for any z T . Therefore, we obtain:
| ϖ ( z ) | = 1 α 0 | S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) | ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q ] ) = | S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , , R L I 0 + δ k * ϖ ( z ) S ˜ * z , 0 , 0 , , 0 | + | S ˜ * z , 0 , 0 , , 0 | × ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ ( r ) | d r d m d q ] ) [ υ ˜ * ( 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) ) ϖ + S ˜ * ] ( 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) x * ( r ) d r d m d q ] ) υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) ϖ + S ˜ * M * x *
for all z T . After some simplifications, we reach the following inequality
r * S ˜ * M * x * 1 υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + + 1 Γ ( δ k * + 1 ) M * x * .
In view of Inequality (13), the impossibility of Condition (b) of Lemma 2 is deduced. Hence, we have ϖ A 1 * ϖ A 2 * ϖ . Eventually, we manage to verify that there exists a fixed point for G * and the sequential hybrid multi-order integro-differential inclusion BVP (1) and (2) involves a solution. □
In what follows, we carry out our procedure to derive other existence criteria for the Caputo sequential nonhybrid multi-order inclusion BVP (3) by applying two novel pure analytical theorems.
Definition 2.
The function ϖ AC R ( T ) is defined to be a solution of the given Caputo sequential nonhybrid multi-order inclusion BVP (3) if κ ^ L R 1 ( T ) exists via κ ^ O ˜ z , ϖ ( z ) for all z T (a.e.), which satisfies the multi-order integro-derivative boundary conditions:
ϖ ( 0 ) = 0 , ϖ ( 0 ) = 0 , R L I 0 + θ 1 * ϖ ( 1 ) + R L I 0 + θ 2 * ϖ ( σ ) = 0 ,
and
ϖ ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ]
for any z T .
In addition, for each ϖ W , we build the following collection of selections for O ˜ :
( SEL ) O ˜ , ϖ = κ ^ L 1 ( T ) : κ ^ ( z ) O ˜ z , ϖ ( z )
for almost all z T . From here onwards, we consider P : W P ( W ) by
P ( ϖ ) = Φ W : there is κ ^ ( SEL ) O ˜ , ϖ so that Φ ( z ) = b ( z ) for any z T ,
where
b ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ] .
Theorem 2.
Let O ˜ : T × W P c m p ( W ) be compact. Suppose that all six hypotheses hold:
(HPS5)
An integrable operator O ˜ is bounded, and also, O ˜ ( · , ϖ ) : T P c m p ( W ) is a measurable set for every ϖ W ;
(HPS6)
There are ψ Ψ and β C R 0 ( T ) provided that for any z T and ϖ , ϖ W , we have
PH d W O ˜ ( z , ϖ ) , O ˜ ( z , ϖ ) ψ | ϖ ϖ | β ( z ) M * β ,
where sup z T | β ( z ) | = β and M * is displayed in (14);
(HPS7)
ε ˜ : W × W R exists so that ε ˜ ϖ , ϖ 0 for each ϖ , ϖ W ;
(HPS8)
Let { ϖ n } n 1 W go to ϖ and ε ˜ ϖ n ( z ) , ϖ n + 1 ( z ) 0 for any z T . Then, a subsequence { ϖ n l } l 1 of { ϖ n } exists such that ε ˜ ϖ n l ( z ) , ϖ ( z ) 0 for any z T and l 1 ;
(HPS9)
There exist two elements ϖ 0 W and Φ P ( ϖ 0 ) such that ε ˜ ϖ 0 ( z ) , Φ ( z ) 0 for all z T , in which P : W P ( W ) is the same operator illustrated by (15);
(HPS10)
For every ϖ W and Φ P ( ϖ ) along with ε ˜ ϖ ( z ) , Φ ( z ) 0 , there exists b P ( ϖ ) provided that ε ˜ Φ ( z ) , b ( z ) 0 for each z T .
Then, the Caputo sequential nonhybrid multi-order inclusion BVP (3) has a solution.
Proof. 
For the same reason as before, it is explicit that a solution of the given Caputo sequential nonhybrid multi-order FBVP (3) is a fixed point of P : W P ( W ) defined by (15). By ( HPS 5 ) , the measurability of z O ˜ z , ϖ ( z ) is an obvious fact, and so, it is closed-valued for any ϖ W . Therefore, we realize that O ˜ involves a measurable selection and ( SEL ) O ˜ , ϖ is nonempty. Subsequently, we attempt to verify that P ( ϖ ) is closed in W for each ϖ W . To implement the process, we regard a sequence { ϖ n } n 1 of P ( ϖ ) having the property ϖ n ϖ . For each n, we choose κ ^ n ( SEL ) O ˜ , ϖ so that
ϖ n ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ n ( r ) d r d m d q ]
for all z T (a.e.). In light of the compactness of O ˜ , we can pass it into a subsequence to reach { κ ^ n } n 1 , which approaches some κ ^ L 1 ( T ) . Hence, there exists κ ^ ( SEL ) O ˜ , ϖ , and so,
lim n ϖ n ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ ( r ) d r d m d q ] = ϖ ( z )
for any z T . Based on this argument, it follows that ϖ P ( ϖ ) , and so, P is closed-valued. According to the hypothesis of the theorem, O ˜ is compact, so one can easily verify that P ( ϖ ) is bounded for each ϖ W . Here, we want to show that P is α - ψ -contractive. To prove this claim, we introduce α : W 2 [ 0 , ) by α ( ϖ , ϖ ) = 1 if ε ˜ ϖ ( z ) , ϖ ( z ) 0 , and α ( ϖ , ϖ ) = 0 otherwise. Assume that ϖ , ϖ W and Φ 1 P ( ϖ ) are arbitrary. Furthermore, we select κ ^ 1 ( SEL ) O ˜ , ϖ so that
Φ 1 ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d m d q ]
for any z T . Based on the hypothesis (16), we obtain
PH d W O ˜ z , ϖ ( z ) , O ˜ z , ϖ ( z ) β ( z ) ψ | ϖ ( z ) ϖ ( z ) | 1 M * β ,
for each ϖ , ϖ W having the property ε ˜ ϖ ( z ) , ϖ ( z ) 0 for any z T . Hence, b O ˜ z , ϖ ( z ) exists with | κ ^ 1 ( z ) b | ψ | ϖ ( z ) ϖ ( z ) | β ( z ) M * β . Next, we introduce Q * : T P ( W ) given by
Q * ( z ) = b W : | κ ^ 1 ( z ) b | ψ | ϖ ( z ) ϖ ( z ) | β ( z ) M * β
for each z T . Because of the measurability of both κ ^ 1 and ϱ = β ψ | ϖ ϖ | 1 M * β , it follows that the intersection Q * ( · ) O ˜ · , ϖ ( · ) is measurable. In this direction, we select κ ^ 2 , which belongs to O ˜ z , ϖ ( z ) so that
| κ ^ 1 ( z ) κ ^ 2 ( z ) | ψ | ϖ ( z ) ϖ ( z ) | β ( z ) M * β
for all z T . Now, we regard the member Φ 2 P ( ϖ ) by
Φ 2 ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d m d q ]
for any z T . Then, by computing the following estimates, we have
| Φ 1 ( z ) Φ 2 ( z ) | 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) | κ ^ 1 ( r ) κ ^ 2 ( r ) | d r d q + | 1 e μ 2 * z | + μ 2 * z μ 1 * | Ω ^ * | [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ 1 ( r ) κ ^ 2 ( r ) | d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) | κ ^ 1 ( r ) κ ^ 2 ( r ) | d r d m d q ] 1 e μ 2 * μ 1 * μ 2 * Γ ( η * ) + | 1 e μ 2 * | + μ 2 * μ 1 * | Ω ^ * | μ 2 * + e μ 2 * 1 μ 2 * 2 Γ ( θ 1 * ) Γ ( η * ) + σ θ 2 * + η * 1 ( μ 2 * σ + e μ 2 * σ 1 ) μ 2 * 2 Γ ( θ 2 * ) Γ ( η * ) × β ψ ϖ ϖ 1 M * β = M * β ψ ϖ ϖ 1 M * β = ψ ( ϖ ϖ )
for all z T . Therefore, we obtain
Φ 1 Φ 2 = sup z T | Φ 1 ( z ) Φ 2 ( z ) | ψ ( ϖ ϖ )
and so, we obtain
α ( ϖ , ϖ ) PH d W ( P ( ϖ ) , P ( ϖ ) ) ψ ( ϖ ϖ )
for every ϖ , ϖ W . This indicates that P is α - ψ -contractive. Now, by following the proof, we consider ϖ W and ϖ P ( ϖ ) coupled with α ( ϖ , ϖ ) 1 . Then, by taking into account the definition of ε ˜ , we obtain ε ˜ ϖ ( z ) , ϖ ( z ) 0 , and so, there exists b P ( ϖ ) so that ε ˜ ϖ ( z ) , b ( z ) 0 . Hence, α ( ϖ , b ) 1 , and the latter inequality demonstrates that P is α -admissible. To complete the proof process, we consider ϖ 0 W and ϖ P ( ϖ 0 ) so that ε ˜ ϖ 0 ( z ) , ϖ ( z ) 0 for all z. Thus, we obtain α ( ϖ 0 , ϖ ) 1 .
Moreover, consider { ϖ n } n 1 in W via ϖ n ϖ and α ( ϖ n , ϖ n + 1 ) 1 for all n. In this phase, we reach ε ˜ ϖ n ( z ) , ϖ n + 1 ( z ) 0 . Now, by using the assumption ( HPS 8 ) , we figure out that there exists a subsequence { ϖ n l } l 1 of { ϖ n } so that ε ˜ ϖ n l ( z ) , ϖ ( z ) 0 for each z T . As a consequence, α ( ϖ n l , ϖ ) 1 for l 1 , and it is deduced that the space W has the ( C α ) -property. In the final step, by considering Lemma 3, we can easily realize that P includes a fixed point and the sequential nonhybrid multi-order inclusion BVP (3) has at least one solution. □
As a final criterion for the existence of the solution, we review the following theorem based on a new interesting condition attributed to Amini Harandi [32]. On this basis, we shall employ the APPX-endpoint property for P illustrated in (15).
Theorem 3.
Let O ˜ : T × W P c m p ( W ) be compact. Further:
(HPS11)
ψ : [ 0 , ) [ 0 , ) is a nondecreasing map having the upper semi-continuity property via lim inf z ( z ψ ( z ) ) > 0 and ψ ( z ) < z for all z > 0 ;
(HPS12)
O ˜ : T × W P c m p ( W ) is bounded integrable so that O ˜ ( · , ϖ ) : T P c p ( W ) is measurable for every ϖ W ;
(HPS13)
Y C R 0 ( T ) exists such that for all z T and ϖ , ϖ W , we have
PH d W O ˜ ( z , ϖ ) O ˜ ( z , ϖ ) Y ( z ) ψ | ϖ ϖ | 1 M * Y ,
where sup z T | Y ( z ) | = Y and M * is defined by (14);
(HPS14)
P is an operator having the APPX-endpoint property, where P is formulated by (15).
Then, the Caputo sequential nonhybrid multi-order inclusion BVP (3) has at least one solution.
Proof. 
To begin the proof, we want to check the existence of at least one endpoint for given multifunction P : W P ( W ) . To reach this objective, we have to verify that P ( ϖ ) is a closed set for each ϖ W . By ( HPS 12 ) , due to the measurability of z O ˜ z , ϖ ( z ) and also since it is closed-valued for all ϖ W , it follows that O ˜ has a measurable selection and, thus, ( SEL ) O ˜ , ϖ is nonempty for every ϖ W . In this case, similar to the argument of Theorem 2, one can simply see that the subset P ( ϖ ) of W is closed, and so, we omit it. From another angle, we know that P ( ϖ ) is a set having the boundedness property for each ϖ W due to the compactness of O ˜ . Eventually, to end our proof, we must control whether the inequality PH d W ( P ( ϖ ) , P ( ϖ ) ) ψ ( ϖ ϖ ) holds or not. To check this, let ϖ , ϖ W and Φ 1 P ( ϖ ) . Let us select κ ^ 1 ( SEL ) O ˜ , ϖ such that
Φ 1 ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 1 ( r ) d r d m d q ]
for all z T (a.e.). In light of the inequality (17) demonstrated in the assumption ( HPS 13 ) , we know that
PH d W O ˜ ( z , ϖ ) O ˜ ( z , ϖ ) Y ( z ) ψ | ϖ ϖ | 1 M * Y
for any z T ; thus; there exists b * O ˜ z , ϖ ( z ) for which one can write
| κ ^ 1 ( z ) b * | Y ( z ) ψ | ϖ ( z ) ϖ ( z ) | 1 M * Y
for each z T . Now, we introduce R * : T P ( W ) , which is illustrated by
R * ( z ) = b * W : | κ ^ 1 ( z ) b * | Y ( z ) ψ | ϖ ( z ) ϖ ( z ) | 1 M * Y .
For the sake of the measurability of κ ^ 1 and ϱ = Y ψ | ϖ ϖ | 1 M * Y , we can easily realize that R * ( · ) O ˜ ( · , ϖ ( · ) ) is measurable. At this moment, choose the member κ ^ 2 ( z ) O ˜ z , ϖ ( z ) such that
| κ ^ 1 ( z ) κ ^ 2 ( z ) | Y ( z ) ψ | ϖ ( z ) ϖ ( z ) | 1 M * Y
for each z T . Now, we can choose Φ 2 P ( ϖ ) , provided that
Φ 2 ( z ) = 1 μ 1 * 0 z e μ 2 * ( z q ) 0 q ( q r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d q + 1 e μ 2 * z μ 2 * z μ 1 * Ω ^ * [ 0 1 ( 1 q ) θ 1 * 1 Γ ( θ 1 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d m d q + 0 σ ( σ q ) θ 2 * 1 Γ ( θ 2 * ) 0 q e μ 2 * ( q m ) 0 m ( m r ) η * 2 Γ ( η * 1 ) κ ^ 2 ( r ) d r d m d q ]
for all z T . Hence, if we repeat the same process implemented in Theorem 2, then we arrive at inequality
Φ 1 Φ 2 = sup z T | Φ 1 ( z ) Φ 2 ( z ) | M * Y ψ ( ϖ ϖ ) 1 M * Y = ψ ( ϖ ϖ ) .
These findings verify that PH d W ( P ( ϖ ) , P ( ϖ ) ) ψ ( ϖ ϖ ) holds for any ϖ , ϖ W . In addition, from the hypothesis ( HPS 14 ) , we are sure that P is an operator having the APPX-endpoint property, so by referring to Lemma 4, we derive the conclusion, which shows the existence of ϖ * W uniquely such that P ( ϖ * ) = { ϖ * } . Therefore, it is concluded that ϖ * is a solution of the sequential multi-order inclusion FBVP (3). □

4. Examples

This part of the current paper is devoted to supporting the obtained theoretical results by proposing two simulative examples to demonstrate the correctness and the applicability of these outcomes.
Example 1.
In view of the given boundary problem (1) and (2), we here formulate a Caputo sequential hybrid integro-differential inclusion as the following form
3 1000 C D 0 + 2.34 + 5 1000 C D 0 + 1.34 ϖ ( z ) z sin ϖ ( z ) + arcsin ( R L I 0 + 0.71 ϖ ( z ) ) + arctan ( R L I 0 + 0.21 ϖ ( z ) ) 100000 + 37 10000 0 , 0.0002 z 4999 + 1 5000 cos ϖ ( z )
subject to hybrid multi-order integro-derivative boundary conditions
ϖ ( z ) z sin ϖ ( z ) + arcsin ( R L I 0 + 0.71 ϖ ( z ) ) + arctan ( R L I 0 + 0.21 ϖ ( z ) ) 100000 + 37 10000 | z = 0 = 0 , C D 0 + 1 ϖ ( z ) z sin ϖ ( z ) + arcsin ( R L I 0 + 0.71 ϖ ( z ) ) + arctan ( R L I 0 + 0.21 ϖ ( z ) ) 100000 + 37 10000 | z = 0 = 0 , R L I 0 + 1.02 ϖ ( z ) z sin ϖ ( z ) + arcsin ( R L I 0 + 0.71 ϖ ( z ) ) + arctan ( R L I 0 + 0.21 ϖ ( z ) ) 100000 + 37 10000 | z = 1 + R L I 0 + 1.08 ϖ ( z ) z sin ϖ ( z ) + arcsin ( R L I 0 + 0.71 ϖ ( z ) ) + arctan ( R L I 0 + 0.21 ϖ ( z ) ) 100000 + 37 10000 | z = 0.65 = 0 ,
where z T : = [ 0 , 1 ] , η * = 2.34 , η * 1 = 1.34 , μ 1 * = 0.003 , μ 2 * = 0.005 , σ = 0.65 , θ 1 * = 1.02 . θ 2 * = 1.08 , and for k = 2 , we have δ 1 * = 0.71 and δ 2 * = 0.21 . In this case, we obtain Ω ^ * 0.000006 . Here, one can specify the map S ˜ * : T × R 3 R { 0 } defined continuously as follows:
S ˜ * ( z , ϖ 1 ( z ) , ϖ 2 ( z ) , ϖ 3 ( z ) ) = z sin ϖ 1 ( z ) + arcsin ( R L I 0 + 0.71 ϖ 2 ( z ) ) + arctan ( R L I 0 + 0.21 ϖ 3 ( z ) ) 100000 + 37 10000
in which S ˜ * = sup z T | S ˜ * ( z , 0 , 0 , 0 ) | = 0.0037 . We claim that S ˜ * is Lipschitzian. To confirm this claim, for each ϖ , ϖ R , we may write
| S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , R L I 0 + δ 2 * ϖ ( z ) S ˜ * z , ϖ ( z ) , R L I 0 + δ 1 * ϖ ( z ) , R L I 0 + δ 2 * ϖ ( z ) | υ ˜ ( z ) 1 + z δ 1 * Γ ( δ 1 * + 1 ) + z δ 2 * Γ ( δ 2 * + 1 ) | ϖ ( z ) ϖ ( z ) | = z 100000 1 + z 0.71 Γ ( 1.71 ) + z 0.21 Γ ( 1.21 ) | ϖ ( z ) ϖ ( z ) | ,
where υ ˜ ( z ) = z 100000 and υ ˜ * = sup z T | υ ˜ ( z ) | = 0.00001 . Hence, one can observe that the Lipschitz constant of S ˜ * equals
υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + 1 Γ ( δ 2 * + 1 ) = 1 100000 1 + 1 Γ ( 1.71 ) + 1 Γ ( 1.21 ) 0.0000319 > 0 .
In the next stage, we formulate the multifunction O ˜ : T × R P ( R ) in the following framework:
O ˜ z , ϖ ( z ) = 0 , ( 0.0002 z 4999 + 1 5000 ) cos ϖ ( z ) .
As for every v O ˜ z , ϖ ( z ) , we have
| v | max 0 , ( 0.0002 z 4999 + 1 5000 ) cos ϖ ( z ) 0.0002 z 4999 + 0.0002 ,
so we obtain
O ˜ z , ϖ ( z ) = sup { | κ ^ | : κ ^ O ˜ z , ϖ ( z ) } 0.0002 z 4999 + 0.0002 .
By setting x * ( z ) = 0.0002 z 4999 + 0.0002 for any z T , we obtain x * = 0.0004 . Furthermore, by the above obtained values, M * 721419.86984 . In view of the above results and by some direct calculations, we select the positive number r * so that r * > 1.0776 . Then, the following inequality holds:
υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + 1 Γ ( δ 2 * + 1 ) M * x * ( 0.0000319 ) ( 721419.86984 ) ( 0.0004 ) 0.0092053 < 1 2 .
The latter inequality and all the above numerical findings show that Theorem 1 is settled on the present BVP, and so, it is verified that the Caputo sequential hybrid multi-order integro-differential inclusion (18) and (19) has a solution.
Remark 1.
Note that in the above example, we can review other numerical results to check the correctness of the existence criteria of solutions for different values of order η * . Indeed, for every η * [ 2 , 3 ) , we compute the corresponding values of the constant M * and the minimum value of r * . Then, by assuming Δ : = υ ˜ * 1 + 1 Γ ( δ 1 * + 1 ) + 1 Γ ( δ 2 * + 1 ) M * x * , we shall see that an inequality Δ < 0.5 holds in all cases. These numerical values can be observed in Table 1 and Figure 1.
In the last part, the analytical results obtained in Theorem 3 are investigated by the following numerical example.
Example 2.
In this example, we utilize the same given values η * = 2.34 , η * 1 = 1.34 , μ 1 * = 0.003 , μ 2 * = 0.005 , σ = 0.65 , θ 1 * = 1.02 . θ 2 * = 1.08 . According to the given problem (3), we formulate the following Caputo sequential nonhybrid multi-order differential inclusion:
3 1000 C D 0 + 2.34 + 5 1000 C D 0 + 1.34 ϖ ( z ) 0 , ( z 4999 + 1 5000 ) cos ϖ ( z )
furnished with multi-order integro-derivative boundary conditions
ϖ ( 0 ) = 0 , C D 0 + 1 ϖ ( 0 ) = 0 , R L I 0 + 1.02 ϖ ( 1 ) + R L I 0 + 1.08 ϖ ( 0.65 ) = 0
for all z T , where C D 0 + γ stands for the derivative of order γ { 2.34 , 1.34 } of the Caputo type and R L I 0 + λ indicates the RL integral of order λ { 1.02 , 1.08 } . In a similar manner, we have Ω ^ * 0.000006 and M * 2354320.0929 . To investigate the proposed problem precisely, we construct the Banach space as follows W = { ϖ ( z ) : ϖ ( z ) C R ( [ 0 , 1 ] ) } with ϖ W = sup z T | ϖ ( z ) | . Now, we can introduce O ˜ : T × W P ( W ) by
O ˜ z , ϖ ( z ) = 0 , 0.008 ( z + 6 ) 1628 | sin ϖ ( z ) | 1 + | sin ϖ ( z ) |
for each z T . In view of the above set-valued map O ˜ , we reach the function Y C R 0 ( T ) , which is demonstrated by Y ( z ) = 0.008 ( z + 6 ) 814 for any z along with Y = 0.056 814 = 0.0000687 . In addition, the nondecreasing mapping ψ : [ 0 , ) [ 0 , ) with the upper semi-continuity specification is regarded as ψ ( z ) = z 2 for z > 0 . Furthermore, it is evident that lim inf z ( z ψ ( z ) ) > 0 and ψ ( z ) < z for each z > 0 . For both arbitrary elements ϖ , ϖ W , we estimate the following inequality:
PH d W ( O ˜ z , ϖ ( z ) , O ˜ z , ϖ ( z ) ) 0.008 ( z + 6 ) 814 1 2 | ϖ ϖ | = 0.008 ( z + 6 ) 814 ψ | ϖ ϖ | Y ( z ) ψ | ϖ ϖ | 1 M * Y ,
where 1 M * Y 0.006182 . Finally, with due attention to the method implemented in Theorem 3, we regard P : W P ( W ) illustrated by
P ( ϖ ) = Φ W : t h e r e i s κ ^ ( SEL ) O ˜ , ϖ s u c h t h a t Φ ( z ) = b ( z ) f o r a n y z T ,
where
b ( z ) = 1 0.003 0 z e 0.005 ( z q ) 0 q ( q r ) 2.34 2 Γ ( 2.34 1 ) κ ^ ( r ) d r d q + 1 e 0.005 z 0.005 z ( 0.003 ) ( 0.000006 ) [ 0 1 ( 1 q ) 1.02 1 Γ ( 1.02 ) 0 q e 0.005 ( q m ) 0 m ( m r ) 2.34 2 Γ ( 2.34 1 ) κ ^ ( r ) d r d m d q + 0 0.65 ( 0.65 q ) 1.08 1 Γ ( 1.08 ) 0 q e 0.005 ( q m ) 0 m ( m r ) 2.34 2 Γ ( 2.34 1 ) κ ^ ( r ) d r d m d q ] .
Therefore, by referring to Theorem 3 based on a numerical method, it is realized that the Caputo sequential nonhybrid inclusion (20) along with multi-order integro-differential conditions (21) involves a solution. This means that the numerical outcomes of the current example are compatible with all theoretical arguments given in Theorem 3.

5. Conclusions

In this manuscript, we introduced and combined a new configuration of a Caputo sequential inclusion BVP with the hybrid integro-differential inclusion problem in which the boundary conditions are also formulated as the hybrid multi-order integro-derivative conditions. To reach the desired goal of this abstract general fractional model, we derived some theoretical existence results with the help of analytical techniques due to Dhage on the product operators for the given mixed sequential hybrid BVP (1) and (2). Next, based on some conditions on a defined space including the APPX-endpoint property, the ( C α ) -property, and the compactness, we investigated the nonhybrid structure of the suggested BVP (3). Lastly, two examples were designed in this regard. For the next works, one can study the qualitative behaviors and the stability of coupled systems of such sequential hybrid inclusion BVPs by using newly defined fractional operators. Specifically, we can conduct research on the fractal–fractional hybrid models of thermostat or pantograph systems by considering different cases for fractional orders and fractal dimensions.

Author Contributions

Conceptualization, S.E. and S.R.; formal analysis, S.E., S.K.N., B.A., S.R. and J.T.; funding acquisition, J.T.; methodology, S.E., S.K.N., B.A., S.R. and J.T.; software, S.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science, Research and Innovation Fund (NSRF) and King Mongkut’s University of Technology North Bangkok with Contract No. KMUTNB-FF-65-36.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article as no datasets were generated nor analyzed during the current study.

Acknowledgments

The first and fourth authors would like to thank Azarbaijan Shahid Madani University.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FBVPFractional Boundary Value Problem

References

  1. Adjabi, Y.; Samei, M.E.; Matar, M.M.; Alzabut, J. Langevin differential equation in frame of ordinary and Hadamard fractional derivatives under three point boundary conditions. AIMS Math. 2021, 6, 2796–2843. [Google Scholar] [CrossRef]
  2. Mohammadi, H.; Kumar, S.; Rezapour, S.; Etemad, S. A theoretical study of the Caputo-Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control. Chaos Solitons Fractals 2021, 144, 110668. [Google Scholar] [CrossRef]
  3. Thabet, S.T.M.; Abdo, M.S.; Shah, K.; Abdeljawad, T. Abdeljawad, Study of transmission dynamics of COVID-19 mathematical model under ABC fractional order derivative. Res. Phys. 2020, 19, 103507. [Google Scholar]
  4. Thabet, S.T.M.; Etemad, S.; Rezapour, S. On a new structure of the pantograph inclusion problem in the Caputo conformable setting. Bound. Value Probl. 2020, 2020, 171. [Google Scholar] [CrossRef]
  5. Abbas, M.I.; Ragusa, M.A. Solvability of Langevin equations with two Hadamard fractional derivatives via Mittag-Leffler functions. Appl. Anal. 2022, 101, 3231–3245. [Google Scholar] [CrossRef]
  6. Abdo, M.S.; Shah, K.; Wahash, H.A.; Panchal, S.K. On a comprehensive model of the novel coronavirus (COVID-19) under Mittag-Leffler derivative. Chaos Solitons Fractal 2020, 135, 109867. [Google Scholar] [CrossRef]
  7. Amara, A.; Etemad, S.; Rezapour, S. Topological degree theory and Caputo–Hadamard fractional boundary value problems. Adv. Differ. Equ. 2020, 2020, 369. [Google Scholar] [CrossRef]
  8. Boucenna, D.; Boulfoul, A.; Chidouh, A.; Makhlouf, A.B.; Tellab, B. Some results for initial value problem of nonlinear fractional equation in Sobolev space. J. Appl. Math. Comput. 2021, 67, 605–621. [Google Scholar] [CrossRef]
  9. Derbazi, C.; Baitiche, Z.; Abdo, M.S.; Abdeljawad, T. Qualitative analysis of fractional relaxation equation and coupled system with Ψ-Caputo fractional derivative in Banach spaces. AIMS Math. 2020, 6, 2486–2509. [Google Scholar] [CrossRef]
  10. Ntouyas, S.K.; Etemad, S. On the existence of solutions for fractional differential inclusions with sum and integral boundary conditions. Appl. Math. Comput. 2015, 266, 235–243. [Google Scholar] [CrossRef]
  11. Kamenskii, M.; Petrosyan, G.; Wen, C.F. An existence result for a periodic boundary value problem of fractional semilinear differential equations in a Banach space. J. Nonlinear Var. Anal. 2021, 5, 155–177. [Google Scholar]
  12. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley-Interscience, John-Wiley and Sons: New York, NY, USA, 1993. [Google Scholar]
  13. Alsaedi, A.; Ntouyas, S.K.; Agarwal, R.P.; Ahmad, B. On Caputo type sequential fractional differential equations with nonlocal integral boundary conditions. Adv. Differ. Equations 2015, 2015, 33. [Google Scholar] [CrossRef] [Green Version]
  14. Klimek, M. Sequential fractional differential equations with Hadamard derivative. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4689–4697. [Google Scholar] [CrossRef]
  15. Baleanu, D.; Etemad, S.; Pourrazi, S.; Rezapour, S. On the new fractional hybrid boundary value problems with three-point integral hybrid conditions. Adv. Differ. Equations 2019, 2019, 473. [Google Scholar] [CrossRef]
  16. Hilal, K.; Kajouni, A. Boundary value problems for hybrid differential equations with fractional order. Adv. Differ. Equations 2015, 2015, 183. [Google Scholar] [CrossRef] [Green Version]
  17. Mohammadi, H.; Rezapour, S.; Etemad, S.; Baleanu, D. Two sequential fractional hybrid differential inclusions. Adv. Differ. Equations 2020, 2020, 385. [Google Scholar] [CrossRef]
  18. Sun, S.; Zhao, Y.; Han, Z.; Li, Y. The existence of solutions for boundary value problem of fractional hybrid differential equations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 4961–4967. [Google Scholar] [CrossRef]
  19. Sutar, S.T.; Kucche, K.D. On nonlinear hybrid fractional differential equations with Atangana-Baleanu-Caputo derivative. Chaos Solitons Fractals 2021, 143, 110557. [Google Scholar] [CrossRef]
  20. Dhage, B.C. Existence and attractivity theorems for nonlinear hybrid fractional integrodifferential equations with anticipation and retardation. J. Nonlinear Funct. Anal. 2020, 2020, 47. [Google Scholar] [CrossRef]
  21. Dhage, B.C.; Lakshmikantham, V. Basic results on hybrid differential equation. Nonlinear Anal. Hybrid Syst. 2010, 4, 414–424. [Google Scholar] [CrossRef]
  22. Zhao, Y.; Sun, S.; Han, Z.; Li, Q. Theory of fractional hybrid differential equations. Comput. Math. Appl. 2011, 62, 1312–1324. [Google Scholar] [CrossRef] [Green Version]
  23. Ahmad, B.; Ntouyas, S.K.; Tariboon, J. On hybrid Caputo fractional integro-differential inclusions with nonlocal conditions. J. Nonlinear Sci. Appl. 2016, 9, 4235–4246. [Google Scholar] [CrossRef] [Green Version]
  24. Baleanu, D.; Etemad, S.; Rezapour, S. A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions. Bound. Value Probl. 2020, 2020, 64. [Google Scholar] [CrossRef] [Green Version]
  25. Abbas, M.I.; Ragusa, M.A. On the hybrid fractional differential equations with fractional proportional derivatives of a function with respect to a certain function. Symmetry 2021, 13, 264. [Google Scholar] [CrossRef]
  26. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  27. Samko, G.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Longhorne, PA, USA, 1993. [Google Scholar]
  28. Deimling, K. Multi-Valued Differential Equations; Walter de Gruyter: Berlin, Germany, 1992. [Google Scholar]
  29. Aubin, J.; Cellna, A. Differential Inclusions: Set-Valued Maps and Viability Theory; Springer: Berlin/Heidelberg, Germany, 1984. [Google Scholar]
  30. Samet, B.; Vetro, C.; Vetro, P. Fixed point theorems for α-ψ-contractive type mappings. Nonlinear Anal. 2012, 75, 2154–2165. [Google Scholar] [CrossRef] [Green Version]
  31. Mohammadi, B.; Rezapour, S.; Shahzad, N. Some results on fixed points of α-ψ-Ciric generalized multifunctions. Fixed Point Theory Appl. 2013, 2013, 24. [Google Scholar] [CrossRef] [Green Version]
  32. Amini-Harandi, A. Endpoints of set-valued contractions in metric spaces. Nonlinear Anal. 2010, 72, 132–134. [Google Scholar] [CrossRef]
  33. Lasota, A.; Opial, Z. An application of the Kakutani-Ky Fan theorem in the theory of ordinary differential equations. Bull. Acad. Pol. Sci. Set. Sci. Math. Astronom. Phy. 1965, 13, 781–786. [Google Scholar]
  34. Dhage, B.C. Existence results for neutral functional differential inclusions in Banach algebras. Nonlinear Anal. 2006, 64, 1290–1306. [Google Scholar] [CrossRef]
Figure 1. Values of Δ and min r * with respect to different values of η * .
Figure 1. Values of Δ and min r * with respect to different values of η * .
Mathematics 10 02090 g001
Table 1. Numerical values for M * , min r * and Δ corresponding to different values of order η * .
Table 1. Numerical values for M * , min r * and Δ corresponding to different values of order η * .
η * M * r * > Δ < 0.5
2889,292.919611.33130.011347
2.05865,898.416691.29580.011049
2.10841,868.561651.25950.010742
2.15817,328.274601.22240.010429
2.20792,395.794441.18470.010111
2.25767,182.510101.14670.0097892
2.30741,792.858501.10830.0094653
2.34721,419.86984 1.07760.0092053
2.35716,324.283321.06990.0091403
2.40690,867.248161.03160.0088155
2.45665,505.298580.993380.0084918
2.50640,315.167050.955470.0081704
2.55615,366.915630.917950.0078521
2.60590,724.111230.880910.0075376
2.65566,444.028510.844440.0072278
2.70542,577.876030.808610.0069233
2.75519,171.041390.773500.0066246
2.80496,263.351460.739150.0063323
2.85473,889.344100.705620.0060468
2.90452,078.548150.672960.0057685
2.95430,855.768660.641190.0054977
2.99414,314.746630.616440.0052867
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Etemad, S.; Ntouyas, S.K.; Ahmad, B.; Rezapour, S.; Tariboon, J. Sequential Fractional Hybrid Inclusions: A Theoretical Study via Dhage’s Technique and Special Contractions. Mathematics 2022, 10, 2090. https://doi.org/10.3390/math10122090

AMA Style

Etemad S, Ntouyas SK, Ahmad B, Rezapour S, Tariboon J. Sequential Fractional Hybrid Inclusions: A Theoretical Study via Dhage’s Technique and Special Contractions. Mathematics. 2022; 10(12):2090. https://doi.org/10.3390/math10122090

Chicago/Turabian Style

Etemad, Sina, Sotiris K. Ntouyas, Bashir Ahmad, Shahram Rezapour, and Jessada Tariboon. 2022. "Sequential Fractional Hybrid Inclusions: A Theoretical Study via Dhage’s Technique and Special Contractions" Mathematics 10, no. 12: 2090. https://doi.org/10.3390/math10122090

APA Style

Etemad, S., Ntouyas, S. K., Ahmad, B., Rezapour, S., & Tariboon, J. (2022). Sequential Fractional Hybrid Inclusions: A Theoretical Study via Dhage’s Technique and Special Contractions. Mathematics, 10(12), 2090. https://doi.org/10.3390/math10122090

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop