Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Next Article in Journal
Fractional Derivative Model on Physical Fractal Space: Improving Rock Permeability Analysis
Previous Article in Journal
Pavement Crack Detection Using Fractal Dimension and Semi-Supervised Learning
Previous Article in Special Issue
An Application of Multiple Erdélyi–Kober Fractional Integral Operators to Establish New Inequalities Involving a General Class of Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

β–Ulam–Hyers Stability and Existence of Solutions for Non-Instantaneous Impulsive Fractional Integral Equations

1
Department of Mathematics, National Kaohsiung Normal University, Kaohsiung 82444, Taiwan
2
Department of Mathematical Analysis and Numerical Mathematics, Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, Mlynskaá Dolina, 842 48 Bratislava, Slovakia
3
Faculty of Technical Sciences, University of Novi Sad, Trg D. Obradovića 6, 21125 Novi Sad, Serbia
4
Department for Mathematics and Informatics, Faculty of Civil Engineering, Ss. Cyril and Methodius University in Skopje, Partizanski Odredi 24, P.O. Box 560, 1000 Skopje, North Macedonia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(8), 469; https://doi.org/10.3390/fractalfract8080469
Submission received: 26 June 2024 / Revised: 7 August 2024 / Accepted: 9 August 2024 / Published: 12 August 2024
(This article belongs to the Special Issue Fractional Integral Inequalities and Applications, 2nd Edition)

Abstract

:
In this paper, we investigate a class of non-instantaneous impulsive fractional integral equations. Utilizing the Banach contraction mapping principle, we establish the existence and uniqueness of solutions for the considered problem. Additionally, employing Schauder’s fixed-point theorem, we demonstrate the existence of solutions within the framework of β -Banach spaces. Moreover, we examine the β –Ulam–Hyers stability of the solutions, providing insights into the stability behavior under small perturbations. An illustrative example is presented to demonstrate the practical applicability and effectiveness of the theoretical results obtained.

1. Introduction and Preliminaries

The impulsive differential equations describe evolution processes characterized by the fact that at certain moments they experience an abrupt change of state, i.e., these processes are subject to short-term perturbations whose duration is negligible compared with the duration of the whole process (see, e.g., [1,2,3,4,5,6] for more details on the subject). Hence, there is a substantial level of interest in assessing the quality of solutions for these impulsive systems. In a broader context, the classical instantaneous impulses fall short in capturing certain aspects of evolution processes. The establishment of hemodynamic equilibrium in an individual exemplifies this principle, wherein drugs are introduced into the bloodstream and assimilated by the body over a prolonged and continuous duration (see [7] (Example 1.1.3)). This particular scenario depicts a unique manifestation of impulsive behavior, commencing at a defined moment and enduring for a limited duration thereafter. In ecological studies, the population dynamics of species are often subject to various forms of disturbances, such as natural disasters, climate changes, and human interventions [1,7,8]. Traditional ecological models, which rely on instantaneous perturbations, may struggle to represent the gradual and continuous impacts of these disturbances. As a general theory, non-instantaneous impulsive fractional differential equations have significant potential to be applied in these areas to better capture the continuous and memory-like characteristics of drug absorption, ecology, and population dynamics. For example, this approach provides a more accurate representation of how drugs are gradually absorbed into the bloodstream, making it an essential tool for pharmaceutical research and dosage optimization.
We recall the definition of a β -Banach space from [9]. Let X be a vector space over the field K and β ( 0 , 1 ] . A function · β : X [ 0 , + ) is called a β -norm if and only if the following holds:
(i)
w β = 0 if and only if w = 0 ;
(ii)
λ w β = | λ | β w β for all λ K and w X ;
(iii)
v + w β v β + w β .
The pair ( X , · β ) is called a β -normed space. A complete β -normed space is called a β -Banach space. Let I = [ 0 , 1 ] . The space C ( I : R ) signifies the β -Banach space comprising continuous functions w : I R , endowed with the β -norm
w β : = max { | w ( s ) | β : s I , 0 < β 1 } .
The β -Banach space P C ( I : R ) is given by P C ( I : R ) = { f : I R , f C ( ( t i , t i + 1 ] : R ) , f ( t i ) = f ( t i ) a n d f ( t i + ) exists for any i = 0 , 1 , 2 , , m } , where the symbols f ( t i ) and f ( t i + ) denote the left and the right limits of the function f ( t ) at the point t = t i , i = 0 , 1 , 2 , , m , and it is endowed with the P β -norm f P β = sup { | f ( s ) | β : s I , 0 < β 1 } . Moreover, we denote f P C = max { sup s I f ( s + ) , sup s I f ( s ) } , because f P C ( I : R ) . Let w L 1 ( I ) , γ > 0 . Recall that the definition of the Riemann–Liouville fractional integral of w of order γ is given by
J 0 , t γ w ( t ) : = 1 Γ ( γ ) 0 t ( t s ) γ 1 w ( s ) d s .
For a more detailed approach to fractional calculus, interested readers may consult [1,10,11,12] and the references therein.
In light of more generalized theories, E. Hernández and D. O’Regan [13] originally embarked on a study of a novel class of abstract cases involving semilinear impulsive differential modeling that excluded instantaneous impulses; cf. also the research article [14] by M. Pierri, D. O’Regan, and V. Rolnik. Nevertheless, it is noteworthy that drug absorption exhibits a memory effect. Furthermore, fractional calculus, with its memory and hereditary properties, offers a more realistic approach to modeling of the aforementioned ecological systems. As indicated by I. Podlubny [11], fractional-order differential equations can effectively describe complex ecological dynamics by accounting for the lasting effects of past disturbances.
The famous problem for the stability of functional equations originally posed by Ulam [15] was the following:
(Ulam) 
Assume that   Z 1  is a group and  Z 2  is a metric group with the metric  ρ . Given   ε > 0 , does there exist a   δ > 0  such that if a function  U : Z 1 Z 2  satisfies the inequality
ρ ( U ( a b ) , U ( a ) U ( b ) ) < δ
for all   a , b Z 1 , then there exists a homomorphism   P : Z 1 Z 2  with
ρ ( U ( c ) , P ( c ) ) < ε
for all   c Z 1 ?
The problem of Ulam stability [15], as extensively explored by many mathematicians, has garnered considerable attention. For a more comprehensive understanding of this topic, interested readers are encouraged to consult the remarkable monographs by D.H. Hyers et al. [16], S.-M. Jung [9], and M.T. Rassias [17], as well as the recent contributions [18,19,20,21,22,23,24] in normed spaces. The corresponding results in the realm of β -normed spaces can be found in [25,26,27].
In recent studies, significant advancements have been made in the theory and applications of non-instantaneous impulsive fractional differential equations. Kostić et al. [28] studied generalized ρ -almost periodic sequences and established several new results on the existence and uniqueness of solutions for the abstract impulsive Volterra integro-differential equations. In [29], the existence and uniqueness of continuous solutions for a class of the fractional quadratic functional integro-differential equations with nonlocal fractional-order integro-differential condition were studied. Some recent results on the β –Ulam–Hyers stability of fractional differential equations and also on the existence and uniqueness of piecewise continuous solutions on non-instantaneous impulsive fractional differential equations are given in [27]. The results on the solvability of nonlocal fractional integro-differential equations are presented in [30,31], while the impulsive integro–differential equations are considered in [32]. Agarwal et al. extensively covered both basic and advanced theoretical aspects of non-instantaneous impulsive Caputo fractional differential equations in their works [7,33]. Shah et al. investigated the Hyers–Ulam stability and exponential dichotomy in impulsive linear systems of the first order [34], while Luo et al. considered the existence and Hyers–Ulam stability of solutions for fractional differential equations involving time-varying delays and non-instantaneous impulses [35]. Abbas et al. provided unique insights into Ulam-type stability concepts for the Darboux problem of partial functional differential equations with non-instantaneous impulses in Banach spaces [36]. In a related study, Ibrahim examined the generalized Hyers–Ulam stability for certain types of fractional differential equations in complex Banach spaces [37]. Furthermore, Parthasarathy explored the existence and Hyers–Ulam stability of nonlinear impulsive differential equations with nonlocal conditions [38].
Agarwal et al. also provided sufficient conditions for various types of Ulam–Hyers stability in the context of boundary-value problems for scalar nonlinear differential equations with delays and generalized proportional Caputo fractional derivatives [39]. Dhayal and his collaborators discussed the stability and exact controllability of non-instantaneous ϕ -Caputo fractional systems [40]. The Ulam–Hyers stability of nonlinear ψ -Hilfer fractional differential equations has been addressed in several studies [41,42,43].
In another significant contribution, Agarwal et al. obtained uniqueness results using the Banach contraction mapping principle and established existence results with Krasnoselskii’s fixed-point theorem and Leray–Schauder’s nonlinear alternative, while also discussing Ulam–Hyers stability for the considered problem [44]. Qian et al. focused on establishing stability theorems for fractional differential systems with the Riemann–Liouville derivative, covering linear, perturbed, and time-delayed systems [45]. Alam and Shah investigated the Hyers–Ulam stability of coupled implicit fractional integro-differential equations with Riemann–Liouville derivatives under specific conditions [46].
In [27], inspired by the investigation in [14], the following class of non-instantaneous impulsive Caputo fractional differential equations was considered:
D s i , t α 1 C w ( t ) = c w ( t ) + f ( t , w ( t ) ) , s i < t t i + 1 , i = 1 , 2 , , m , c 0 , w ( t ) = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) , t i < t s i , i = 1 , 2 , , m ,
where D s i , t α 1 C is the Caputo fractional derivative of the order α 1 , α 3 ( 0 , 1 ) , α 1 α 3 , with the lower limit s i ; 0 = s 0 < t 1 s 1 t 2 < . . . < t m s m t m + 1 = T are prefixed numbers; f : [ 0 , T ] × R R is continuous; and g i : [ t i , s i ] × R R is continuous for all i = 1 , 2 , , m , a R , J t i , t α 3 g i ( t , w ( t ) ) and J 0 , s i α 1 f ( s i , w ( s i ) ) .
The motivation for investigating problem Equation (1) stems from the above problem and its close connection to the quadratic fractional functional integro-differential equation with a nonlocal fractional integro-differential condition studied in [29], having significant potential for applications. Inspired by the aforementioned papers, this paper will explore the following non-instantaneous impulsive fractional integral problem:
w ( t ) = w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) , 0 < t t 1 , w ( t ) = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) , t i < t s i , i = 1 , 2 , , m , α 1 α 3 , w ( t ) = a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) , s i < t t i + 1 ,
where
v ( t ) : = J 0 , t 1 α 1 f t , v ( t ) · 0 t ( t θ ) α 2 1 Γ ( α 2 ) g θ , 0 θ ( θ ν ) α 1 α 3 1 Γ ( α 1 α 3 ) v ( ν ) d ν d θ
and a , w 0 R , α 1 , α 2 , α 3 , α 4 ( 0 , 1 ] , α 3 , α 4 α 1 , α 1 α 3 , and 0 s i t i + 1 1 for i = 0 , 1 , 2 , , m are prefixed numbers, defining the intervals where non-instantaneous impulses occur, and J 0 , t α is the fractional Riemann–Liouville integral of order α ; f : I × R R , g i : [ t i , s i ] × R R for all i = 1 , 2 , , m , and f 1 : I × R × R R are functions. Note that,
J t i , t α 3 g i ( t , w ( t ) ) = 1 Γ ( α 3 ) t i t ( t s ) α 3 1 g i ( s , w ( s ) ) d s
and
J 0 , s i α 1 f ( s i , w ( s i ) ) = 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 1 f ( s , w ( s ) ) d s
for all i = 1 , 2 , , m .
As far as we are aware, there has been no prior exploration of the existence of a solution or examination of Ulam-type stability for the Equation (1), nor for analogous instances of non-instantaneous impulsive fractional integral equations in β -normed spaces. In this work, we shall employ techniques from fractional differential equations, non–instantaneous impulsive differential equations, and fixed-point theory within β -Banach spaces to investigate our problem under consideration. The provided example emphasizes the importance of the obtained results and potential for application in the aforementioned and many other areas.
The structure of the paper can be outlined as follows. Following an initial presentation of preliminary definitions and results, the paper proceeds to establish results pertaining to the existence and uniqueness of solutions for the class of non-instantaneous impulsive fractional integral equations under consideration. In Section 3, the focus shifts to investigating the β –Ulam–Hyers stability for non-instantaneous impulsive fractional integral equations. Additionally, in Section 4, the paper incorporates two illustrative examples that underscores the significance of the derived results. The paper is ended by summarizing the obtained outcomes and outlining prospective research directions.

