Next Article in Journal
Symmetric Bernstein Polynomial Approach for the System of Volterra Integral Equations on Arbitrary Interval and Its Convergence Analysis
Next Article in Special Issue
Analysis of Tempered Fractional Calculus in Hölder and Orlicz Spaces
Previous Article in Journal
Murine Motion Behavior Recognition Based on DeepLabCut and Convolutional Long Short-Term Memory Network
Previous Article in Special Issue
On Fractional Hybrid Non-Linear Differential Equations Involving Three Mixed Fractional Orders with Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations

1
Department of Mathematics & Statistics, University of Swat, Mingora 19130, Pakistan
2
Department of Mathematics, Texas A & M University Kingsville, Kingsville, TX 78363, USA
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1342; https://doi.org/10.3390/sym14071342
Submission received: 20 May 2022 / Revised: 17 June 2022 / Accepted: 21 June 2022 / Published: 29 June 2022

Abstract

:
In the recent past, multi-term fractional equations have been studied using symmetry methods. In some cases, many practical test problems with some symmetries are provided to demonstrate the authenticity and utility of the used techniques. Fractional-order differential equations can be formulated by using two types of differential operators: single-term and multi-term differential operators. Boundary value problems with single- as well as multi-term differential operators have been extensively studied, but several multi-term fractional differential equations still need to be formulated, and examination should be done with symmetry or any other feasible techniques. Therefore, the purpose of the present research work is the formulation and study of a new couple system of multi-term fractional differential equations with delay, as well as supplementation with nonlocal boundary conditions. After model formulation, the existence of a solution and the uniqueness conditions will be developed, utilizing fixed point theory and functional analysis. Moreover, results related to Ulam’s and other types of functional stability will be explored, and an example is carried out to illustrate the findings of the work.

1. Introduction

The study of fractional-order mathematical models, which govern real-world problems appearing in various disciplines of engineering, has been extensively improved by the methods of fractional calculus. The main reason for the importance of the selection of fractional-order differential and integral operators is the nonlocal behavior of these operators. Theoretical and applied aspects of fractional-order differential equations have been extensively explored. Agarwal et al. in [1] used Pettis integral and nonreflexive Banach spaces to establish results related to the existence of Abel’s integral equation and Cauchy type problem. A new aspect known as practical stability has been studied by the authors in [2]. This new notion of stability enables us to study the behavior of a couple of solutions when both initial values and intervals are different. Positive solutions for a class of boundary value problem, equipped with p-Laplacian operator and multipoint, are studied with the help of a fixed-point theorem of mixed-type monotone operators [3]. To see the effect of initial conditions on the fractional derivative, a model is formulated in [4], which also discusses the delay differential equation for the existence of a solution.
Besides the theoretical aspects of fractional differential equations, some applications of this emerging field can also be found in economics, ecology, immune systems, pathology, and chaos theory [5,6,7,8,9,10]. In this connection, various researchers have been inspired by the great popularity of fractional differential equations (FDEs) and used the theory to some extent. Numerical investigation of Whitham–Broer–Kaup equations was carried out by Ali et al. in [11], utilizing Laplace transform together with Adomian decomposition method. For more exploration of the numerical solution of partial and ordinary FDEs, one can refer to [12,13,14,15]. More information regarding the theoretical and practical development of FDE can be found in [16,17,18,19,20,21] and the references therein.
Another important aspect of FDE is the investigation of the well posedness of a model formulated via fractional differential operator. In this connection, the existence and uniqueness of a solution studied by various techniques, such as monotone iterative technique, were discussed in [22]; topological degree theory was the concern in [23]; a method of successive approximation was developed in [24]; fixed point index theory was used in [25], and tools of fixed point theory (FPT) were used in [26]. For more details, see the recent work cited as [27,28,29,30].
Initial and boundary values problems (BVP) of ordinary and fractional-order differential equations have many applications in the modeling of biological problems and dynamical systems. Nevertheless, many physical, chemical, and other processes depend not only on boundary points but on the interior of the domain. Therefore, to analyze the underlying problem, we need a novel boundary and initial conditions called the nonlocal conditions. These special types of conditions relate initial and boundary values of the solution to their values in the interior points of the domain [31]. For further details, [32,33,34] can enrich the reader’s understanding regarding nonlocal conditions.
The first FDE supplemented the classical boundary condition with two points and was studied by Benchohra et al. in [35] using results from nonlinear analysis and FPT. The authors first converted the underlying problem to the corresponding integral equation and then used Schaefer’s fixed point theorem for the existence of a solution. Furthermore, the solution uniqueness result was also developed by utilizing the Banach contraction principle. To be precise, the model under consideration is given as follows:
c D α 1 V ( t ) = E ( t , V ( t ) ) , 0 t T , a 0 V ( 0 ) + a 1 V ( 1 ) = c ,
where a 0 , a 1 , c , and T are arbitrary constants, with the given condition that the sum of a 0 and a 1 must be a non-zero real number together with T > 0 . The fractional differential operator used in the above-mentioned model is the Caputo operator, which is symbolized as c D α 1 , where α 1 is any real number from the semi-open interval (0,1]. Furthermore, the operator E is a given continuous function in the corresponding domain. The authors also stated that the problem under consideration takes the form of an initial value problem for a 0 = 1 and a 1 = 0 , while the terminal value problem can be obtained by fixing a 0 = 0 and a 1 = 1 . . An appropriate example was also given to justify the obtained results.
The imposing nonlocal conditions on the unknown function of second-order FDE were investigated in [36]. The author derived the desired result by using Krasnoselskii’s fixed point theorem and Banach contraction to demonstrate the existence of a solution to the following model:
c D α 1 V ( t ) = H ( t , V ( t ) ) , 0 t 1 , a 0 V ( 0 ) a 1 V ( 0 ) = a 2 V ( σ 1 ) , b 0 V ( 1 ) + b 1 V ( 1 ) = b 2 V ( σ 2 ) .
Here, the nonlinear function H : [ 0 , 1 ] × R R is a given continuous function, while a i , b i ( i = 1 , 2 , 3 ) are taken from the set of all real numbers. The points σ 1 and σ 2 are taken from the interior of the domain, and the considered model uses a Caputo fractional derivative of order α 1 , which is mathematically symbolized as c D α 1 , where α 1 ( 0 , 2 ] . Moreover, the arbitrary constant should not satisfy the following condition: ( a 2 + a 3 σ 1 ) ( b 1 b 3 ) = ( b 1 + b 2 b 3 σ 2 ) ( a 1 a 3 ) .
Another example is related to geophysical morphodynamics. This type of investigation is associated with the dynamical aspect of coastlines and sandbanks [37]. Based on these applications, researchers recently studied mathematical aspects and different types of stabilities for a coupled system of fractional-order differential equations.
In real-world problems, there exist several phenomena that can be modeled and studied by a coupled system of differential equations, rather than other types of differential equations. Among various other applications, binary fluid convection is an example of fractional calculus in physics [38]. In binary fluid convection, two ideal infinitely long plates are considered as a channel for fluid flow.
In this regard, Ahmad et al. [39] investigated the following coupled system of pantograph differential equations:
c D η V 1 ( t ) = F ( t , V 1 ( t ) , V 1 ( λ t ) , V 2 ( t ) ) , c D β V 2 ( t ) = K ( t , V 2 ( t ) , V 2 ( λ t ) , V 1 ( t ) ) , η , β ( 0 , 1 ] , 0 < λ < 1 , t [ 0 , 1 ] , h 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , h 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) .
Here, the domain of the unknown function is unit interval [0,1], while the nonlinear functions F , K : [ 0 , 1 ] × ( , ) [ 3 ] ( , ) are continuous. The given continuous function h 1 , h 2 : C ( [ 0 , 1 ] , R ) ( , ) , which is involved in the auxiliary conditions, makes the considered problem non-local, and the arbitrary constants a j , b j ( j = 1 , 2 , 3 ) will be chosen from the set of real numbers. The aforementioned system of FDEs utilized a Caputo derivative of order η and β , symbolized as c D η and c D β , respectively. The arbitrary constants from the set of real numbers can be chosen such that a j b j c j 0 for j=1,2. The solution’s existence and uniqueness were obtained by using FPT and results from non-linear analysis. Furthermore, after deriving the integral form of the solution, the considered problem Equation (1) was studied with the help of an operator. Banach contraction principle and Shafer’s fixed point theorems were considered in developing the desired results. Moreover, the results for Ulam–Hyres (UH) and generalized Ulam–Hyres (GUH) stability together with appropriate examples were also considered.
The study of differential equations equipped with multi-term differential operators of fractional orders has been given very little attention. Exploration of the functional stability of such kinds of multi-term differential operators has been given negligible attention. Therefore, motivated by the existence results discussed above and applications of multi-term FED, we formulate a new model here. This novel model consists of a coupled system of FDEs containing more than one fractional derivative of first order supplemented with non-local boundary conditions. The novel aspects of our newly formulated model are as follows: To the best of our knowledge no contribution exists regarding the formulation of coupled systems of FDEs, equipped with multi-term fractional differential operators. Therefore, the formulation of a coupled system with n-fractional differential multi-term differential operators is the first novel aspect of our work. As mentioned, less attention is given to the functional stability of multi-term FDEs, so our analysis of the four types of UH stability is our second novel aspect. The final novel aspect is the incorporation of nonlocal conditions and proportional type delay. The model under consideration is given as follow:
i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) , λ 1 , λ 2 , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) , α i , β i ( 0 , 1 ] , for i = 1 , 2 , 3 , , n g 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , g 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) .
In the considered problem Equation (2), standard Caputo fractional-order derivatives of order α i , β i ( 0 , 1 ] , symbolized as c D α i , c D β i ( i = 1 , 2 , n ) , are utilized with the conditions that α 1 > α i and β 1 > β i for i = 2, 3, ..., n. The auxiliary conditions imposed on the given problem are nonlocal due to the points ζ 1 and ζ 2 being chosen from the interior of the domain. The arbitrary constants η i , σ i , a j , b j , c j , ( j = 1 , 2 ) , ( i = 1 , 2 , , n ) can be selected from the set of real numbers, with conditions that η 1 0 , σ 1 0 , a 1 b 1 + c 1 , and a 2 b 2 + c 2 . The non-linear functions f 1 , f 2 : [ 0 , 1 ] × R × R × R R and g 1 , g 2 : C ( [ 0 , 1 ] , R ) × R R are continuous, where 1 R is any real number greater than zero.
The scope of work for studying the model represented in Equation (2) is given as follows: The tools of FPT and functional analysis will be used to obtain the conditions for solution existence and uniqueness. In this regard, the problem under consideration will be converted into an integral equation and then to a fixed point problem. To achieve our goals, we will impose certain assumptions on the given functions and arbitrary constants. Furthermore, Banach’s contraction principle and Krasnoselskii’s fixed point theorem will be used in the investigation of uniqueness and at least one solution, respectively. Moreover, the solution of the problem under consideration has interesting behavior under certain assumptions imposed on the given data. In addition to the above, results regarding Ulam’s type stability, such as UH, GUH, UH Rassias, and GUH Rassias, will be derived. For verification of the obtained results, we give an appropriate example.
The structure of the rest of the paper is as follows. Basic supporting results will be given in Section 2. The results that rely on Krasnoselskii’s fixed point theorem and Banach contraction principle will be provided in Section 3. Stability-related results are explored in Section 4. An elaborating example is provided at the end of the obtained results.

