Next Article in Journal
A Hierarchy of Interactions between Pathogenic Virus and Vertebrate Host
Next Article in Special Issue
Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses
Previous Article in Journal
Exact Solutions for Vector Phase-Matching Conditions in Nonlinear Uniaxial Crystals
Previous Article in Special Issue
The Fractional Hilbert Transform of Generalized Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions

by
Muath Awadalla
1,*,
Muthaiah Subramanian
2,
Kinda Abuasbeh
1 and
Murugesan Manigandan
3
1
Department of Mathematics and Statistics, College of Science, King Faisal University, Hafuf, Al Ahsa 31982, Saudi Arabia
2
Department of Mathematics, KPR Institute of Engineering and Technology, Coimbatore 641407, Tamil Nadu, India
3
Department of Mathematics, Sri Ramakrishna Mission Vidyalaya College of Arts and Science, Coimbatore 641020, Tamil Nadu, India
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(11), 2273; https://doi.org/10.3390/sym14112273
Submission received: 10 October 2022 / Revised: 21 October 2022 / Accepted: 26 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus II)

Abstract

:
In this study, we examine the existence and Hyers–Ulam stability of a coupled system of generalized Liouville–Caputo fractional order differential equations with integral boundary conditions and a connection to Katugampola integrals. In the first and third theorems, the Leray–Schauder alternative and Krasnoselskii’s fixed point theorem are used to demonstrate the existence of a solution. The Banach fixed point theorem’s concept of contraction mapping is used in the second theorem to emphasise the analysis of uniqueness, and the results for Hyers–Ulam stability are established in the next theorem. We establish the stability of Ulam–Hyers using conventional functional analysis. Finally, examples are used to support the results. When a generalized Liouville–Caputo ( ρ ) parameter is modified, asymmetric results are obtained. This study presents novel results that significantly contribute to the literature on this topic.

1. Introduction

We consider the nonlinear coupled fractional differential equations with generalized Liouville–Caputo derivatives
C ρ D 0 + ξ p ( τ ) = f ( τ , p ( τ ) , q ( τ ) ) , τ G : = [ 0 , T ] , C ρ D 0 + ζ q ( τ ) = g ( τ , p ( τ ) , q ( τ ) ) , τ G : = [ 0 , T ] ,
enhanced with boundary conditions which are defined by:
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) = ϵ ρ 1 ς Γ ( ς ) 0 ϖ θ ρ 1 ( ϖ ρ θ ρ ) 1 ς q ( θ ) d θ , q ( T ) = π ρ I 0 + ϱ p ( σ ) = π ρ 1 ϱ Γ ( ϱ ) 0 σ θ ρ 1 ( σ ρ θ ρ ) 1 ϱ p ( θ ) d θ , 0 < σ < ϖ < T ,
where C ρ D 0 + ξ , C ρ D 0 + ζ are the Liouville–Caputo-type generalized fractional derivative of order 1 < ξ , ζ 2 , C ρ I 0 + ς , C ρ I 0 + ϱ are the generalized fractional integral of order (Katugampola type) ϱ , ς > 0 , ρ > 0 , f , g : G × R × R R are continuous functions, ϵ , π R . The strip conditions states that the value of the unknown function at the right end point τ = T of the given interval is proportional to the values of the unknown function on the strips of varying lengths. When ρ = 1 , the generalized Liouville–Caputo equation is changed to the Caputo sense, which leads to asymmetric results. In a similar way, when ρ = 1 , the Katugampola integrals are changed to Riemann-Liouville integrals, which leads to cases that are not symmetric. To the best of our knowledge, the stability analysis of boundary value problems (BVPs) is still in its early stages. This paper’s primary contribution is to study existence and Ulam-Hyers stability analysis. In addition, we demonstrate the problem (1)–(2) employed by Leray–Schauder, Banach and Krasnoselskii’s fixed point theorems to prove the existence and uniqueness of solutions. The system (1) is the well-known fractional-order coupled logistic system [1]:
D α u ( τ ) = r 1 u ( τ ) r 1 k 1 u ( τ ) ( u ( τ ) + v ( τ ) ) , τ I , D β v ( τ ) = r 2 v ( τ ) r 2 k 2 v ( τ ) ( v ( τ ) + u ( τ ) ) ,
and the Lotka–Volterra prey-predator system [1]:
D α u ( τ ) = u ( τ ) ( a u ( τ ) E γ v ( τ ) ) , τ I , D β v ( τ ) = v ( τ ) ( b + γ E v ( τ ) β E ) .
We now provide some recent results related to our problem (1)–(2). In [2], the authors discussed the existence results for coupled system of fractional differential equations Riemann–Liouville derivatives
D 0 + α 1 ( D 0 + β 1 x ( t ) ) + f ( t , x ( t ) , y ( t ) ) , t [ 0 , 1 ] , D 0 + α 2 ( D 0 + β 2 y ( t ) ) + f ( t , x ( t ) , y ( t ) ) , t [ 0 , 1 ] ,
with the Riemann–Stieltjes integral boundary conditions:
D 0 + β 1 x ( 0 ) = 0 , x ( 0 ) = 0 , D 0 + β 2 y ( 0 ) = 0 , y ( 0 ) = 0 , x ( 1 ) = γ 1 I 0 + δ 1 y ( ξ ) + i = 1 p 0 1 y ( τ ) d H i ( τ ) , y ( 1 ) = γ 2 I 0 + δ 2 x ( η ) + j = 1 q 0 1 x ( τ ) d K i ( τ ) ,
where α 1 is in the interval ( 0 , 1 ) , β 1 is in the interval ( 1 , 2 ) , α 2 is in the interval ( 0 , 1 ] , β 2 is in the interval ( 1 , 2 ] , p , q N , and γ 1 , γ 2 , δ 1 , δ 2 > 0 , 0 < ξ , η < 1 K j ( t ) , j = 1 , , q , H i ( t ) , i = 1 , , p are bounded variation functions. Both function f and g are nonlinear. They used several theorems from fixed point index theory to prove the main results. In [3], the authors investigated existence of solutions for coupled system of fractional differential equations with Hilfer derivatives
( H D 0 + α 1 , β 1 x ) ( t ) + λ 1 ( H D 0 + α 1 1 , β 1 x ) ( t ) = f ( t , x ( t ) , R ( δ q , , δ 1 ) x ( t ) , y ( t ) ) , t [ 0 , T ] , ( H D 0 + α 2 , β 2 y ) ( t ) + λ 2 ( H D 0 + α 2 1 , β 2 y ) ( t ) = f ( t , x ( t ) , y ( t ) , R ( ζ q , , ζ 1 ) y ( t ) ) , t [ 0 , T ] ,
with Riemann–Liouville and Hadamard-type iterated integral boundary conditions:
x ( 0 ) = 0 , y ( 0 ) = 0 , x ( T ) = i = 1 m ϵ i R ( μ ρ , , μ 1 ) y ( η i ) η i ( 0 , T ) , y ( T ) = j = 1 n θ j R ( ν ρ , , ν 1 ) x ( ξ j ) ξ i ( 0 , T ) ,
where H D α l , β l is the Hilfer fractional derivative operator of order α l with parameters β l , l 1 , 2 , 1 < α l < 2 , 0 β l 1 , λ 1 , λ 2 , ϵ i , θ j R { 0 } , i = 1 , 2 , , m , j = 1 , 2 , , n , f , g : [ 0 , T ] × R × R × R × R are nonlinear continuous functions and R ( ϕ τ , , ϕ 1 ) , ϕ r { δ , ζ , μ , ν } , r { q , p , ρ | q , p , ρ N } , involves the iterated Riemann–Liouville and Hadamard fractional integral operators. They used several theorems from fixed point index theory to prove the main results. Numerous scientific and engineering phenomena are mathematically modelled using fractional order differential and integral operators. The main benefit of adopting these operators is their nonlocality, which enables the description of the materials and processes involved in the history of the phenomenon. As a result, compared to their integer-order counterparts, fractional-order models are more precise and informative. As a result of the extensive use of fractional calculus techniques in a range of real-world occurrences, such as those described in the texts cited [4,5,6,7,8] numerous researchers developed this significant branch of mathematical study. In recent years, a lot of research has been done on fractional differential equations with different boundary conditions. Nonlocal nonlinear fractional-order boundary value problems, in particular, have attracted a lot of attention (BVPs). The idea of nonlocal circumstances, which help to describe physical processes occurring inside the confines of a specific domain, was originally introduced in the work of Bitsadze and Samarski [9]. It is challenging to defend the assumption of a circular cross-section in computational fluid dynamics investigations of blood flow problems because to the changing shape of a blood vessel throughout the vessel. To solve that problem, integral boundary conditions have been developed. In addition, the ill-posed parabolic backward problems are solved under integral boundary conditions. Integral boundary conditions are also essential in mathematical models of bacterial self-regularization, as shown in [10]. Fractional order differential equations, as well as inclusions including Riemann–Liouville, Liouville–Caputo (Caputo), and Hadamard-type derivatives, among others, have all been included in the literature on the topic recently. For some recent works on the topic, we point the reader to several papers [11,12,13,14,15] and the references listed therein. The use of fractional differential systems in mathematical representations of physical and engineering processes has drawn considerable interest. See [16,17,18,19,20,21,22] for additional details on the theoretical evolution of such systems. The following is the remainder of the article: Section 2 introduces some fundamental definitions, lemmas, and theorems that support our main results. For the existence and uniqueness of solutions to the given system (1) and (2), we use various conditions and some standard fixed-point theorems in Section 3. Section 4 discusses the Ulam–Hyers stability of the given system (1) and (2) under certain conditions. In Section 6, examples are provided to demonstrate the main results. Finally, the consequences of existence, uniqueness, and stability for the problem (1) and (75) are provided.

