Next Article in Journal
Hydroxyl Spectroscopy of Laboratory Air Laser-Ignition
Next Article in Special Issue
Analysis of Sequential Caputo Fractional Differential Equations versus Non-Sequential Caputo Fractional Differential Equations with Applications
Previous Article in Journal
A Possible Explanation of the Proton Radius Puzzle Based on the Second Flavor of Muonic Hydrogen Atoms
Previous Article in Special Issue
Green’s Functions for a Fractional Boundary Value Problem with Three Terms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of a Nonlinear Coupled (k, ψ)–Hilfer Fractional Differential System with Coupled (k, ψ)–Riemann–Liouville Fractional Integral Boundary Conditions

1
Department of Mathematics, Miyaneh Branch, Islamic Azad University, Miyaneh 5315836511, 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
Intelligent and Nonlinear Dynamic Innovations, Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Author to whom correspondence should be addressed.
Foundations 2022, 2(4), 918-933; https://doi.org/10.3390/foundations2040063
Submission received: 25 September 2022 / Revised: 13 October 2022 / Accepted: 14 October 2022 / Published: 18 October 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
This paper is concerned with the existence of solutions for a new boundary value problem of nonlinear coupled ( k , ψ ) –Hilfer fractional differential equations subject to coupled ( k , ψ ) –Riemann–Liouville fractional integral boundary conditions. We prove two existence results by applying the Leray–Schauder alternative, and Krasnosel’skiĭ’s fixed-point theorem under different criteria, while the third result, concerning the uniqueness of solutions for the given problem, relies on the Banach’s contraction mapping principle. Examples are included for illustrating the abstract results.

1. Introduction

Fractional-order boundary value problems involving classical, nonlocal, multi-point, and integral boundary conditions have attracted significant interest in the recent years, as such problems appear in diverse disciplines of engineering and applied sciences. For theoretical details of the topic, see the monographs [1,2,3]. Recently, some existence results for nonlocal integro-multipoint boundary value problems involving Hilfer-type fractional integro-differential equations and inclusions were presented in [4]. In [5], the authors analyzed a multipoint boundary value problem involving a Caputo–Hadamard fractional derivative operator, while a nonlinear Caputo–Hadamard fractional initial value problem was studied in [6].
Great interest has also been shown in studying the coupled systems of fractional-order differential equations, as such systems constitute the mathematical models of many physical phenomena occurring in bio-engineering [7], fractional dynamics [8], and financial economics [9], etc. Let us now dwell on some recent works on boundary value problems involving coupled fractional differential systems. In [10,11], the authors proved some interesting results for coupled systems of sequential fractional differential equations involving a ψ –Hilfer fractional derivative. Boundary value problems for Hilfer-type sequential fractional differential equations and inclusions with Riemann–Stieltjes integral multi-strip boundary conditions were studied in [12].
It is interesting to notice that there do exist different definitions of fractional derivative operators such as Riemann–Liouville, Caputo–Liouville, Hadamard, Erdélyi–Kober, ψ –Riemann–Liouville [1], Hilfer [13], k–Riemann–Liouville [14], ψ –Hilfer [15], ( k , ψ ) –Riemann–Liouville [14] and ( k , ψ ) –Hilfer [16], to name a few.
The topic of ( k , ψ ) –Hilfer-type fractional boundary value problems constitutes an interesting and useful area of research, as a ( k , ψ ) –Hilfer-type fractional derivative operator generalizes many well-known fractional derivative operators; for details, see [17] and the references cited therein. Recently, in [18], the authors studied the existence of solutions for a nonlocal nonlinear ( k , ψ ) –Hilfer-type fractional boundary value problem. For some more results on ( k , ψ ) –Hilfer fractional boundary value problems, see [19]. In [20], the authors studied a coupled system of nonlinear fractional differential equations involving the ( k , φ ) –Hilfer fractional derivative operators with nonlocal multipoint boundary conditions:
k , H D α ¯ , β ¯ ; φ w ( θ ) = f ( θ , w ( θ ) , z ( θ ) ) , θ ( c , d ] , k , H D α 1 , β 1 ; φ z ( θ ) = f 1 ( θ , w ( θ ) , z ( θ ) ) , θ ( c , d ] , w ( c ) = 0 , w ( d ) = i = 1 m λ i z ( ξ i ) , z ( c ) = 0 , z ( d ) = j = 1 k μ j w ( η j ) ,
where k , H D α ¯ , β ¯ ; φ , k , H D α 1 , β 1 ; φ denote the ( k , φ ) –Hilfer fractional derivative operators of orders α ¯ , α 1 , 1 < α ¯ , α 1 < 2 and parameters β ¯ , β 1 , 0 β ¯ , β 1 1 , respectively; f , f 1 : [ c , d ] × R R are continuous functions, λ i , μ j R ; and a < ξ i , η j < b , i = 1 , 2 , , m , j = 1 , 2 , , k .
Exploring the literature on ( k , ψ ) –Hilfer fractional differential systems, it is found that the only work dealing with such a coupled system was discussed in [20]. In order to develop this new area of research, we consider a nonlinear coupled system of ( k , ψ ) –Hilfer fractional differential equations of different orders in ( 1 , 2 ] , complemented with coupled ( k , ψ ) –Riemann-Liouville fractional integral boundary conditions. In precise terms, we investigate the existence of solutions for the following problem:
k , H D α ^ , β ^ ; ψ r ^ ( t ) = f ( t , r ^ ( t ) , z ^ ( t ) ) , t ( a , b ] , k , H D p ^ , q ^ ; ψ z ^ ( t ) = f ˜ ( t , r ^ ( t ) , z ^ ( t ) ) , t ( a , b ] , r ^ ( a ) = 0 , r ^ ( b ) = λ z ^ ( ξ ) + μ k I v , ψ z ^ ( σ ) , z ^ ( a ) = 0 , z ^ ( b ) = ς r ^ ( u ) + θ k I w , ψ r ^ ( τ ) ,
where k , H D α ^ , β ^ ; ψ , k , H D p ^ , q ^ ; ψ denote the ( k , ψ ) –Hilfer fractional derivative operators of orders α ^ , p ^ , 1 < α ^ , p ^ < 2 and parameters β ^ , q ^ , 0 β ^ , q ^ 1 , respectively, f , f ˜ : [ a , b ] × R R are continuous functions, λ , μ , ς , θ R , k I v , ψ , k I w , ψ are the ( k , ψ ) –Riemann–Liouville fractional integrals of order v , w > 0 , and a < ξ , σ , u , τ < b . The existence and uniqueness results for the nonlinear system (1) will be established with the aid of the Leray–Schauder alternative, Krasnosel’skiĭ’s fixed-point theorem and Banach’s contraction mapping principle. The tools of the fixed point theory provide a decent platform from which to develop the existence theory for initial and boundary value problems once such problems are convertible to fixed-point problems. In this context, the present strategy is preferred over the other existing methods such as Lie symmetry, quasilinearization, etc.
Here, we emphasize that the problem (1) formulated and investigated in the present work is novel and the existence and uniqueness results derived for this problem enrich the literature on nonlocal integral boundary value problems involving ( k , ψ ) –Hilfer fractional differential and ( k , ψ ) –Riemann–Liouville integral operators of different orders. For more details, see the Conclusion section.
The structure of this paper is as follows. In Section 2, we recall some definitions and prove an auxiliary result which plays a fundamental role in transforming the given nonlinear system into a fixed-point problem. The main results are proved in Section 3. In Section 4, we provide numerical examples for illustration of our theoretical results.

2. Preliminaries