2. Preliminaries

This section of research is devoted to basic results and definitions of FPT and fractional calculus, which will be needed for investigation of the main work.
Definition 1
([17]). “The integral of fractional-order α of a function V ( t ) L [ 0 , d ] is denoted by I α V ( t ) and defined as
I α V ( t ) = 0 t V ( χ ) Γ ( α ) ( t χ ) 1 α d χ ,
provided that the right-hand side is pointwise convergent“.
Definition 2
([17]). “The fractional-order Caputo derivative for a function V ( t ) L 1 ( [ 0 , d ] , R + ) on the interval [ 0 , d ] is defined as
c D α V ( t ) = 0 t V n ( χ ) Γ ( n α ) ( t χ ) α + 1 n d χ ,
provided that the right-hand side is pointwise convergent. Where n = α , and α is defined to be the smallest integer equal or greater than α”.
Theorem 1
([17]). “The solution of FDE
c D α V ( t ) = 0 , w h e r e n 1 < α n ,
is given by
V ( t ) = A 1 + A 2 t + A 3 t 2 + A 4 t 3 + . . . . . . . . . . A n t n 1 ,
where A i R for i = 1 , 2 , , n ”.
Lemma 1
([17]). “The relation between fractional-order integral and derivative is given as
I α [ c D α V ( t ) ] = A 1 + A 2 t + A 3 t 2 + A 4 t 3 + . . . . . . . . . . A n t n 1 + V ( t ) ,
where A i R for i = 1 , 2 , , n ”.
Theorem 2
([40]). “ Assume that H is a non-empty, convex, bounded, and closed subset of a Banach space X . Let J 1 and J 1 be two operators, provided that J 1 U 1 + J 2 U 2 H whenever U , U 2 H , J 1 is continuous and compact, and J 2 is a contraction map. Then, we get U H , provided that U = J 1 U + J 2 U ”.

3. Main Result