2. Preliminaries

For our research, we recall some preliminary definitions of generalized Liouville–Caputo fractional derivatives and Katugampola fractional integrals.
The space of all complex-valued Lebesgue measurable functions ϕ on ( c , d ) equipped with the norm is denoted by Z b q ( c , d ) :
| | ϕ | | Z b q = c d | z b ϕ ( z ) | q d z z 1 q < , b R , 1 q .
Let L 1 ( c , d ) represent the space of all Lebesgue measurable functions φ on ( c , d ) endowed with the norm:
φ L 1 = c d | φ ( z ) | d z < .
We further recall that A C n ( E , R ) = { p : E R : p , p , , p ( n 1 ) C ( E , R ) a n d p ( n 1 ) is absolutely continuous. For 0 ϵ < 1 , we define C ϵ , ρ ( E , R ) = { f : E R : ( τ ρ a ρ ) ϵ f ( τ ) C ( E , R ) endowed with the norm f C ϵ , ρ = ( τ ρ a ρ ) ϵ f ( τ ) C . Moreover, we define the class of functions f that have absolute continuous δ n 1 derivative, denoted by A C γ n ( E , R ) , as follows: A C γ n ( E , R ) = { f : E R : γ n 1 f AC ( E , R ) , γ = τ 1 ρ d d τ } , which is equipped with the norm f C γ , ϵ n = k = 0 n 1 γ k f C + γ n f C ϵ , ρ is defined by
C γ , ϵ n ( E , R ) = f : E R : γ n 1 f C ( E , R ) , γ n f C ϵ , ρ ( E , R ) , γ = τ 1 ρ d d τ .
Notice that C γ , 0 n = C γ n . We define space P = { p ( τ ) : p ( τ ) C ( E , R ) } equipped with the norm | | p | | = sup { | p ( τ ) | , τ E } - this is a Banach space. Furthermore Q = { q ( τ ) : q ( τ ) C ( E , R ) } equipped with the norm is | | q | | = sup { | q ( τ ) | , τ E } is a Banach space. Then the product space ( P × Q , | | ( p , q ) | | ) is also a Banach space with norm | | ( p , q ) | | = | | p | | + | | q | | .
Definition 1
([23]). The left and right-sided generalized fractional integrals (GFIs) of f Z b q ( c , d ) of order ξ > 0 and ρ > 0 for < c < τ < d < , are defined as follows:
( ρ I c + ξ f ) ( τ ) = ρ 1 ξ Γ ( ξ ) c τ θ ρ 1 ( τ ρ θ ρ ) 1 ξ f ( θ ) d θ ,
( ρ I d ζ f ) ( τ ) = ρ 1 ξ Γ ( ξ ) τ d θ ρ 1 ( θ ρ τ ρ ) 1 ξ f ( θ ) d θ .
Definition 2
([24]). The generalized fractional derivatives (GFDs) which are associated with GFIs (7) and (8) for 0 c < τ < d < , are defined as follows:
( ρ D c + ξ f ) ( τ ) = τ 1 ρ d d τ n ( ρ I c + n ξ f ) ( τ ) = ρ ξ n + 1 Γ ( n ξ ) τ 1 ρ d d τ n c τ θ ρ 1 ( τ ρ θ ρ ) ξ n + 1 f ( θ ) d θ ,
( ρ D d ξ f ) ( τ ) = τ 1 ρ d d τ n ( ρ I d n ξ f ) ( τ ) = ρ ξ n + 1 Γ ( n ξ ) τ 1 ρ d d τ n τ d θ ρ 1 ( τ ρ θ ρ ) ξ n + 1 f ( θ ) d θ ,
if the integrals exist.
Definition 3
([25]). The above GFDs define the left and right-sided generalized Liouville–Caputo type fractional derivatives of f A C γ n [ c , d ] of order ξ 0
C ρ D c + ξ f ( z ) = ρ D c + ξ f ( τ ) k = 0 n 1 γ k f ( c ) k ! τ ρ c ρ ρ k ( z ) , γ = z 1 ρ d d z ,
C ρ D d ξ f ( z ) = ρ D d ξ f ( τ ) k = 0 n 1 ( 1 ) k γ k f ( d ) k ! d ρ τ ρ ρ k ( z ) , γ = z 1 ρ d d z ,
when n = [ ξ ] + 1 .
Lemma 1
([25]). Let ξ 0 , n = [ ξ ] + 1 and f A C γ n [ c , d ] , where 0 < c < d < . Then,
1.
if ξ N
C ρ D c + ξ f ( τ ) = 1 Γ ( n ξ ) c τ τ ρ θ ρ ρ n ξ 1 ( γ n f ) ( θ ) d θ θ 1 ρ = ρ I c + n ξ ( γ n f ) ( τ ) ,
C ρ D d ξ f ( τ ) = 1 Γ ( n ξ ) τ d θ ρ τ ρ ρ n ξ 1 ( 1 ) n ( γ n f ) ( θ ) d θ θ 1 ρ = ρ I d n ξ ( γ n f ) ( τ ) .
2.
if ξ N
C ρ D c + ξ f = γ n f , C ρ D d ξ f = ( 1 ) n γ n f .
Lemma 2
([25]). Let f A C γ n [ c , d ] or C γ n [ c , d ] and ξ R . Then,
ρ I c + ξ C ρ D c + ξ f ( z ) = f ( z ) k = 0 n 1 γ k f ( c ) k ! z ρ c ρ ρ k ,
ρ I d ξ C ρ D d ξ f ( z ) = f ( z ) k = 0 n 1 ( 1 ) k γ k f ( d ) k ! d ρ z ρ ρ k .
In particular, for 0 < ξ 1 , we have
ρ I c + ξ C ρ D c + ξ f ( z ) = f ( z ) f ( c ) , ρ I d ξ C ρ D d ξ f ( z ) = f ( z ) f ( d ) .
We introduce the following notations for computational ease:
E 1 = ϵ ϖ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) , E 2 = π σ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) , E ^ = T ρ ρ ,
G = E ^ 2 E 1 E 2 0 ,
δ ( τ ) = τ ρ ρ G .
Next, we are proving a lemma, which is vital in converting the given problem to a fixed-point problem.
Lemma 3.
Given the functions f ^ , g ^ C ( 0 , T ) L ( 0 , T ) , p , q A C γ 2 ( E ) and Λ 0 . Then the solution of the coupled BVP:
C ρ D 0 + ξ p ( τ ) = f ^ ( τ ) , τ E : = [ 0 , T ] , C ρ D 0 + ζ q ( τ ) = g ^ ( τ ) , τ E : = [ 0 , T ] , p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) , q ( T ) = π ρ I 0 + ϱ p ( σ ) 0 < σ < ϖ < T ,
is given by
p ( τ ) = ρ I 0 + ξ f ^ ( τ ) + δ ( τ ) E ^ ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + E 1 π ρ I 0 + ξ + ϱ f ^ ( σ ) ρ I 0 + ζ g ^ ( T )
and
q ( τ ) = ρ I 0 + ζ g ^ ( τ ) + δ ( τ ) E ^ π ρ I 0 + ξ + ϱ f ^ ( σ ) ρ I 0 + ζ g ^ ( T ) + E 2 ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) .
Proof. 
When ρ I 0 + ξ , ρ I 0 + ζ are applied to the FDEs in (19) and Lemma 2 is used, the solution of the FDEs in (19) for τ E is
p ( τ ) = ρ I 0 + ξ f ^ ( τ ) + a 1 + a 2 τ ρ ρ = ρ 1 ξ Γ ( ξ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) ξ 1 f ^ ( θ ) d θ + a 1 + a 2 τ ρ ρ ,
q ( τ ) = ρ I 0 + ζ g ^ ( τ ) + b 1 + b 2 τ ρ ρ = ρ 1 ζ Γ ( ζ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) ζ 1 g ^ ( θ ) d θ + b 1 + b 2 τ ρ ρ ,
respectively, for some a 1 , a 2 , b 1 , b 2 R . Making use of the boundary conditions p ( 0 ) = q ( 0 ) = 0 in (22) and (23) respectively, we get a 1 = b 1 = 0 . Next, we obtain by using the generalized integral operators ρ I 0 + ξ , ρ I 0 + ζ (22) and (23) respectively,
ρ I 0 + ϱ p ( τ ) = ρ I 0 + ξ + ϱ f ^ ( τ ) + a 1 τ ρ ϱ ρ ϱ Γ ( ϱ + 1 ) + a 2 τ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) ,
ρ I 0 + ς q ( τ ) = ρ I 0 + ζ + ς g ^ ( τ ) + b 1 τ ρ ς ρ ς Γ ( ς + 1 ) + b 2 τ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) ,
which, when combined with the boundary conditions p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) , q ( T ) = π ρ I 0 + ϱ p ( σ ) , gives the following results:
ρ I 0 + ξ f ^ ( T ) + a 1 + a 2 T ρ ρ = ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) + b 1 ϵ ϖ ρ ς ρ ς Γ ( ς + 1 ) + b 2 ϵ ϖ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) ,
ρ I 0 + ζ g ^ ( T ) + b 1 + b 2 T ρ ρ = π ρ I 0 + ξ + ϱ f ^ ( σ ) + a 1 π σ ρ ϱ ρ ϱ Γ ( ϱ + 1 ) + a 2 π σ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) .
Next, we obtain
a 2 E ^ b 2 E 1 = ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) ,
b 2 E ^ a 2 E 2 = π ρ I 0 + ξ + ϱ f ^ ( σ ) ρ I 0 + ζ g ^ ( T ) ,
by employing the notations (16) in (26) and (27) respectively. We find that when we solve the system of Equations (28) and (29) for a 2 and b 2 ,
a 2 = 1 G E ^ ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + E 1 π ρ I 0 + ξ + ϱ f ^ ( σ ) ρ I 0 + ζ g ^ ( T ) ,
b 2 = 1 G E 2 ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + E ^ π ρ I 0 + ξ + ϱ f ^ ( σ ) ρ I 0 + ζ g ^ ( T ) .
Substituting the values of a 1 , a 2 , b 1 , b 2 in ( 22 ) and ( 23 ) respectively, we get the solution for the BVP (19). □