Let us first recall the definitions of a ( k , ψ ) –Riemann–Liouville fractional integral and ( k , ψ ) –Hilfer fractional derivative operators.
Definition 1
([21]). Let ψ : [ a , b ] R be an increasing function with ψ ( t ) 0 for all t [ a , b ] . Then, the ( k , ψ ) –Riemann–Liouville fractional integral of order α > 0 ( α R ) of a function h L 1 ( [ a , b ] , R ) is given by
k I a + α ; ψ h ( t ) = 1 k Γ k ( α ) a t ψ ( u ) ( ψ ( t ) ψ ( u ) ) α k 1 h ( u ) d u , k > 0 .
Definition 2
([16]). Let α , k R + = ( 0 , ) , β [ 0 , 1 ] , ψ C n ( [ a , b ] , R ) , ψ ( t ) 0 , t [ a , b ] and h C n ( [ a , b ] , R ) . Then, the ( k , ψ ) –Hilfer fractional derivative of the function h of order α and type β , is defined by
k , H D α , β ; ψ h ( t ) = k I a + β ( n k α ) ; ψ k ψ ( t ) d d t n k I a + ( 1 β ) ( n k α ) ; ψ h ( t ) , n = α k .
Now, we prove an auxiliary lemma that plays a key role in the study of the nonlinear system (1).
Lemma 1.
Assume that a < b , k > 0 , 1 < α ^ , p ^ 2 , β ^ , q ^ [ 0 , 1 ] , ϑ k = α ^ + β ^ ( 2 k α ^ ) , η k = p ^ + q ^ ( 2 k p ^ ) , h , h ˜ C 2 ( [ a , b ] , R ) and
A : = A 1 A 4 A 2 A 3 0 .
Then, the pair ( r ^ , z ^ ) is a solution of the nonlocal ( k , ψ ) –Hilfer fractional system
k , H D α ^ , β ^ ; ψ r ^ ( t ) = h ( t ) , t ( a , b ] , k , H D p ^ , q ^ ; ψ z ^ ( t ) = h ˜ ( t ) , t ( a , b ] , r ^ ( a ) = 0 , r ^ ( b ) = λ z ^ ( ξ ) + μ k I v , ψ z ^ ( σ ) , z ^ ( a ) = 0 , z ^ ( b ) = ς r ^ ( u ) + θ k I w , ψ r ^ ( τ ) ,
if and only if
r ^ ( t ) = k I α ^ , ψ h ( t ) + ( ψ ( t ) ψ ( a ) ) ϑ k k 1 A Γ k ( ϑ k ) [ A 4 λ k I p ^ , ψ h ˜ ( ξ ) + μ k I p ^ + v , ψ h ˜ ( σ ) k I α ^ , ψ h ( b ) + A 2 ς k I α ^ , ψ h ( u ) + θ k I α ^ + w , ψ h ( τ ) k I p ^ , ψ h ˜ ( b ) ] ,
and
z ^ ( t ) = k I p ^ , ψ h ˜ ( t ) + ( ψ ( t ) ψ ( a ) ) η k k 1 A Γ k ( η k ) [ A 1 ς k I α ^ , ψ h ( u ) + θ k I α ^ + w , ψ h ( τ ) k I p ^ , ψ h ˜ ( b ) + A 3 λ k I p ^ , ψ h ˜ ( ξ ) + μ k I p ^ + v , ψ h ˜ ( σ ) k I α ^ , ψ h ( b ) ] ,
where
A 1 = ( ψ ( b ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) , A 2 = λ ( ψ ( ξ ) ψ ( a ) ) η k k 1 Γ k ( η k ) + μ ( ψ ( σ ) ψ ( a ) ) η k + v k 1 Γ k ( η k + v ) , A 3 = ς ( ψ ( u ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) + θ ( ψ ( τ ) ψ ( a ) ) ϑ k + w k 1 Γ k ( ϑ k + w ) , A 4 = ( ψ ( b ) ψ ( a ) ) η k k 1 Γ k ( η k ) .
Proof. 
Let ( r ^ , z ^ ) be a solution of the system (5). As argued in [18], operating k I α ^ , ψ and k I p ^ , ψ on both sides of the first and second equations in (5), respectively, we obtain
r ^ ( t ) = k I α ^ , ψ h ( t ) + c 0 ( ψ ( b ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) + c 1 ( ψ ( t ) ψ ( a ) ) ϑ k k 2 Γ k ( ϑ k k ) , z ^ ( t ) = k I p ^ , ψ h ˜ ( t ) + d 0 ( ψ ( t ) ψ ( a ) ) η k k 1 Γ k ( η k ) + d 1 ( ψ ( t ) ψ ( a ) ) η k k 2 Γ k ( η k k ) ,
where c 0 , c 1 , d 0 and d 1 are unknown arbitrary constants. Using the conditions r ^ ( a ) = 0 and z ^ ( a ) = 0 in (9), we obtain c 1 = 0 and d 1 = 0 , since ϑ k k 2 < 0 , η k k 2 < 0 . Therefore, the equations in (9) take the form:
r ^ ( t ) = k I α ^ , ψ h ( t ) + c 0 ( ψ ( b ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) , z ^ ( t ) = k I p ^ , ψ h ˜ ( t ) + d 0 ( ψ ( t ) ψ ( a ) ) η k k 1 Γ k ( η k ) ,
which, on inserting in the boundary conditions: r ^ ( b ) = λ z ^ ( ξ ) + μ k I v , ψ z ^ ( σ ) and z ^ ( b ) = ς r ^ ( u ) + θ k I w , ψ r ^ ( τ ) , yields
k I α ^ , ψ h ( b ) + c 0 ( ψ ( b ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) = λ k I p ^ , ψ h ˜ ( ξ ) + λ d 0 ( ψ ( ξ ) ψ ( a ) ) η k k 1 Γ k ( η k ) + μ k I p ^ + v , ψ h ˜ ( σ ) + d 0 μ ( ψ ( σ ) ψ ( a ) ) η k + v k 1 Γ k ( η k + v ) ,
and
k I p ^ , ψ h ˜ ( b ) + d 0 ( ψ ( b ) ψ ( a ) ) ϑ k k 1 Γ k ( η k ) = ς k I α ^ , ψ h ( u ) + τ c 0 ( ψ ( u ) ψ ( a ) ) ϑ k k 1 Γ k ( ϑ k ) + θ k I α ^ + w , ψ h ( τ ) + θ c 0 ( ψ ( τ ) ψ ( a ) ) ϑ k + w k 1 Γ k ( ϑ k + w ) .
Solving the system (11) and (12) for c 0 and d 0 together with the notation in (4) and (8), we find that
c 0 = 1 A [ A 4 λ k I p ^ , ψ h ˜ ( ξ ) + μ k I p ^ + v , ψ h ˜ ( σ ) k I α ^ , ψ h ( b ) + A 2 ς k I α ^ , ψ h ( u ) + θ k I α ^ + w , ψ h ( τ ) k I p ^ , ψ h ˜ ( b ) ] , d 0 = 1 A [ A 1 ς k I α ^ , ψ h ( u ) + θ k I α ^ + w , ψ h ( τ ) k I p ^ , ψ h ˜ ( b ) + A 3 λ k I p ^ , ψ h ˜ ( ξ ) + μ k I p ^ + v , ψ h ˜ ( σ ) k I α ^ , ψ h ( b ) ] .
Replacing c 0 and d 0 by their above values in (10), we obtain (6) and (7). The converse is obtained by direct calculation and this completes the proof. □
In the following, we denote by X = C ( [ a , b ] , R ) the Banach space of all continuous real-valued functions on [ a , b ] equipped with the norm r ^ = max { | r ^ ( t ) | ; t [ a , b ] } . The space ( X × X , ( r ^ , z ^ ) ) is a Banach space endowed with norm ( r ^ , z ^ ) = r ^ + z ^ .

3. Existence and Uniqueness Results

In the following, for t [ a , b ] , we set f r ^ , z ^ ( t ) = f ( t , r ^ ( t ) , z ^ ( t ) ) and f ˜ r ^ , z ^ ( t ) = f ˜ ( t , r ^ ( t ) , z ^ ( t ) ) . By Lemma 1, we define an operator T : X × X X × X associated with the nonlinear system (1) as
T ( r ^ , z ^ ) ( t ) = T 1 ( r ^ , z ^ ) ( t ) T 2 ( r ^ , z ^ ) ( t ) ,
where
T 1 ( r ^ , z ^ ) ( t ) = k I α ^ , ψ f r ^ , z ˇ ( t ) + ( ψ ( t ) ψ ( a ) ) ϑ k k 1 A Γ k ( ϑ k ) [ A 4 λ k I p ^ , ψ f ˜ r ^ , z ^ ( ξ ) + μ k I p ^ + v , ψ f ˜ r ^ , z ^ ( σ ) k I α ^ , ψ f r ^ , z ^ ( b ) + A 2 ς k I α ^ , ψ f r ^ , z ^ ( u ) + θ k I α ^ + w , ψ f r ^ , z ^ ( τ ) k I p ^ , ψ f ˜ r , z ( b ) ] ,
and
T 2 ( r ^ , z ^ ) ( t ) = k I p ^ , ψ f ˜ r ^ , z ^ ( t ) + ( ψ ( t ) ψ ( a ) ) η k k 1 A Γ k ( η k ) [ A 1 ς k I α ^ , ψ f r ^ , z ^ ( u ) + θ k I α ^ + w , ψ h ( τ ) k I p ^ , ψ f ˜ r ^ , z ^ ( b ) + A 3 λ k I p ^ , ψ f ˜ r ^ , z ^ ( ξ ) + μ k I p ^ + v , ψ f ˜ r ^ , z ^ ( σ ) k I α ^ , ψ f r ^ , z ^ ( b ) ] .
Observe that the fixed-point problem T ( r ^ , z ^ ) = ( r ^ , z ^ ) is equivalent to the nonlinear system (1) and, hence, the existence of a fixed point for the operator T will imply the existence of a solution of the nonlinear system (1).
For simplicity, we set
Q 1 = ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α + w ^ k Γ k ( α ^ + w + k ) ) ] , Q 2 = ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) + ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) ] , Q 3 = ( ψ ( t ) ψ ( a ) ) η k k 1 | A | Γ k ( η k ) [ A 1 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k | A | Γ k ( α ^ + w + k ) ) + A 3 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) ] , Q 4 = ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + ( ψ ( b ) ψ ( a ) ) η k k 1 | A | Γ k ( ϑ k ) [ A 1 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + A 3 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) ] ,
and
Q 1 * = Q 1 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) , Q 4 * = Q 4 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) .