This section is devoted to the main results associated with our newly formulated model.
Theorem 3.
For N 1 , N 2 C ( J , R ) , the solution of system of linear MFDDEs,
i = 1 n σ i c D α i V 1 ( t ) = N 1 ( t ) , i = 1 n η i c D β i V 2 ( t ) = N 2 ( t ) , α i , β i ( 0 , 1 ] , for i = 1 , 2 , 3 , , n , t [ 0 , 1 ] , g 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , g 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) ,
is given by
V 1 ( t ) = d 1 g 1 ( V 1 ) d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 N 1 ( U ) d U + i = 2 n d 1 b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 N 1 ( U ) d U + i = 2 n d 1 c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 N 1 ( U ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U , V 2 ( t ) = d 2 g 2 ( V 2 ) d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 N 2 ( U ) d U + i = 2 n d 2 b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 N 2 ( U ) d U + i = 2 n d 2 c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 N 2 ( U ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U ,
where d 1 = 1 a 1 + b 1 + c 1 , d 2 = 1 a 2 + b 2 + c 2 , a 1 b 1 c 1 0 , and a 2 b 2 c 2 0 .
Proof. 
Applying fractional-order integral I α 1 , I β 1 on Equation (3), and in view of Lemma (1),
V 1 ( t ) = A 1 + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 N 1 ( U ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U , V 2 ( t ) = B 1 + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 N 2 ( U ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U .
Now, by using the subsidiary condition Equation (3) in Equation (5),
A 1 = g 1 ( V 1 ) 1 a 1 + b 1 + c 1 1 a 1 + b 1 + c 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 N 1 ( U ) d U + 1 a 1 + b 1 + c 1 i = 2 n b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U 1 a 1 + b 1 + c 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 N 1 ( U ) d U + 1 a 1 + b 1 + c 1 i = 2 n c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U , B 1 = g 2 ( V 2 ) 1 a 2 + b 2 + c 2 1 a 2 + b 2 + c 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 N 2 ( U ) d U + 1 a 2 + b 2 + c 2 i = 2 n b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U 1 a 2 + b 2 + c 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 N 2 ( U ) d U + 1 a 2 + b 2 + c 2 i = 2 n c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U .
One can obtain the desired solution Equation (4) by using Equation (6) in (5) and assuming that d 1 = 1 a 1 + b 1 + c 1 0 , d 2 = 1 a 2 + b 2 + c 2 0 . □
Corollary 1.
In view of Theorem (3), the solution of the given system of MFDDEs
i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) , λ 1 , λ 2 [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) , α i , β i ( 0 , 1 ] , for i = 1 , 2 , 3 , , n , t [ 0 , 1 ] , g 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , g 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) ,
is given by
V 1 ( t ) = d 1 g 1 ( V 1 ) d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U ,
V 2 ( t ) = d 2 g 2 ( V 2 ) d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U .
where d 1 = 1 a 1 + b 1 + c 1 , d 2 = 1 a 2 + b 2 + c 2 , a 1 b 1 + c 1 , and a 2 b 2 + c 2 .

Assumptions and Existence Results

There are several approaches for studying the existence of a solution and stability. These include iterative methods, topological degree theory, fixed point theory, and so on. In each approach, we need to impose some auxiliary information in the form of hypotheses. Here, we use operator theory, fixed point theorems, and some results from functional analysis. Therefore, the required solution of the underlying MFDDEs Equation (2) will basically be expressed in the form of an operator equation and will also provide some subsidiary assumptions for the proposed model.
Let X = C (J,R) denote the subspace of all continuous functions defined on interval J, and let us formulate an operator T : X × X X × X by the following:
T V 1 ( t ) , V 2 ( t ) = T 1 V 1 ( t ) , V 2 ( t ) , T 2 V 1 ( t ) , V 2 ( t ) , where X × X , | | ( V 1 , V 2 ) | | is Banach space with norm given as | | ( V 1 , V 2 ) | | = sup t J { | V 1 ( t ) | } + sup t J { | V 2 ( t ) | } = | | V 1 | | + | | V 2 | | for every V 1 , V 2 X and
T 1 V 1 ( t ) , V 2 ( t ) = d 1 g 1 ( V 1 ) + i = 2 n d 1 b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U ,
T 2 V 1 ( t ) , V 2 ( t ) = d 2 g 2 ( V 2 ) + i = 2 n d 2 b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U .
We consider the following assumptions, which will be needed for further progress in this work.
H1. For V 1 , V 1 X , there exist L f 1 , L f 1 λ , L f 1 c 0 , such that,
| f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) | L f 1 | | V 1 V 1 | | + L f 1 λ | | V 1 ( λ t ) V 1 ( λ 1 t ) | | + L f 1 c | | V 2 V 2 | | .
H2. For V 1 , V 1 X , there exist L f 1 , L f 1 λ , L f 1 c 0 , such that
| f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) | L f 2 | | V 2 V 2 | | + L f 2 λ | | V 2 ( λ 2 t ) V 2 ( λ 2 t ) | | + L f 2 c | | V 1 V 1 | | .
H3. For V j , V j ( j = 1 , 2 ) X , there exist L g j ( j = 1 , 2 ) > 0 , such that
| g j ( V j ( t ) ) g j ( V j ( t ) ) | L g j | | V j V j | | .
H4. For any V 1 X , there exist H f 1 a , H f 1 b . H f 1 c : C ( J , R + ) , such that
| f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) ) , V 2 | H f 1 a ( t ) + H f 1 b ( t ) | V 1 | + H f 1 c ( t ) | V 2 | .
H5. For any V 2 X , there exist H f 2 a , H f 2 b . H f 2 c : C ( J , R + ) , such that
| f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) ) , V 1 | H f 2 a ( t ) + H f 2 b ( t ) | V 1 | + H f 2 c ( t ) | V 2 | .
H6. For any V 2 X , there exist H g j ( j = 1 , 2 ) C ( J , R + ) , such that
| g j ( V j ) | H g j ( t ) .
Remark 1.
To move the calculations from tedious operations, the following notations will be used
L 1 = | d 1 | L g 1 + i = 2 n | d 1 | | b 1 | ζ 1 α 1 α i + | d 1 | c 1 | + 1 | σ i | | σ 1 | Γ ( α 1 α i + 1 )
+ | d 1 | | b 1 | ζ 1 α 1 + | d 1 | c 1 | + 1 L f 1 + L f 1 λ + L f 1 c | σ 1 | Γ ( α 1 + 1 ) ,
L 2 = | d 2 | L g 2 + i = 2 n | d 2 | | b 2 | ζ 2 β 1 β i + | d 2 | c 2 | + 1 | η i | | η 1 | Γ ( β 1 β i + 1 )
+ | d 2 | | b 2 | ζ 2 β 1 + | d 2 | c 2 | + 1 L f 2 + L f 2 λ + L f 2 c | η 1 | Γ ( β 1 + 1 ) ,
L 1 a = i = 2 n | σ i | | σ 1 | Γ ( α 1 α i + 1 ) + L f 1 + L f 1 λ + L f 1 c | σ 1 | Γ ( α 1 + 1 ) ,
L 2 a = i = 2 n | η i | | η 1 | Γ ( β 1 β i + 1 ) + L f 2 + L f 2 λ + L f 2 c | η 1 | Γ ( β 1 + 1 ) ,
R 1 > | d 1 | | | H g 1 ( t ) | | + | | H f 1 a ( t ) | | | σ 1 | Γ ( α 1 + 1 ) ( | d 1 | | b 1 | σ 1 α 1 + | d 1 | | c 1 | + 1 ) 1 | | H f 1 b ( t ) | | + | | H f 1 c ( t ) | | | σ 1 | Γ ( α 1 + 1 ) ( | d 1 | | b 1 | σ 1 α 1 + | d 1 | | c 1 | + 1 ) i = 2 n | σ i | ( | d 1 | | b 1 | σ 1 α 1 α i + | d 1 | | c 1 | 1 ) | σ 1 | Γ ( α 1 α i + 1 ) ,
R 2 > | d 2 | | | H g 2 ( t ) | | + | | H f 2 a ( t ) | | | η 1 | Γ ( β 1 + 1 ) ( | d 2 | | b 2 | ζ 1 β 1 + | d 2 | | c 2 | + 1 ) 1 | | H f 2 b ( t ) | | + | | H f 2 c ( t ) | | | η 1 | Γ ( β 1 + 1 ) ( | d 2 | | b 2 | ζ 1 β 1 + | d 2 | | c 2 | + 1 ) i = 2 n | η i | ( | d 2 | | b 2 | ζ 1 β 1 β i + | d 2 | | c 1 | 1 ) | η 1 | Γ ( β 1 β i + 1 ) .
Theorem 4.
Consider that (H1)–(H3) holds and L < 1 ; then, the operator T has at most one fixed point, where L = max { L 1 , L 2 } and L 1 , L 2 are defined by Equations (7) and (9), respectively.
Proof. 
Consider that V j ( t ) , V j ( t ) ( j = 1 , 2 ) X ,
| T 1 V 1 ( t ) , V 2 ( t ) T 1 V 1 ( t ) , V 2 ( t ) | | d 1 | | g 1 ( V 1 ) g 1 ( V 1 ) | + | d 1 | | b 1 | | σ 1 | Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | b 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 | V 1 ( U ) V 1 ( U ) | d U + | d 1 | | c 1 | | σ 1 | Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | c 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 | V 1 ( U ) V 1 ( U ) | d U + 1 | σ 1 | Γ ( α 1 ) 0 t ( t U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | σ i | | σ 1 | Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 | V 1 ( U ) V 1 ( U ) | d U .
Now, by making use of (H1) and (H3) in the inequality Equation (15), we obtain the following:
| T 1 V 1 ( t ) , V 2 ( t ) T 1 V 1 ( t ) , V 2 ( t ) | | d 1 | | g 1 ( V 1 ) g 1 ( V 1 ) | + | d 1 | | b 1 | L f 1 + L f 1 λ | | V 1 V 1 | | + L f 1 c | | V 2 V 2 | | | σ 1 | Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 d U + i = 2 n | d 1 | | b 1 | | σ i | | V 1 V 1 | | σ 1 | Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 d U + | d 1 | | c 1 | L f 1 + L f 1 λ | | V 1 V 1 | | + L f 1 c | | V 2 V 2 | | | σ 1 | Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 d U + i = 2 n | d 1 | | c 1 | | σ i | | V 1 V 1 | | σ 1 | Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 d U + L f 1 + L f 1 λ | | V 1 V 1 | | + L f 1 c | | V 2 V 2 | | | σ 1 | Γ ( α 1 ) 0 t ( t U ) α 1 1 d U + i = 2 n | σ i | | V 1 V 1 | | σ 1 | Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 d U .
Evaluating the integral involved in Equation (16) and using Equation (7), we get
| | T 1 V 1 ( t ) , V 2 ( t ) T 1 V 1 ( t ) , V 2 ( t ) | | L 1 | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | .
By similar calculations, one can infer the following expression:
| | T 2 V 1 ( t ) , V 2 ( t ) T 2 V 1 ( t ) , V 2 ( t ) | | L 2 | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | .
Now, from Equations (17) and (18), we obtain
| | T V 1 ( t ) , V 2 ( t ) T V 1 ( t ) , V 2 ( t ) | | = | | T 1 V 1 ( t ) , V 2 ( t ) T 1 V 1 ( t ) , V 2 ( t ) | | + | | T 2 V 1 ( t ) , V 2 ( t ) T 2 V 1 ( t ) , V 2 ( t ) | | , L 1 | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | + L 2 | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | , L | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | .
Therefore, the use of the Banach contraction principle implies that the operator T possesses a fixed point. Consequently, our proposed model Equation (2) has a solution. Furthermore, this solution is unique. □
Theorem 5.
The considered model Equation (2) possesses at least one solution if (H1) and (H1) and (H4)–(H6) hold together with L a < 1 , where L a = max ( L 1 a , L 2 a ) and L 1 a , L 2 a are defined by Equations (11) and (12), respectively.
Proof. 
In order to prove the existence of at least one solution, we define the operators F , G : X × X X × X given by F V 1 ( t ) , V 2 ( t ) = F 1 V 1 ( t ) , V 2 ( t ) , F 2 V 1 ( t ) , V 2 ( t ) and G V 1 ( t ) , V 2 ( t ) = G 1 V 1 ( t ) , V 2 ( t ) , G 2 V 1 ( t ) , V 2 ( t ) , where
F 1 V 1 ( t ) , V 2 ( t ) = d 1 g 1 ( V 1 ) + i = 2 n d 1 b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U ,
F 2 V 1 ( t ) , V 2 ( t ) ( t ) = d 2 g 2 ( V 2 ) + i = 2 n d 2 b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U , G 1 V 1 ( t ) , V 2 ( t ) = 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U , G 2 V 1 ( t ) , V 2 ( t ) = 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U .
Moreover, we construct a ball H = { ( V 1 , V 2 ) X × X : | | ( V 1 , V 2 ) | | R } , with positive radius R > max ( R 1 , R 2 ) , where R 1 and R 2 are defined by Equations (13) and (14). For simplicity, we divide the proof into various steps.
Step 1: We claim in this step that F V 1 ( t ) , V 2 ( t ) + G V 1 ( t ) , V 2 ( t ) H X × X for every V 1 ( t ) , V 2 ( t ) , V 1 ( t ) , V 2 ( t ) H . To do this, we proceed as follows:
| F 1 V 1 ( t ) , V 2 ( t ) + G 1 V 1 ( t ) , V 2 ( t ) | | d 1 | g 1 ( V 1 ) | + | d 1 | | b 1 | | σ 1 | Γ α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | b 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 | V 1 ( U ) | d U + | d 1 | | c 1 | | σ 1 | Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | c 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 | V 1 ( U ) | d U + 1 | σ 1 | Γ ( α 1 ) 0 t ( t U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n 1 | σ i | | σ 1 | Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 | V 1 ( U ) | d U .
Now, by making use of (H4) and (H6) in the inequality Equation (19), we get
| | F 1 V 1 ( t ) , V 2 ( t ) + G 1 V 1 ( t ) , V 2 ( t ) | | | d 1 | | | H g 1 ( t ) | | + | | H f 1 a ( t ) | | + | | H f 1 b ( t ) | | | | V 1 | | + | | H f 1 c ( t ) | | | | V 2 | | | σ 1 | Γ ( α 1 + 1 ) ( | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 ) + i = 2 n | σ i | | | V 1 | | | σ 1 | Γ ( α 1 α i + 1 ) ( | d 1 | | b 1 | ζ 1 α 1 α i + | d 1 | | c 1 | + 1 ) .
or
| | F 1 V 1 ( t ) , V 2 ( t ) + G 1 V 1 ( t ) , V 2 ( t ) | | < R 1 .
Similarly, one can obtain the following inequality for F 2 :
| | F 2 V 1 ( t ) , V 2 ( t ) + G 2 V 1 ( t ) , V 2 ( t ) | | < R 2 .
Hence, from Equations (20) and (21), one can infer that F V 1 ( t ) , V 2 ( t ) + G V 1 ( t ) , V 2 ( t ) H .
Step 2: In this step, we claim that F is uniformly bounded. For verification, we proceed as follows:
Let V 1 , V 2 X . Then, we have
| F 1 V 1 ( t ) , V 2 ( t ) | | d 1 | g 1 ( V 1 ) | + d 1 | b 1 | | σ 1 | Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | b 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 | V 1 ( U ) | d U + | d 1 | | c 1 | | σ 1 | Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 | f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | c 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 | V 1 ( U ) | d U .
Now, by making use of (H4) and (H6) in the inequality Equation (22), we get
| | F 1 V 1 ( t ) , V 2 ( t ) | | | d 1 | | | H g 1 ( t ) | | + | d 1 | | | H f 1 a ( t ) | | + | | H f 1 b ( t ) | | | | V 1 | | + | | H f 1 c ( t ) | | | | V 2 | | | σ 1 | Γ ( α 1 + 1 ) ( | b 1 | ζ 1 α 1 + | c 1 | ) + i = 2 n | d 1 | | σ i | | | V 1 | | | σ 1 | Γ ( α 1 α i + 1 ) ( | b 1 | ζ 1 α 1 α i + | c 1 | ) .
Similarly, one can obtain the following inequality for F 2 :
| | F 2 V 1 ( t ) , V 2 ( t ) | | | d 2 | | | H g 2 ( t ) | | + | d 2 | | | H f 2 a ( t ) | | + | | H f 2 b ( t ) | | | | V 2 | | + | | H f 2 c ( t ) | | | | V 1 | | | η 1 | Γ ( β 1 + 1 ) ( | b 2 | ζ 1 β 1 + | c 2 | ) + i = 2 n | d 2 | | η i | | | V 2 | | | η 1 | Γ ( β 1 β i + 1 ) ( | b 2 | ζ 1 β 1 β i + | c 2 | ) .
Therefore, from Equations (23) and (24), one can get that F is uniformly bounded.
Step 3: Now, we claim that F is continuous. For the proof, consider a sequence ( V 1 n , V 2 n ) H , which converges to ( V 1 , V 2 ) . Now, we need to prove the relation that < F ( V 1 n , V 2 n ) > converges to < F ( V 1 , V 2 ) > , as n . Consider the following:
| F 1 V 1 n ( t ) , V 2 n ( t ) F 1 V 1 ( t ) , V 2 ( t ) | | d 1 | | g 1 ( V 1 n ) g 1 ( V 1 ) | + | d 1 | | b 1 | | σ 1 | Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 | f 1 ( U , V 1 n ( t U , V 1 n ( λ 1 U ) , V 2 n ( U ) ) f 1 ( U , V 1 ( t U , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | b 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 | V 1 n ( U ) V 1 ( U ) | d U + | d 1 | | c 1 | | σ 1 | Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 | f 1 ( U , V 1 n ( t U , V 1 n ( λ 1 U ) , V 2 n ( U ) ) f 1 ( U , V 1 ( t U , V 1 ( λ 1 U ) , V 2 ( U ) ) | d U + i = 2 n | d 1 | | c 1 | | σ i | | σ 1 | Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 | V 1 ( U ) V 1 ( U ) | d U .
Hence, by Lebsegue’s dominated convergent theorem, we have F 1 ( V 1 n , V 2 n ) F 1 ( V 1 ( t ) , V 2 ( t ) ) 0 , as n .
In the same way, the operator F 2 can be proved continuous. Furthermore, from the continuity of F 1 and F 2 , we infer that F is continuous.
Step 4: In this step, it will be proved that bounded sets will map to equicontinuous sets under F . Consider for t 1 , t 2 J , t 1 < t 2 and ( V 1 , V 2 ) H X × X ,
| | F 1 V 1 ( t 2 ) , V 2 ( t 2 ) F 1 V 1 ( t 1 ) , V 2 ( t 1 ) | | | d 1 | | g 1 ( V 1 ( t 2 ) ) g 1 ( V 1 ( t 1 ) ) | .
Clearly, as t 1 t 2 , we have | | F 1 V 1 ( t 2 ) , V 2 ( t 2 ) F 1 V 1 ( t 1 ) , V 2 ( t 1 ) | | = 0 . By a similar procedure, it can be shown that bounded sets map to the set of equicontinous functions under the operator F 2 . Moreover, the aforementioned result also holds for F .
Step 5:The contraction G can be derived in a similar manner as for Theorem 4.
Thus, all requirements of Krasnoselskii’s fixed point theorem are fulfilled; therefore, the model of fractional-order differential equations has at least one solution. □

4. Stability Analysis

The results carried out in this part of paper are specific to the functional stability of the model under consideration in Equation (2). Four types of functional-stability-related results will be explored in the context of fractional calculus. Initially, definitions of these kinds of stability results will be provided. Later on, based on some auxiliary hypotheses, results related to the existence of functional stability will be explored.
Definition 3.
The Ulam–Hyres stability for the model Equation (2) can be achieved if one can find B = max ( B 1 , B 2 ) ( c o n s t a n t ) > 0 such that for each solution ( V 1 , V 2 ) X × X of the following differential inequality and ϵ = max ( ϵ 1 , ϵ 2 ) > 0 ,
| i = 1 n σ i c D α i V 1 ( t ) f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) | ϵ 1 , t [ 0 , 1 ] , | i = 1 n η i c D β i V 2 ( t ) f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) | ϵ 1 , t [ 0 , 1 ] .
and a unique solution ( V 1 , V 2 ) X × X of the given problem Equation (2) such that | ( V 1 , V 2 ) ( V 1 , V 2 ) | B ϵ , while the solution satisfies the definition of generalized Ulam–Hyers (GUH) stability, if a positive function K : ( 0 , ) ( 0 , ) can be found with the condition K ( 0 ) = 0 such that | ( V 1 , V 2 ) ( V 1 , V 2 ) | B K ( t ) .
Definition 4.
The problem Equation (2)’s solution satisfies the UH Rassias stability criteria, with regard to a continuous function χ = max ( χ 1 , χ 2 ) X × X if we have B = max ( B 1 , B 2 ) (positive constant) > 0 , and ϵ = max ( ϵ 1 , ϵ 2 ) > 0 , for each solution ( V 1 , V 2 ) X × X of the following differential inequality:
| i = 1 n σ i c D α i V 1 ( t ) f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) | χ 1 ( t ) ϵ 1 , t [ 0 , 1 ] , | i = 1 n η i c D β i V 2 ( t ) f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) | χ 2 ( t ) ϵ 2 , t [ 0 , 1 ] .
and a unique solution ( V 1 , V 2 ) X × X for the given problem Equation (2) such that | ( V 1 , V 2 ) ( V 1 , V 2 ) | B χ ( t ) ϵ .
Definition 5.
The solution of the considered problem is GUH Rassias stable, with regard to continuous function χ = max ( χ 1 , χ 2 ) X × X if we have B = max ( B 1 , B 2 ) (positive constant) > 0 , and for each solution ( V 1 , V 2 ) X × X of the following differential inequality
| i = 1 n σ i c D α i V 1 ( t ) f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) | χ 1 ( t ) , t [ 0 , 1 ] , | i = 1 n η i c D β i V 2 ( t ) f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) | χ 2 ( t ) , t [ 0 , 1 ] .
and a unique solution ( V 1 , V 2 ) X × X of the given problem Equation (2) such that | ( V 1 , V 2 ) ( V 1 , V 2 ) | B χ ( t ) .
Remark 2.
The solution of the inequality Equation (25) is ( V 1 , V 2 ) X × X , iff one can find functions, ξ 1 , ξ 2 X such that
( i ) | ξ 1 ( t ) | ϵ 1 , | ξ 2 ( t ) | ϵ 2 , t J ,
and
( i i ) i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) + ξ 1 ( t ) , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) + ξ 2 ( t ) , t [ 0 , 1 ] .
Remark 3.
The solution of the inequality Equation (26) is ( V 1 , V 2 ) X × X , iff one can find functions ξ 1 , ξ 2 X such that
( i ) | ξ 1 ( t ) | χ 1 ( t ) ϵ 1 , | ξ 2 ( t ) | χ 2 ( t ) ϵ v , t J ,
and
( i i ) i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) + ξ 1 ( t ) , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) + ξ 2 ( t ) , t [ 0 , 1 ] .
Remark 4.
The solution of the inequality Equation (27) is ( V 1 , V 2 ) X × X , iff one can find functions, ξ 1 , ξ 2 X such that
( i ) | ξ 1 ( t ) | χ 1 ( t ) , | ξ 2 ( t ) | χ 2 ( t ) , t J ,
and
( i i ) i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) + ξ 1 ( t ) , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) + ξ 2 ( t ) , t [ 0 , 1 ] .
Lemma 2.
Consider that ( V 1 , V 2 ) X × X is a solution of system of MFDDEs
i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) + ξ 1 ( t ) , λ 1 , λ 2 , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) + ξ 2 ( t ) , α i , β i ( 0 , 1 ] , for i = 1 , 2 , 3 , , n , g 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , g 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) .
that satisfies the following relations:
| | V 1 ( t ) T 1 ( V 1 , V 2 ) | | | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 ϵ 1 | σ 1 | Γ ( α 1 + 1 ) , | | V 2 ( t ) T 2 ( V 1 , V 2 ) | | | d 2 | | b 2 | ζ 2 β 1 + | d 2 | | c 2 | + 1 ϵ 2 | η 1 | Γ ( β 1 + 1 ) .
Proof. 
Suppose ( V 1 , V 2 ) X × X to be the solution of the problem Equation (28). Then, by corollary, (1), we have
V 1 ( t ) = d 1 g 1 ( V 1 ) d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 b 1 σ i σ 1 Γ ( α 1 α i ) 0 ζ 1 ( ζ 1 U ) α 1 α i 1 V 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U + i = 2 n d 1 c 1 σ i σ 1 Γ ( α 1 α i ) 0 1 ( 1 U ) α 1 α i 1 V 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 f 1 ( U , V 1 ( U ) , V 1 ( λ 1 U ) , V 2 ( U ) ) d U i = 2 n σ i σ 1 Γ ( α 1 α i ) 0 t ( t U ) α 1 α i 1 V 1 ( U ) d U d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 ξ 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 ξ 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 ξ 1 ( U ) d U ,
V 2 ( t ) = d 2 g 2 ( V 2 ) d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 b 2 η i η 1 Γ ( β 1 β i ) 0 ζ 2 ( ζ 2 U ) β 1 β i 1 V 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U + i = 2 n d 2 c 2 η i η 1 Γ ( β 1 β i ) 0 1 ( 1 U ) β 1 β i 1 V 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 f 2 ( U , V 2 ( U ) , V 2 ( λ 2 U ) , V 1 ( U ) ) d U i = 2 n η i η 1 Γ ( β 1 β i ) 0 t ( t U ) β 1 β i 1 V 2 ( U ) d U d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 ξ 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 ξ 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 ξ 2 ( U ) d U .
For simplicity, let us introduce T 1 ( V 1 , V 2 ) and T 2 ( V 1 , V 2 ) in Equations (29) and (30). Then, we obtain
V 1 ( t ) T 1 ( V 1 , V 2 ) = d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 + ξ 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 ξ 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 ξ 1 ( U ) d U , V 2 ( t ) T 2 ( V 1 , V 2 ) = d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 ξ 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 ξ 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 ξ 2 ( U ) d U .
By taking the absolute and using Remark (2), we have
| | V 1 ( t ) T 1 ( V 1 , V 2 ) | | | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 ϵ 1 | σ 1 | Γ ( α 1 + 1 ) , | | V 2 ( t ) T 2 ( V 1 , V 2 ) | | | d 2 | | b 2 | ζ 2 β 1 + | d 2 | | c 2 | + 1 ϵ 2 | η 1 | Γ ( β 1 + 1 ) .
This proves the required result. □
Theorem 6.
Under the assumptions (H1)–(H3), the fractional-order model Equation (2) is UH stable and GUH stable if L < 1 , where L = max ( L 1 , L 2 ) and L 1 , L 2 are defined by Equations ( 7 ) a n d ( 9 ) , respectively.
Proof. 
Given any solution ( V 1 , V 2 ) × X , and the unique solution ( V 1 , V 2 ) of the given problem Equation (2), the following holds:
| | ( V 1 , V 2 ) ( V 1 , V 2 ) | | = | | ( V 1 , V 2 ) T ( V 1 , V 2 ) | | , | | V 1 T 1 ( V 1 , V 2 ) | | + | | T 1 ( V 1 , V 2 ) T 1 ( V 1 , V 2 ) | | + | | V 2 T 2 ( V 1 , V 2 ) | | + | | T 2 ( V 1 , V 2 ) T 2 ( V 1 , V 2 ) | | .
Using Theorem 4 and Lemma 2, we have
| | ( V 1 , V 2 ) ( V 1 , V 2 ) | | | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 ϵ 1 | σ 1 | Γ ( α 1 + 1 ) + | d 2 | | b 2 | ζ 2 β 1 + | d 2 | | c 2 | + 1 ϵ 2 | η 1 | Γ ( β 1 + 1 ) + L | | ( V 1 , V 2 ) ( V 1 , V 2 ) | | , [ | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 ϵ 1 | σ 1 | Γ ( α 1 + 1 ) + | d 2 | | b 2 | ζ 2 β 1 + | d 2 | | c 2 | + 1 ϵ 2 | η 1 | Γ ( β 1 + 1 ) ] 1 1 L , B ϵ .
Let B = max | d 1 | | b 1 | ζ 1 α 1 + | d 1 | | c 1 | + 1 | σ 1 | Γ ( α 1 + 1 ) ( 1 L ) , | d 2 | | b 2 | ζ 2 β 1 + | d 2 | | c 2 | + 1 | η 1 | Γ ( β 1 + 1 ) ( 1 L ) ; then, the solution of the considered problem Equation (2) is UH stable. Furthermore, we can set K ( ϵ ) = ϵ , ; then, the considered problem Equation (2) is GUH stable. □
To prove the next stability result, we need the following assumption, given as (H7). For any non-decreasing function χ 1 , χ 2 X , there exist the positive constants G 1 and G 2 , such that
1 Γ ( α ) 0 t ( t U ) α 1 1 χ 1 ( U ) d U G 1 χ 1 ( t ) , 1 Γ ( β ) 0 t ( t U ) β 1 1 χ 2 ( U ) d U G 2 χ 2 ( t ) .
Lemma 3.
Consider that ( V 1 , V 2 ) X × X is a solution of a system of MFDDEs,
i = 1 n σ i c D α i V 1 ( t ) = f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) + ξ 1 ( t ) , λ 1 , λ 2 , t [ 0 , 1 ] , i = 1 n η i c D β i V 2 ( t ) = f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) + ξ 2 ( t ) , α i , β i ( 0 , 1 ] , for i = 1 , 2 , 3 , , n , g 1 ( V 1 ) = a 1 V 1 ( 0 ) b 1 V 1 ( ζ 1 ) c 1 V 1 ( 1 ) , g 2 ( V 2 ) = a 2 V 2 ( 0 ) b 2 V 2 ( ζ 2 ) c 2 V 2 ( 1 ) .
that satisfies the following relations,
| | V 1 ( t ) T 1 ( V 1 , V 2 ) | | | d 1 | | b 1 | + | d 1 | | c 1 | + 1 G 1 χ 1 ( t ) ϵ 1 | σ 1 | Γ ( α 1 + 1 ) , | | V 2 ( t ) T 2 ( V 1 , V 2 ) | | | d 2 | | b 2 | + | d 2 | | c 2 | + 1 G 2 χ 2 ( t ) ϵ 2 | η 1 | Γ ( β 1 + 1 ) .
Proof. 
Suppose a solution ( V 1 , V 2 ) X × X of the problem Equation (31). Then, by Corollary (1), we have
V 1 ( t ) T 1 ( V 1 , V 2 ) = d 1 b 1 σ 1 Γ ( α 1 ) 0 ζ 1 ( ζ 1 U ) α 1 1 + ξ 1 ( U ) d U d 1 c 1 σ 1 Γ ( α 1 ) 0 1 ( 1 U ) α 1 1 ξ 1 ( U ) d U + 1 σ 1 Γ ( α 1 ) 0 t ( t U ) α 1 1 ξ 1 ( U ) d U , V 2 ( t ) T 2 ( V 1 , V 2 ) = d 2 b 2 η 1 Γ ( β 1 ) 0 ζ 2 ( ζ 2 U ) β 1 1 ξ 2 ( U ) d U d 2 c 2 η 1 Γ ( β 1 ) 0 1 ( 1 U ) β 1 1 ξ 2 ( U ) d U + 1 η 1 Γ ( β 1 ) 0 t ( t U ) β 1 1 ξ 2 ( U ) d U .
By taking the absolute and using Remark (3), we have
| | V 1 ( t ) T 1 ( V 1 , V 2 ) | | | d 1 | | b 1 | + | d 1 | | c 1 | + 1 G 1 χ 1 ( t ) ϵ 1 | σ 1 | Γ ( α 1 + 1 ) , | | V 2 ( t ) T 2 ( V 1 , V 2 ) | | | d 2 | | b 2 | + | d 2 | | c 2 | + 1 G 2 χ 2 ( t ) ϵ 2 | η 1 | Γ ( β 1 + 1 ) .
This proves the required result. □
Theorem 7.
If the assumptions (H1)–(H3) and (H7) are supposed, then the problem Equation (2) is UHR stable and GUHR stable, if L < 1 , where L = max ( L 1 , L 2 ) and L 1 , L 2 are given by Equations ( 7 ) and ( 9 ) , respectively.
Proof. 
Suppose a given general solution U C J , R ) , as well as a unique solution U , for the model Equation (2). Then,
| | ( V 1 , V 2 ) ( V 1 , V 2 ) | | | | V 1 T 1 ( V 1 , V 2 ) | | + | | T 1 ( V 1 , V 2 ) T 1 ( V 1 , V 2 ) | | + | | V 2 T 2 ( V 1 , V 2 ) | | + | | T 2 ( V 1 , V 2 ) T 2 ( V 1 , V 2 ) | | .
Using Theorem (4) and Lemma (3), we have
| | ( V 1 , V 2 ) ( V 1 , V 2 ) | | [ | d 1 | | b 1 | + | d 1 | | c 1 | + 1 G 1 χ 1 ( t ) ϵ 1 | σ 1 | Γ ( α 1 + 1 ) + | d 2 | | b 2 | + | d 2 | | c 2 | + 1 G 2 χ 2 ( t ) ϵ 2 | η 1 | Γ ( β 1 + 1 ) ] 1 1 L , B χ ( t ) ϵ .
Let B = max | d 1 | | b 1 | + | d 1 | | c 1 | + 1 G 1 | σ 1 | Γ ( α 1 + 1 ) ( 1 L ) , | d 2 | | b 2 | + | d 2 | | c 2 | + 1 G 2 | η 1 | Γ ( β 1 + 1 ) ( 1 L ) ; then, the solution of the considered problem Equation (2) is UHR stable. Furthermore, set ϵ = 1 , ; then, the proposed problem Equation (2) is GUHR stable. □
To see the applications of the obtained results, an illustrative example is provided in the next part.
Example 1.
Taking some specific functions in the model Equation (2), we construct the following MFDDEs:
i = 1 100 14 3 i c D 3 5 i V 1 ( t ) = e sin ( t ) 14 V 1 ( t ) + log ( 1 + sin ( t ) ) 71 V 1 ( t 5 ) c o h ( e t ) 34 V 2 ( t ) + cosh ( t 3 ) 9 e t t 2 , i = 1 100 33 8 i c D 3 7 i V 2 ( t ) = tan ( x ) 45 V 2 ( t ) + c o s 1 ( t 0.5 ) ) 13 V 2 ( t 7 ) c o h ( t 2 + 1 ) 53 V 1 ( t ) + cosh ( t 3 ) 9 e t t 2 , 3 V 1 18 + | V 1 | = 1 14 V 1 ( 0 ) V 1 ( ζ 1 ) 2 50 V 1 ( 1 ) , t [ 0 , 1 ] e t V 2 12 ( t 2 + 5 ) + | V 2 | = 2 15 V 2 ( 0 ) 3 4 V 2 ( ζ 2 ) 1 100 V 2 ( 1 ) .
Here, n = 100 , α i = 3 5 i , β i = 3 7 i , a 1 = 1 14 , a 2 = 2 15 , b 1 = 1 ; b 2 = 3 4 , c 1 = 2 50 , c 2 = 1 100 , σ i = 14 3 i , η i = 33 8 i , ζ 1 = 1 24 , ζ 2 = 2 31 , λ 1 = 1 5 , λ 2 = 1 7 ,
f 1 ( t , V 1 ( t ) , V 1 ( λ 1 t ) , V 2 ( t ) ) = e sin ( t ) 14 V 1 ( t ) + log ( 1 + sin ( t ) ) 71 V 1 ( t 5 ) c o h ( e t ) 34 V 2 ( t ) + cosh ( t 3 ) 9 e t t 2 , f 2 ( t , V 2 ( t ) , V 2 ( λ 2 t ) , V 1 ( t ) ) = tan ( x ) 45 V 2 ( t ) + c o s 1 ( t 0.5 ) ) 13 V 2 ( t 7 ) c o h ( t 2 + 1 ) 53 V 1 ( t ) + cosh ( t 3 ) 9 e t t 2 , g 1 ( V 1 ) = 3 V 1 18 + | V 1 | , and g 2 ( V 2 ) = e t V 2 12 ( t 2 + 5 ) + | V 2 | .
Now, the reader can easily obtain from Equation (33) that L f 1 = 1 14 , L f 1 λ = log ( 1 + sin ( 1 ) ) 71 , L f 1 c = cosh ( e ) 34 , L f 2 = tan ( 1 ) 45 , L f 2 λ = 1.0471975512 ; L f 2 c = cosh ( 2 ) 53 , L g 1 = 3 18 , L g 2 = e 65 , consequently L = 0.9539 < 1 . Therefore, by Theorem (4), the problem Equation (32) has a unique solution. Set χ 1 ( t ) = t + 1 , and χ 2 ( t ) = 2 t 2 + 3 t [ 0 , 1 ] . Then, we can get the following:
1 Γ ( 3 5 ) 0 t ( t U ) 3 5 1 ( U + 1 ) d U 1.12 ( t + 1 ) , 1 Γ ( 3 5 ) 0 t ( t U ) 3 5 1 ( 2 U 2 + 3 ) d U 1.12 ( t + 1 ) 1.13 ( 2 t 2 + 3 )
As a condition, (H7) is satisfied with G 1 = 1.12 and G 2 = 1.13 . Therefore, all requirements of Theorem (7) are fulfilled, so this leads to the fact that the solution of the problem Equation (2) is UHR stable with respect to a continuous function χ ( t ) = 2 t 2 + 3 together with constant B = 13.1184 . Furthermore, the solution of the aforementioned problem is GUHR stable. In addition to this, the problem Equation (2) is UH and GUH stable because of Theorem (6).