3. Existence Results for the Problem (1) and (2)

As a result of Lemma 3, we define an operator Δ : P × Q P × Q by
Δ ( p , q ) ( τ ) = ( Δ 1 ( p , q ) ( τ ) , Δ 2 ( p , q ) ( τ ) ) ,
where
Δ 1 ( p , q ) ( τ ) = ρ I 0 + ξ f ( τ , p ( τ ) , q ( τ ) ) + δ ( τ ) [ E ^ ϵ ρ I 0 + ζ + ς g ( ϖ , p ( ϖ ) , q ( ϖ ) ) ρ I 0 + ξ f ( T , p ( T ) , q ( T ) ) + E 1 π ρ I 0 + ξ + ϱ f ( σ , p ( σ ) , q ( σ ) ) ρ I 0 + ζ g ( T , p ( T ) , q ( T ) ) ] ,
Δ 2 ( p , q ) ( τ ) = ρ I 0 + ζ g ( τ , p ( τ ) , q ( τ ) ) + δ ( τ ) [ E ^ π ρ I 0 + ξ + ϱ f ( σ , p ( σ ) , q ( σ ) ) ρ I 0 + ζ g ( T , p ( T ) , q ( T ) ) + E 2 ϵ ρ I 0 + ζ + ς g ( ϖ , p ( ϖ ) , q ( ϖ ) ) ρ I 0 + ξ f ( T , p ( T ) , q ( T ) ) ] .
For brevity’s sake, we’ll use the following notations:
J 1 = T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ,
K 1 = | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) ,
J 2 = | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ,
K 2 = T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) ,
Φ = min { 1 [ ψ 1 ( J 1 + J 2 ) + ψ 1 ^ ( K 1 + K 2 ) ] , 1 [ ψ 2 ( J 1 + J 2 ) + ψ 2 ^ ( K 1 + K 2 ) ] } .
Theorem 1.
Assume that f , g : E × R × R R are continuous functions satisfying the condition: ( A 1 ) there exists constants ψ m , ψ m ^ 0 ( m = 1 , 2 ) and ψ 0 , ψ 0 ^ > 0 such that
| f ( τ , o 1 , o 2 ) | ψ 0 + ψ 1 | o 1 | + ψ 2 | o 2 | , | g ( τ , o 1 , o 2 ) | ψ 0 ^ + ψ 1 ^ | o 1 | + ψ 2 ^ | o 2 | , o m R , m = 1 , 2 .
If ψ 1 ( J 1 + J 2 ) + ψ 1 ^ ( K 1 + K 2 ) < 1 , ψ 2 ( J 1 + J 2 ) + ψ 2 ^ ( K 1 + K 2 ) < 1 . Thenat least one solution for the BVP (1) and (2) on E , where J 1 , K 1 , J 2 , K 2 are given by (35)–(38) respectively.
Proof. 
We define operator Δ : P × Q P × Q as being completely continuous in the first step. The continuity of the functions f and g implies that the operators Δ 1 and Δ 2 are continuous. As a result, the operator Δ is continuous. Let Ψ P × Q be a bounded set to demonstrate the uniformly bounded operator Δ . Then N 1 ^ and N 2 ^ are positive constants such that | f ( τ , p ( τ ) , q ( τ ) ) | N 1 ^ , | g ( τ , p ( τ ) , q ( τ ) ) | N 2 ^ , ( p , q ) Ψ . Then we have
| Δ 1 ( p , q ) ( τ ) | ρ I 0 + ξ | f ( τ , p ( τ ) , q ( τ ) ) | + | δ ( τ ) | [ | E ^ | | ϵ | ρ I 0 + ζ + ς | g ( ϖ , p ( ϖ ) , q ( ϖ ) ) | + ρ I 0 + ξ | f ( T , p ( T ) , q ( T ) ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | f ( σ , p ( σ ) , q ( σ ) ) | + ρ I 0 + ζ | g ( T , p ( T ) , q ( T ) ) | ] N ^ 1 | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + N 2 ^ | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) | δ | ,
when taking the norm and using (35) and (36), that yields for ( p , q ) Ψ ,
| | Δ 1 ( p , q ) | | J 1 N ^ 1 + K 1 N ^ 2 .
Likewise, we obtain
| | Δ 2 ( p , q ) | | N 2 ^ | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + N 1 ^ | δ | | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) J 2 N ^ 1 + K 2 N ^ 2 ,
using (37) and (38). Based on the inequalities (40) and (41), we can conclude that Δ 1 and Δ 2 are uniformly bounded, which indicates that the operator Δ is uniformly bounded. Next, we show that Δ is equicontinuous. Let τ 1 , τ 2 E with τ 1 < τ 2 . Then we have
| Δ 1 ( p , q ) ( τ 2 ) Δ 1 ( p , q ) ( τ 1 ) | | ρ I 0 + ξ f ( τ 2 , p ( τ 2 ) , q ( τ 2 ) ) ρ I 0 + ξ f ( τ 1 , p ( τ 1 ) , q ( τ 1 ) ) | + | δ ( τ 2 ) δ ( τ 1 ) | [ E ^ | ϵ | ρ I 0 + ζ + ς | g ( ϖ , p ( ϖ ) , q ( ϖ ) ) | + ρ I 0 + ξ | f ( T , p ( T ) , q ( T ) ) | + E 1 | π | ρ I 0 + ξ + ϱ | f ( σ , p ( σ ) , q ( σ ) ) | + ρ I 0 + ζ | g ( T , p ( T ) , q ( T ) ) | ] ρ 1 ξ N 1 ^ Γ ( ξ ) 0 τ 1 θ ρ 1 ( τ 2 ρ θ ρ ) 1 ξ θ ρ 1 ( τ 1 ρ θ ρ ) 1 ξ d θ + τ 1 τ 2 θ ρ 1 ( τ 2 ρ θ ρ ) 1 ξ d θ + | δ ( τ 2 ) δ ( τ 1 ) | | E ^ | N 2 ^ | ϵ | ϖ ρ ζ + ς ρ ζ + ς Γ ( ζ + ς + 1 ) + N 1 ^ T ρ ξ ρ ξ Γ ( ξ + 1 ) + | δ ( τ 2 ) δ ( τ 1 ) | | E 1 | N 1 ^ | π | σ ρ ξ + ϱ ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + N 2 ^ T ρ ζ ρ ζ Γ ( ζ + 1 ) 0 a s τ 2 τ 1 .
independent of ( p , q ) with respect to | f ( τ , p ( τ 1 ) , q ( τ 1 ) ) | N 1 ^ and | g ( τ , p ( τ 1 ) , q ( τ 1 ) ) | N 2 ^ . Similarly, we can express | Δ 2 ( p , q ) ( τ 2 ) Δ 2 ( p , q ) ( τ 1 ) | 0 as τ 2 τ 1 independent of ( p , q ) in terms of the boundedness of f and g. As a result of the equicontinuity of Δ 1 and Δ 2 , operator Δ is equicontinuous. As a result of the Arzela–Ascoli theorem, the operator is compact. Finally, we demonstrate that the set Π ( Δ ) = { ( p , q ) P × Q : λ Δ ( p , q ) ; 0 < λ < 1 } is bounded. Let ( p , q ) Π ( Δ ) .Then ( p , q ) = λ Δ ( p , q ) . For any τ E , we have p ( τ ) = λ Δ 1 ( p , q ) ( τ ) , q ( τ ) = λ Δ 2 ( p , q ) ( τ ) . By utilizing ( A 1 ) in (33), we obtain
| p ( τ ) | ρ I 0 + ξ ( ψ 0 , ψ 1 | p ( τ ) | , ψ 2 | q ( τ ) | ) + | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς ( ψ 0 ^ + ψ 1 ^ | p ( ϖ ) | + ψ 2 ^ | q ( ϖ ) | ) + ρ I 0 + ξ ( ψ 0 + ψ 0 | p ( T ) | + ψ 2 | q ( T ) | ) + | E 1 | | π | ρ I 0 + ξ + ϱ ( ψ 0 + ψ 1 | p ( σ ) | + ψ 2 | q ( σ ) | ) + ρ I 0 + ζ ( ψ 0 ^ + ψ 1 ^ | p ( T ) | + ψ 2 ^ | q ( T ) | ) ) ,
which results when taking the norm for τ E ,
| | p | | ( ψ 0 + ψ 1 | | p | | + ψ 2 | | q | | ) J 1 + ( ψ 0 ^ + ψ 1 ^ | | p | | + ψ 2 ^ | | q | | ) K 1 .
Similarly, we are capable of obtaining that
| | q | | ( ψ 0 ^ + ψ 1 ^ | | p | | + ψ 2 ^ | | q | | ) K 2 + ( ψ 0 + ψ 1 | | p | | + ψ 2 | | q | | ) J 2 .
From (43) and (44), we get
| | p | | + | | q | | = ψ 0 ( J 1 + J 2 ) + ψ 0 ^ ( K 1 + K 2 ) + | | p | | ψ 1 ( J 1 + J 2 ) + ψ 1 ^ ( K 1 + K 2 ) + | | q | | ψ 1 ( J 1 + J 2 ) + ψ 1 ^ ( K 1 + K 2 ) ,
which results, with | | ( p , q ) | | = | | p | | + | | q | | ,
| | ( p , q ) | | ψ 0 ( J 1 + J 2 ) + ψ 0 ^ ( K 1 + K 2 ) Φ .
As a result, Π ( Δ ) is bounded. Thus, the nonlinear alternative of Leray–Schauder [26] is valid and the operator Δ has at least one fixed point. It implies that the BVP (1) and (2) contain at least one solution on E . □
Theorem 2.
Assume that f , g : E × R × R R are continuous functions satisfying the condition: ( A 2 ) there exists constants ϕ m , ϕ m ^ 0 ( m = 1 , 2 ) such that
| f ( τ , o 1 , o 2 ) f ( τ , o ^ 1 , o ^ 2 ) | ϕ 1 | o 1 o ^ 1 | + ϕ 2 | o 2 o ^ 2 | , | g ( τ , o 1 , o 2 ) g ( τ , o ^ 1 , o ^ 2 ) | ϕ 1 ^ | o 1 o ^ 1 | + ϕ 2 ^ | o 2 o ^ 2 | , o m , o ^ m R , m = 1 , 2 .
Furthermore, there exist S 1 , S 2 > 0 such that | f ( τ , 0 , 0 ) | S 1 , | g ( τ , 0 , 0 ) | S 2 , Then, given that
( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) < 1 ,
the BVP (1) and (2) has a unique solution on E , where J 1 , K 1 , J 2 , K 2 are given by (35)–(38) respectively.
Proof. 
Let us fix φ ( J 1 + J 2 ) S 1 + ( K 1 + K 2 ) S 2 1 ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) and demonstrate that Δ B φ B φ when operator Δ is given by (32) and B φ = { ( p , q ) P × Q : | | ( p , q ) | | φ } . For ( p , q ) B φ , τ E
| f ( τ , p ( τ ) , q ( τ ) ) | ϕ 1 | p ( τ ) | + ϕ 2 | q ( τ ) | + S 1 ϕ 1 | | p | | + ϕ 2 | | q | | + S 1 ,
and
| g ( τ , p ( τ ) , q ( τ ) ) | ϕ 1 ^ | p ( τ ) | + ϕ 2 ^ | q ( τ ) | + S 2 ϕ 1 ^ | | p | | + ϕ 2 ^ | | q | | + S 2 .
This guides to
| Δ 1 ( p , q ) ( τ ) | ρ I 0 + ξ | f ( τ , p ( τ ) , q ( τ ) ) f ( τ , 0 , 0 ) | + | f ( τ , 0 , 0 ) | + | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς g [ ( ϖ , p ( ϖ ) , q ( ϖ ) ) g ( ϖ , 0 , 0 ) | + | g ( ϖ , 0 , 0 ) | ] + ρ I 0 + ξ f [ ( T , p ( T ) , q ( T ) ) f ( T , 0 , 0 ) | + | f ( T , 0 , 0 ) ] | + | E 1 | | π | ρ I 0 + ξ + ϱ f [ f ( σ , p ( σ ) , q ( σ ) ) f ( σ , 0 , 0 ) | + | f ( σ , 0 , 0 ) | ] + ρ I 0 + ζ [ | g ( T , p ( T ) , q ( T ) ) g ( T , 0 , 0 ) | + | g ( T , 0 , 0 ) | ] )
( ϕ 1 | | p | | + ϕ 2 | | q | | + S 1 ) T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + ( ϕ 1 ^ | | p | | + ϕ 2 ^ | | q | | + S 2 ) | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 )
| | Δ 1 ( p , q ) | | ( ϕ 1 | | p | | + ϕ 2 | | q | | + S 1 ) J 1 + ( ϕ 1 ^ | | p | | + ϕ 2 ^ | | q | | + S 2 ) K 1 .
Similarly, we obtain
| Δ 2 ( p , q ) ( τ ) | ( ϕ 1 ^ | | p | | + ϕ 2 ^ | | q | | + S 2 ) T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + ( ϕ 1 | | p | | + ϕ 2 | | q | | + S 1 ) | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) | | Δ 2 ( p , q ) | | ( ϕ 1 ^ | | p | | + ϕ 2 ^ | | q | | + S 2 ) K 2 + ( ϕ 1 | | p | | + ϕ 2 | | q | | + S 1 ) J 2 .
As a result, (47) and (48) follow | | Δ ( p , q ) | | φ , and thus Δ B φ B φ . Now, for ( p 1 , q 1 ) , ( p 2 , q 2 ) P × Q and any τ E , we get
| Δ 1 ( p 1 , q 1 ) ( τ ) Δ 1 ( p 2 , q 2 ) ( τ ) | ρ I 0 + ξ | f ( τ , p 1 ( τ ) , q 1 ( τ ) ) f ( τ , p 2 ( τ ) , q 2 ( τ ) ) | + | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς | g ( ϖ , p 1 ( ϖ ) , q 1 ( ϖ ) ) g ( ϖ , p 2 ( ϖ ) , q 2 ( ϖ ) ) | + ρ I 0 + ξ | f ( T , p 1 ( T ) , q 1 ( T ) ) f ( T , p 2 ( T ) , q 2 ( T ) ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | f ( σ , p 1 ( σ ) , q 1 ( σ ) ) f ( σ , p 2 ( σ ) , q 2 ( σ ) ) | + ρ I 0 + ζ | g ( T , p 1 ( T ) , q 1 ( T ) ) g ( T , p 2 ( T ) , q 2 ( T ) ) | )
( ϕ 1 | | p 1 p 2 | | + ϕ 2 | | q 1 q 2 | | ) T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + ( ϕ 1 ^ | | p 1 p 2 | | + ϕ 2 ^ | | q 1 q 2 | | ) | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 )
( J 1 ( ϕ 1 + ϕ 2 ) + K 1 ( ϕ 1 ^ + ϕ 2 ^ ) ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) .
Similarly, we obtain
| Δ 2 ( p 1 , q 1 ) ( τ ) Δ 2 ( p 2 , q 2 ) ( τ ) | ( ϕ 1 ^ | | p 1 p 2 | | + ϕ 2 ^ | | q 1 q 2 | | ) T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + ( ϕ 1 | | p 1 p 2 | | + ϕ 2 | | q 1 q 2 | | ) | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ( J 2 ( ϕ 1 + ϕ 2 ) + K 2 ( ϕ 1 ^ + ϕ 2 ^ ) ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) ) .
Thus we obtain
| | Δ 1 ( p 1 , q 1 ) ( τ ) Δ 1 ( p 2 , q 2 ) ( τ ) | | ( J 1 ( ϕ 1 + ϕ 2 ) + K 1 ( ϕ 1 ^ + ϕ 2 ^ ) ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) .
In a similar manner,
| | Δ 2 ( p 1 , q 1 ) ( τ ) Δ 2 ( p 2 , q 2 ) ( τ ) | | ( J 2 ( ϕ 1 + ϕ 2 ) + K 2 ( ϕ 1 ^ + ϕ 2 ^ ) ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) .
Hence, using (49) and (50) we can get
| | Δ ( p 1 , q 1 ) ( τ ) Δ ( p 2 , q 2 ) ( τ ) | | ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) .
As a consequence of condition ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) < 1 , Δ is a contraction operator. As an outcome of the Banach fixed point theorem, we can conclude that operator has a unique fixed point, which is the unique solution of the problem (1), and (2). □
For brevity’s sake, we’ll use the following notations:
Ω ^ 1 = J 1 T ρ ξ ρ ξ Γ ( ξ + 1 ) + K 1 ,
Ω ^ 2 = J 2 T ρ ζ ρ ζ Γ ( ζ + 1 ) + K 2 .
Theorem 3.
Assume that f , g : E × R × R R are continuous functions satisfying the assumption ( A 2 ) in Theorem 2. Furthermore, there exist positive constants U 1 , U 2 such that τ E and r i R , i = 1 , 2 .
| f ( τ , r 1 , r 2 ) | U 1 , | g ( τ , r 1 , r 2 ) | U 2 .
If
T ρ ξ ( ϕ 1 + ϕ 2 ) ρ ξ Γ ( ξ + 1 ) + T ρ ζ ( ϕ 1 ^ + ϕ 2 ^ ) ρ ζ Γ ( ζ + 1 ) < 1 ,
then the BVP (1), and (2) has at least one solution on E .
Proof. 
Let us define a closed ball B φ = { ( p , q ) P × Q : | | ( p , q ) | | φ } as follows and split Δ 1 , Δ 2 as:
Δ 1 , 1 ( p , q ) ( τ ) = δ ( τ ) ( E ^ ϵ ρ I 0 + ζ + ς g ( ϖ , p ( ϖ ) , q ( ϖ ) ) ρ I 0 + ξ f ( T , p ( T ) , q ( T ) )
+ E 1 π ρ I 0 + ξ + ϱ f ( σ , p ( σ ) , q ( σ ) ) ρ I 0 + ζ g ( T , p ( T ) , q ( T ) ) ) ,
Δ 1 , 1 ( p , q ) ( τ ) = ρ I 0 + ξ f ( τ , p ( τ ) , q ( τ ) ) ,
Δ 2 , 1 ( p , q ) ( τ ) = δ ( τ ) ( E ^ π ρ I 0 + ξ + ϱ f ( σ , p ( σ ) , q ( σ ) ) ρ I 0 + ζ g ( T , p ( T ) , q ( T ) )
+ E 2 ϵ ρ I 0 + ζ + ς g ( ϖ , p ( ϖ ) , q ( ϖ ) ) ρ I 0 + ξ f ( T , p ( T ) , q ( T ) ) ) ,
Δ 2 , 2 ( p , q ) ( τ ) = ρ I 0 + ξ g ( τ , p ( τ ) , q ( τ ) ) .
In the Banach space P × Q , Δ 1 ( p , q ) ( τ ) = Δ 1 , 1 ( p , q ) ( τ ) + Δ 1 , 2 ( p , q ) ( τ ) , and Δ 2 ( p , q ) ( τ ) = Δ 2 , 1 ( p , q ) ( τ ) + Δ 2 , 2 ( p , q ) ( τ ) on B φ are closed, bounded and convex subsets of P × Q . Let us fix φ max { J 1 U 1 + K 1 U 2 , J 2 U 1 + K 2 U 2 } and show that Δ B φ B φ to verify Krasnoselskii’s theorem [27] condition (i), If we choose p = ( p 1 , p 2 ) , q = ( q 1 , q 2 ) B φ , and utilizing condition (53), we obtain
| Δ 1 , 1 ( p , q ) ( τ ) + Δ 1 , 2 ( p , q ) ( τ ) | ρ I 0 + ξ | f ( τ , p ( τ ) , q ( τ ) ) | + | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς | g ( ϖ , p ( ϖ ) , q ( ϖ ) ) | + ρ I 0 + ξ | f ( T , p ( T ) , q ( T ) ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | f ( σ , p ( σ ) , q ( σ ) ) | + ρ I 0 + ζ | g ( T , p ( T ) , q ( T ) ) | ) U 1 T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + U 2 | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) U 1 J 1 + U 2 K 1 φ .
In a similar manner, we can find that
| Δ 2 , 1 ( p , q ) ( τ ) + Δ 2 , 2 ( p , q ) ( τ ) | U 1 J 2 + U 2 K 2 φ .
Clearly the above two inequalities lead to the fact that Δ 1 ( p , q ) + Δ 2 ( p , q ) B φ . Thus, we define operator ( Δ 1 , 2 , Δ 2 , 2 ) as a contraction-satisfying condition (iii) of Krasnoselskii’s theorem [27]. For ( p 1 , q 1 ) , ( p 2 , q 2 ) B φ , we have
| Δ 1 , 2 ( p 1 , q 1 ) ( τ ) Δ 1 , 2 ( p 2 , q 2 ) ( τ ) | ρ 1 ξ Γ ( ξ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) 1 ξ × | f ( θ , p 1 ( θ ) , q 1 ( θ ) ) f ( θ , p 2 ( θ ) , q 2 ( θ ) ) | d θ T ρ ξ ρ ξ Γ ( ξ + 1 ) ( ϕ 1 | | p 1 p 2 | | + ϕ 2 | | q 1 q 2 | | )
and
| Δ 2 , 1 ( p 1 , q 1 ) ( τ ) Δ 2 , 1 ( p 2 , q 2 ) ( τ ) | ρ 1 ζ Γ ( ζ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) 1 ζ × | g ( θ , p 1 ( θ ) , q 1 ( θ ) ) g ( θ , p 2 ( θ ) , q 2 ( θ ) ) | d θ T ρ ζ ρ ζ Γ ( ζ + 1 ) ( ϕ 1 ^ | | p 1 p 2 | | + ϕ 2 ^ | | q 1 q 2 | | ) .
As a result (59) and (60),
| ( Δ 1 , 2 , Δ 2 , 2 ) ( p 1 , q 1 ) ( τ ) ( Δ 1 , 2 , Δ 2 , 2 ) ( p 2 , q 2 ) ( τ ) | T ρ ξ ( ϕ 1 + ϕ 2 ) ρ ξ Γ ( ξ + 1 ) + T ρ ζ ( ϕ 1 ^ + ϕ 2 ^ ) ρ ζ Γ ( ζ + 1 ) ( | | p 1 p 2 | | + | | q 1 q 2 | | ) ,
is a contraction by (54). Therefore, condition (iii) of the Theorem is satisfied. Following that, we can establish that the operator ( Δ 1 , 1 , Δ 2 , 1 ) satisfies the Krasnoselskii theorem’s [27] condition (ii). We can infer the continuous existence of the ( Δ 1 , 1 , Δ 2 , 1 ) operator by examining the continuity of the f , g functions. For each ( p , q ) B φ we have
| Δ 1 , 1 ( p , q ) ( τ ) | | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς | g ( ϖ , p ( ϖ ) , q ( ϖ ) ) | + ρ I 0 + ξ | f ( T , p ( T ) , q ( T ) ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | f ( σ , p ( σ ) , q ( σ ) ) | + ρ I 0 + ζ | g ( T , p ( T ) , q ( T ) ) | )
U 1 T ρ ξ ( | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + U 2 | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) = Ω ^ 1 ,
| Δ 2 , 1 ( p , q ) ( τ ) | U 2 T ρ ζ ( | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + U 1 | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) = Ω ^ 2 ,
which leads to
| | ( Δ 1 , 1 , Δ 2 , 1 ) ( p , q ) | | Ω ^ 1 + Ω 2 ^ .
From the above inequalities, the set ( Δ 1 , 1 , Δ 2 , 1 ) B φ is uniformly bounded. The following step will demonstrate that the set ( Δ 1 , 1 , Δ 2 , 1 ) B φ is equicontinuous. For τ 1 , τ 2 E with τ 1 < τ 2 and for any ( p , q ) B φ we get
| Δ 1 , 1 ( p , q ) ( τ 2 ) Δ 1 , 1 ( p , q ) ( τ 1 ) |
| δ ( τ 2 ) δ ( τ 1 ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς | g ( ω , p ( ω ) , q ( ω ) ) | + ρ I 0 + ξ | f ( T , p ( T ) , q ( T ) ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | f ( σ , p ( σ ) , q ( σ ) ) | + ρ I 0 + ζ | g ( T , p ( T ) , q ( T ) ) | ) | δ ( τ 2 ) δ ( τ 1 ) ( U 1 T ρ ξ ( | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + U 2 | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) ) .
Likewise, we obtain
| Δ 2 , 1 ( p , q ) ( τ 2 ) Δ 2 , 1 ( p , q ) ( τ 1 ) | | δ ( τ 2 ) δ ( τ 1 ) ( U 2 T ρ ζ ( | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + U 1 | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ) .
Therefore | ( Δ 1 , 1 , Δ 2 , 1 ( τ 2 ) ) ( Δ 1 , 1 , Δ 2 , 1 ( τ 1 ) ) | 0 as τ 2 τ 1 independent of ( p , q ) B φ Thus the set ( Δ 1 , 1 , Δ 2 , 1 ) B φ is equicontinuous. As an outcome, the Arzela–Ascoli theorem implies that the operator ( Δ 1 , 1 , Δ 2 , 1 ) is compact on B φ . Krasnoselskii’s theorem [27] statement leads us to the conclusion that the problem (1) and (2) has at least one solution on E . □

4. Example

Consider the following Liouville–Caputo type generalized FDEs coupled system:
C 3 4 D 0 + 5 4 p ( τ ) = f ( τ , p ( τ ) , q ( τ ) ) , τ E : = [ 0 , 1 ] , C 3 4 D 0 + 31 20 q ( τ ) = g ( τ , p ( τ ) , q ( τ ) ) , τ E : = [ 0 , 1 ] ,
supplemented with boundary conditions:
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( 1 ) = 1 6 3 4 I 13 20 q ( 7 10 ) , q ( 1 ) = 1 7 3 4 I 17 20 p ( 1 2 ) ,
where ξ = 5 4 , ζ = 31 20 , ρ = 3 4 , T = 1 , ϵ = 1 6 , ϖ = 7 10 , π = 1 7 , σ = 1 2 , ς = 13 20 , ϱ = 17 20 and
f ( τ , p ( τ ) , q ( τ ) ) = ( 1 + τ ) 30 | p ( τ ) | 1 + | p ( τ ) | + 1 3 cos ( q ( τ ) ) + 3 τ ,
g ( τ , p ( τ ) , q ( τ ) ) = e τ 25 τ + 1 5 + 1 6 cos ( p ( τ ) ) + | q ( τ ) | 1 + | q ( τ ) | .
With ψ 0 = 1 10 , ψ 1 = 1 30 , ψ 2 = 1 90 , ψ 0 ^ = 1 125 , ψ 1 ^ = 1 25 , and ψ 2 ^ = 1 150 , the functions f and g clearly satisfy the ( A 1 ) condition. Next, we find that ( J 1 ) = 2.5370237266984113 , ( K 1 ) = 0.17111607453629377 , J 2 = 0.0906406939922634 , K 2 = 2.274156747108814 , J i , K i ( i = 1 , 2 ) are respectively given by (35),(36),(37) and (38), based on the data available. Thus ψ 1 ( J 1 + J 2 ) + ψ 1 ^ ( K 1 + K 2 ) 0.18539972688882678 < 1 , ψ 2 ( J 1 + J 2 ) + ψ 2 ^ ( K 1 + K 2 ) 0.04549809015197488 < 1 , all the conditions of Theorem 1 are satisfied, and there is at least one solution for problem (61) and (62) on [ 0 , 1 ] with f and g given by (63) and (64) respectively.
In addition, we’ll use
f ( τ , p ( τ ) , q ( τ ) ) = τ 3 + 3 4 ( τ + 16 ) + | p ( τ ) | 1 + | p ( τ ) | + 2 75 cos ( q ( τ ) ) ,
g ( τ , p ( τ ) , q ( τ ) ) = ( 1 + e τ ) 4 + 19 400 cos ( p ( τ ) ) + 1 60 | q ( τ ) | 1 + | q ( τ ) | ,
to demonstrate Theorem 2. It is simple to demonstrate that f and g are continuous and satisfy the assumption ( A 2 ) with ϕ 1 = 3 64 , ϕ 2 = 2 75 , ϕ 1 ^ = 19 400 and ϕ 2 ^ = 1 60 . All the assumptions of Theorem 2 are also satisfied with ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) 0.35014782699385444 < 1 . As a result, Theorem 2 holds true, and the problem (61) and (62) with f and g given by (65) and (66) respectively, has a unique solution on [0,1].

5. Ulam–Hyers Stability Results for the Problem (1) and (2)

The U–H stability of the solutions to the BVP (1) and (2) will be discussed in this section using the integral representation of their solutions defined by
p ( τ ) = Δ 1 ( p , q ) ( τ ) , q ( τ ) = Δ 2 ( p , q ) ( τ ) ,
where Δ 1 and Δ 2 are given by (33) and (34). Consider the following definitions of nonlinear operators
H 1 , H 2 C ( E , R ) × C ( E , R ) C ( E , R ) ,
C ρ D 0 + ξ p ( τ ) f ( τ , p ( τ ) , q ( τ ) ) = H 1 ( p , q ) ( τ ) , τ E , C ρ D 0 + ζ q ( τ ) g ( τ , p ( τ ) , q ( τ ) ) = H 1 ( p , q ) ( τ ) , τ E .
It considered the following inequalities for some λ 1 ^ , λ ^ 2 > 0 :
| | H 1 ( p , q ) | | λ 1 ^ , | | H 2 ( p , q ) | | λ 2 ^ .
Definition 4.
The coupled system (1) and (2) is said to be U–H stable if V 1 , V 2 > 0 and there exists a unique solution ( p , q ) C ( E , R ) of a problem (1) and (2) with
| | ( p , q ) ( p * , q * ) | | V 1 λ 1 ^ + V 2 λ 2 ^ ,
( p , q ) C ( E , R ) of inequality (68).
Theorem 4.
Assume that ( A 2 ) holds. Then the problem (1) and (2) is U–H stable.
Proof. 
Let ( p , q ) C ( E , R ) × C ( E , R ) be the (1)–(2) solution of the problem that satisfies (33) and (34). Let ( p , q ) be any solution that meets the condition (68):
C ρ D 0 + ξ p ( τ ) = f ( τ , p ( τ ) , q ( τ ) ) + H 1 ( p , q ) ( τ ) , τ E , C ρ D 0 + ζ q ( τ ) = g ( τ , p ( τ ) , q ( τ ) ) + H 1 ( p , q ) ( τ ) , τ E ,
so,
p * ( τ ) = Δ 1 ( p * , q * ) ( τ ) + ρ I 0 + ξ H 1 ( p , q ) ( τ ) + δ ( τ ) ( E ^ ϵ ρ I 0 + ζ + ς H 2 ( p , q ) ( ϖ ) ρ I 0 + ξ H 1 ( p , q ) ( T ) + E 1 π ρ I 0 + ξ + ϱ H 1 ( p , q ) ( σ ) ρ I 0 + ζ H 2 ( p , q ) ( T ) ) .
It follows that
| Δ 1 ( p * , q * ) ( τ ) p * ( τ ) | ρ I 0 + ξ | H 1 ( p , q ) ( τ ) | + | δ ( τ ) | ( | E ^ | | ϵ | ρ I 0 + ζ + ς | H 2 ( p , q ) ( ϖ ) | + ρ I 0 + ξ | H 1 ( p , q ) ( T ) | + | E 1 | | π | ρ I 0 + ξ + ϱ | H 1 ( p , q ) ( σ ) | + ρ I 0 + ζ | H 2 ( p , q ) ( T ) | ) λ ^ 1 T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) + λ ^ 2 | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) J 1 λ 1 ^ + K 1 λ 2 ^ .
Similarly, we obtain
| Δ 2 ( p * , q * ) ( τ ) q * ( τ ) | λ 2 ^ T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) + λ 1 ^ | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | σ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) J 2 λ ^ 1 + K 2 λ ^ 2 ,
where J 1 , K 1 , J 2 , and K 2 are defined in (35)–(38), respectively. As an outcome, we deduce from operator Δ ’s fixed-point property, which is defined by (33) and (34),
| p ( τ ) p * ( τ ) | = | p ( τ ) Δ 1 ( p * , q * ) ( τ ) + Δ 1 ( p * , q * ) ( τ ) p * ( τ ) | | Δ 1 ( p , q ) ( τ ) Δ 1 ( p * , q * ) ( τ ) | + | Δ 1 ( p * , q * ) ( τ ) p * ( τ ) | ( ( J 1 ϕ 1 + K 1 ϕ 1 ^ ) + ( J 1 ϕ 2 + K 1 ϕ 2 ^ ) ) | | ( p , q ) ( p * , q * ) | | + J 1 λ 1 ^ + K 1 λ 2 ^ .
| q ( τ ) q * ( τ ) | = | q ( τ ) Δ 2 ( p * , q * ) ( τ ) + Δ 2 ( p * , q * ) ( τ ) q * ( τ ) | | Δ 2 ( p , q ) ( τ ) Δ 2 ( p * , q * ) ( τ ) | + | Δ 2 ( p * , q * ) ( τ ) q * ( τ ) | ( ( J 2 ϕ 1 + K 2 ϕ 1 ^ ) + ( J 2 ϕ 2 + K 2 ϕ 2 ^ ) ) | | ( p , q ) ( p * , q * ) | | + J 2 λ 1 ^ + K 2 λ 2 ^ .
From the above Equations (69) and (70) it follows that
| | ( p , q ) ( p * , q * ) | | ( J 1 + J 2 ) λ 1 ^ + ( K 1 + K 2 ) λ 2 ^ + ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) | | ( p , q ) ( p * , q * ) | | .
| | ( p , q ) ( p * , q * ) | | ( J 1 + J 2 ) λ 1 ^ + ( K 1 + K 2 ) λ 2 ^ 1 ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) V 1 λ 1 ^ + V 2 λ 2 ^ ,
with
V 1 = J 1 + J 2 1 ( ( J 1 + | J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + | K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) ,
V 2 = K 1 + K 2 1 ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) .
Hence, the problem (1)–(2) is U–H stable. □