3.1. Existence Results

We prove our first existence result by applying the Leray–Schauder alternative [22].
Theorem 1.
Suppose that A 0 ( A is given by (4)) and f , f ˜ : [ a , b ] × R R are two continuous functions such that, for all t [ a , b ] and r ^ i , z ^ i R , i = 1 , 2 , we have
| f ( t , r ^ 1 , z ^ 1 ) | k ^ 0 + k ^ 1 | r ^ 1 | + k ^ 2 | z ^ 1 | , | f ˜ ( t , r ^ 2 , z ^ 2 ) | v ^ 0 + v ^ 1 | r ^ 2 | + v ^ 2 | z ^ 2 | ,
where k ^ i , v ^ i , i = 1 , 2 are real positive constants with k ^ 0 , v ^ 0 > 0 . Moreover, it is assumed that
( Q 1 + Q 3 ) k ^ 1 + ( Q 2 + Q 4 ) v ^ 1 < 1 a n d ( Q 1 + Q 3 ) k ^ 2 + ( Q 2 + Q 4 ) v ^ 2 < 1 ,
where Q i , i = 1 , 2 , 3 , 4 are given in (15). Then, the nonlinear system (1) has at least one solution on [ a , b ] .
Proof. 
Observe that continuity of the functions f and f ˜ ensures the continuity of the operator T . Next, we establish that the operator T is uniformly bounded.
Let M be a bounded set of X × X . Then, we can find positive constants L 1 and L 2 such that
| f ( t , r ^ ( t ) , z ^ ( t ) | L 1 , | f ˜ ( t , r ^ ( t ) , z ^ ( t ) | L 2 , ( r ^ , z ^ ) M .
Thus, for all ( r ^ , z ^ ) M , we have
| T 1 ( r ^ , z ^ ) ( t ) | { ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α + w ^ k Γ k ( α ^ + w + k ) ) } L 1 + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) ] L 2 .
Hence, we obtain
T 1 ( r ^ , z ^ ) Q 1 L 1 + Q 2 L 2 .
Analogously, we have
T 2 ( r ^ , z ^ ) Q 3 L 1 + Q 4 L 2 .
In view of the foregoing inequalities, it follows that
T ( r ^ , z ^ ) T 1 ( r ^ , z ^ ) + T 2 ( r ^ , z ^ ) ( Q 1 + Q 3 ) L 1 + ( Q 2 + Q 4 ) L 2 .
This shows that the operator T is uniformly bounded. To prove that the operator T is equicontinuous, we take t 1 , t 2 [ a , b ] with t 1 < t 2 . Then, we have
| T 1 ( r ^ ( t 2 ) , z ^ ( t 2 ) ) T 1 ( r ^ ( t 1 ) , z ^ ( t 1 ) ) | 1 Γ k ( α ^ ) | a t 1 ψ ( s ) [ ( ψ ( t 2 ) ψ ( s ) ) α k 1 ( ψ ( t 1 ) ψ ( s ) ) α ^ k 1 ] f ( s , r ^ ( s ) , z ^ ( s ) ) d s + t 1 t 2 ψ ( s ) ( ψ ( t 2 ) ψ ( s ) ) α ^ k 1 f ( s , r ^ ( s ) , z ^ ( s ) ) d s | + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( ϑ k ) [ A 4 ( | λ | k I p ^ , ψ | f ˜ ( ξ , r ^ ( ξ ) , z ^ ( ξ ) ) | + | μ | k I p ^ + v , ψ | f ˜ ( σ , r ^ ( σ ) , z ^ ( σ ) ) | + k I p ^ , ψ | f ( b , r ^ ( b ) , z ^ ( b ) ) | ) + A 2 ( | ς | k I α ^ , ψ | f ( u , r ^ ( u ) , z ^ ( u ) ) + | θ | k I α ^ + w , ψ | f ( τ , r ^ ( τ ) , z ^ ( τ ) ) | + k I p ^ , ψ | f ˜ ( b , r ^ ( b ) , z ^ ( b ) ) | ) ] L 1 Γ k ( α ^ + k ) [ 2 ( ψ ( t 2 ) ψ ( t 1 ) ) α ^ k + | ( ψ ( t 2 ) ψ ( c ) ) α ^ k ( ψ ( t 1 ) ψ ( c ) ) α ^ k | ] + ( ψ ( t 2 ) ψ ( a ) ) θ k k 1 ( ψ ( t 1 ) ψ ( a ) ) θ k k 1 | Λ | Γ k ( ϑ k ) [ A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) L 2 + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) L 2 + ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) L 1 ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) L 1 + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k Γ k ( α ^ + w + k ) L 1 + ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) L 2 ) ] .
The right hand of the above inequality is independent of ( r ^ , z ^ ) and tends to zero as t 2 t 1 0 . Hence, T 1 ( r ^ , z ^ ) is equicontinuous. Similarly, it can be shown that T 2 ( r ^ , z ^ ) is equicontinuous. Accordingly, the operator T ( r ^ , z ^ ) is completely continuous.
Finally, we show the boundedness of the set D = { ( r ^ , z ^ ) X × X : ( r ^ , z ^ ) = ν T ( r ^ , z ^ ) , 0 ν 1 } . Let ( r ^ , z ^ ) D . Then, ( r ^ , z ^ ) = ν T ( r ^ , z ^ ) and for all t [ a , b ] , we have
r ^ ( t ) = ν T 1 ( r ^ , z ^ ) ( t ) , z ( t ) = ν T 2 ( r ^ , z ^ ) ( t ) .
Thus, we obtain
| r ^ ( t ) | { ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α + w ^ k Γ k ( α ^ + w + k ) ) ] } ( k ^ 0 + k ^ 1 r ^ + k ^ 2 z ^ ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) { A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) } ( v ^ 0 + v ^ 1 r ^ + v ^ 2 z ^ ) ,
and
| z ^ ( t ) | ( ψ ( b ) ψ ( a ) ) η k k 1 | A | Γ k ( η k ) [ A 1 ( ( | ς | ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k | A | Γ k ( α ^ + w + k ) ) + A 3 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) ] ( k ^ 0 + k ^ 1 r ^ + k ^ 2 z ^ ) + { ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + ( ψ ( t ) ψ ( a ) ) η k k 1 | A | Γ k ( ϑ k ) [ A 1 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + A 3 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) ] } ( v ^ 0 + v ^ 1 r ^ + v ^ 2 z ^ ) .
Thus, we have
r ^ Q 1 ( k ^ 0 + k ^ 1 r ^ + k ^ 2 z ) + Q 2 ( ( v ^ 0 + v ^ 1 r ^ + v ^ 2 z )
and
z ^ Q 3 ( k ^ 0 + k 1 r ^ + k 2 z ) + Q 4 ( ( v ^ 0 + v ^ 1 r ^ + v ^ 2 z ) .
Hence, we obtain
r ^ + z ^ ( Q 1 + Q 3 ) k ^ 0 + ( Q 2 + Q 4 ) v ^ 0 + [ ( Q 1 + Q 3 ) k ^ 1 + ( Q 2 + Q 4 ) v ^ 1 ] r ^ + [ ( Q 1 + Q 3 ) k ^ 2 + ( Q 2 + Q 4 ) v ^ 2 ] z ^ ,
which implies that
( r ^ , z ^ ) ( Q 1 + Q 3 ) k ^ 0 + ( Q 2 + Q 4 ) v ^ 0 M 0 ,
where
M 0 = min { 1 [ ( Q 1 + Q 3 ) k ^ 1 + ( Q 2 + Q 4 ) v ^ 1 ] , 1 [ ( Q 1 + Q 3 ) k ^ 2 + ( Q 2 + Q 4 ) v ^ 2 ] } .
Consequently, using the Leray–Schauder alternative, it follows that the operator T has at least one fixed point, which is a solution of the nonlinear system (1) on [ a , b ] . □
In the following result, we apply Krasnosel’skiĭ’s fixed-point theorem [23] to prove our second existence result.
Theorem 2.
Assume that:
(1) There exist positive real constants m ^ i , n ^ i , i = 1 , 2 such that, f , f ˜ : [ a , b ] × R × R R are continuous and satisfy the following conditions for all t [ a , b ] and r ^ i , z ^ i R , i = 1 , 2 :
| f ( t , r ^ 1 , z ^ 1 ) f ( t , r ^ 2 , z ^ 2 ) | m ^ 1 | r ^ 1 r ^ 2 | + m ^ 2 | z ^ 1 z ^ 2 | ,
| f ˜ ( t , r ^ 1 , z ^ 1 ) f ˜ ( t , r ^ 2 , z ^ 2 ) | n ^ 1 | r ^ 1 r ^ 2 | + n ^ 2 | z ^ 1 z ^ 2 | .
(2) There exist P and Q C ( [ a , b ] , R + ) satisfying
| f ( t , r ^ , z ^ ) | P ( t ) , | f ˜ ( t , r ^ , z ^ ) | Q ( t ) , f o r e a c h ( t , r ^ , z ^ ) [ a , b ] × R × R .
Then, the nonlinear system (1) has at least one solution on [ a , b ] , provided that
[ Q 1 * + Q 3 ] ( m ^ 1 + m ^ 2 ) + [ Q 2 + Q 4 * ] ( n ^ 1 + n ^ 2 ) < 1 ,
where Q i , i = 1 , 2 , 3 , 4 are given in (15).
Proof. 
First, we decompose the operator T into four operators T 1 , 1 , T 1 , 2 , T 2 , 1 and T 2 , 2 as
T 1 , 1 ( r ^ , z ^ ) ( t ) = k I α ^ , ψ f ( t , r ^ ( t ) , z ^ ( t ) ) , t [ a , b ] , T 1 , 2 ( r ^ , z ^ ) ( t ) = ( ψ ( t ) ψ ( a ) ) ϑ k k 1 A Γ k ( ϑ k ) [ A 4 ( λ k I p ^ , ψ f ˜ ( ξ , r ^ ( ξ ) , z ^ ( ξ ) ) + μ k I p ^ + v , ψ f ˜ ( σ , r ^ ( σ ) , z ^ ( σ ) ) k I α ^ , ψ f ( b , r ^ ( b ) , z ^ ( b ) ) ) + A 2 ( ς k I α ^ , ψ f ( u , r ^ ( u ) , z ^ ( u ) ) ) + θ k I α ^ + w , ψ h ( τ , r ^ ( τ ) , z ^ ( τ ) k I p ^ , ψ f ˜ ( b , r ^ ( b ) , z ^ ( b ) ) ) ] , t [ a , b ] , T 2 , 1 ( r ^ , z ^ ) ( t ) = k I p ^ , ψ f ˜ ( t , r ^ ( t ) , z ^ ( t ) ) , t [ a , b ] , T 2 , 2 ( r ^ , z ^ ) ( t ) = ( ψ ( t ) ψ ( a ) ) η k k 1 A Γ k ( η k ) [ A 1 ( ς k I α ^ , ψ f ( u , r ^ ( u ) , z ^ ( u ) ) + θ k I α ^ + w , ψ f ( τ , r ^ ( τ ) , z ^ ( τ ) ) k I p ^ , ψ f ˜ ( b , r ^ ( b ) , z ^ ( b ) ) ) + A 3 ( λ k I p ^ , ψ f ˜ ( ξ , r ^ ( ξ ) , z ^ ( ξ ) ) + μ k I p ^ + v , ψ f ˜ ( σ , r ^ ( σ ) , z ^ ( σ ) ) k I α ^ , ψ f ( b , r ^ ( b ) , z ^ ( b ) ) ) ] , t [ a , b ] .
Notice that T 1 = T 1 , 1 + T 1 , 2 and T 2 = T 2 , 1 + T 2 , 2 . Assume that B ρ ^ = { ( r ^ , z ^ ) X × X : ( r ^ , z ^ ) ρ ^ } is a ball, in which ρ ^ ( Q 1 + Q 3 ) P + ( Q 2 + Q 4 ) Q . Similarly, as in Theorem 1, we have
| T 1 , 1 ( r ^ 1 , r ^ 2 ) ( t ) + T 1 , 2 ( z ^ 1 , z ^ 2 ) ( t ) | Q 1 P + Q 2 Q ,
and
| T 1 , 1 ( r ^ 1 , r ^ 2 ) ( t ) + T 2 , 2 ( z ^ 1 , z ^ 2 ) ( t ) | Q 3 P + Q 4 Q .
Consequently, we obtain
T 1 ( r ^ 1 , r ^ 2 ) + T 2 ( z ^ 1 , z ^ 2 ) ( Q 1 + Q 3 ) P + ( Q 2 + Q 4 ) Q < ρ ^ .
Thus, T 1 ( r ^ 1 , r ^ 2 ) + T 2 ( z ^ 1 , z ^ 2 ) B ρ ^ . Next, it will be shown that the operator ( T 1 , 2 , T 2 , 2 ) is a contraction mapping. For ( r ^ 1 , z ^ 1 ) , ( r ^ 2 , z ^ 2 ) B ρ ^ , we have
| T 1 , 2 ( r ^ 1 , r ^ 2 ) ( t ) T 1 , 2 ( z ^ 1 , z ^ 2 ) ( t ) | ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k Γ k ( α ^ + w + k ) ) ] ( m ^ 1 r ^ 1 z ^ 1 + m ^ 2 r ^ 2 z ^ 2 ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) ] ( n ^ 1 r ^ 1 z ^ 1 + n ^ 2 r ^ 2 z ^ 2 ) = Q 1 * ( m ^ 1 r ^ 1 z ^ 1 + m ^ 2 r ^ 2 z ^ 2 ) + Q 2 ( n ^ 1 r ^ 1 z ^ 1 + n ^ 2 r ^ 2 z ^ 2 ) = [ Q 1 * m ^ 1 + Q 2 n ^ 1 ] r ^ 1 z ^ 1 + [ Q 1 * m ^ 2 + Q 2 n ^ 2 ] r ^ 2 z ^ 2 .
In a similar fashion, one can find that
| T 2 , 2 ( r ^ 1 , r ^ 2 ) ( t ) T 2 , 2 ( z ^ 1 , z ^ 2 ) ( t ) | [ Q 3 m ^ 1 + Q 4 * n ^ 1 ] r ^ 1 z ^ 1 + [ Q 3 m ^ 2 + Q 4 * n ^ 2 ] r ^ 2 z ^ 2 .
From (18) and (19), we obtain
( T 1 , 2 , T 2 , 2 ) ( r ^ 1 , z ^ 1 ) ( T 1 , 2 , T 2 , 2 ) ( z ^ 1 , z ^ 2 ) [ Q 1 * + Q 3 ] ( m ^ 1 + m ^ 2 ) + [ Q 2 + Q 4 * ] ( n ^ 1 + n ^ 2 ) ( r ^ 1 z ^ 1 + r ^ 2 z ^ 2 ) .
which shows that the operator ( T 1 , 2 , T 2 , 1 ) is a contraction owing to the condition (17). On the other hand, in view of the continuity property of f and f ˜ , we conclude that the operator ( T 1 , 1 , T 2 , 1 ) is continuous. In addition, we have
T 1 , 1 ( r ^ , z ^ ) ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) P a n d T 2 , 1 ( r , z ) ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) Q .
Consequently, we obtain
( T 1 , 1 , T 2 , 1 ) ( r ^ , z ^ ) ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) P + ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) Q .
Hence, ( T 1 , 1 , T 2 , 1 ) B ρ is uniformly bounded. In the next step, it will be shown that the set ( T 1 , 1 , T 2 , 1 ) B ρ is equicontinuous. For t 1 , t 2 [ a , b ] with t 1 < t 2 and for all ( r ^ , z ^ ) B ρ ^ , we obtain
| T 1 , 1 ( r ^ , z ^ ) ( t 2 ) T 1 , 1 ( r ^ , z ^ ) ( t 1 ) 1 Γ k ( α ^ ) | a t 1 ψ ( s ) [ ( ψ ( t 2 ) ψ ( s ) ) α ^ k 1 ( ψ ( t 1 ) ψ ( s ) ) α k 1 ] f ( s , r ^ ( s ) , z ^ ( s ) ) d s + t 1 t 2 ψ ( s ) ( ψ ( t 2 ) ψ ( s ) ) α ^ k 1 f ( s , r ^ ( s ) , z ^ ( s ) ) d s | P Γ k ( α ^ + k ) [ 2 ( ψ ( t 2 ) ψ ( t 1 ) ) α ^ k + | ( ψ ( t 2 ) ψ ( a ) ) α ^ k ( ψ ( t 1 ) ψ ( a ) ) α ^ k | ] .
The right hand side of the above inequality is independent of ( r ^ , z ^ ) B ρ ^ and tends to zero as t 1 t 2 . Analogously, we have that | ( T 2 , 1 ( r ^ , z ^ ) ( t 2 ) T 2 , 1 ( r ^ , z ^ ) ( t 1 ) | 0 as t 1 t 2 , independent of ( r ^ , z ^ ) B ρ ^ .
Thus, | ( T 1 , 1 , T 2 , 1 ) ( r ^ , z ^ ) ( t 2 ) ( T 1 , 1 , T 2 , 1 ) ( r ^ , z ^ ) ( t 1 ) | 0 as t 1 t 2 . Therefore, ( T 1 , 1 , T 2 , 1 ) is equicontinuous. Now, applying the Arzelá–Ascoli theorem, we deduce that the operator ( T 1 , 1 , T 2 , 1 ) is compact on B ρ ^ . Thus, the hypotheses of Krasnosel’skiĭ’s fixed-point theorem are satisfied, and, hence, there exists at least one solution of the nonlinear system (1) on [ a , b ] .

3.2. Existence of a Unique Solution

In this subsection, we present a result concerning the uniqueness of solutions for the nonlinear system (1), which is derived with the aid of the Banach fixed-point theorem.
Theorem 3.
Suppose that A 0 ( A is given by (4)) and ( H 1 ) holds. Moreover, it is assumed that
( Q 1 + Q 3 ) ( m ^ 1 + m ^ 2 ) + ( Q 2 + Q 4 ) ( n ^ 1 + n ^ 2 ) < 1 ,
where Q i , i = 1 , 2 , 3 , 4 are defined in (15). Then, the nonlinear system (1) has a unique solution on [ a , b ] .
Proof. 
Letting sup t [ a , b ] | f ( t , 0 , 0 ) | = N < and sup t [ a , b ] | f ^ ( t , 0 , 0 ) | = N 1 , we fix
r ( Q 1 + Q 3 ) N + ( Q 2 + Q 4 ) N 1 1 [ ( Q 1 + Q 3 ) ( m ^ 1 + m ^ 2 ) + ( Q 2 + Q 4 ) ( n ^ 1 + n ^ 2 ) ] ,
and show that T B r B r , where B r = { ( r ^ , z ^ ) X × X : ( r ^ , z ^ ) r } . For ( r ^ , z ^ ) B r , we have
| T 1 ( r ^ , z ^ ) ( t ) | k I α ^ , ψ [ | f ( t , r ^ ( t ) , z ^ ( t ) ) f ( t , 0 , 0 ) | + | f ( t , 0 , 0 ) | ] + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( k I p ^ , ψ [ | f ^ ( ξ , r ^ ( ξ ) , z ^ ( ξ ) ) f ˜ ( ξ , 0 , 0 ) | + | f ˜ ( ξ , 0 , 0 ) | ] + | μ | k I p ^ + v , ψ [ | f ( σ , r ^ ( σ ) , z ^ ( σ ) ) f ( σ , 0 , 0 ) | + | f ( σ , 0 , 0 ) | ] + k I α ^ , ψ [ | f ( b , r ^ ( b ) , z ^ ( b ) ) f ( b , 0 , 0 ) | + | f ( b , 0 , 0 ) | ] ) + A 2 ( | ς | k I α ^ , ψ [ | f ( u , r ^ ( u ) , z ^ ( u ) f ( u , 0 , 0 ) | + | f ( u , 0 , 0 ) | ] + | θ | k I α ^ , ψ | f ( τ , r ^ ( τ ) , z ^ ( τ ) ) f ( τ , 0 , 0 ) | + | f ( τ , 0 , 0 ) | ] + k I p ^ , ψ [ [ f ˜ ( b , r ^ ( b ) , z ^ ( b ) ) f ˜ ( b , 0 , 0 ) | + | f ˜ ( b , 0 , 0 ) | ] ) ] k I α ^ , ψ [ m ^ 1 r ^ + m ^ 2 z ^ + N ] ( b ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | k I p ^ , ψ [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] ( ξ ) + | μ | k I p ^ + v , ψ [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] ( σ ) + k I p ^ , ψ [ m ^ 1 r ^ + m ^ 2 z ^ + N ] ( b ) + A 2 ( | ς | k I α ^ , ψ [ m ^ 1 r ^ + m ^ 2 z ^ + N ] ( u ) + | θ | k I α + w ^ , ψ [ m ^ 1 r ^ + m ^ 2 z ^ + N ] ( τ ) + k I p ^ , ψ [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] ( b ) ) ] ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] + ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) [ m ^ 1 r ^ + m ^ 2 z ^ + N ] ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k Γ k ( α ^ + w + k ) [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] ) ]
= { ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α + w ^ k Γ k ( α ^ + w + k ) ) ] } [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + { A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) } [ n ^ 1 r ^ + A n ^ 2 z ^ + N 1 ] = Q 1 [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + Q 2 [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] = ( Q 1 m ^ 1 + Q 2 n ^ 1 ) r ^ + ( Q 1 m ^ 2 + Q 2 n ^ 2 ) z ^ + Q 1 N + Q 2 N 1 Q 1 m ^ 1 + Q 2 n ^ 1 + Q 1 m ^ 2 + Q 2 n ^ 2 ) r + Q 1 N + Q 2 N 1 .
Analogously, we obtain
| T 2 ( r ^ , z ^ ) ( t ) | { ( ψ ( b ) ψ ( a ) ) η k k 1 | A | Γ k ( η k ) [ A 1 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α ^ + w k | A | Γ k ( α ^ + w + k ) ) + A 3 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) } [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + { ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + ( ψ ( b ) ψ ( a ) ) η k k 1 | A | Γ k ( ϑ k ) [ A 1 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + A 3 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) ] } [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] = Q 3 [ m ^ 1 r ^ + m ^ 2 z ^ + N ] + Q 4 [ n ^ 1 r ^ + n ^ 2 z ^ + N 1 ] = ( Q 3 m ^ 1 + Q 4 n ^ 1 ) r ^ + ( Q 3 m ^ 2 + Q 4 n ^ 2 ) z ^ + Q 3 N + Q 4 N 1 Q 3 m ^ 1 + Q 4 n ^ 1 + Q 3 m ^ 2 + Q 4 n ^ 2 ) r + Q 3 N + Q 4 N 1 .
Accordingly, we have
T ( r ^ , z ^ ) = T 1 ( r ^ , z ^ ) + T 2 ( r ^ , z ^ ) [ ( Q 1 + Q 3 ) ) ( m ^ 1 + m ^ 2 ) + ( Q 2 + Q 4 ) ( n ^ 1 + n ^ 2 ) ] r + ( Q 1 + Q 3 ) N + ( Q 2 + Q 4 ) N 1 r ,
which implies that T ( B r ) B r . On the other hand, for ( r ^ 2 , z ^ 2 ) , ( r ^ 1 , z ^ 1 ) X × X and t [ a , b ] , we have
| T 1 ( r 2 , z 2 ) ( t ) T 1 ( r 1 , z 1 ) ( t ) | k I α ^ , ψ | f ( t , r ^ 2 ( t ) , z ^ 2 ( t ) ) f ( t , r ^ 1 ( t ) , z ^ 1 ( t ) ) | + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( | λ | k I p ^ , ψ | f ˜ ( ξ , r ^ 2 ( ξ ) , z ^ 2 ( ξ ) ) f ˜ ( ξ , r ^ 1 ( ξ ) , z ^ 1 ( ξ ) ) | + | μ | k I p ^ + v , ψ | f ˜ ( σ , r ^ 2 ( σ ) , z ^ 2 ( σ ) ) f ˜ ( σ , r ^ 1 ( σ ) , z ^ 1 ( σ ) ) | + k I p ^ , ψ | f ( b , r ^ 2 ( b ) , z ^ 2 ( b ) ) f ( b , r ^ 1 ( b ) , z ^ 1 ( b ) ) | ) + A 2 ( | ς | k I α ^ , ψ | f ( u , r ^ 2 ( u ) , z ^ 2 ( u ) ) f ( u , r ^ 1 ( u ) , z ^ 1 ( u ) ) | + | θ | k I α ^ + w , ψ | f ( τ , r ^ 2 ( τ ) , z ^ 2 ( τ ) ) f ( τ , r ^ 1 ( τ ) , z ^ 1 ( τ ) ) | + k I p ^ , ψ | f ˜ ( b , r ^ 2 ( b ) , z ^ 2 ( b ) ) f ˜ ( b , r ^ 1 ( b ) , z ^ 1 ( b ) ) | ) ] { ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ A 4 ( ψ ( b ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + A 2 ( | ς | ( ψ ( u ) ψ ( a ) ) α ^ k Γ k ( α ^ + k ) + | θ | ( ψ ( τ ) ψ ( a ) ) α + w ^ k Γ k ( α ^ + w + k ) ) ] } ( m ^ 1 r ^ 2 r ^ 1 + m ^ 2 z ^ 2 z ^ 1 ) + ( ψ ( b ) ψ ( a ) ) ϑ k k 1 | A | Γ k ( ϑ k ) [ { A 4 ( | λ | ( ψ ( ξ ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) + | μ | ( ψ ( σ ) ψ ( a ) ) p ^ + v k Γ k ( p ^ + v + k ) ) + A 2 ( ψ ( b ) ψ ( a ) ) p ^ k Γ k ( p ^ + k ) } ] ( n ^ 1 r ^ 2 r ^ 1 + n ^ 2 z ^ 2 z ^ 1 ) = ( Q 1 m ^ 1 + Q 2 n ^ 1 ) ( r ^ 2 r ^ 1 ) + ( Q 1 m ^ 2 + Q 2 n ^ 2 ) ( z ^ 2 z ^ 1 ) .
Hence, we obtain
T 1 ( r ^ 2 , z ^ 2 ) ( t ) T 1 ( r ^ 1 , z ^ 1 ) ( ( Q 1 m ^ 1 + Q 2 n ^ 1 + Q 1 m ^ 2 + Q 2 n ^ 2 ) [ r ^ 2 r ^ 1 + z ^ 2 z ^ 1 ] .
Similarly, one can obtain
T 2 ( r ^ 2 , z ^ 2 ) ( t ) T 2 ( r ^ 1 , z ^ 1 ) ( ( Q 3 m ^ 1 + Q 4 n ^ 1 + Q 3 m ^ 2 + Q 4 n ^ 2 ) [ r ^ 2 r ^ 1 + z ^ 2 z ^ 1 ] .
From (22) and (23), we deduce that
T ( r ^ 2 , z ^ 2 ) ( t ) T ( r ^ 1 , z ^ 1 ) [ ( Q 1 + Q 3 ) ( m ^ 1 + m ^ 2 ) + ( Q 2 + Q 4 ) ( n ^ 1 + n ^ 2 ) ] ( r ^ 2 r ^ 1 + z ^ 2 z ^ 1 ) ,
which, by the condition (21), implies that T is a contraction. Thus, we deduce by the conclusion of Banach’s contraction mapping principle that T has a unique fixed point which corresponds to a unique solution of the nonlinear system (1). □

4. Examples

Consider the following six-point boundary value problems for a ( k , ϕ ) –Hilfer fractional system after fixing the coefficients in problem (1):
7 6 , H D 4 3 , 3 4 ; 2 + log t r ^ ( t ) = f ( t , r ^ ( t ) , z ^ ( t ) ) , 1 4 < t < 11 4 , 7 6 , H D 5 3 , 1 4 ; 2 + log t z ^ ( t ) = f ˜ ( t , r ^ ( t ) , z ^ ( t ) ) , 1 4 < t < 11 4 , r ^ 1 4 = 0 , r ^ 11 4 = 1 π 2 z ^ 3 4 + 7 39 7 6 I 5 3 ; 2 + log t z ^ 7 4 , z ^ 1 4 = 0 , z ^ 11 4 = 9 37 r ^ 5 4 + 1 e 2 7 6 I 3 5 ; 2 + log t r ^ 9 4 .
Here k = 7 / 6 , α ^ = 4 / 3 , p ^ = 5 / 3 , β ^ = 3 / 4 , q ^ = 1 / 4 , v = 5 / 3 , w = 3 / 5 , ψ ( t ) = 2 + log t , λ = 1 / π 2 , μ = 7 / 39 , ς = 9 / 37 , θ = 1 / e 2 , a = 1 / 4 , b = 11 / 4 , ξ = 3 / 4 , σ = 7 / 4 , u = 5 / 4 , τ = 9 / 4 . Using the given data, we find that ϑ k = 25 / 12 , η k = 11 / 6 , Γ k ( ϑ k ) 1.047114902 , Γ k ( η k ) 0.9726277286 , Γ k ( ϑ k + w ) 1.425589394 , Γ k ( η k + v ) 2.722222222 , Γ k ( α ^ + k ) 1.275021091 , Γ k ( p ^ + k ) 1.577652779 , Γ k ( α ^ + w + k ) 1.928101260 , Γ k ( p ^ + v + k ) 7.812151797 , A 1 1.898644621 , A 2 0.3595879078 , A 3 0.6017756924 , A 4 1.694724814 ( A 1 , A 2 , A 3 , and A 4 are given in (8)), A 3.001288890 ( A is given by (4)), Q 1 4.549161327 , Q 2 0.7467564627 , Q 3 1.353804071 , Q 4 4.484708590 ( Q i , i = 1 , 2 , 3 , 4 are defined in (15)), Q 1 * 2.418206215 , Q 4 * 2.273635184 ( Q 1 * , Q 4 * are given in (16)).
(a) Illustration of Theorem 1.
Let the functions f, f ˜ : [ 1 / 4 , 11 / 4 ] × R × R R be defined by
f ( t , r ^ , z ^ ) = 1 + cos 2 ( r ^ z ^ ) 2 π t + e z ^ 2 sin r ^ 4 t + 11 + | z ^ | 15 e ( 4 t 1 ) 2 13 ( 1 + z ^ 14 ) ,
f ˜ ( t , r ^ , z ^ ) = 4 t ( r ^ z ^ ) 2 1 + ( r ^ z ^ ) 2 + r ^ 12 11 ( 1 + | r ^ | 11 ) + z ^ sin 8 ( r ^ ) 4 t + 9 .
It is easy to find that | f ( t , r ^ , z ^ ) | ( 4 / π ) + ( 1 / 12 ) | r ^ | + ( 1 / 13 ) | z ^ | and | f ˜ ( t , r ^ , z ^ ) | 11 + ( 1 / 11 ) | r ^ | + ( 1 / 10 ) | z ^ | . Setting k ^ 0 = 4 / π , k ^ 1 = 1 / 12 , k ^ 2 = 1 / 13 , v ^ 0 = 11 , v ^ 1 = 1 / 11 , v ^ 2 = 1 / 10 , we obtain ( Q 1 + Q 3 ) k ^ 1 + ( Q 2 + Q 4 ) v ^ 1 0.9675015153 < 1 and ( Q 1 + Q 3 ) k ^ 2 + ( Q 2 + Q 4 ) v ^ 2 0.9772207667 < 1 . Therefore, we deduce by Theorem 1 that problem (25) with functions f, f ˜ given by (26) and (27), respectively, has at least one solution on [ 1 / 4 , 11 / 4 ] .
(b) Illustration of Theorem 2.
Define two nonlinear Lipschitzian functions f, f ˜ : [ 1 / 4 , 11 / 4 ] × R × R R by
f ( t , r ^ , z ^ ) = 1 3 sin 2 π t + | r ^ | e ( 4 t 1 ) 16 ( 1 + | r ^ | ) + 1 14 tan 1 z ^ ,
f ˜ ( t , r ^ , z ^ ) = 1 5 cos 2 π t + 1 15 sin r ^ + | z ^ | 52 t ( 1 + | z ^ | ) .
Obviously, we have
| f ( t , r ^ , z ^ ) | 1 3 sin 2 π t + 1 16 e ( 4 t 1 ) + π 28 , | f ˜ ( t , r ^ , z ^ ) | 1 5 cos 2 π t + 1 52 t + 1 15 ,
and
| f ( t , r ^ 1 , z ^ 1 ) f ( t , r ^ 2 , z ^ 2 ) | 1 16 | r ^ 1 r ^ 2 | + 1 14 | z ^ 1 z ^ 2 | ,
| f ˜ ( t , r ^ 1 , z ^ 1 ) f ˜ ( t , r ^ 2 , z ^ 2 ) | 1 15 | r ^ 1 r ^ 2 | + 1 13 | z ^ 1 z ^ 2 | .
Setting m ^ 1 = 1 / 16 , m ^ 2 = 1 / 14 , n ^ 1 = 1 / 15 , n ^ 2 = 1 / 13 , we find that [ Q 1 * + Q 3 ] ( m ^ 1 + m ^ 2 ) + [ Q 2 + Q 4 * ] ( n ^ 1 + n ^ 2 ) 0.9388772111 < 1 . Therefore, all the assumptions of Theorem 2 are fulfilled and, hence, the ( k , ϕ ) –Hilfer fractional system (25) with functions f, f ˜ given by (28) and (29), respectively, has at least one solution on [ 1 / 4 , 11 / 4 ] .
(c) Illustration of Theorem 3.
Let two nonlinear Lipschitzian unbounded functions f, f ˜ : [ 1 / 4 , 11 / 4 ] × R × R R be given by
f ( t , r ^ , z ^ ) = e ( 4 t 1 ) 2 8 ( 3 t + 5 ) r ^ 2 + 2 | r ^ | 1 + | r ^ | + sin | z ^ | 4 ( 4 t + 5 ) + 1 3 ,
f ˜ ( t , r ^ , z ^ ) = sin 2 t tan 1 ( r ^ ) 4 ( 2 t + 5 ) + cos 2 ( π t ) 28 ( 2 t + 1 ) z ^ 2 + 2 | z ^ | 1 + | z ^ | + 1 4 .
Clearly
| f ( t , r ^ 1 , z ^ 1 ) f ( t , r ^ 2 , z ^ 2 ) | 1 23 | r ^ 1 r ^ 2 | + 1 24 | z ^ 1 z ^ 2 | ,
and
| f ˜ ( t , r ^ 1 , z ^ 1 ) f ˜ ( t , r ^ 2 , z ^ 2 ) | 1 22 | r ^ 1 r ^ 2 | + 1 21 | z ^ 1 z ^ 2 | .
Setting m ^ 1 = 1 / 23 , m ^ 2 = 1 / 24 , n ^ 1 = 1 / 22 and n ^ 2 = 1 / 21 , we have that ( Q 1 + Q 3 ) ( m ^ 1 + m ^ 2 ) + ( Q 2 + Q 4 ) ( n ^ 1 + n ^ 2 ) 0.9895188106 < 1 . Thus, the conclusion of Theorem 3 applies and, hence, the problem (25) with functions f, f ˜ given by (30) and (31), respectively, has a unique solution on the interval [ 1 / 4 , 11 / 4 ] .

5. Conclusions

In this paper, we presented sufficient criteria for the existence and uniqueness of solutions for a coupled system of nonlinear ( k , ψ ) –Hilfer fractional differential equations complemented with coupled ( k , ψ ) –Riemann–Liouville fractional integral boundary conditions. We proved the desired results for the given problem by applying the Leray–Schauder alternative, Krasnosel’skiĭ’s fixed-point theorem and Banach’s contraction mapping principle. As a special case, we obtained the new results for a coupled system of nonlinear ( k , ψ ) –Hilfer fractional differential equations equipped with
(i)
four-point nonlocal coupled boundary conditions:
r ^ ( a ) = 0 , r ^ ( b ) = λ z ^ ( ξ ) , z ^ ( a ) = 0 , z ^ ( b ) = ς r ^ ( u ) ,
by taking μ = 0 and θ = 0 in the results of the paper; and
(ii)
purely coupled ( k , ψ ) –Riemann–Liouville fractional integral boundary conditions:
r ^ ( a ) = 0 , r ^ ( b ) = μ k I v , ψ z ^ ( σ ) , z ^ ( a ) = 0 , z ^ ( b ) = θ k I w , ψ r ^ ( τ ) ,
by fixing λ = 0 , ς = 0 in the present results.
Comparing the present results with the related work on the topic, it is worthwhile to mention that the authors studied the ( k , ψ ) –Hilfer fractional differential equations and inclusions equipped with the different boundary conditions in the articles [17,18,19], while a system of ( k , φ ) –Hilfer fractional differential equations with nonlocal boundary conditions was investigated in [20]. Thus, the work established in this paper is not only new in the given configuration but also produces some new results as special cases. Therefore, our work is a significant contribution to the existing literature on the boundary value problems of ( k , ψ ) –Hilfer fractional differential systems. In future, we plan to investigate coupled systems of nonlinear sequential ( k , ψ ) –Hilfer fractional differential equations and inclusions equipped with different kinds of boundary conditions.

Author Contributions

Conceptualization, S.K.N.; methodology, A.S., S.K.N., B.A. and J.T.; validation, A.S., S.K.N., B.A. and J.T.; formal analysis, A.S., S.K.N., B.A. and J.T. writing—original draft preparation, A.S., S.K.N., B.A. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the reviewers for their constructive remarks on their work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of the Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  2. Ahmad, B.; Ntouyas, S.K. Nonlocal Nonlinear Fractional-Order Boundary Value Problems; World Scientific: Singapore, 2021. [Google Scholar]
  3. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  4. Nuchpong, C.; Ntouyas, S.K.; Tariboon, J. Boundary value problems of Hilfer-type fractional integro-differential equations and inclusions with nonlocal integro-multipoint boundary conditions. Open Math. 2020, 18, 1879–1894. [Google Scholar] [CrossRef]
  5. Subramanian, M.; Gopal, T.N. Analysis of boundary value problem with multi-point conditions involving Caputo-Hadamard fractional derivative. Proyecciones 2020, 39, 155–1575. [Google Scholar] [CrossRef]
  6. Belbali, H.; Benbachir, M.; Etemad, S.; Park, C.; Rezapour, S. Existence theory and generalized Mittag-Leffler stability for a nonlinear Caputo-Hadamard FIVP via the Lyapunov method. AIMS Math. 2022, 7, 14419–14433. [Google Scholar] [CrossRef]
  7. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Danbury, CT, USA, 2006. [Google Scholar]
  8. Zaslavsky, G.M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2005. [Google Scholar]
  9. Fallahgoul, H.A.; Focardi, S.M.; Fabozzi, F.J. Fractional Calculus and Fractional Processes with Applications to Financial Economics. Theory and Application; Elsevier/Academic Press: London, UK, 2017. [Google Scholar]
  10. Almalahi, M.A.; Abdo, M.S.; Panchal, S.K. Existence and Ulam-Hyers stability results of a coupled system of ψ–Hilfer sequential fractional differential equations. Results Appl. Math. 2021, 10, 100142. [Google Scholar] [CrossRef]
  11. Wongcharoen, A.; Ntouyas, S.K.; Wongsantisuk, P.; Tariboon, J. Existence results for a nonlocal coupled system of sequential fractional differential equations involving ψ–Hilfer fractional derivatives. Adv. Math. Phys. 2021, 2021, 5554619. [Google Scholar] [CrossRef]
  12. Nuchpong, C.; Ntouyas, S.K.; Samadi, A.; Tariboon, J. Boundary value problems for Hilfer type sequential fractional differential equations and inclusions involving Riemann-Stieltjes integral multi-strip boundary conditions. Adv. Differ. Equ. 2021, 2021, 268. [Google Scholar] [CrossRef]
  13. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  14. Dorrego, G.A. An alternative definition for the k–Riemann-Liouville fractional derivative. Appl. Math. Sci. 2015, 9, 481–491. [Google Scholar] [CrossRef] [Green Version]
  15. Vanterler da C. Sousa, J.; Capelas de Oliveira, E. On the ψ–Hilfer fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 60, 72–91. [Google Scholar] [CrossRef]
  16. Kucche, K.D.; Mali, A.D. On the nonlinear (k,ψ)-Hilfer fractional differential equations. Chaos Solitons Fractals 2021, 152, 111335. [Google Scholar] [CrossRef]
  17. Ntouyas, S.K.; Ahmad, B.; Tariboon, J.; Alhodaly, M.S. Nonlocal integro-multi-point (k,ψ)-Hilfer type fractional boundary value problems. Mathematics 2022, 10, 2357. [Google Scholar] [CrossRef]
  18. Ntouyas, S.K.; Ahmad, B.; Tariboon, J. Nonlocal boundary value problems for (k,ψ)-Hilfer fractional differential equations and inclusions. Foundations 2022, 2, 681–696. [Google Scholar] [CrossRef]
  19. Ntouyas, S.K.; Ahmad, B.; Nuchpong, C.; Tariboon, J. On (k,ψ)-Hilfer fractional differential equations and inclusions with mixed (k,ψ)-derivative and integral boundary conditions. Axioms 2022, 11, 403. [Google Scholar] [CrossRef]
  20. Samadi, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal coupled system for (k,φ)-Hilfer fractional differential equations. Fractal Fract. 2022, 6, 234. [Google Scholar] [CrossRef]
  21. Kwun, Y.C.; Farid, G.; Nazeer, W.; Ullah, S.; Kang, S.M. Generalized Riemann-Liouville k-fractional integrals associated with Ostrowski type inequalities and error bounds of Hadamard inequalities. IEEE Access 2018, 6, 64946–64953. [Google Scholar] [CrossRef]
  22. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003. [Google Scholar]
  23. Krasnosel’skiĭ, M.A. Two remarks on the method of successive approximations. Uspekhi 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

Samadi, A.; Ntouyas, S.K.; Ahmad, B.; Tariboon, J. Investigation of a Nonlinear Coupled (k, ψ)–Hilfer Fractional Differential System with Coupled (k, ψ)–Riemann–Liouville Fractional Integral Boundary Conditions. Foundations 2022, 2, 918-933. https://doi.org/10.3390/foundations2040063

AMA Style

Samadi A, Ntouyas SK, Ahmad B, Tariboon J. Investigation of a Nonlinear Coupled (k, ψ)–Hilfer Fractional Differential System with Coupled (k, ψ)–Riemann–Liouville Fractional Integral Boundary Conditions. Foundations. 2022; 2(4):918-933. https://doi.org/10.3390/foundations2040063

Chicago/Turabian Style

Samadi, Ayub, Sotiris K. Ntouyas, Bashir Ahmad, and Jessada Tariboon. 2022. "Investigation of a Nonlinear Coupled (k, ψ)–Hilfer Fractional Differential System with Coupled (k, ψ)–Riemann–Liouville Fractional Integral Boundary Conditions" Foundations 2, no. 4: 918-933. https://doi.org/10.3390/foundations2040063

Article Metrics

Back to TopTop