5. Conclusions

Findings of the research are described here. We derived some adequate results regarding the existence and functional stability of a class of nonlocal BVP of coupled systems of MFDEs with proportional type delay term. We have used the results from functional analysis and FPT to derive the conditions for the solution’s existence and uniqueness. Furthermore, results related to UH stabilities were also investigated. The results reached in this manuscript generalize some results from the existing literature. In addition, the outcomes have been illustrated through a proper example. Moreover, the special case of the considered problem Equation (2) can be obtained by fixing the parameters, η i = σ i = 0 , for i = 2 , 3 , , n ) , which was investigated in [39]. As future work, we recommend studying the present problem, as well as other similar expected problems, using symmetry methods.

Author Contributions

Conceptualization, G.u.R., D.A. and R.P.A.; methodology, G.u.R. and D.A.; software, G.u.R. and D.A.; validation, G.u.R. and R.P.A.; formal analysis, G.u.R. and D.A.; investigation, G.u.R., D.A. and R.P.A.; resources, G.u.R. and D.A.; data curation, D.A. and R.P.A.; writing—G.u.R. and D.A.; writing—review and editing, G.u.R., D.A. and R.P.A.; visualization, G.u.R., D.A. and R.P.A.; supervision, G.u.R. and R.P.A.; project administration, G.u.R. and R.P.A.; funding acquisition, R.P.A. 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.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agarwal, R.P.; Lupulescu, V.; O’Regan, D.; Rahman, G. Fractional calculus and fractional differential equations in nonreflexive Banach spaces. Commun. Nonlinear Sci. Numer. Simul. 2005, 20, 59–73. [Google Scholar] [CrossRef]
  2. Agarwal, R.P.; O’Regan, D.; Hristova, S.; Ciceke, M. Practical stability with respect to initial time difference for Caputo fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2017, 42, 106–120. [Google Scholar] [CrossRef]
  3. Jong, K.; Choi, H.; Ri, Y. Existence of positive solutions of a class of multi-point boundary value problems for p-Laplacian fractional differential equations with singular source terms. Commun. Nonlinear Sci. Numer. Simul. 2019, 72, 272–281. [Google Scholar] [CrossRef]
  4. Garrappa, R.; Kaslik, E. On initial conditions for fractional delay differential equations. Commun. Nonlinear Sci. Numer. Simul. 2020, 90, 105359. [Google Scholar] [CrossRef]
  5. Fallahgoul, H.; Focardi, S.; Fabozzi, F. Fractional Calculus and Fractional Processes with Applications to Financial Economics: Theory and Application; Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  6. Javidi, M.; Ahmad, B. Dynamic analysis of time fractional order phytoplankton–toxic phytoplankton–zooplankton system. Ecol. Model. 2015, 318, 8–18. [Google Scholar] [CrossRef]
  7. Ding, Y.; Wang, Z.; Ye, H. Optimal control of a fractional-order HIV-immune system with memory. IEEE Trans. Control Syst. Technol. 2011, 20, 763–769. [Google Scholar] [CrossRef]
  8. Carvalho, A.; Pinto, C.M. A delay fractional order model for the co-infection of malaria and HIV/AIDS. Int. J. Dynam. Control 2017, 5, 168–186. [Google Scholar] [CrossRef]
  9. Xu, Y.; Li, W. Finite-time synchronization of fractional-order complex-valued coupled systems. Phys. Stat. Mech. Its Appl. 2020, 549, 123903. [Google Scholar] [CrossRef]
  10. Zhang, F.; Chen, G.; Li, C.; Kurths, J. Chaos synchronization in fractional differential systems. Philos. Trans. Soc. 2013, 371, 1990. [Google Scholar] [CrossRef]
  11. Ali, A.; Shah, K.; Khan, R.A. Numerical treatment for traveling wave solutions of fractional Whitham-Broer-Kaup equations. Alex. Eng. J. 2018, 57, 1991–1998. [Google Scholar] [CrossRef]
  12. Chen, Y.; Sun, Y.; Liu, L. Numerical solution of fractional partial differential equations with variable coefficients using generalized fractional-order Legendre functions. Appl. Math. Comput. 2014, 244, 847–858. [Google Scholar] [CrossRef]
  13. Torres-Hernandez, A.; Brambila-Paz, F.; Torres-Martínez, C. Numerical solution using radial basis functions for multidimensional fractional partial differential equations of type Black–Scholes. Comput. Appl. Math. 2021, 40, 245. [Google Scholar] [CrossRef]
  14. Rehman, M.; Baleanu, D.; Alzabut, J.; Ismail, M.; Saeed, U. Green-Haar wavelets method for generalized fractional differential equations. Adv. Diff. Equ. 2020, 2020, 515. [Google Scholar] [CrossRef]
  15. Izadi, M.; Cattani, C. Generalized Bessel polynomial for multi-order fractional differential equations. Symmetry 2020, 12, 1260. [Google Scholar] [CrossRef]
  16. Tuan, H.T.; Thai, H.D.; Garrappa, R. An analysis of solutions to fractional neutral differential equations with delay. Commun. Nonlinear Sci. Numer. Simul. 2021, 10, 105854. [Google Scholar] [CrossRef]
  17. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  18. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integral and Derivatives (Theory and Applications); Taylor & Francis: Oxfordshire, UK, 1993. [Google Scholar]
  19. Zhou, Y.; Wang, J.; Zhang, L. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2016. [Google Scholar]
  20. Abbas, S.; Benchohra, M.; Graef, J.R.; Henderson, J. Implicit Differential and Integral Equations: Existence and Stability; De Gruyter: London, UK, 2018. [Google Scholar]
  21. Hilfer, R. Applications of Fractional Calculus in Physics; World Science Publishing: River Edge, NJ, USA, 2000. [Google Scholar]
  22. Haq, F.; Akram, M.; Shah, K.; Rahman, G. Study of new monotone iterative technique for a class of arbitrary order differential equations. Comput. Methods Differ. Equ. 2020, 8, 639–647. [Google Scholar]
  23. Shah, K.; Hussain, W. Investigating a class of nonlinear fractional differential equations and its Hyers-Ulam stability by means of topological degree theory. Numer. Funct. Anal. Optim. 2019, 40, 1355–1372. [Google Scholar] [CrossRef]
  24. Vu, H.; Van Hoa, N. Hyers-Ulam stability of fuzzy fractional Volterra integral equations with the kernel ψ-function via successive approximation method. Fuzzy Sets Syst. 2021, 419, 67–98. [Google Scholar] [CrossRef]
  25. Liu, X.; Liu, L.; Wu, Y. Existence of positive solutions for a singular nonlinear fractional differential equation with integral boundary conditions involving fractional derivatives. Bound. Value Probl. 2018, 2018, 24. [Google Scholar] [CrossRef] [Green Version]
  26. Shah, K.; Khalil, H.; Khan, R.A. Investigation of positive solution to a coupled system of impulsive boundary value problems for nonlinear fractional order differential equations. Chaos Solitons Fract. 2015, 77, 240–246. [Google Scholar] [CrossRef]
  27. Ali, G.; Shah, K.; Abdeljawad, T.; Khan, H.; Rahman, G.U.; Khan, A. On existence and stability results to a class of boundary value problems under Mittag-Leffler power law. Adv. Differ. Equ. 2020, 2020, 407. [Google Scholar] [CrossRef]
  28. Ali, G.; Shah, K.; Rahman, G.U. Existence of solution to a class of fractional delay differential equation under multi-points boundary conditions. Arab. J. Basic Appl. Sci. 2020, 27, 471–479. [Google Scholar] [CrossRef]
  29. Faizullah, F.; Bux, M.; Rana, M.A.; Rahman, G.U. Existence and stability of solutions to non-linear neutral stochastic functional differential equations in the framework of G-Brownian motion. Adv. Differ. Equ. 2017, 2017, 350. [Google Scholar] [CrossRef]
  30. Ahmad, I.; Shah, K.; Rahman, G.U.; Baleanu, D. Stability analysis for a nonlinear coupled system of fractional hybrid delay differential equations. Math. Methods Appl. Sci. 2020, 43, 8669–8682. [Google Scholar] [CrossRef]
  31. Ahmad, B.; Alsaedia, A.; Alsulamia, M.; Ntouyasa, S.K. A study of a coupled system of nonlinear second-order ordinary differential equations with nonlocal integral multi-strip boundary conditions on an arbitrary domain. J. Comput. Anal. Appl. 2021, 29, 215–235. [Google Scholar]
  32. Ali, A.; Shah, K.; Li, Y. Topological degree theory and Ulam’s stability analysis of a boundary value problem of fractional differential equations. In Frontiers in Functional Equations and Analytic Inequalities; Springer: Berlin/Heidelberg, Germany, 2019; pp. 73–92. [Google Scholar]
  33. Ahmed, I.; Kuman, P.; Shah, K.; Borisut, P.; Sitthithakerngkiet, K.; Demba, M.A. Stability results for implicit fractional pantograph differential equations via ψ-Hilfer fractional derivative with a nonlocal Riemann–Liouville fractional integral condition. Mathematics 2020, 8, 94. [Google Scholar] [CrossRef] [Green Version]
  34. Shah, K.; Ullah, A.; Nieto, J.J. Study of fractional order impulsive evolution problem under nonlocal Cauchy conditions. Math. Methods Appl. Sci. 2021, 44, 8516–8527. [Google Scholar] [CrossRef]
  35. Benchohra, M.; Hamani, S.; Ntouyas, S.K. Boundary value problems for differential equations with fractional order. Surv. Math. Appl. 2008, 3, 2391–2396. [Google Scholar]
  36. Ahmad, B. Existence results for multi-point nonlinear boundary value problems for fractional differential equations. Mem. Differ. Equ. Math. Phys. 2010, 49, 83–94. [Google Scholar]
  37. Komarova, N.L.; Newell, A.C. Nonlinear dynamics of sand banks and sand waves. J. Fluid Mech. 2000, 415, 285–321. [Google Scholar] [CrossRef] [Green Version]
  38. Riecke, H. Self-trapping of traveling-wave pulses in binary mixture convection. Phys. Rev. Lett. 1992, 68, 301. [Google Scholar] [CrossRef]
  39. Ahmad, I.; Nieto, J.J.; Rahman, G.U.; Shah, K. Existence and stability for fractional order pantograph equations with nonlocal conditions. Electron. J. Differ. Equ. 2020, 2020, 1–16. [Google Scholar]
  40. Baleanu, D.; Etemad, S.; Mohammadi, H.; Rezapour, S. A novel modeling of boundary value problems on the glucose graph. Commun. Nonlinear Sci. Numer. Simul. 2021, 100, 105844. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ahmad, D.; Agarwal, R.P.; ur Rahman, G. Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations. Symmetry 2022, 14, 1342. https://doi.org/10.3390/sym14071342

AMA Style

Ahmad D, Agarwal RP, ur Rahman G. Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations. Symmetry. 2022; 14(7):1342. https://doi.org/10.3390/sym14071342

Chicago/Turabian Style

Ahmad, Dildar, Ravi P. Agarwal, and Ghaus ur Rahman. 2022. "Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations" Symmetry 14, no. 7: 1342. https://doi.org/10.3390/sym14071342

APA Style

Ahmad, D., Agarwal, R. P., & ur Rahman, G. (2022). Formulation, Solution’s Existence, and Stability Analysis for Multi-Term System of Fractional-Order Differential Equations. Symmetry, 14(7), 1342. https://doi.org/10.3390/sym14071342

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