6. Example

Consider the following Liouville–Caputo type generalized FDEs coupled system:
C 19 20 D 0 + 5 4 p ( τ ) = τ 2 + 1 5 ( τ + 25 ) | p ( τ ) | 1 + | p ( τ ) | + 3 80 cos ( q ( τ ) ) , τ [ 0 , 1 ] , C 19 20 D 0 + 31 20 q ( τ ) = τ 5 + 17 300 cos ( p ( τ ) ) + 1 70 | q ( τ ) | 1 + | q ( τ ) | , τ [ 0 , 1 ] ,
supplemented with boundary conditions:
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( 1 ) = 5 6 19 20 I 13 20 q ( 9 20 ) , q ( 1 ) = 6 7 19 20 I 17 20 p ( 13 20 ) ,
where ξ = 5 4 , ζ = 31 20 , ρ = 19 20 , T = 1 , ϵ = 5 6 , ϖ = 9 20 , π = 6 7 , σ = 13 20 , ς = 13 20 , ϱ = 17 20 and
| f ( τ , p 1 ( τ ) , q 1 ( τ ) ) f ( τ , p 2 ( τ ) , q 2 ( τ ) ) | = 1 125 | p 1 ( τ ) p 2 ( τ ) | + 3 80 | q 1 ( τ ) q 2 ( τ ) | ,
| g ( τ , p 1 ( τ ) , q 1 ( τ ) ) g ( τ , p 2 ( τ ) , q 2 ( τ ) ) | = 17 300 | p 1 ( τ ) p 2 ( τ ) | + 1 70 | q 1 ( τ ) q 2 ( τ ) | .
With ϕ 1 = 1 125 , ϕ 2 = 3 80 , ϕ 1 ^ = 17 300 , and ϕ 2 ^ = 1 70 , the functions f and g clearly satisfy the ( A 2 ) condition. Next, we find that ( J 1 ) = 1.9529307397739033 , ( K 1 ) = 0.21135021378560123 , J 2 = 0.42682560046779994 , K 2 = 1.6225052940838325 , J i , K i ( i = 1 , 2 ) are respectively given by (35),(36),(37) and (38), based on the data available. Thus ( ( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) ) 0.2383953280869716 < 1 , all the conditions of Theorem 5.2 are satisfied, and there is a unique solution for problem (71) and (72) on [ 0 , 1 ] , which is stable for Ulam–Hyers, with f and g given by (73) and (74) respectively.