2. The Existence and Uniqueness of Solutions

In this study, we employ a set of conditions that will be specified subsequently. These conditions can be broadly classified into three groups: (C1)–(C3) represent Lipschitz conditions, (C4)–(C6) denote boundedness conditions, and (C7)–(C8) are technical conditions, which are closely related to the considered problem. The conditions (C1)–(C6) are widely recognized in the literature devoted to the existence and uniqueness of the solutions of various types of differential equations.
(C1)
The function f : I × R R is measurable in t I for all y R and there is a constant L f > 0 , such that | f ( t , y 1 ) f ( t , y 2 ) | L f | y 1 y 2 | for all t I and all y 1 , y 2 R ;
(C2)
The function g i : [ t i , s i ] × R R is measurable with respect to t [ t i , s i ] , i = 1 , 2 , , m for all y R and there are constants L g i > 0 , such that | g i ( t , y 1 ) g i ( t , y 2 ) | L g i | y 1 y 2 | for all t [ t i , s i ] , i = 1 , 2 , , m and all y 1 , y 2 R ;
(C3)
Suppose f 1 : I × R × R R is continuous with the respect to t I , for all ( y 1 , z 1 ) R × R and there is a constant L f 1 > 0 , such that | f 1 ( t , y 1 , z 1 ) f 1 ( t , y 2 , z 2 ) | L f 1 | y 1 y 2 | + | z 1 z 2 | for all t I and all y 1 , y 2 , z 1 , z 2 R ;
(C4)
Suppose f : I × R R is measurable with respect to t I for all y R and continuous with respect to y R for all t I . Additionally, there exists a bounded measurable function a f : I R and a constant b f , such that | f ( t , y ) | | a f ( t ) | + b f | y | a f * + b f | y | , where a f * = sup t I J 0 , t 1 α 1 | a f ( t ) | ;
(C5)
Suppose g : I × R R is measurable with respect to t I for all y R , continuous with respect to y R for all t I and g i : [ t i , s i ] × R R is measurable in t [ t i , s i ] , i = 1 , 2 , , m for all y R and continuous with respect to y R for all t [ t i , s i ] , i = 1 , 2 , , m . Additionally, there exists a bounded measurable function a g : I R and a real number b g and there are bounded measurable functions a g i : [ t i , s i ] R and a constant b g i , such that | g ( t , y ) | | a g ( t ) | + b g | y | a g * + b g | y | , where a g * = sup t I J 0 , t α 2 | a g ( t ) | and | g i ( t , y ) | | a g i ( t ) | + b g i | y | for all i = 1 , 2 , , m ;
(C6)
Suppose f 1 : I × R × R R is measurable with respect to t I for all ( y 1 , y 2 ) R × R and continuous with respect to y 1 , y 2 for all t I . Additionally, there exists a bounded measurable function a f 1 : I R and a constant b f 1 > 0 , such that | f 1 ( t , y 1 , y 2 ) | | a f 1 ( t ) | + b f 1 ( | y 1 | + | y 2 | ) , for all t I , y 1 , y 2 R , where sup t I 0 t | a f 1 ( s ) | d s N ;
(C7)
The number r 1 is a positive solution of the equation
a f * + b f a g * Γ ( 2 α 1 ) 1 r 1 + b f b g r 1 2 Γ ( 2 α 1 ) Γ ( α 1 + α 2 α 3 + 1 ) = 0 ;
(C8)
It holds b f 1 < 1 .
Employing the Banach contraction principle and the Schauder fixed-point theorem, we will present two theorems concerning the existence and uniqueness of the solutions of the non-instantaneous impulsive fractional integral problem Equation (1).
Theorem 1. 
Let (C1)–(C3) be satisfied and let 0 < ρ < 1 , where
ρ = max 1 i m L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β + L f s i α 1 Γ ( α 1 + 1 ) β + ( L f 1 ) β .
Then, the Equation (1) has a unique solution.
Proof. 
We define the operator R : P C ( I : R ) P C ( I : R ) by
( R w ) ( t ) : = w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) , 0 < t t 1 , a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) , t i < t s i , i = 1 , 2 , , m , a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) , s i < t t i + 1 , i = 1 , 2 , , m ,
where
v ( t ) = J 0 , t 1 α 1 f t , v ( t ) · 0 t ( t θ ) α 2 1 Γ ( α 2 ) g θ , 0 θ ( θ ν ) ( α 1 α 3 ) 1 Γ ( α 1 α 3 ) v ( ν ) d ν d θ .
It is obvious that the mapping R is well defined. Hence, we have the following cases:
Case 1. 
Let w 1 , w 2 P C ( I : R ) and 0 < t t 1 . Then,
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | = | w 0 0 1 f 1 s , w 1 ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) w 0 + 0 1 f 1 s , w 2 ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | 0 1 | f 1 s , w 1 ( s ) , J 0 , t α 1 α 4 v ( s ) h s , w 2 ( s ) , J 0 , t α 1 α 4 v ( s ) | d s 0 1 L f 1 | w 1 ( s ) w 2 ( s ) | d s L f 1 w 1 w 2 P C .
Now, since
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | β ( L f 1 ) β w 1 w 2 P β ,
we obtain
R w 1 R w 2 P β L f 1 β w 1 w 2 P β .
Case 2. 
Let w 1 , w 2 P C ( I : R ) and t i < t s i , i = 1 , 2 , , m . Since
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | 1 Γ ( α 3 ) t i t ( t s ) α 3 1 | g i ( s , w 1 ( s ) ) g i ( s , w 2 ( s ) ) | d s + 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 1 | f ( s , w 1 ( s ) ) f ( s , w 2 ( s ) ) | d s L g i Γ ( α 3 ) t i t ( t s ) α 3 1 | w 1 ( s ) w 2 ( s ) | d s + L f Γ ( α 1 ) 0 s i ( s i s ) α 1 1 | w 1 ( s ) w 2 ( s ) | d s L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) + L f s i α 1 Γ ( α 1 + 1 ) w 1 w 2 P C ,
we obtain
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | β L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β + L f s i α 1 Γ ( α 1 + 1 ) β w 1 w 2 P β .
Hence,
R w 1 R w 2 L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β + L f s i α 1 Γ ( α 1 + 1 ) β w 1 w 2 P β .
Case 3. 
Let w 1 , w 2 P C ( I : R ) and s i < t t i + 1 , i = 1 , 2 , , m . Since
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | 1 Γ ( α 3 ) t i t ( t s ) α 3 1 | g i ( s , w 1 ( s ) ) g i ( s , w 2 ( s ) ) | d s + 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 1 | f ( s , w 1 ( s ) ) f ( s , w 2 ( s ) ) | d s + 0 1 | f 1 s , w 1 ( s ) , J 0 , t α 1 α 4 v ( s ) f 1 s , w 2 ( s ) , J 0 , t α 1 α 4 v ( s ) | d s L g i Γ ( α 3 ) t i t ( t s ) α 3 1 | w 1 ( s ) w 2 ( s ) | d s + L f Γ ( α 1 ) 0 s i ( s i s ) α 1 1 | w 1 ( s ) w 2 ( s ) | d s + 0 1 L f 1 | w 1 ( s ) w 2 ( s ) | d s L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) + L f s i α 1 Γ ( α 1 + 1 ) + L f 1 w 1 w 2 P C ,
we have
| ( R w 1 ) ( t ) ( R w 2 ) ( t ) | β L g i ( s i t i ) α 3 Γ ( α 3 ) β + L f s i α 1 Γ ( α 1 + 1 ) β + L f 1 β w 1 w 2 P β ,
which implies
R w 1 R w 2 P β L g i ( s i t i ) α 3 Γ ( α 3 ) β + L f s i α 1 Γ ( α 1 + 1 ) β + L f 1 β w 1 w 2 P β .
Let ρ : = max 1 i m { L g i ( s i t i ) α 3 Γ ( α 3 ) β + L f s i α 1 Γ ( α 1 + 1 ) β + L f 1 β } . Accordingly, by cases 1, 2, and 3, we have
R w 1 R w 2 P β ρ w 1 w 2 P β ,
which shows that R is a contraction mapping on P C ( I : R ) . Therefore, by the Banach contraction principle, we can prove that Equation (1) has a unique solution. □
Theorem 2. 
Let (C4)–(C8) be fulfilled. Then, the Equation (1) possesses at least one solution.
Proof. 
Define the closed ball B r 2 by
B r 2 : = { w P C ( I : R ) : w r 2 } ,
where
max 1 i m | a g i ( s ) | = a g i * * ,
max 1 i m | a f ( s ) | = a f * *
and
r 2 : = max 1 i m { | w 0 | + N + b f 1 r 1 Γ ( α 1 α 4 + 1 ) + r 1 Γ ( α 1 + 1 ) 1 b f 1 , | a | + ( s i t i ) α 3 a g i * * Γ ( α 3 + 1 ) + s i α 1 a f * * Γ ( α 1 + 1 ) 1 ( s i t i ) α 3 b g i Γ ( α 3 + 1 ) s i α 1 b f Γ ( α 1 + 1 ) , | a | + ( s i t i ) α 3 a g i * * Γ ( α 3 + 1 ) + s i α 1 a f * * Γ ( α 1 + 1 ) + N + b f 1 r 1 Γ ( α 1 α 4 + 1 ) + r 1 Γ ( α 1 + 1 ) 1 b f 1 ( s i t i ) α 3 b g i Γ ( α 3 + 1 ) s i α 1 b f Γ ( α 1 + 1 ) } .
We claim that there exists r 1 > 0 , such that | v ( t ) | r 1 for all t I . To calculate r 1 , we have
| v ( t ) | = J 0 , t 1 α 1 | f t , v ( t ) · J 0 , t α 2 g t , J 0 , t α 1 α 3 v ( t ) | J 0 , t 1 α 1 | a f ( t ) | + b f | v ( t ) · J 0 , t α 2 g t , J 0 , t α 1 α 3 v ( t ) | J 0 , t 1 α 1 | a f ( t ) | + b f v ( a g * + b g J 0 , t α 1 + α 2 α 3 v ) J 0 , t 1 α 1 | a f ( t ) | + b f r 1 a g * + r 1 b g Γ ( α 1 + α 2 α 3 + 1 ) ,
so
r 1 = a f * + 1 Γ ( 2 α 1 ) r 1 b f a g * + b f b g r 1 2 Γ ( α 1 + α 2 α 3 + 1 ) .
We define the operator S by
( S w ) ( t ) : = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) , t i < t s i , i = 1 , 2 , , m , w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) , 0 < t t 1 , a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) , s i < t t i + 1 , i = 1 , 2 , , m ,
where
v ( t ) : = J 0 , t 1 α 1 f t , v ( t ) · 0 t ( t θ ) α 2 1 Γ ( α 2 ) g θ , 0 θ ( θ ν ) ( α 1 α 3 ) 1 Γ ( α 1 α 3 ) v ( ν ) d ν d θ .
Let w B r 2 . First, we are going to prove that S : B r 2 B r 2 and { S w } are uniformly bounded on B r 2 . We use (C4)–(C6) when considering the following cases:
Case 1. 
Let 0 < t t 1 . Then, we have
| S w ( t ) | = | w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) | | w 0 | + 0 1 | f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) | d s + J 0 , t α 1 | v ( t ) | | w 0 | + 0 1 | a f 1 ( s ) | + b f 1 | w ( s ) | + J 0 , t α 1 α 4 | v ( s ) | d s + J 0 , t α 1 | v ( t ) | | w 0 | + 0 1 | a f 1 ( s ) | d s + b f 1 0 1 | w ( s ) | d s + b f 1 v Γ ( α 1 α 4 + 1 ) 0 1 d s + v Γ ( α 1 + 1 ) | w 0 | + N + b f 1 r 2 + b f 1 r 1 Γ ( α 1 α 4 + 1 ) + r 1 Γ ( α 1 + 1 ) r 2 .
Case 2. 
Let t i < t s i , i = 1 , 2 , , m . Then,
| S w ( t ) | = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) | a | + t i t ( t s ) α 3 1 Γ ( α 3 ) | g i ( s , w ( s ) ) | d s + 0 s i ( s i s ) α 1 1 Γ ( α 1 ) | f ( s , w ( s ) ) | d s | a | + t i t ( t s ) α 3 1 Γ ( α 3 ) | a g i ( s ) | + b g i | w ( s ) | d s + 0 s i ( s i s ) α 1 1 Γ ( α 1 ) | a f ( s ) | + b f | w ( s ) | d s | a | + ( s i t i ) α 3 Γ ( α 3 + 1 ) ( a g i * * + b g i r 2 ) + s i α 1 Γ ( α 1 + 1 ) a f * * + b f r 2 r 2 ,
where max 1 i m | a g i ( s ) | = a g i * * and max 1 i m | a f ( s ) | = a f * * .
Case 3. 
Let s i < t t i + 1 , for i = 1 , 2 , , m . Then, we have
| S w ( t ) | = | a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) | | a | + t i t ( t s ) α 3 1 Γ ( α 3 ) | g i ( s , w ( s ) ) | d s + 0 s i ( s i s ) α 1 1 Γ ( α 1 ) | f ( s i , w ( s i ) ) | d s + 0 1 | f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) | d s + J s i , t α 1 | v ( t ) | | a | + ( s i t i ) α 3 Γ ( α 3 + 1 ) ( a g i * * + b g i r 2 ) + s i α 1 Γ ( α 1 + 1 ) a f * * + b f r 2 + 0 1 | a f 1 ( s ) | d s + b f 1 0 1 | w ( s ) | d s + b f 1 0 1 J 0 , t α 1 α 4 | v ( s ) | d s + J s i , t α 1 | v ( t ) | | a | + ( s i t i ) α 3 Γ ( α 3 + 1 ) ( a g i * * + b g i r 2 ) + s i α 1 Γ ( α 1 + 1 ) a f * * + b f r 2 + N + b f 1 r 2 + b f 1 r 1 Γ ( α 1 α 4 + 1 ) + r 1 Γ ( α 1 + 1 ) r 2 .
We now prove that S is continuous. Let ( w n ) B r 2 and w n w , when n . By the assumptions in the statement of the theorem, we note that the conditions for Lebesgue dominated convergence theorem in the fractional calculations in the corresponding domains are fulfilled. Then, we have the following cases:
Case (a). 
Let 0 < t t 1 . Then,
S w n ( t ) = w 0 0 1 f 1 s , w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) .
By using the Lebesgue dominated convergence theorem, we obtain
lim n S w n ( t ) = lim n w 0 0 1 f 1 s , w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) = w 0 0 1 f 1 s , lim n w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) = w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) = S w ( t ) .
Case (b). 
Let t i < t s i , i = 1 , 2 , , m . So,
S w n ( t ) = a + J t i , t α 3 g i ( t , w n ( t ) ) J 0 , s i α 1 f ( s i , w n ( s i ) ) .
Therefore, it follows that
lim n S w n ( t ) = lim n a + J t i , t α 3 g i ( t , w n ( t ) ) J 0 , s i α 1 f ( s i , w n ( s i ) ) = a + 1 Γ ( α 3 ) t i t ( t s ) α 3 1 g i ( s , lim n w n ( s ) ) d s 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 f ( s , lim n w n ( s ) ) d s = a + 1 Γ ( α 3 ) t i t ( t s ) α 3 1 g i ( s , w ( s ) ) d s 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 f ( s , w ( s ) ) d s = S w ( t ) .
Case (c). 
Let s i < t t i + 1 , i = 1 , 2 , , m . Since
S w n ( t ) = a + J t i , s i α 3 g i ( s i , w n ( s i ) ) J 0 , s i α 1 f ( s i , w n ( s i ) ) 0 1 f 1 s , w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) ,
we acquire
lim n S w n ( t ) = lim n ( a + J t i , s i α 3 g i ( s i , w n ( s i ) ) J 0 , s i α 1 f ( s i , w n ( s i ) ) 0 1 f 1 s , w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) ) = a + J t i , s i α 3 g i ( s i , lim n w n ( s i ) ) J 0 , s i α 1 f ( s i , lim n w n ( s i ) ) 0 1 f 1 s , lim n w n ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) = a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) = S w ( t ) .
Hence, by cases (a), (b), and (c), we show that the operator S is continuous.
Next, we are going to prove that the family { S w } is equicontinuous and the operator S is relatively compact. Let w B r 2 and let t , t I , t > t and | t t | δ , for δ > 0 . Hence, we need to consider the following possible cases:
Case (i). 
Let t , t ( 0 , t 1 ] . Then, we have
| S w ( t ) S w ( t ) | = | w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v ( t ) w 0 + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s = 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s + 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s = t t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s + 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s .
By t t , we obtain that S w ( t ) S w ( t ) .
Case (ii). 
Let t , t ( t i , s i ] , i = 1 , 2 , , m . Then,
| S w ( t ) S w ( t ) | = | a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) a + J t i , t α 3 g i ( t , w ( t ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) | 1 Γ ( α 3 ) t i t ( t s ) α 3 1 | g i ( s , w ( s ) ) | d s 1 Γ ( α 3 ) t i t ( t s ) α 3 1 | g i ( s , w ( s ) ) | d s = 1 Γ ( α 3 ) t t ( t s ) α 3 1 | g i ( s , w ( s ) ) | d s 0 as t t .
Case (iii). 
Let t , t ( s i , t i + 1 ] , i = 1 , 2 , , m . Hence, we have
| S w ( t ) S w ( t ) | = | a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v ( t ) a J t i , s i α 3 g i ( s i , w ( s i ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s = t t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s + 0 t | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | | f s , v ( s ) · J 0 , t α 2 g s , J 0 , s α 1 α 3 v ( s ) | d s .
Keeping in mind that f is a continuous function with respect to w for all t I (by (C4)), we find that S w ( t ) S w ( t ) when t t .
Finally, by applying the Schauder fixed-point theorem, we conclude the existence of at least one solution of Equation (1). □

3. β –Ulam–Hyers Stability of Solutions

Let 0 < β < 1 and ε > 0 . Let us consider the following inequality:
| w ( t ) w 0 + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | ε , 0 < t t 1 | w ( t ) a J t i , s i α 3 g i ( s i , w ( s i ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | ε , s i < t t i + 1 , i = 0 , 1 , 2 , m | w ( t ) a J t i , t α 3 g i ( t , w ( t ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) | ε , t i < t s i , i = 1 , 2 , , m ,
where v ( t ) is defined as in Equation (1). The definition provided bears a resemblance to the definition of solution stability in the context of β –Ulam–Hyers, as outlined in [27].
Definition 1. 
The Equation (1) is β–Ulam–Hyers-stable if there exists a positive real number C α 1 , α 2 , α 3 , f , g i , f 1 > 0 , such that for each ε > 0 and for each solution w P C ( I : R ) of the inequality (2) there exists a solution u P C ( I : R ) of (1) whereby
| u ( t ) w ( t ) | β C α 1 , α 2 , α 3 , f , g i , f 1 ε β
for all t I .
In terms of Ulam’s stability, the mapping and homomorphism are defined on the respective intervals. Therefore, we consider several cases.
Remark 1. 
It is worth noting that w P C ( I : R ) is a solution of Equation (2) and is equivalent to saying that there exists a real number K > 0 , such that
(i) 
| K | ε ;
(ii) 
w ( t ) = w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v K ( t ) , 0 < t t 1 ;
(iii) 
w ( t ) = a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s
+ J s i , t α 1 v K ( t ) , s i < t t i + 1 ;
(iv) 
w ( t ) = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) + K , t i < t s i , i = 1 , 2 , , m , where
v K ( t ) = J 0 , t 1 α 1 f t , v ( t ) · 0 t ( t θ ) α 2 1 Γ ( α 2 ) g θ , 0 θ ( θ ν ) α 1 α 3 1 Γ ( α 1 α 3 ) v ( ν ) d ν d θ + K .
We continue by stating the following result:
Proposition 1. 
Let w ( t ) be a solution of inequality Equation (2). Then, w ( t ) is a solution of the following inequality:
| w ( t ) w 0 + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | ε t 1 α 1 α 1 π , 0 < t t 1 | w ( t ) a J t i , t α 1 g i ( t , w ( t ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) | ε , t i < t s i , i = 1 , 2 , , m | w ( t ) a J t i , s i α 3 g i ( s i , w ( s i ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | ε ( 1 + ( t i + 1 s i ) α 1 α 1 π ) , s i < t t i + 1 , i = 1 , 2 , , m .
Proof. 
We have
w ( t ) = w 0 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J 0 , t α 1 v K ( t ) for 0 < t t 1 ,
w ( t ) = a + J t i , t α 3 g i ( t , w ( t ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) + K for t i < t s i , i = 1 , 2 , , m
and
w ( t ) = a + J t i , s i α 3 g i ( s i , w ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s + J s i , t α 1 v K ( t ) for s i < t t i + 1 , i = 1 , 2 , , m ,
where
v K ( t ) = J 0 , t 1 α 1 f t , v ( t ) · 0 t ( t θ ) α 2 1 Γ ( α 2 ) g θ , 0 θ ( θ ν ) α 1 α 3 1 Γ ( α 1 α 3 ) v ( ν ) d ν d θ + K .
Let 0 < t t 1 . Then, we have
w ( t ) w 0 + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) = J 0 , t α 1 1 Γ ( 1 α 1 ) 0 t ( t s ) α 1 | K | d s J 0 , t α 1 | K | Γ ( 1 α 1 ) = 1 Γ ( 1 α 1 ) · | K | Γ ( α 1 ) · 0 t ( t s ) α 1 1 d s | K | t 1 α 1 sin ( π α 1 ) π α 1 ε t 1 α 1 π α 1 .
For t i < t s i , i = 1 , 2 , , m , we obtain
| w ( t ) a J t i , t α 3 g i ( t , w ( t ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) | | K | ε .
Finally, for s i < t t i + 1 , i = 1 , 2 , , m , we acquire
| w ( t ) a J t i , s i α 3 g i ( s i , w ( s i ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | | K | + J s i , t α 1 J s i , t 1 α 1 | K | = | K | + J s i , t α 1 1 Γ ( 1 α 1 ) s i t ( t s ) α 1 | K | d s | K | + 1 Γ ( 1 α 1 ) · | K | Γ ( α 1 ) s i t ( t s ) α 1 1 d s | K | + 1 Γ ( 1 α 1 ) · | K | Γ ( α 1 ) · 1 α 1 ( t i + 1 s i ) α 1 | K | ( 1 + ( t i + 1 s i ) α 1 · sin ( π α 1 ) α 1 π ) ε ( 1 + ( t i + 1 s i ) α 1 α 1 π ) .
The proof is completed. □
We now present the following β –Ulam–Hyers stability result, which is one of the main results of this paper.
Theorem 3. 
Let the conditions from Theorem 1 be satisfied. Suppose that
1 L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β L f s i α 1 Γ ( α 1 + 1 ) β L f 1 β > 0
and 1 C 1 β L f 1 β > 0 hold, where | u 0 w 0 | β C 1 u w P β and w ( 0 ) = u ( 0 ) . Then, Equation (1) is β–Ulam–Hyers-stable with respect to ε.
Proof. 
Let us denote by w ( · ) the unique solution of Equation (1). Let u P C ( I : R ) be a solution of the inequality Equation (2). By virtue of Proposition 1, we conclude the following:
(1).
For 0 < t t 1 ,
| u ( t ) u ( 0 ) J 0 , t α 1 v ( t ) | ε t 1 α 1 α 1 π .
(2).
For t i < t s i , i = 1 , 2 , , m ,
| u ( t ) a J t i , t α 3 g i ( t , u ( t ) ) + J 0 , s i α 1 f ( s i , u ( s i ) ) | ε .
(3).
For s i < t t i + 1 , i = 1 , 2 , , m ,
| u ( t ) a J t i , t α 3 g i ( t , u ( t ) ) J 0 , s i α 1 f ( s i , u ( s i ) ) + 0 1 f 1 s , u ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | ε 1 + ( t i + 1 s i ) α 1 α 1 π .
So, we need to verify the following cases:
Case 1. 
Let 0 < t t 1 . Then,
| u ( t ) w ( t ) | β = | u ( t ) w 0 + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | β | u ( t ) u 0 + 0 1 f 1 s , u ( s ) , J 0 , t α 1 α 4 v ( s ) d s J 0 , t α 1 v ( t ) | β + | u 0 w 0 | + | 0 1 f 1 s , u ( s ) , J 0 , t α 1 α 4 v ( s ) f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s | β ε t 1 α 1 α 1 π β + C 1 + L f 1 0 1 d s β u w P β = ε t 1 α 1 α 1 π β + C 1 + L f 1 β u w P β ,
where | u 0 w 0 | β C 1 u w P β with C 1 ( 0 , 1 ) . Accordingly, we obtain
1 C 1 β L f 1 β | u ( t ) w ( t ) | β t 1 α 1 α 1 π β · ε β ,
which implies
| u ( t ) w ( t ) | β C α 1 , α 2 , α 3 , f , g i , f 1 · ε β ,
where
C α 1 , α 2 , α 3 , f , g i , f 1 = t 1 α 1 α 1 π β 1 C 1 β L f 1 β .
Case 2. 
Let t i < t s i , i = 1 , 2 , , m . Since
| u ( t ) w ( t ) | β = | u ( t ) a J t i , t α 3 g i ( t , w ( t ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) | β | u ( t ) a J 0 , t α 3 g i ( t , u ( t ) ) + J 0 , s i α 1 f ( s i , u ( s i ) ) | β + | J t i , t α 3 g i ( t , u ( t ) ) J t i , t α 3 g i ( t , w ( t ) ) | β + | J 0 , s i α 1 f ( s i , u ( s i ) ) J 0 , s i α 1 f ( s i , u ( s i ) ) | β ε β + 1 Γ ( α 3 ) t i t ( t s ) α 3 1 | g i ( t , u ( t ) ) g i ( t , w ( t ) ) | d s β + 1 Γ ( α 1 ) 0 s i ( s i s ) α 1 1 | f ( s i , u ( s i ) ) f ( s i , w ( s i ) ) | d s β ε β + L f s i α 1 Γ ( α 1 + 1 ) β + L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β u w P β ,
we derive
1 L f s i α 1 Γ ( α 1 + 1 ) β L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β u w P β ε β .
Hence, we obtain
| u ( t ) w ( t ) | β C α 1 , α 2 , α 3 , f , g i , f 1 · ε β , t i < t s i , i = 1 , 2 , , m ,
where
C α 1 , α 2 , α 3 , f , g i , f 1 = 1 1 L f s i α 1 Γ ( α 1 + 1 ) β L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β .
Case 3. 
Let s i < t t i + 1 , i = 1 , 2 , , m . Then,
| u ( t ) w ( t ) | β = | u ( t ) a J t i , s i α 3 g i ( s i , w ( s i ) ) + J 0 , s i α 1 f ( s i , w ( s i ) ) + 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | β | u ( t ) a J t i , s i α 3 g i ( s i , u ( s i ) ) + J 0 , s i α 1 f ( s i , u ( s i ) ) + 0 1 f 1 s , u ( s ) , J 0 , t α 1 α 4 v ( s ) d s J s i , t α 1 v ( t ) | β + | J t i , s i α 3 g i ( s i , u ( s i ) ) J t i , s i α 3 g i ( s i , w ( s i ) ) | β + | J 0 , s i α 1 f ( s i , u ( s i ) ) J 0 , s i α 1 f ( s i , w ( s i ) ) | β + | 0 1 f 1 s , u ( s ) , J 0 , t α 1 α 4 v ( s ) d s 0 1 f 1 s , w ( s ) , J 0 , t α 1 α 4 v ( s ) d s | β ε β 1 + ( t i + 1 s i ) α 1 α 1 π β + ( L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β + L f s i α 1 Γ ( α 1 + 1 ) β + L f 1 β ) u w P β .
Hence,
1 L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β L f s i α 1 Γ ( α 1 + 1 ) β L f 1 β | u ( t ) w ( t ) | β ε β 1 + ( t i + s i ) α 1 α 1 π β ,
implying that
| u ( t ) w ( t ) | β C α 1 , α 2 , α 3 , f , g i , f 1 · ε β , s i < t t i + 1 , i = 1 , 2 , , m ,
where
C α 1 , α 2 , α 3 , f , g i , f 1 = 1 + ( t i + 1 s i ) α 1 α 1 π β 1 L g i ( s i t i ) α 3 Γ ( α 3 + 1 ) β L f s i α 1 Γ ( α 1 + 1 ) β L f 1 β .
Finally, keeping in mind the considered three cases, we showed in each case the existence of a constant C α 1 , α 2 , α 3 , f , g i , f 1 > 0 , such that | u ( t ) w ( t ) | β C α 1 , α 2 , α 3 , f , g i , f 1 · ε β for all t I , i.e., Equation (1) is β –Ulam–Hyers-stable, with respect to ε . □

4. Nontrivial Example of Application of Theorem 1

In this section, we give a nontrivial example illustrating the application of Theorem 1.
Example 1. 
Let us consider the following equation:
w ( t ) = 1 0 1 s 2 3 + 1 2 w ( s ) + J 0 , t 1 2 v ( s ) d s + J 0 , t 3 4 v ( t ) , 0 < t 1 2 , w ( t ) = a + 1 Γ ( 1 5 ) 1 2 t ( t s ) 4 5 | w ( s ) | 16 ( 1 + | w ( s ) | ) d s 1 Γ ( 3 4 ) 0 1 ( 1 s ) 1 4 arctan ( s 2 + w ( s ) ) 8 + e s + s 2 d s , 1 2 < t 1 , w ( 0 ) = 1 0 1 s 2 3 + 1 2 w ( s ) + J 0 , t 1 2 v ( s ) d s ,
where
v ( s ) = J 0 , s 1 4 ( 1 8 + e s + s 2 arctan ( s 2 + v ( s ) · 0 s ( s θ ) 1 2 Γ 1 2 · ( θ 3 + 1 2 0 θ ( θ ν ) 9 20 Γ 11 20 v ( ν ) d ν ) d θ ) )
and a R .
  • Here, we set
f ( t , s ) = 1 8 + e t + t 2 arctan ( t 2 + s ) ,
g 1 ( t , s ) = | s | 16 ( 1 + | s | ) ,
g ( t , s ) = t 3 + 1 2 s ,
f 1 ( t , s , p ) = t 2 3 + 1 2 ( s + p ) ,
α 1 = 3 4 , α 2 = 1 2 , α 3 = 1 5 , α 4 = 1 2 , w 0 = 1 , and a R . Moreover, I = [ 0 , 1 ] , where 0 = t 0 = s 0 < t 1 < s 1 = 1 , with t 1 = 1 2 and s 1 = 1 .
It is clear that (C1)–(C3) hold, with L f = 1 9 , L g 1 = 1 16 , L g = 1 2 , and L f 1 = 1 2 . Since
ρ = 1 16 · 1 2 1 5 Γ 6 5 β + 1 9 · 1 3 4 Γ 7 4 β + ( 1 2 ) β < 1 ,
for β > 0.688 , by Theorem 1, the considered non-instantaneous impulsive fractional integral problem has a unique solution.
Let u P C ( I : R ) be a solution of inequality Equation (2). Put ε = 1 . Then, there exists a constant K R , such that | K | ε = 1 and
u ( t ) = 1 0 1 s 2 3 + 1 2 u ( s ) + J 0 , t 1 2 v ( s ) d s + J 0 , t 3 4 v K ( t ) , 0 < t 1 2 ,
where
v K ( s ) = J 0 , s 1 4 ( 1 8 + e s + s 2 arctan ( s 2 + v ( s ) · 0 s ( t θ ) 1 2 Γ 1 2 · ( θ 3 + 1 2 0 θ ( θ ν ) 9 20 Γ 11 20 v ( ν ) d ν ) d θ ) + K )
and
u ( t ) = a + 1 Γ 1 5 1 2 t ( t s ) 4 5 | u ( s ) | 16 ( 1 + | u ( s ) | ) d s 1 Γ 3 4 0 1 3 4 s 1 4 arctan ( s 2 + u ( s ) ) 8 + e s + s 2 d s + K , 1 2 < t 1 .
From the previous discussion,
w ( t ) = 1 0 1 s 2 3 + 1 2 w ( s ) + J 0 , t 1 2 v ( s ) d s + J 0 , t 3 4 v ( t ) , 0 < t 1 2 , w ( t ) = a + 1 Γ ( 1 5 ) 1 2 t ( t s ) 4 5 | w ( s ) | 16 ( 1 + | w ( s ) | ) d s 1 Γ ( 3 4 ) 0 1 ( 1 s ) 1 4 arctan ( s 2 + w ( s ) ) 8 + e s + s 2 d s , 1 2 < t 1 w ( 0 ) = u ( 0 ) ,
has a unique solution for β = 9 4 . Now, for 0 < t 1 2 , we obtain
| u ( t ) w ( t ) | 9 4 1 2 3 4 3 π 4 9 4 1 1 0 . 9 9 4 1 2 9 4 · ε 9 4 = 54.3934 · 1 9 4 = 54.3934 .
For 1 2 < t 1 , we have
| u ( t ) w ( t ) | 9 4 1 1 1 9 · 1 3 4 Γ 7 4 9 4 1 16 · 1 2 1 5 Γ 6 5 9 4 · ε 9 4 = 1.0105 · 1 9 4 = 1.0105 .
Hence, we obtain
| u ( t ) w ( t ) | 9 4 54.3934 = 54.3934 · 1 9 4 ,
for all t [ 0 , 1 ] , implying that the considered problem is 9 4 –Ulam–Hyers-stable with respect to ε = 1 .
Example 2. 
Consider the following equation:
w ( t ) = 1 0 1 s 3 + 1 4 w ( s ) + J 0 , t 1 6 v ( s ) d s + J 0 , t 2 3 v ( t ) , 0 < t 1 2 , w ( t ) = a + 1 Γ ( 1 3 ) 1 2 t ( t s ) 1 2 w ( s ) 2 8 ( 1 + w ( s ) 2 ) d s 1 Γ ( 2 3 ) 0 1 1 3 ( 1 s ) 1 3 e s ( sin ( s 2 + w ( s ) ) ) d s , 1 2 < t 1 , w ( 0 ) = 1 0 1 s 3 + 1 4 w ( s ) + J 0 , t 1 6 v ( s ) d s ,
where
v ( s ) = J 0 , s 1 3 ( 1 3 e s sin s 2 + v ( s ) · 0 s ( s θ ) 4 5 Γ 1 5 · θ 3 + + 1 4 0 θ ( θ ν ) 2 3 Γ 1 3 v ( ν ) d ν d θ )
and a R .
  • Now,
f ( t , s ) = 1 3 e t sin ( t 2 + s ) ,
g 1 ( t , s ) = f r a c s 2 8 ( 1 + s 2 ) ,
g ( t , s ) = t 3 + 1 4 s ,
f 1 ( t , s , p ) = t 3 + 1 4 ( s + p ) .
Also, α 1 = 2 3 , α 2 = 1 5 , α 3 = 1 3 , α 4 = 1 2 , w 0 = 1 , I = [ 0 , 1 ] , and 0 = t 0 = s 0 < t 1 < s 1 = 1 , with t 1 = 1 2 and s 1 = 1 . We have L f = 1 3 , L g 1 = 1 8 , L g = 1 4 , and L f 1 = 1 4 , hence the conditions (C1)–(C3) hold and for β > 0.717 , ρ < 1 . So, the conditions of Theorem 1 are fulfilled, hence the considered problem has a unique solution.
Like in the previous example, the considered problem has a unique solution for β = 3 2 . For 0 < t 1 2 , we have
| u ( t ) w ( t ) | 3 2 1 2 2 3 3 π 4 3 2 1 1 0 . 9 3 2 1 4 3 2 · ε 3 2 = 7.7867 · 1 3 2 = 7.7867 .
For 1 2 < t 1 , we have
| u ( t ) w ( t ) | 3 2 1 1 1 3 · 1 1 3 Γ 5 3 3 2 1 8 · 1 2 1 3 Γ 1 3 3 2 · ε 3 2 = 1.3012 · 1 3 2 = 1.3012 .
Hence, we obtain
| u ( t ) w ( t ) | 3 2 7.7867 = 7.7867 · 1 3 2 ,
for all t [ 0 , 1 ] , implying that the considered problem is 3 2 –Ulam–Hyers-stable with respect to ε = 1 .

5. Conclusions

In this research, we focused on investigating a class of non-instantaneous impulsive fractional integral equations. We have examined the existence and uniqueness of solutions and the β –Ulam–Hyers stability of solutions for this class of integral impulsive equations. Finally, we provided examples to illustrate the validity and effectiveness of our results. Our findings have unveiled the phenomenon of β –Ulam–Hyers stability among these solutions, providing a valuable perspective for both researchers and practitioners in their pursuit of dependable and resilient mathematical models in future studies.
It remains an open question whether the imposed conditions on Lipschitz continuity or boundedness can be replaced with weaker ones, like the ones proposed in [47], while still obtaining the same results. This issue warrants further investigation. Additionally, we propose exploring the semigroup approach (see, for example, [48] and references therein), which may provide new tools and perspectives for further generalizations.

Author Contributions

Conceptualization, W.-S.D., M.F., M.K. and D.V.; methodology, W.-S.D., M.F., M.K. and D.V.; software, W.-S.D., M.F., M.K. and D.V.; validation, W.-S.D., M.F., M.K. and D.V.; formal analysis, W.-S.D., M.F., M.K. and D.V.; investigation, W.-S.D., M.F., M.K. and D.V.; writing—original draft preparation, W.-S.D., M.F., M.K. and D.V.; writing—review and editing, W.-S.D., M.F., M.K. and D.V.; visualization, W.-S.D., M.F., M.K. and D.V.; supervision, W.-S.D., M.F., M.K. and D.V.; project administration, W.-S.D., M.F., M.K. and D.V. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is partially supported by Grant No. NSTC 113-2115-M-017-004 of the National Science and Technology Council of the Republic of China. The second author is partially supported by the Slovak Research and Development Agency under the Contract no. APVV-23-0039, and the Slovak Grant Agency VEGA No. 1/0084/23 and No. 2/0062/24. The third author and the fourth author are partially supported by grant 174024 of Ministry of Science and Technological Development, Republic of Serbia and Bilateral project between MANU and SANU.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors wish to express their thanks to the anonymous referees for their valuable suggestions and comments.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bainov, D.; Simeonov, P. Impulsive Differential Equations: Periodic Solutions and Applications; Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Halanay, A.; Wexler, D. Qualitative Theory of Impulse Systems; Mir: Moscow, Russia, 1971. [Google Scholar]
  3. Hristova, S.G.; Bainov, D.D. Integral surfaces for hyperbolic ordinary differential equations with impulse effect. COMPEL 1995, 4, 1–18. [Google Scholar] [CrossRef]
  4. Lakshmikantham, V.; Bainov, D.D.; Simeonov, P.S. Theory of Impulsive Differential Equations; World Scientific Publishing Co. Pte. Ltd.: Singapore, 1989. [Google Scholar]
  5. Li, X.; Martin, B.; Wang, C. Impulsive differential equations: Periodic solutions and applications. Automatica 2015, 52, 173–178. [Google Scholar] [CrossRef]
  6. Samoilenko, A.M.; Perestyuk, N.A. Impulsive Differential Equations; World Scientific: Singapore, 1995. [Google Scholar]
  7. Agarwal, R.P.; Hristova, S.; O’Regan, D. Non-Instantaneous Impulses in Differential Equations; Springer: New York, NY, USA, 2017. [Google Scholar]
  8. Bainov, D.; Simeonov, P.S. Systems with Impulse Effect: Stability, Theory, and Applications; Ellis Horwood: Chichester, UK, 1989. [Google Scholar]
  9. Jung, S.-M. Hyers–Ulam–Rassias Stability of Functional Equations in Nonlinear Analysis; Springer: New York, NY, USA, 2011. [Google Scholar]
  10. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  11. Podlubny, I. Fractional Differential Equations. In Mathematics in Science and Engineering; Academic Press: San Diego, CA, USA, 1999; Volume 198. [Google Scholar]
  12. Tarasov, V.E. Fractional Dynamic: Application of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  13. Hernández, E.; O’Regan, D. On a new class of abstract impulsive differential equations. Proc. Am. Math. Soc. 2013, 141, 1641–1649. [Google Scholar] [CrossRef]
  14. Pierri, M.; O’Regan, D.; Rolnik, V. Existence of solutions for semi-linear abstract differential equations with not instantaneous impulses. Appl. Math. Comput. 2013, 219, 6743–6749. [Google Scholar] [CrossRef]
  15. Ulam, S.M. A Collection of Mathematical Problems; Interscience: New York, NY, USA, 1960. [Google Scholar]
  16. Hyers, D.H.; Isac, G.; Rassias, T.M. Stability of Functional Equations in Several Variables; Birkhäuser: Basel, Switzerland, 1998. [Google Scholar]
  17. Rassias, T.M. On the stability of linear mappings in Banach spaces. Proc. Am. Math. Soc. 1978, 72, 297–300. [Google Scholar] [CrossRef]
  18. András, S.; Mészáros, A.R. Ulam-Hyers stability of dynamic equations on time scales via Picard operators. Appl. Math. Comput. 2013, 219, 4856–4864. [Google Scholar] [CrossRef]
  19. Burger, M.; Ozawa, N.; Thom, A. On Ulam stability. Isr. J. Math. 2013, 193, 109–129. [Google Scholar] [CrossRef]
  20. Du, W.-S. A generalization of Diaz-Margolis’s fixed point theorem and its application to the stability of generalized Volterra integral equations. J. Inequal. Appl. 2015, 2015, 407. [Google Scholar] [CrossRef]
  21. Jung, S.-M.; Kim, T.S.; Lee, K.S. A fixed point approach to the stability of quadratic functional equation. Bull. Korean Math. Soc. 2006, 43, 531–541. [Google Scholar] [CrossRef]
  22. Lungu, N.; Popa, D. Hyers-Ulam stability of a first order partial differential equation. J. Math. Anal. Appl. 2011, 381, 530–537. [Google Scholar]
  23. Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpath. J. Math. 2010, 26, 103–107. [Google Scholar]
  24. Wang, J.; Fečkan, M.; Zhou, Y. Ulam’s type stability of impulsive ordinary differential equations. J. Math. Anal. Appl. 2012, 395, 258–264. [Google Scholar] [CrossRef]
  25. Ciepliński, K. Stability of multi-additive mappings in β-Banach spaces. Nonlinear Anal. TMA 2012, 75, 4205–4212. [Google Scholar] [CrossRef]
  26. Xu, T.Z. On the stability of multi-Jensen mappings in β-normed spaces. Appl. Math. Lett. 2012, 25, 1866–1870. [Google Scholar] [CrossRef]
  27. Yu, X. Existence and β-Ulam-Hyers stability for a class of fractional differential equations with non-instantaneous impulses. Adv. Differ. Equ. 2015, 104, 13. [Google Scholar] [CrossRef]
  28. Kostić, M.; Chaouchi, B.; Du, W.-S.; Velinov, D. Generalized ρ-Almost Periodic Sequences and Applications. Fractal Fract. 2023, 7, 410. [Google Scholar] [CrossRef]
  29. El-Sayed, A.M.A.; Alhmali, A.A.A.; Hamdallah, E.M.A. Analysis of fractional-order quadratic functional integro-differential equation with nonlocal fractional-order integro-differential condition. Axioms 2023, 12, 788. [Google Scholar] [CrossRef]
  30. Aicha, S.; Merad, A. Solvability of nonlinear fractional integro-differential equation with nonlocal condition. Arab. J. Math. Sci. 2020, 29, 172–190. [Google Scholar] [CrossRef]
  31. Subashini, R.; Jothimani, K.; Nisar, K.S.; Ravichandran, C. New results on nonlocal functional integrodifferential equations via Hilfer fractional derivative. Alex. Eng. J. 2020, 59, 2891–2899. [Google Scholar] [CrossRef]
  32. Hong, S.H.; Qiu, Z. Existence of solutions on nth order impulsive integro–differential equations in Banach spaces. Electron. J. Qual. Theory Differ. Equ. 2008, 22, 1–11. [Google Scholar] [CrossRef]
  33. Agarwal, R.P.; Hristova, S.; O’Regan, D. Non-instantaneous impulses in Caputo fractional differential equations. Fract. Calc. Appl. 2017, 20, 595–622. [Google Scholar] [CrossRef]
  34. Shah, D.; Riaz, U.; Zada, A. Exponential and Hyers-Ulam stability of impulsive linear system of first order. Differ. Equ. Appl. 2023, 15, 1–11. [Google Scholar] [CrossRef]
  35. Luo, D.; Luo, Z. Existence and Hyers-Ulam stability results for a class of fractional order delay differential equations with non-instantaneous impulses. Math. Slovaca 2020, 70, 1231–1248. [Google Scholar] [CrossRef]
  36. Abbas, S.; Benchohra, M. Stability results for fractional differential equations with state-dependant delay and not instantaneous impulses. Math. Slovaca 2017, 67, 875–894. [Google Scholar] [CrossRef]
  37. Ibrahim, R.W. Generalized Ulam-Hyers stability for fractional differential equations. Internat. J. Math. 2012, 23, 1250056. [Google Scholar] [CrossRef]
  38. Parthasarathy, C. Existence and Hyers-Ulam stability of nonlinear impulsive differential equations with nonlocal conditions. Electron. J. Math. Anal. Appl. 2016, 4, 106–115. [Google Scholar]
  39. Agarwal, R.P.; Hristova, S. Ulam-Type Stability for a Boundary-Value Problem for Multi-Term Delay Fractional Differential Equations of Caputo Type. Axioms 2022, 11, 742. [Google Scholar] [CrossRef]
  40. Dhayal, R.; Malik, M.; Nisar, K.S. Results on non-instantaneous impulsive ϕ-Caputo fractional differential systems: Stability and controllability. Differ. Equ. Appl. 2024, 16, 113–134. [Google Scholar]
  41. Abdo, M.S.; Thabet, S.T.M.; Ahmad, B. The existence and Ulam-Hyers stability results for ψ-Hilfer fractional integrodifferential equations. J. Pseudo-Differ. Oper. Appl. 2020, 11, 1757–1780. [Google Scholar] [CrossRef]
  42. Luo, D.; Shah, K.; Luo, Z. On the novel Ulam-Hyers stability for a class of nonlinear ψ-Hilfer fractional differential equation with time-varying delays. Mediterr. J. Math. 2019, 16, 112. [Google Scholar] [CrossRef]
  43. Sousa, J.V.d.C.; Capelas de Oliveira, E. On the Ulam-Hyers-Rassias stability for nonlinear fractional differential equations using the ψ-Hilfer operator. J. Fixed Point Theory Appl. 2018, 20, 96. [Google Scholar] [CrossRef]
  44. Agarwal, R.P.; Assolami, A.; Alsaedi, A.; Ahmad, B. Existence Results and Ulam-Hyers Stability for a Fully Coupled System of Nonlinear Sequential Hilfer Fractional Differential Equations and Integro-Multistrip-Multipoint Boundary Conditions. Qual. Theory Dyn. Syst. 2022, 21, 125. [Google Scholar] [CrossRef]
  45. Qian, D.; Li, C.; Agarwal, R.P.; Wong, P.J.Y. Stability analysis of fractional differential system with Riemann-Liouville derivative. Math. Comput. Model. 2010, 52, 861–874. [Google Scholar] [CrossRef]
  46. Alam, M.; Shah, D. Hyers-Ulam stability of coupled implicit fractional integro-differential equations with Riemann-Liouville derivatives. Chaos Solit. Fractals 2021, 150, 111122. [Google Scholar] [CrossRef]
  47. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993. [Google Scholar]
  48. Kukushkin, M.V. Abstract fractional calculus for m-accretive operators. Int. J. Appl. Math. 2021, 34, 1–41. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Du, W.-S.; Fečkan, M.; Kostić, M.; Velinov, D. β–Ulam–Hyers Stability and Existence of Solutions for Non-Instantaneous Impulsive Fractional Integral Equations. Fractal Fract. 2024, 8, 469. https://doi.org/10.3390/fractalfract8080469

AMA Style

Du W-S, Fečkan M, Kostić M, Velinov D. β–Ulam–Hyers Stability and Existence of Solutions for Non-Instantaneous Impulsive Fractional Integral Equations. Fractal and Fractional. 2024; 8(8):469. https://doi.org/10.3390/fractalfract8080469

Chicago/Turabian Style

Du, Wei-Shih, Michal Fečkan, Marko Kostić, and Daniel Velinov. 2024. "β–Ulam–Hyers Stability and Existence of Solutions for Non-Instantaneous Impulsive Fractional Integral Equations" Fractal and Fractional 8, no. 8: 469. https://doi.org/10.3390/fractalfract8080469

Article Metrics

Back to TopTop