7. Existence Results for the Problem (1) and (75)

Furthermore, we are investigating the system (1) under the following conditions:
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) = ϵ ρ 1 ς Γ ( ς ) 0 ϖ θ ρ 1 ( ϖ ρ θ ρ ) 1 ς q ( θ ) d θ , q ( T ) = π ρ I 0 + ϱ p ( ϖ ) = π ρ 1 ϱ Γ ( ϱ ) 0 ϖ θ ρ 1 ( ϖ ρ θ ρ ) 1 ϱ p ( θ ) d θ , 0 < ϖ < T .
Bear in mind that the conditions (2) contain strips of varying lengths, whereas the one in (75) contains only one strip of the same length ( 0 , ϖ ) . We introduce the following notations for computational ease:
E 1 = ϵ ϖ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) , E 2 = π ϖ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) , E ^ = T ρ ρ ,
G = E ^ 2 E 1 E 2 0 ,
δ ( τ ) = ( τ ρ ρ G ) .
Lemma 4.
Given the functions f ^ , g ^ C ( 0 , T ) L ( 0 , T ) , p , q A C γ 2 ( E ) and Λ 0 . Then the solution of the coupled BVP:
C ρ D 0 + ξ p ( τ ) = f ^ ( τ ) , τ E : = [ 0 , T ] , C ρ D 0 + ζ q ( τ ) = g ^ ( τ ) , τ E : = [ 0 , T ] , p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) , q ( T ) = π ρ I 0 + ϱ p ( ϖ ) , 0 < ϖ < T ,
is given by
p ( τ ) = ρ I 0 + ξ f ^ ( τ ) + δ ( τ ) ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + π ρ I 0 + ξ + ς f ^ ( ϖ ) ρ I 0 + ζ g ^ ( T )
and
q ( τ ) = ρ I 0 + ζ g ^ ( τ ) + δ ( τ ) π ρ I 0 + ζ + ς f ^ ( ϖ ) ρ I 0 + ζ g ^ ( T ) + ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) .
Proof. 
When ρ I 0 + ξ , ρ I 0 + ζ are applied to the FDEs in (79) and Lemma 4 is used the solution of the FDEs in (79) for τ E is
p ( τ ) = ρ I 0 + ξ f ^ ( τ ) + a 1 + a 2 τ ρ ρ = ρ 1 ξ Γ ( ξ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) ξ 1 f ^ ( θ ) d θ + a 1 + a 2 τ ρ ρ ,
q ( τ ) = ρ I 0 + ζ g ^ ( τ ) + b 1 + b 2 τ ρ ρ = ρ 1 ζ Γ ( ζ ) 0 τ θ ρ 1 ( τ ρ θ ρ ) ζ 1 g ^ ( θ ) d θ + b 1 + b 2 τ ρ ρ ,
respectively, for some a 1 , a 2 , b 1 , b 2 R . Making use of the boundary conditions p ( 0 ) = q ( 0 ) = 0 in (82) and (83) respectively, we get a 1 = b 1 = 0 . We obtain by using the generalized integral operators ρ I 0 + ϱ , ρ I 0 + ζ (82) and (83) respectively,
ρ I 0 + ϱ p ( τ ) = ρ I 0 + ξ + ϱ f ^ ( τ ) + a 1 τ ρ ϱ ρ ϱ Γ ( ϱ + 1 ) + a 2 τ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) ,
ρ I 0 + ς q ( τ ) = ρ I 0 + ζ + ς g ^ ( τ ) + b 1 τ ρ ς ρ ς Γ ( ς + 1 ) + b 2 τ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) ,
which, when combined with the boundary conditions p ( T ) = ϵ ρ I 0 + ς q ( ϖ ) , q ( T ) = π ρ I 0 + ϱ p ( ϖ ) , gives the following results:
ρ I 0 + ξ f ^ ( T ) + a 1 + a 2 T ρ ρ = ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) + b 1 ϵ ϖ ρ ς ρ ς Γ ( ς + 1 ) + b 2 ϵ ϖ ρ ( ς + 1 ) ρ ς + 1 Γ ( ς + 2 ) ,
ρ I 0 + ζ g ^ ( T ) + b 1 + b 2 T ρ ρ = π ρ I 0 + ξ + ϱ f ^ ( ϖ ) + a 1 π ϖ ρ ϱ ρ ϱ Γ ( ϱ + 1 ) + a 2 π ϖ ρ ( ϱ + 1 ) ρ ϱ + 1 Γ ( ϱ + 2 ) .
Next, we obtain
a 2 E ^ b 2 E 1 = ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) ,
b 2 E ^ a 2 E 2 = π ρ I 0 + ξ + ϱ f ^ ( ϖ ) ρ I 0 + ζ g ^ ( T ) ,
by employing the notations (76)–(78) in (86) and (87) respectively. We find that when we solve the system of Equations (88) and (89) for a 2 and b 2 ,
a 2 = 1 G E ^ ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + E 1 π ρ I 0 + ξ + ϱ f ^ ( ϖ ) ρ I 0 + ζ g ^ ( T ) ,
b 2 = 1 G E 2 ϵ ρ I 0 + ζ + ς g ^ ( ϖ ) ρ I 0 + ξ f ^ ( T ) + E ^ π ρ I 0 + ξ + ϱ f ^ ( ϖ ) ρ I 0 + ζ g ^ ( T ) .
Substituting the values of a 1 , a 2 , b 1 , b 2 in (82) and (83) respectively, we get the solution for (79). □
For brevity’s sake, we’ll use the following notations:
J 1 = T ρ ξ ( 1 + | δ | | E ^ | ) ρ ξ Γ ( ξ + 1 ) + | δ | | π | | E 1 | ϖ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ,
K 1 = | δ | | E 1 | T ρ ζ ρ ζ Γ ( ζ + 1 ) + | E ^ | | ϵ | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) ,
J 2 = | δ | T ρ ξ | E 2 | ρ ξ Γ ( ξ + 1 ) + | π | | E ^ | ϖ ρ ( ξ + ϱ ) ρ ξ + ϱ Γ ( ξ + ϱ + 1 ) ,
K 2 = T ρ ζ ( 1 + | δ | | E ^ | ) ρ ζ Γ ( ζ + 1 ) + | δ | | ϵ | | E 2 | ϖ ρ ( ζ + ς ) ρ ζ + ς Γ ( ζ + ς + 1 ) .
To finish up, we will go over the results of existence, uniqueness, and Ulam–Hyers stability for problems (1) and (75), respectively. For reasons that are similar to those in Section 3, Section 4, Section 5 and Section 6, we are not providing the proof.
Corollary 1.
Assume that f , g : E × R × R R are continuous functions satisfying the condition: ( A 1 ) there exists constants ψ m , ψ m ^ 0 ( m = 1 , 2 ) and ψ 0 , ψ 0 ^ > 0 such that
| f ( τ , o 1 , o 2 ) | ψ 0 + ψ 1 | o 1 | + ψ 2 | o 2 | , | g ( τ , o 1 , o 2 ) | ψ 0 ^ + ψ 1 ^ | o 1 | + ψ 2 ^ | o 2 | , o m R , m = 1 , 2 .
If ψ 1 ( J 1 ^ + J 2 ^ ) + ψ 1 ^ ( K 1 ^ + K 2 ^ ) < 1 , ψ 2 ( J 1 ^ + J 2 ^ ) + ψ 2 ^ ( K 1 ^ + K 2 ^ ) < 1 . Then at least one solution for the BVP (1) and (75) on E , where J 1 ^ , K 1 ^ , J 2 ^ , K 2 ^ are given by (92)–(95) respectively.
Corollary 2.
Assume that f , g : E × R × R R are continuous functions satisfying the condition: ( A 2 ) there exists constants ϕ m , ϕ m ^ 0 ( m = 1 , 2 ) such that
| f ( τ , o 1 , o 2 ) f ( τ , o ^ 1 , o ^ 2 ) | ϕ 1 | o 1 o ^ 1 | + ϕ 2 | o 2 o ^ 2 | , | g ( τ , o 1 , o 2 ) g ( τ , o ^ 1 , o ^ 2 ) | ϕ 1 ^ | o 1 o ^ 1 | + ϕ 2 ^ | o 2 o ^ 2 | , o m , o ^ m R , m = 1 , 2 .
Moreover, there exist S 1 , S 2 > 0 such that | f ( τ , 0 , 0 ) | S 1 , | f ( τ , 0 , 0 ) | S 2 , Then, given that
( J 1 + J 2 ) ( ϕ 1 + ϕ 2 ) + ( K 1 + K 2 ) ( ϕ 1 ^ + ϕ 2 ^ ) < 1 ,
the BVP (1) and (75) has a unique solution on E , where J 1 ^ , K 1 ^ , J 2 ^ , K 2 ^ are given by (92)–(95) respectively.
Corollary 3.
Assume that f , g : E × R × R R are continuous functions satisfying the assumption ( A 2 ) in Theorem 2. Further more, there exist positive constants U 1 , U 2 such that τ E and r i R , i = 1 , 2 .
| f ( τ , r 1 , r 2 ) | U 1 , | g ( τ , r 1 , r 2 ) | U 2 .
If
T ρ ξ ( ϕ 1 + ϕ 2 ) ρ ξ Γ ( ξ + 1 ) + T ρ ζ ( ϕ 1 ^ + ϕ 2 ^ ) ρ ζ Γ ( ζ + 1 ) < 1 ,
then the BVP (1), and (75) has at least one solution on E .
Corollary 4.
Assume that (A2) holds. Then the problem (1) and (75) is Ulam–Hyers stable.

8. Asymmetric Cases

Remark 1.
If ρ = 1 , the problem (1) generalized Liouville–Caputo type reduces to the classical Caputo form.
C D 0 + ξ p ( τ ) = f ( τ , p ( τ ) , q ( τ ) ) , τ G : = [ 0 , T ] , C D 0 + ζ q ( τ ) = g ( τ , p ( τ ) , q ( τ ) ) , τ G : = [ 0 , T ] .
Remark 2.
If ρ = 1 in the boundary conditions (2) and (75) generalized Riemann–Liouville integral boundary conditions reduces to the Riemann–Liouville integral conditions respectively.
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ I 0 + ς q ( ϖ ) = ϵ Γ ( ς ) 0 ϖ ( ϖ θ ) ς 1 q ( θ ) d θ , q ( T ) = π I 0 + ϱ p ( σ ) = π Γ ( ϱ ) 0 σ ( σ θ ) ϱ 1 p ( θ ) d θ , 0 < σ < ϖ < T ,
and
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ I 0 + ς q ( ϖ ) = ϵ Γ ( ς ) 0 ϖ ( ϖ θ ) ς 1 q ( θ ) d θ , q ( T ) = π I 0 + ϱ p ( ϖ ) = π Γ ( ϱ ) 0 ϖ ( ϖ θ ) ϖ 1 p ( θ ) d θ , 0 < ϖ < T .
Remark 3.
If ρ = 1 and ς = ϱ = 1 in the boundary conditions (2) and (75) generalized Riemann–Liouville integral boundary conditions reduces to the classical integral conditions respectively.
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ 0 ϖ q ( θ ) d θ , q ( T ) = π 0 σ p ( θ ) d θ 0 < σ < ϖ < T
and
p ( 0 ) = 0 , q ( 0 ) = 0 , p ( T ) = ϵ 0 ϖ q ( θ ) d θ , q ( T ) = π 0 ϖ p ( θ ) d θ 0 < ϖ < T .

9. Conclusions

This paper employs coupled nonlinear generalized Liouville–Caputo fractional differential equations and Katugampola fractional integral operators to solve a novel class of boundary value problems. Applying the techniques of fixed-point theory to discover the existence criterion for solutions is efficient. While the second outcome provides a sufficient criterion to establish the problem’s unique solution, the first and third results define various criteria for the presence of solutions to the given problem. In the fourth section, the Hyers–Ulam stability of the solution was determined. In the remarks, we have shown the asymmetric cases of the assigned problem. Moreover, the form of the solution in these kinds of remarks can be used to study the positive solution and its asymmetry in more depth. We conclude that our results are novel and can be viewed as an expansion of the qualitative analysis of fractional differential equations. Our results are novel in this configuration and add to the literature on nonlinear coupled generalized Liouville–Caputo fractional differential equations with nonlocal boundary conditions utilizing Katugampola-type integral operators. Future research could focus on various conceptions of stability and existence in relation to a Lotka–Volterra prey-predator system/coupled logistic system.

Author Contributions

Conceptualization, M.S.; formal analysis, M.A. and M.S.; methodology, M.A, M.S., K.A. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Grant No.793].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Britton, N.F. Essential Mathematical Biology; Springer: London, UK, 2003. [Google Scholar]
  2. Ma, Y.; Ji, D. Existence of Solutions to a System of Riemann-Liouville Fractional Differential Equations with Coupled Riemann-Stieltjes Integrals Boundary Conditions. Fractal Fract. 2022, 6, 543. [Google Scholar] [CrossRef]
  3. Theswan, S.; Ntouyas, S.K.; Ahmad, B.; Tariboon, J. Existence Results for Nonlinear Coupled Hilfer Fractional Differential Equations with Nonlocal Riemann–Liouville and Hadamard-Type Iterated Integral Boundary Conditions. Symmetry 2022, 14, 1948. [Google Scholar] [CrossRef]
  4. Klafter, J.; Lim, S.; Metzler, R. Fractional Dynamics: Recent Advances; World Scientific: Singapore, 2012. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  6. Valerio, D.; Machado, J.T.; Kiryakova, V. Some pioneers of the applications of fractional calculus. Fract. Calc. Appl. Anal. 2014, 17, 552–578. [Google Scholar] [CrossRef] [Green Version]
  7. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  8. Kilbas, A.A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science Limited, North-Holland Mathematics Studies: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  9. Bitsadze, A.; Samarskii, A. On some simple generalizations of linear elliptic boundary problems. Soviet Math. Dokl. 1969, 10, 398–400. [Google Scholar]
  10. Ciegis, R.; Bugajev, A. Numerical approximation of one model of bacterial self-organization. Nonlinear Anal. Model. Control. 2012, 17, 253–270. [Google Scholar] [CrossRef] [Green Version]
  11. Subramanian, M.; Alzabut, J.; Baleanu, D.; Samei, M.E.; Zada, A. Existence, uniqueness and stability analysis of a coupled fractional-order differential systems involving hadamard derivatives and associated with multi-point boundary conditions. Adv. Differ. Equ. 2021, 2021, 1–46. [Google Scholar] [CrossRef]
  12. Rahmani, A.; Du, W.S.; Khalladi, M.T.; Kostić, M.; Velinov, D. Proportional Caputo Fractional Differential Inclusions in Banach Spaces. Symmetry 2022, 14, 1941. [Google Scholar] [CrossRef]
  13. Tudorache, A.; Luca, R. Positive Solutions for a Fractional Differential Equation with Sequential Derivatives and Nonlocal Boundary Conditions. Symmetry 2022, 14, 1779. [Google Scholar] [CrossRef]
  14. Ahmad, B.; Alghanmi, M.; Alsaedi, A.; Nieto, J.J. Existence and uniqueness results for a nonlinear coupled system involving caputo fractional derivatives with a new kind of coupled boundary conditions. Appl. Math. Lett. 2021, 116, 107018. [Google Scholar] [CrossRef]
  15. Alsaedi, A.; Alghanmi, M.; Ahmad, B.; Ntouyas, S.K. Generalized liouville–caputo fractional differential equations and inclusions with nonlocal generalized fractional integral and multipoint boundary conditions. Symmetry 2018, 10, 667. [Google Scholar] [CrossRef]
  16. Boutiara, A.; Etemad, S.; Alzabut, J.; Hussain, A.; Subramanian, M.; Rezapour, S. On a nonlinear sequential four-point fractional q-difference equation involving q-integral operators in boundary conditions along with stability criteria. Adv. Differ. Equ. 2021, 2021, 1–23. [Google Scholar] [CrossRef]
  17. Baleanu, D.; Alzabut, J.; Jonnalagadda, J.; Adjabi, Y.; Matar, M. A coupled system of generalized sturm–liouville problems and langevin fractional differential equations in the framework of nonlocal and nonsingular derivatives. Adv. Differ. Equ. 2020, 2020, 1–30. [Google Scholar] [CrossRef]
  18. Muthaiah, S.; Baleanu, D. Existence of solutions for nonlinear fractional differential equations and inclusions depending on lower-order fractional derivatives. Axioms 2020, 9, 44. [Google Scholar] [CrossRef]
  19. Saeed, A.M.; Abdo, M.S.; Jeelani, M.B. Existence and Ulam–Hyers stability of a fractional-order coupled system in the frame of generalized Hilfer derivatives. Mathematics 2021, 9, 2543. [Google Scholar] [CrossRef]
  20. Ahmad, D.; Agarwal, R.P.; Rahman, G.U.R. Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations. Symmetry 2022, 14, 1342. [Google Scholar] [CrossRef]
  21. Samadi, A.; Ntouyas, S.K.; Tariboon, J. On a nonlocal coupled system of Hilfer generalized proportional fractional differential equations. Symmetry 2022, 14, 738. [Google Scholar] [CrossRef]
  22. Awadalla, M.; Abuasbeh, K.; Subramanian, M.; Manigandan, M. On a System of ψ-Caputo Hybrid Fractional Differential Equations with Dirichlet Boundary Conditions. Mathematics 2022, 10, 1681. [Google Scholar] [CrossRef]
  23. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef] [Green Version]
  24. Katugampola, U.N. A new approach to generalized fractional derivatives. arXiv 2011, arXiv:1106.0965. [Google Scholar]
  25. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their caputo modification. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. [Google Scholar] [CrossRef]
  26. Granas, A.; Dugundji, J. Fixed Point Theory; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar]
  27. Krasnoselskiı, M. Two remarks on the method of successive approximations, uspehi mat. Nauk 1955, 10, 123–127. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Awadalla, M.; Subramanian, M.; Abuasbeh, K.; Manigandan, M. On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions. Symmetry 2022, 14, 2273. https://doi.org/10.3390/sym14112273

AMA Style

Awadalla M, Subramanian M, Abuasbeh K, Manigandan M. On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions. Symmetry. 2022; 14(11):2273. https://doi.org/10.3390/sym14112273

Chicago/Turabian Style

Awadalla, Muath, Muthaiah Subramanian, Kinda Abuasbeh, and Murugesan Manigandan. 2022. "On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions" Symmetry 14, no. 11: 2273. https://doi.org/10.3390/sym14112273

APA Style

Awadalla, M., Subramanian, M., Abuasbeh, K., & Manigandan, M. (2022). On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions. Symmetry, 14(11), 2273. https://doi.org/10.3390/sym14112273

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