Next Article in Journal
Optimizing Variational Quantum Neural Networks Based on Collective Intelligence
Previous Article in Journal
A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them

by
Ravi P. Agarwal
1,2,
Snezhana Hristova
3,* and
Donal O’Regan
4
1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Department of Mathematics and Systems Engineering, Florida Institute of Technology, Melbourne, FL 32901, USA
3
Faculty of Mathematics and Informatics, Plovdiv University, Tzar Asen 24, 4000 Plovdiv, Bulgaria
4
School of Mathematical and Statistical Sciences, University of Galway, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(11), 1626; https://doi.org/10.3390/math12111626
Submission received: 2 May 2024 / Revised: 17 May 2024 / Accepted: 20 May 2024 / Published: 22 May 2024

Abstract

:
Ulam type stability is an important property studied for different types of differential equations. When this type of stability is applied to boundary value problems, there are some misunderstandings in the literature. In connection with this, initially, we give a brief overview of the basic ideas of the application of Ulam type stability to initial value problems. We provide several examples with simulations to illustrate the main points in the application. Then, we focus on some misunderstandings in the application of Ulam stability to boundary value problems. We suggest a new way to avoid these misunderstandings and how to keep the main idea of Ulam type stability when it is applied to boundary value problems of differential equations. We present one possible way to connect both the solutions of the given problem and the solutions of the corresponding inequality. In addition, we provide several examples with simulations to illustrate the ideas for boundary value problems and we also show the necessity of the new way of applying the Ulam type stability. To illustrate the theoretical application of the suggested idea to Ulam type stability, we consider a linear boundary value problem for nonlinear impulsive fractional differential equations with the Caputo fractional derivative with respect to another function and piecewise-constant variable order. We define the Ulam–Hyers stability and obtain sufficient conditions on a finite interval. As partial cases, integral presentations of the solutions of boundary value problems for various types of fractional differential equations are obtained and their Ulam type stability is studied.

1. Introduction

Ulam type stability, sometimes known as Hyers–Ulam stability or Hyers–Ulam–Rassias stability, has been defined and studied by many authors for functional equations (see the formulation given in [1,2] and some recent studies [3,4] and the books [5,6,7]). This type of stability is also applied to various problems for different types of differential equations; it has applications in optimization, biology, economics, etc., and it is a special type of data dependence of solutions. The study in this area has grown and is a central subject in the mathematical analysis area. In this paper, we consider and discuss only its application to differential equations.
The main goal of this paper is to present Ulam type stability for boundary value problems (BVP). We point out the basic misunderstandings in some published paper when this type of stability is considered, and we give one of the possible ways to avoid these misunderstandings and to study Ulam type stability for boundary value problems in a similar way as in the case of initial value problems. A novel approach involving the introduction of a parameter to mitigate these issues is proposed. The presented method is illustrated with specific examples involving a general linear boundary value. To be more precise, we present our ideas on differential equations with a tempered Caputo fractional derivative with respect to another function and with impulses. Note the problem with non-instantaneous impulses is considered in [8] but for the initial value problem and with the generalized proportional Caputo fractional derivative. In this paper, because of the main goal, we study a boundary value problem. Based on some ideas about non-instantaneous impulses in fractional differential equations presented in [8], we consider a differential equation with the more general case of a Caputo type fractional derivative. We present an integral presentation of the solution and prove the existence of the solution of the given BVP. Then, we investigate the Ulam type stability for the studied BVP and sufficient conditions for the Ulam–Hyers stability are derived.
Note the proposed ideas could be applied to study Ulam type stability for various types of boundary value problems for other types of differential equations, such as ordinary differential equations, fractional differential equations with any type of Caputo type fractional derivative, integro-differential fractional equations, and partial differential equations with time-fractional derivatives.
The proposed methodology could be successfully applied to problems in engineering and mathematical physics. For example, in [9], the equilibrium shape (or configuration or state) of the Elastica is described by a boundary value problem for a second-order ordinary differential equation, and since the stability is important in this model, one could investigate the applicability of Ulam type stability. In [10], a generalized two-dimensional space-time reaction–diffusion model is applied to model the calcium dynamics that firmly orchestrates exchanges of calcium flux through intracellular/extracellular sources of calcium to investigate cellular activities and calcium homeostasis. If an appropriate boundary condition is set up in the model, then the suggested ideas could be applied to study Ulam type stability. Therefore, the suggested methodology could broaden the scope of the study and enhance its practical relevance not only to mathematicians but also to numerical analysts, biologists, physicists, and engineers.
The main contributions of the article to the field of stability theory of differential equations, particularly in fractional differential equations, can be summarized as follows:
-
In Section 2, we provide a brief overview of the concepts of Ulam type stability and its applications to various problems such as:
*
Initial value problems to any type of differential equation with various types of derivatives, such as an ordinary derivative or a partial derivative of integer order, any type of fractional derivatives—we briefly make a review and give an algorithm for the application of Ulam type stability;
*
Boundary value problems to any type of differential equation with various types of derivatives—we emphasize the basic misunderstandings of the application of Ulam type stability;
*
A new methodology—we give one of the possible ways to avoid the mentioned misunderstandings by introducing a parameter in the boundary condition;
*
All of the above theoretical explanations are illustrated with appropriate examples.
-
In Section 3, we apply the proposed new methodology to study Ulam type stability for a linear boundary value problem for a differential equation with instantaneous impulses and the Caputo fractional derivative with respect to another function.
-
To summarize and to emphasize some possible applications of the suggested approach, we finish the paper with a conclusion.

2. Brief Overview of the Concepts of Ulam Type Stability

Without loss of generality, we assume the initial time point is zero.
We recall the basic algorithm in the study of Ulam type stability. The steps of the algorithm are described below:
1.
Statement of the problem: Define the main problem of study which consists of two parts:
-
The differential equation D t x ( t ) = F ( t , x ) , t [ 0 , b ] , where x : ( 0 , b ] R n is the unknown function, and D t is an applied derivative. Note b , n 1 , the applied derivative could be an ordinary derivative or a fractional derivative, and the function F could have more than two arguments (it could contain any type of integral or derivative);
-
The initial condition (IVP), the boundary condition (BVP), or both IVP and BVP.
2.
Integral presentation: Obtain the equivalent integral presentation (integral equation) of the solution of the problem defined in Step 1.
3.
Existence of solution: Define an operator based on the integral equation obtained in Step 2, and prove the existence (and uniqueness) of a fixed point (or points), i.e., prove the existence of a solution of the problem defined in Step 1.
4.
Definition of Ulam type stability: Define in an appropriate way the differential inequality deeply connected with the problem given in Step 1. Based on this inequality, define the Ulam type stability (US) of the problem given in Step 1.
5.
Proof of US: Based on Step 2, Step 3, and Step 4, prove the Ulam type stability of the solutions of the studied problem.
We first discuss the application of Ulam type stability to some types of problems for differential equations.

2.1. Initial Value Problems

First, we discuss the application of Ulam type stability to initial value problems of differential equations.
Consider the differential equation (DE)
D t x ( t ) = F ( t , x ) , t [ 0 , b ] ,
with the initial condition
x ( 0 ) = x 0 ,
where x 0 R n , and D t x ( t ) is an ordinary derivative or any type of Caputo fractional derivative.
To focus on Ulam type stability, we skip both Step 2 and Step 3, assuming that for any initial value x 0 R n , the initial value problem (IVP) (1), (2) has a solution on [ 0 , b ] .
Remark 1.
Note that any solution of the IVP (1), (2) depends significantly on the initial value x 0 , i.e., x 0 R n could be considered as a parameter.
Let ε > 0 . Consider the differential inequality (DI)
| D t ν ( t ) F ( t , ν ( t ) ) | ε , t [ 0 , b ] ,
where | · | is an appropriate norm in R n .
We now recall the definition for Ulam type stability based on the classical papers [11,12]:
Definition 1.
The IVP (1), (2) is Ulam–Hyers stable if there exists a real number C F > 0 such that for each ε > 0 and for each solution ν : [ 0 , b ] R n of the DI (3), there exists a solution x : [ 0 , b ] R n of IVP (1), (2) with | ν ( t ) x ( t ) | C F ε for t [ 0 , b ] .
The main idea in the practical application of Ulam type stability in this case is that for any solution ν : [ 0 , b ] R n of the differential inequality (3), we consider a solution x : [ 0 , b ] R n of the initial value problem (1), (2). As it is mentioned in Remark 1, this solution x ( · ) depends significantly on the initial value x 0 . In the application of US, we choose the initial value x 0 = ν ( 0 ) (it exists because of the corresponding Step 3, or the assumption in our case).
Remark 2.
The main point in the application of US to IVPs is the chosen initial value of the given IVP which is equal to the initial value of the chosen arbitrary solution of the corresponding differential inequality.
We illustrate the main idea of the application of Ulam type stability to initial value problems for different types of differential equations.
Example 1
(Ordinary derivative). Consider the simple scalar ordinary differential equation
x = x , t [ 0 , 1 ] ,
with the initial condition
x ( 0 ) = x 0 , x 0 R .
Its solution is x ( t ) = x 0 e t .
Let ε > 0 be an arbitrary number and consider the differential inequality
| ν ( t ) ν ( t ) | < ϵ .
We could choose the solution ν ( t ) = ϵ ( t + 1 ) , t ( 0 , 1 ] .
Now, how do we choose the solution of the differential equation which depends on the initial condition?
Case 1.1. The initial value does not depend on the solution of the inequality. Consider the solution of the given differential equation with the initial value x 0 = 2 , i.e., the solution is x ( t ) = 2 e t , and then, the difference | x ( t ) ν ( t ) | = | 2 e t ϵ ( t + 1 ) | depends significantly on ϵ, and we cannot find a constant C F > 0 independent from ϵ such that | x ( t ) ν ( t ) | C F ϵ (see Figure 1).
Case 1.2. The initial value depends on the solution of the inequality. Consider the solution of the given differential equation with x 0 = ν ( 0 ) = ϵ , i.e., x ( t ) = ϵ e t . Then, the inequality | x ( t ) ν ( t ) | = ϵ | t + 1 e t | ϵ , t ( 0 , 1 ] holds, i.e., according to Definition 1, the solution of the considered initial value problem is Ulam–Hyers stable with a constant C F = 1 .
Therefore, the choice of the initial value x 0 = ν ( 0 ) of the solution of the IVP is very important for US.
Example 2
(Caputo fractional derivative). Consider the Caputo fractional derivative
D t 0.3 x ( t ) = 1 Γ ( 1 0.3 ) 0 t ( t s ) 0.3 x ( s ) d s , α ( 0 , 1 )
as the derivative D t x ( t ) in (1), i.e., consider the scalar fractional differential equation
D t 0.3 x ( t ) = x , t ( 0 , 1 ] ,
with the initial condition
x ( 0 ) = x 0 , x 0 R .
Its solution is x ( t ) = x 0 E 0.3 ( t 0.3 ) , t [ 0 , 1 ] , where E 0.3 ( · ) is the Mittag-Leffler function with one parameter.
Let ε > 0 be an arbitrary number, and consider the differential inequality
| ν ( t ) D t 0.3 ν ( t ) | < ϵ , t [ 0 , 1 ] .
We have the inequality | t + 1 D t 0.3 ( t + 1 ) | = | t + 1 t 1 0.3 ( 1 0.3 ) Γ ( 1.3 ) | < 1 , t ( 0 , 1 ) (see Figure 2). Then, we could choose the solution of the differential inequality ν ( t ) = ϵ ( t + 1 ) , t [ 0 , 1 ] .
Case 2.1. The initial value does not depend on the solution of the inequality. Consider the solution of the given differential equation with the initial value x 0 = 2 , i.e., the solution is x ( t ) = 2 e t , and then the difference | x ( t ) ν ( t ) | = | 2 E 0.3 ( t 0.3 ) ϵ ( t + 1 ) | depends significantly on ϵ, and we cannot find a constant C F > 0 independent from ϵ such that | x ( t ) ν ( t ) C F ϵ (see Figure 3).
Case 2.2. The initial value depends on the solution of the inequality. Consider the solution of the given differential equation with x 0 = ν ( 0 ) = ϵ , i.e., x ( t ) = ϵ E 0.3 ( t 0.3 ) . Then, the inequality | x ( t ) ν ( t ) | = ϵ | t + 1 E 0.3 ( t 0.3 ) | C F ϵ , t ( 0 , 1 ] holds with the constant C F = | 2 E 0.3 ( 1 ) | > 0 , i.e., according to Definition 1, the solution of the considered initial value problem is Ulam–Hyers stable (see Figure 4).
Ulam type stability is studied for IVPs for various kinds of differential equations such as ordinary differential equations [11,12,13], ordinary differential equation with impulses [14,15,16,17], and fractional differential equations (FDE) (for example, ref. [18] for the Caputo fractional derivative, ref. [19] for the variable-order Caputo hybrid FDE, ref. [20] for the Caputo type fuzzy FDE with time delays, ref. [21] for the Caputo fractional difference equation, ref. [22] for the Caputo–Fabrizio FDE, ref. [23] for the linear FDE, ref. [24] for the impulsive FDE, ref. [25] for the FDE with the Hilfer–Katugampola operator and impulses, ref. [26] for the Riemann–Liouville FDE with a delay, ref. [27] for the functional FDE, ref. [28] for the tempered Caputo fractional derivative of variable order) (these are only some of the many published papers about initial value problems for various types of differential equations).

2.2. Boundary Value Problems

We discuss the application of Ulam type stability to boundary value problems for various types of differential equations with different types of derivatives. Note this problem has been studied by many authors in many published papers.

2.2.1. Some Misunderstandings in Study Ulam Stability for Boundary Value Problems

Consider the differential Equation (1) with the boundary condition
G ( x ( 0 ) , x ( b ) ) = 0 ,
where the function G could be a linear function or a nonlinear function, and it could have more than two arguments.
Remark 3.
If we apply Step 2 and Step 3 to the BVP (1), (5), then we prove the existence (uniqueness) of the solution which depends only on the time variable (compare with the case of the IVP and Remark 1).
Let ε > 0 . Consider the differential inequality (3).
We recall the definition for Ulam type stability used by many authors in the literature, which is based on the classical papers [11,12] for initial value problems.
Definition 2.
The BVP (1), (5) is Ulam–Hyers stable if there exists a real number C F > 0 such that for each ε > 0 and for each solution ν : [ 0 , b ] R n of the DI (3), there exists a solution x : [ 0 , b ] R n of (1), (5) with | ν ( t ) x ( t ) | C F ε for t [ 0 , b ] .
We focus on two main points in the above definition:
P1.
According to Remark 3, any solution x ( · ) of BVP (1), (5) depends only on the time variable, and if Step 2 is applied, then this solution exists (and is unique). Thus, Definition 2 means that for any solution ν ( · ) of DI (3), the inequality | ν ( t ) x ( t ) | C F ε , t [ 0 , b ] holds, where x ( · ) is the unique solution of BVP (1), (5).
P2.
Any solution ν ( · ) of DI (3) might not satisfy the boundary condition (5), and the integral presentation obtained in Step 2 for the solution of BVP (1), (5) might not be true for ν ( · ) .
We illustrate both above-mentioned points (P1) and (P2) with examples.
Example 3
(Ordinary derivative). Let the derivative D t in differential Equation (1) be an ordinary derivative, n = 1 , b = 1 , F ( t , x ) = 1.5 x , and G ( u , v ) = u + e 1.5 v 2 , i.e., consider the scalar linear BVP
x ( t ) = 1.5 x ( t ) , t ( 0 , 1 ] , x ( 0 ) + e 1.5 x ( 1 ) = 2 .
The BVP (6) has a unique solution x ( t ) = e 1.5 t , t ( 0 , 1 ] (see P1).
Let ε > 0 be an arbitrary number. Consider the differential inequality
| ν ( t ) 1.5 ν ( t ) | < ε , t [ 0 , 1 ] .
Choose a solution of DI (7), for example, ν ( t ) = ε t 2 , t [ 0 , 1 ] . Note the solution ν ( t ) does not satisfy the boundary condition ν ( 0 ) + e 1.5 ν ( 1 ) = 2 (see P2).
Note the unique solution of BVP (6) does not depend on any solution of DI (7).
The difference | x ( t ) ν ( t ) | = | e 1.5 t ε t 2 | depends significantly on ϵ. This difference is unbounded with respect to ε. Therefore, there does not exist a constant C > 0 such that | ν ( t ) x ( t ) | = | ε t 2 e 1.5 t | C ε for any ε and t [ 0 , 1 ] , i.e., according to Definition 2, the BVP (6) is not Ulam–Hyers stable.
Choose another solution of DI (7), for example, ν ( t ) = ε e 1.6 t , t [ 0 , 1 ] . This solution does not satisfy the boundary condition ν ( 0 ) + e 1.5 ν ( 1 ) = 2 (see P2).
To apply Ulam type stability according to Definition 2, we have to compare the new solution ν ( t ) = ε e 1.6 t , t [ 0 , 1 ] of DI (7) with the same unique solution x ( t ) = e 1.5 t , t ( 0 , 1 ] of BVP (6).
The difference | x ( t ) ν ( t ) | = | 1 ε e 0.1 t | e 1.5 t depends significantly on ϵ. This difference is unbounded with respect to ε. Therefore, there does not exist a constant C > 0 such that | ν ( t ) x ( t ) | = | ε e 0.1 t 1 | e 1.5 t C ε for any ε and t [ 0 , 1 ] , i.e., according to Definition 2, the BVP (6) is not Ulam–Hyers stable.
Example 4
(Fractional derivative). Let the derivative D t in differential Equation (1) be Caputo fractional derivative (4), n = 1 , b = 1 , F ( t , x ) = 1.5 x , and G ( u , v ) = u + v E α ( 1.5 ) 2 , i.e., consider the scalar linear BVP
D t α x ( t ) = 1.5 x ( t ) , t ( 0 , 1 ] , x ( 0 ) + x ( 1 ) E α ( 1.5 ) = 2 .
The BVP (8) has a unique solution x ( t ) = E α ( 1.5 t α ) , t ( 0 , 1 ] (see P1).
Let ε > 0 be an arbitrary number. Consider the differential inequality
| D t α ν ( t ) 1.5 ν ( t ) | < ε , t [ 0 , 1 ] .
Choose a solution of DI (9), for example, ν ( t ) = ε t 2 , t [ 0 , 1 ] (see Figure 5).
Note the solution ν ( t ) does not satisfy the boundary condition ν ( 0 ) + ν ( 1 ) E α ( 1.5 ) = 2 (see P2).
The difference | x ( t ) ν ( t ) | = | E α ( 1.5 t α ) ε t 2 | depends significantly on ϵ (see Figure 6). This difference is unbounded with respect to ε. Therefore, there does not exist a constant C > 0 such that | ν ( t ) x ( t ) | = | ε t 2 e 1.5 t | C ε for any ε and t [ 0 , 1 ] , i.e., according to Definition 2, the BVP (8) is not Ulam–Hyers stable.
There are some recently published papers about Ulam type stability with the above-mentioned misunderstandings (P1) and (P2) with proofs assuming that any solution of the corresponding differential inequality satisfies the given boundary condition and applying the integral presentation of the solution of the given BVP to any solution of the differential inequality [29,30,31,32,33,34,35,36,37].

2.2.2. How to Avoid the Misunderstandings

Now, let us recall the situation about the initial value problem. The initial value of the solution of the differential equation has to depend on the initial value of the chosen solution of the differential inequality (see Case 1.2, Example 1 and Case 2.2, Example 2).
To apply the same idea as in the case of the IVP (see Section 2.1), we have to connect the boundary value condition of the given problem (6) with the values of the arbitrary chosen solution of the differential inequality (7). As it is mentioned in Remark 1, we need to have a parameter in the boundary value condition, and later, we have to choose this parameter in an appropriate way to depend on the chosen solution of the differential inequality.
We suggest a possible way to avoid the above misunderstandings and to connect the solution of the given BVP and the arbitrary solution of the corresponding differential inequality.
In the general case, we change the boundary condition (5) to a boundary condition with a parameter
g ( x ( 0 ) , x ( b ) , θ ) = 0 ,
where θ R is a parameter.
Then, we prove that the BVP (1), (10) has a solution for any parameter θ R .
Remark 4.
Any solution of BVP (1), (10) depends on both the time variable t and the parameter θ.
For an arbitrary ε > 0 , we choose an arbitrary solution ν ( t ) of the differential inequality (3), and we change Definition 2 for the US of BVP (1), (5) to the following definition of the US of BVP (1), (10):
Definition 3.
The BVP (1), (10) is Ulam–Hyers stable if there exists a real number C F > 0 such that for each ε > 0 and for each solution ν : [ 0 , b ] R n of the DI (3), there exists a solution x : [ 0 , b ] R n of (1), (10) with | ν ( t ) x ( t ) | C F ε for t [ 0 , b ] .
In a practical application of US, for any ε , we choose a solution ν ( · ) of the DI (3), choose the parameter θ to depend on ν ( 0 ) and ν ( b ) , i.e., Θ = ξ ( ν ( 0 ) , ν ( b ) ) , and consider the differential Equation (1) with the function g ( x ( 0 ) , x ( b ) , ξ ( ν ( 0 ) , ν ( b ) ) ) in the boundary condition (10). Consider the solution x ( t ) , t [ 0 , b ] of BVP (1), (10) with the above-defined boundary value condition and study the difference between this solution x ( t ) and the chosen solution ν ( t ) of the differential inequality (3). Now, the solution x ( t ) depends on the solution ν ( · ) of DI(3). Any change in the solution of DI (3) leads to a change in the solution of BVP (1), (10).
Example 5.
Consider again the simple differential equation
x ( t ) = 1.5 x ( t ) , t ( 0 , 1 ] ,
but change the boundary value condition x ( 0 ) + e 1.5 x ( 1 ) = 2 to a boundary value condition with a parameter, i.e., consider the boundary condition
x ( 0 ) + e 1.5 x ( 1 ) = Θ ,
where Θ R is a parameter.
The new BVP (11), (12) has a unique solution x ( t ) = 0.5 Θ e 1.5 t , t [ 0 , 1 ] for any parameter Θ R .
Then, for an arbitrary solution ν ( t ) of DI (7), we choose the parameter Θ = ν ( 0 ) + e 1.5 ν ( 1 ) .
For example, if ν ( t ) = ε t 2 , then Θ = e 1.5 ε .
Then, the solution of BVP (11), (12) is x ( t ) = 0.5 ε e 1.5 ( t 1 ) , t [ 0 , 1 ] and | ν ( t ) x ( t ) | = ε | t 2 0.5 e 1.5 ( t 1 ) | 0.5 ε for any ε > 0 and t [ 0 , 1 ] . Therefore, according to Definition 3, the BVP (11), (12) is Ulam–Hyers stable with a constant C = 0.5 (compare with Example 3).
Remark 5.
Some authors change Definition 2 concerning the choice of the solution of the differential inequality (3) to a solution of the differential inequality (3) satisfying the boundary condition (5). However, this misses the meaning of Ulam stability (compare with the IVP). In addition, it is more difficult to obtain a solution of the differential inequality (3) and the boundary condition (5) than to obtain any solution of the inequality (3).
Remark 6.
Note the idea of considering a parameter in the boundary condition was suggested and applied to a BVP for a nonlinear generalized proportional Caputo fractional differential equations in [38] and to a BVP for a system of multi-term delay fractional differential equations of Caputo type in [39].

3. Ulam Stability of BVP for Impulsive ψ -Caputo Fractional Differential Equation

To illustrate the ideas in Section 2.2, we use the Caputo fractional derivative with respect to another function (called also the ψ -Caputo fractional derivative), impulses, and boundary value conditions to the differential equation. To focus on the application of the new idea and the involvement of a parameter in the boundary value condition, we consider a linear boundary value condition.
Ulam type stability is studied for BVPs for differential equations with a ψ -Caputo fractional derivative in [40]. Unfortunately, the authors use similar misunderstandings as mentioned in Section 2.2.1. They choose a solution of the corresponding fractional differential inequality and assume this solution satisfies the given boundary condition (see (18) [40]).
In this paper, we suggest a way to apply the US to BVPs for differential equations with a ψ -Caputo fractional derivative, and more generally, we consider impulses in the equation and consider the fractional derivative of piecewise-constant order.

3.1. Preliminary Results

3.1.1. Some Results from Fractional Calculus

Let c < d < be fixed numbers and ψ : [ c , d ] R be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ c , d ] .
We now give the definitions for fractional differintegral with respect to a given function.
Definition 4
([41]). Let δ > 0 . The Riemann–Liouville fractional integral with respect to the function ψ (RLI) is defined by
I ψ ( t ) δ c υ ( t ) = 1 Γ ( δ ) c t ψ ( t ) ψ ( s ) δ 1 ψ ( s ) υ ( s ) d s , t ( c , d ] .
Definition 5
([41]). Let δ ( 0 , 1 ) . The Caputo fractional derivative with respect to the function ψ ( t ) (CFD) is defined by
D ψ ( t ) δ c C υ ( t ) = I ψ ( t ) 1 δ c 1 ψ ( t ) d d t υ ( t ) = 1 Γ ( 1 δ ) c t ( ψ ( t ) ψ ( s ) ) δ ν ( s ) d s , t ( c , d ] .
Remark 7.
Note the RLI given in Definition 4 and the CFD given in Definition 5 are also called the ψ-fractional integral and the ψ-Caputo fractional derivative, respectively.
Lemma 1
([41]). Let δ ( 0 , 1 ) . Then, we have
( I δ c ( D δ c C ν ) ) ( t ) = ν ( t ) ν ( c ) , t ( c , d ] ,
and
( D δ c C ( I δ c ν ) ) ( τ ) = ν ( τ ) , τ ( c , d ] .
Lemma 2
([41]). The solution of the IVP for the ψ-Caputo fractional differential equation ( D δ c C ν ) ( t ) = F ( t ) for t ( c , d ] and ν ( c ) = V 0 , where δ ( 0 , 1 ) , F C ( [ c , d ] , R ) , is ν ( t ) = V 0 + ( I δ c F ) ( t ) , t [ c , d ] .
We now present some results for impulsive ψ -Caputo fractional differential equations on the interval [ 0 , b ] .

3.1.2. Preliminary Results for Scalar Linear Impulsive ψ -Caputo Fractional Differential Equations

Let ψ : [ 0 , b ] R be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , b ] .
Note there are different approaches in the statement of the fractional differential equation with impulses. We use fractional derivatives with changeable lower limit at the impulsive points.
Let the sequences { τ k } k = 0 p + 1 , 0 = τ 0 < τ k < τ k + 1 < τ p + 1 = b , k = 1 , 2 , , p , and { δ k } k = 0 p , δ k ( 0 , 1 ) , be given.
In the paper, we use the notations ν ( τ k 0 ) = lim t τ k , t < τ k ν ( t ) = ν ( τ k ) and ν ( τ k + 0 ) = lim t τ k , t > τ k ν ( t ) , k = 1 , 2 , , p , for a function ν : [ 0 , b ] R .
Define
P C [ 0 , b ] = { ν : j = 0 p ( τ k , τ k + 1 ] R : ν C ( ( τ k , τ k + 1 ] , R ) , k = 0 , 1 , 2 , , p , ν ( τ k ) < , ν ( τ k + 0 ) < , k = 1 , 2 , , p }
with a norm | | ν | | = sup t j = 0 p ( τ k , τ k + 1 ] | ν ( t ) | .
Consider the following scalar linear impulsive ψ - Caputo fractional differential equation:
( D δ k τ k C ν ) ( t ) = F ( t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , 2 , , p , ν ( τ k + 0 ) = μ k ν ( τ k ) + γ k , k = 1 , 2 , , p ,
with the initial condition
ν ( 0 ) = ν 0 .
We give an explicit formula for the solution of IVP (17), (18).
Lemma 3.
IVP (17), (18) where F C ( [ 0 , b ] , R ) , μ k , γ k R , k = 1 , 2 , , p , has a solution
ν ( t ) = ν 0 j = 1 k μ j + Q ( F , k , t ) , k = 0 , 1 , , p ,
where
Q ( F , k , t ) = j = 0 k 1 m = j + 1 k μ m ( I δ j τ j F ) ( τ j + 1 ) + j = 1 k m = j + 1 k μ m γ j + ( I δ k τ k F ) ( t ) ,
( I δ k τ k F ) ( t ) is defined by (13) with c = τ k , δ = δ k , d = τ k + 1 , and
( I δ j τ j F ) ( τ j + 1 ) = 1 Γ ( δ j ) τ j τ j + 1 ψ ( τ j + 1 ) ψ ( s ) δ j 1 ψ ( s ) F ( s ) d s , j = 0 , 1 , 2 , , k 1 .
Proof. 
We use induction.
For t [ 0 , τ 1 ] , by direct application of Lemma 2 with c = 0 , d = τ 1 , δ = δ 0 , V 0 = ν 0 , we obtain (19).
Let t ( τ 1 , τ 2 ] . Then, from Lemma 2 with c = τ 1 , d = τ 2 , δ = δ 1 , and V 0 = ν ( τ 1 + 0 ) = μ 1 ν ( t 1 ) + γ 1 , we obtain
ν ( t ) = μ 1 ν 0 + μ 1 ( I δ 0 τ 0 F ) ( τ 1 ) + γ 1 + ( I δ 1 τ 1 F ) ( t ) , t ( τ 1 , τ 2 ] .
Let t ( τ 2 , τ 3 ] . Then, from Lemma 2 with c = τ 2 , d = τ 3 , δ = δ 2 , and V 0 = ν ( τ 2 + 0 ) = μ 2 ν ( t 2 ) + γ 2 , we obtain
ν ( t ) = μ 2 μ 1 ν 0 + μ 2 μ 1 ( I δ 0 τ 0 F ) ( τ 1 ) + μ 2 ( I δ 1 τ 1 F ) ( τ 2 ) + μ 2 γ 1 + γ 2 + ( I δ 2 τ 2 F ) ( t ) .
Continuing this process inductively, we obtain (19). □
Consider the impulsive ψ -Caputo fractional differential Equation (17) with the linear boundary condition
A ν ( 0 ) + B ν ( b ) = θ ,
where A , B R are given constants, θ R is a parameter.
We obtain the explicit form of the solution of BVP (17), (23).
Lemma 4.
Let F C ( [ 0 , b ] , R ) ; μ k , γ k R , k = 1 , 2 , , p , A , B R : A + B j = 1 p μ j 0 , are given constants, and θ R is a parameter. Then, the solution of BVP (17), (23) is
ν ( t ) = θ B Q ( F , p , b ) A + B j = 1 p μ j j = 1 k μ j + Q ( F , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p ,
where Q ( F , k , t ) is defined by (20).
Proof. 
According to Lemma 3, the solution of (17) satisfies (19) for t = b , i.e.,
ν ( b ) = ν ( 0 ) j = 1 p μ j + Q ( F , p , b ) .
From the boundary condition (23), we obtain A ν ( 0 ) + B ν ( b ) = ν ( 0 ) A + B j = 1 p μ j + B Q ( F , p , b ) = θ , i.e., ν ( 0 ) = θ B Q ( F , p , b ) A + B j = 1 p μ j .
From (19), we obtain (24). □

3.2. Application of the Algorithm for Ulam Type Stability to BVPs

We apply the algorithm presented in Section 2 for the practical application of Ulam type stability for BVPs.
  • Step 1. Statement of the problem
Consider the boundary value problem (BVP) for the nonlinear impulsive differential equation with the ψ -Caputo fractional derivative in the form
( D ψ δ k τ k C ν ) ( t ) = f ( t , ν ( t ) ) for t ( 0 , b ] , t t k , k = 1 , 2 , , p , ν ( τ k + 0 ) = μ k ν ( τ k ) + γ k for k = 1 , 2 , , p ,
with the boundary value condition
A x ( 0 ) + B x ( b ) ) = θ ,
where θ R is a parameter. Assume the following assumptions hold:
(A1).
The function ψ : [ 0 , b ] R is smooth and increasing with ψ ( t ) > 0 almost everywhere in [ 0 , b ] .
(A2).
The sequence { τ k } k = 0 p + 1 , 0 = τ 0 < τ k < τ k + 1 < τ p + 1 = b , k = 1 , 2 , , p .
(A3).
The sequence { δ k } k = 0 p , δ k ( 0 , 1 ) .
(A4).
The constants μ k , γ k R , A , B R are such that A + B j = 1 p μ j 0 .
(A5).
The function f C ( [ 0 , b ] × R , R ) .
(A6).
There exists a constant K > 0 such that | f ( t , x ) f ( t , y ) | K | x y | for t [ 0 , b ] and x , y R .
  • Step 2. Integral presentation.
Based on Lemma 4, we introduce the following definition:
Definition 6.
The mild solution of BVP (26), (27) is the solution of the integral equation
ν ( t ) = θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k μ j + L ( f , ν , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p ,
where ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) = 1 Γ ( δ k ) τ k t ψ ( t ) ψ ( s ) δ k 1 ψ ( s ) f ( s , ν ( s ) ) ) d s ,
( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) = 1 Γ ( δ j ) τ j τ j + 1 ψ ( τ j + 1 ) ψ ( s ) δ j 1 ψ ( s ) f ( s , ν ( s ) ) ) d s , f o r j = 0 , 1 , 2 , , k 1 ,
and
L ( f , ν , k , t ) = j = 0 k 1 m = j + 1 k μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + j = 1 k m = j + 1 k μ m γ j + ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) .
Theorem 1.
Let conditions (A1)–(A5) be satisfied and θ R be a parameter. The mild solution of BVP (26), (27) is a solution of BVP (26), (27) and vice versa for any value of θ R .
Proof. 
Let the function ν P C ( [ 0 , b ] , R ) be a mild solution of BVP (26), (27) (see Definition 6).
Note the expressions
j = 0 k 1 m = j + 1 k μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) ,
j = 1 k m = j + 1 k μ m γ j , k = 1 , 2 , , p ,
and
θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k μ j
do not depend on the time variable t, and their Caputo fractional derivatives with respect to the function ψ are zero.
Then, for any k = 1 , 2 , , p , applying Lemma 1 and (16) with c = τ k , d = τ k + 1 , δ = δ k , we obtain
( D δ k τ k C L ( f , ν , k , · ) ) ( t ) = D δ k τ k C ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) = f ( t , ν ( t ) ) .
Take a fractional derivative ( D δ k τ k C ( · ) ) ( t ) , t ( τ k , τ k + 1 ] of both sides of (28), use (30), and obtain ( D δ k τ k C ν ) ( t ) = f ( t , ν ( t ) ) for τ ( t k , t k + 1 ] , k = 0 , 1 , , p , i.e., the function ν ( · ) satisfies the first equation of (26).
Let t = τ k , k = 1 , 2 , , p . From (28), for t = τ k + 0 ( τ k , τ k + 1 ) and t k 0 ( t k 1 , t k ] , applying ( I δ k τ k f ( · , ν ( · ) ) ) ( τ k ) = 0 , we have
ν ( t k + 0 ) = θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k μ j + L ( f , ν , k , τ k ) = μ k θ B L ( ν , p , b ) A + B j = 1 p μ j j = 1 k 1 μ j + j = 0 k 1 μ k m = j + 1 k 1 μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + j = 1 k m = j + 1 k μ m γ j + ( I δ k τ k f ( · , ν ( · ) ) ) ( τ k ) = μ k θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k 1 μ j + j = 0 k 1 μ k m = j + 1 k 1 μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + j = 1 k 1 μ k m = j + 1 k 1 μ m γ j + m = k + 1 k μ m γ k = μ k θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k 1 μ j + L ( f , ν , k 1 , τ k ) + γ k = μ k ν ( τ k ) + γ k .
The equalities in (31) imply that the second equation of (26) is satisfied for all k = 1 , 2 , , p .
From (28), using L ( f , ν , 0 , 0 ) = 0 , we obtain
A ν ( 0 ) + B ν ( b ) = A θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 0 μ j + B θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 p μ j + B L ( f , ν , p , b ) = θ B L ( f , ν , p , b ) + B L ( f , ν , p , b ) = Θ ,
i.e., the function ν ( · ) satisfies the boundary value condition (27).
Therefore, the mild solution ν ( · ) is a solution of BVP (26), (27).
Let the function ν P C ( [ 0 , b ] be a solution of BVP (26), (27). Taking an integral ( I δ τ k ( · ) ) ( t ) on both sides of the first equation of (26) for t ( τ k , τ k + 1 ] and using Lemma 1, similar to the proof of Lemma 3, we obtain the integral equality (28). □
  • Step 3. Existence of a solution.
Define the fractional integral operator Ω : P C [ 0 , b ] P C [ 0 , b ] by
( Ω ν ) ( τ ) = θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k μ j + L ( f , ν , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p .
Theorem 2.
Let conditions (A1)–(A6) be satisfied, θ R , and the inequality
K M P < 1
hold, where
M = max { 1 Γ ( δ j + 1 ) ψ ( τ j + 1 ) ψ ( τ j ) δ j : j = 0 , 1 , 2 , , p } ,
and
P = max k = 0 , 1 , 2 , , p j = 0 k m = j + 1 k | μ m | .
Then, the operator Ω has a fixed point.
Proof. 
Let ν , ν P C [ 0 , b ] . Note 1 Γ ( δ j ) τ j τ ψ ( τ ) ψ ( s ) δ j 1 ψ ( s ) d s = ψ ( τ j + 1 ) ψ ( τ j ) δ j Γ ( δ j + 1 ) for τ ( τ j , τ j + 1 ] , j = 0 , 1 , 2 , , p , and we obtain for t ( τ k , τ k + 1 ] , k = 0 , 1 , , p that
| ( Ω ν ) ( t ) | ( Ω ν ) ( t ) | | L ( f , ν , k , t ) L ( f , ν , k , t ) | j = 0 k 1 m = j + 1 k | μ m | ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) ν ν j = 0 k 1 m = j + 1 k | μ m | K 1 Γ ( δ j + 1 ) ψ ( τ j + 1 ) ψ ( τ j ) δ j + 1 Γ ( δ k + 1 ) ψ ( τ k + 1 ) ψ ( τ j ) δ k K ν ν = K | | ν ν | | j = 0 k m = j + 1 k | μ m | 1 Γ ( δ j + 1 ) ψ ( τ j + 1 ) ψ ( τ j ) δ j K ν ν j = 0 p m = j + 1 p | μ m | 1 Γ ( δ j + 1 ) ψ ( τ j + 1 ) ψ ( τ j ) δ j .
Therefore,
| | ( Ω ν ) ( Ω ν ) | | K | | ν ν | | j = 0 p m = j + 1 p | μ m | 1 Γ ( δ j + 1 ) ψ ( τ j + 1 ) ψ ( τ j ) δ j .
Then, the claim of Theorem 2 follows from the Banach construction principle. □
Corollary 1.
Let all the conditions of Theorem 2 be satisfied. Then, the BVP (26), (27) has a solution for any θ R , and it satisfies the integral Equation (28).
The proof follows from Theorem 2 and Theorem 1.
  • Step 4. Definition of Ulam type stability
Let ε > 0 and consider the following impulsive ψ -Caputo fractional differential inequalities (FDI):
( D δ k τ k C y ) ( t ) f ( t , y ( t ) ) ε , t ( τ k , τ k + 1 ] , k = 0 , 1 , 2 , , p , | y ( τ k + 0 ) μ k y ( τ k ) γ k | ε , k = 1 , 2 , , p .
Note if ν P C [ 0 , b ] is a solution of (38), then there exist a function g C [ 0 , b ] : | g ( t ) | ε , t [ 0 , b ] and constants η k : | η k | ε , k = 1 , 2 , , p , such that ν ( t ) is a solution of the following impulsive fractional differential equation
( D δ k τ k C y ) ( t ) = f ( t , y ( t ) ) + g ( t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , 2 , , p , y ( τ k + 0 ) = μ k y ( τ k ) + γ k + η k , k = 1 , 2 , , p .
Remark 8.
Note that the solution of the impulsive fractional differential inequalities (38) might not satisfy the boundary condition of (27). A similar comment concerns any solution of the impulsive fractional differential Equation (39).
Definition 7.
The BVP (26), (27) is Ulam–Hyers stable if there exists a real number C f > 0 such that for each ε > 0 and for each solution y P C [ 0 , b ] of the FDI (38), there exists a solution ν P C [ 0 , b ] of (26), (27) with θ = A y ( 0 ) + B y ( b ) such that | | ν y | | C f ε .
  • Step 5. Proof of US.
Theorem 3.
Let the assumptions (A1)–(A6) hold and
K M P T < 1 ,
where P is defined by (36), M is defined by (35), Q = max k = 0 , 1 , 2 , , p m = 1 k | μ m | , and T = 1 + | B | Q A + B j = 1 p μ j .
Then, the BVP (26), (27) is US.
Proof. 
Let ε > 0 and y P C [ 0 , b ] be a solution of the fractional differential inequality (38). Note it satisfies the impulsive differential Equation (39). Let θ = A y ( 0 ) + B y ( b ) . Then, the conditions (A5) and (A6) are satisfied for the function G ( t , v ) = f ( t , v ) + g ( t ) and according to Corollary 1, the impulsive differential Equation (39) with a boundary condition A y ( 0 ) + B y ( b ) = θ has a solution satisfying the integral equality
y ( t ) = θ B L ( f + g , y , p , b ) A + B j = 1 p μ j j = 1 k μ j + L ( f + g , y , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p .
If inequality (40) holds, then the inequality (34) is fulfilled, and all conditions of Theorem 2 are satisfied. According to Corollary 1, there exists a solution ν P C ( [ 0 , b ] ) of (26), (27) for the chosen parameter θ , and this solution satisfies the integral equality (28).
Then,
| y ( t ) ν ( t ) | L ( f + g , y , p , b ) L ( f , ν , p , b ) B j = 1 k μ j A + B j = 1 p μ j + L ( f + g , y , k , t ) L ( f , ν , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p .
For any k = 0 , 1 , 2 , , p and t ( τ k , τ k + 1 ] , we apply the inequalities
( I δ k τ k | g | ) ( t ) ε Γ ( 1 + δ k ) ( ψ ( t ) ψ ( τ k ) ) δ k ,
( I δ j τ j | g | ) ( t j + 1 ) ε Γ ( 1 + δ j ) ( ψ ( τ j + 1 ) ψ ( τ j ) ) δ j ,
| f ( s , y ( s ) ) f ( s , ν ( s ) ) | K | y ( s ) ν ( s ) | K | | y ν | | , s ( τ k , t ] ,
and
( I δ j τ j | f ( · , y ( · ) ) f ( · , ν ( · ) ) | ) ( τ j + 1 ) K | | y ν | | Γ ( 1 + δ j ) ( ψ ( τ j + 1 ) ψ ( τ j ) ) δ j
and we obtain
L ( f + g , y , k , t ) L ( f , ν , k , t ) j = 0 k 1 m = j + 1 k μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) ( I δ j τ j f ( · , y ( · ) ) ) ( τ j + 1 ) + j = 0 k 1 m = j + 1 k μ m ( ( I δ j τ j g ) ( τ j + 1 ) + ( I δ k τ k g ) ( t ) + ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) ( I δ k τ k f ( · , y ( · ) ) ) ( t ) j = 0 k 1 m = j + 1 k | μ m | K | | y ν | | Γ ( 1 + δ j ) ( ψ ( τ j + 1 ) ψ ( τ j ) ) δ j + j = 0 k 1 m = j + 1 k | μ m | ε Γ ( 1 + δ j ) ( ψ ( τ j + 1 ) ψ ( τ j ) ) δ j + ε Γ ( 1 + δ k ) ( ψ ( t ) ψ ( τ k ) ) δ k + K | | y ν | | Γ ( 1 + δ k ) ( ψ ( τ k + 1 ) ψ ( τ k ) ) δ k K M | | y ν | | j = 0 k m = j + 1 k | μ m | + ε M j = 0 k m = j + 1 k | μ m | K M P | | y ν | | + ε M P .
From (42) and (43), we obtain for t ( τ k , τ k + 1 ] , k = 0 , 1 , , p ,
| y ( t ) ν ( t ) | K M P | | y ν | | + ε M P 1 + B j = 1 k μ j A + B j = 1 p μ j K M P T | | y ν | | + ε M P T .
Therefore,
| | y ν | | ε M P T 1 K M P T ,
or the BVP (26), (27) is Ulam–Hyers stable with a constant C f = M P T 1 K M P T . □

3.3. Partial Cases

From the obtained result concerning US, we can obtain as partial cases several results for various problems.

3.3.1. Boundary Value Problem for a Fractional Differential Equation without Any Impulses

Consider the fractional differential equation
( D ψ δ 0 C ν ) ( t ) = f ( t , ν ( t ) ) f o r t ( 0 , b ] ,
with the boundary condition (27), i.e., μ k = 1 , γ k = 0 , k = 1 , 2 , , p in (26).
The integral Equation (28) reduces to the following integral equation for the solution of BVP (46), (27):
ν ( t ) = θ B ( I δ 0 f ( · , ν ( · ) ) ( b ) A + B + ( I δ 0 f ( · , ν ( · ) ) ) ( t ) , t ( 0 , b ] ,
where ( I δ 0 f ( · , ν ( · ) ) ) ( t ) = 1 Γ ( δ ) 0 t ψ ( t ) ψ ( s ) δ 1 ψ ( s ) f ( s , ν ( s ) ) d s .
Then, Theorem 3 reduces to the following result:
Theorem 4.
Let the assumptions (A1), (A5), (A6) hold, δ ( 0 , 1 ) , and
K ( ψ ( b ) ψ ( 0 ) ) δ 1 + | B | | A + B | < Γ ( 1 + δ ) .
Then, the BVP (46), (27) is Ulam–Hyers stable.

3.3.2. Boundary Value Problem for an Impulsive Caputo Fractional Differential Equation

Let the function ψ ( t ) t , t [ 0 , b ] . Then, the ψ -Caputo fractional derivative given by (14) reduces to the Caputo fractional derivative D δ c C υ ( t ) .
Consider the impulsive Caputo fractional differential equation
( D δ k 0 C ν ) ( t ) = f ( t , ν ( t ) ) for t ( τ k , τ k + 1 ] , k = 0 , 1 , 2 , , p , ν ( τ k + 0 ) = μ k ν ( τ k ) + γ k , k = 1 , 2 , , p ,
with the boundary condition (27).
The integral Equation (28) reduces to the following integral equation for the solution of BVP (49), (27):
ν ( t ) = θ B L ( f , ν , p , b ) A + B j = 1 p μ j j = 1 k μ j + L ( f , ν , k , t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p ,
where
( I δ k τ k f ( · , ν ( · ) ) ) ( t ) = 1 Γ ( δ k ) τ k t t s δ k 1 f ( s , ν ( s ) ) ) d s ,
( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) = 1 Γ ( δ j ) τ j τ j + 1 t s δ j 1 f ( s , ν ( s ) ) ) d s , j = 0 , 1 , 2 , , k 1 ,
and
L ( f , ν , k , t ) = j = 0 k 1 m = j + 1 k μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + j = 1 k m = j + 1 k μ m γ j + ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) .
Then, Theorem 3 reduces to the following result:
Theorem 5.
Let the assumptions (A1), (A5), (A6) hold, δ ( 0 , 1 ) , and
K M P T ,
where M = max { τ j + 1 τ j Γ ( δ j + 1 ) : j = 0 , 1 , 2 , , p } ; P is defined by (36), T is defined in Theorem 3.
Then, the BVP (46), (27) is Ulam–Hyers stable.

3.3.3. Initial Value Problem for an Impulsive Fractional Differential Equation

Consider the impulsive fractional differential Equation (26) with the initial condition
ν ( 0 ) = θ ,
i.e., A = 1 , B = 0 in (27).
The integral Equation (28) reduces to the following integral equation for the solution of IVP (26), (52):
ν ( t ) = θ j = 1 k μ j + j = 1 k m = j + 1 k μ m γ j + j = 0 k 1 m = j + 1 k μ m ( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) + ( I δ k τ k f ( · , ν ( · ) ) ) ( t ) , t ( τ k , τ k + 1 ] , k = 0 , 1 , , p ,
where
( I δ k τ k f ( · , ν ( · ) ) ) ( t ) = 1 Γ ( δ k ) τ k t ψ ( t ) ψ ( s ) δ k 1 ψ ( s ) f ( s , ν ( s ) ) ) d s ,
and
( I δ j τ j f ( · , ν ( · ) ) ) ( τ j + 1 ) = 1 Γ ( δ j ) τ j τ j + 1 ψ ( τ j + 1 ) ψ ( s ) δ j 1 ψ ( s ) f ( s , ν ( s ) ) ) d s , for j = 0 , 1 , , k 1 .
Then, Theorem 3 reduces to the following result:
Theorem 6.
Let the assumptions (A1), (A2), (A3), (A5), (A6) hold, μ k , γ k R , k = 1 , 2 , , p and
K M P < 1 ,
where P is defined by (36), M is defined by (35).
Then, the IVP (26), (52) is Ulam–Hyers stable.
Remark 9.
The US for the IVP for an impulsive Caputo fractional differential equation (with ψ ( t ) t in (26)) is studied in [24] but in the case where the lower limit of the fractional derivative is fixed at the initial time.

3.3.4. Initial Value Problem for a Fractional Differential Equation without Any Impulses

Consider the fractional differential Equation (46) with the initial condition (52), i.e., A = 1 , B = 0 in (27) and μ k = 1 , γ = 0 , k = 1 , 2 , , p in (26).
The integral Equation (28) reduces to the following integral equation for the solution of IVP (46), (52) (compare with [41])
ν ( t ) = θ + 1 Γ ( δ ) 0 t ψ ( t ) ψ ( s ) δ 1 ψ ( s ) f ( s , ν ( s ) ) ) d s , t ( 0 , b ] .
Then, Theorem 3 reduces to the following result:
Theorem 7.
Let the assumptions (A1), (A5), (A6) hold, and K ( ψ ( b ) ψ ( 0 ) ) δ < Γ ( 1 + δ ) .
Then, the IVP (46), (52) is Ulam–Hyers stable.
Remark 10.
In the case where ψ ( t ) t in Equation (46), the partial case of the IVP for a Caputo fractional differential equation is studied in [18].

3.4. Application

Example 6.
Let b = 1 , p = 4 , τ 1 = 0.1 , τ 2 = 0.5 , τ 3 = 0.7 , τ 4 = 0.9 , μ 1 = 2 , μ 2 = 0.5 , μ 3 = 3 , μ 4 = 1 3 , δ = 0.3 , ψ ( t ) = t t + 1 , δ k = 0.3 , and k = 0 , 1 , 2 , 3 .
Consider the following scalar impulsive ψ-Caputo fractional differential equation:
( D ψ 0.3 τ k C ν ) ( t ) = t 8 t + 1 e | ν ( t ) | , t ( τ k , τ k + 1 ] , k = 0 , 1 , 2 , 3 , 4 , ν ( τ k + 0 ) = μ k ν ( τ k ) + γ k , k = 1 , 2 , 3 , 4 ,
with boundary condition
ν ( 0 ) + 0.2 ν ( 1 ) = θ ,
where θ R is a parameter.
The conditions (A1)–(A3) are satisfied. Also, k = 1 4 μ k = 1 and A + B = 1.2 0 , i.e., condition (A4) holds.
In this case, f ( t , y ) = t 8 t + 1 e | y | , i.e., condition (A6) is satisfied with K = 0.125 .
We have
M = max { 0.5427 , 0.7284 , 0.5192 , 0.4837 } = 0.7284 ,
P = max { 1 , 1 + | μ 1 | , 1 + | μ 2 | + | μ 1 | | μ 2 | , 1 + | μ 3 | + | μ 3 | | μ 2 | + | μ 3 | | μ 2 | | μ 1 | } = 6.5 ,
and K M P = 0.5966 < 1 , i.e., all the conditions of Theorem 2 are satisfied; therefore, BVP (56), (57) has a solution.
Also, Q = max { 1 , | μ 1 | , | μ 1 | | μ 2 | , | μ 1 | | μ 2 | | μ 3 | , | μ 1 | | μ 2 | | μ 3 | | μ 4 | } = 3 , | A + B | = 1.2 , T = 1.5 and K M P T = 0.894839 < 1 , i.e., all the conditions of Theorem 3 are satisfied, and the solution of BVP (56), (57) is Ulam–Hyers stable.

4. Conclusions

The main goal of this paper was to focus on the application of Ulam type stability to boundary value problems for any type of differential equations and to point out the misunderstandings appearing in the literature. We suggested a possible way to avoid these misunderstandings by involving a parameter. We illustrated the suggested idea on a general linear boundary value problem for a nonlinear impulsive differential equation with a Caputo fractional derivative with respect to another function and a piecewise-constant variable order. This involved a parameter in the boundary condition. The existence and uniqueness of the solution depending on this parameter was proved, the Ulam–Hyers stability was defined, and sufficient conditions were obtained. The proposed ideas could be applied to study Ulam type stability for various types of boundary value problems for other types of differential equations, such as ordinary differential equations, fractional differential equation with any type of Caputo type fractional derivative, integro-differential fractional equations, and partial differential equations with time-fractional derivatives.
The suggested new way of the application of Ulam type stability to boundary value problems and the idea of the connection between the solutions of the given problem and the solutions of the corresponding inequality will give researchers a tool for a correct application of Ulam type stability to various types of boundary value problems.

Author Contributions

R.P.A., S.H. and D.O.; methodology, R.P.A., S.H. and D.O.; formal analysis, R.P.A., S.H. and D.O.; investigation, R.P.A., S.H. and D.O.; writing—original draft preparation, R.P.A., S.H. and D.O.; writing—review and editing, R.P.A., S.H. and D.O. All authors have read and agreed to the published version of the manuscript.

Funding

S.H. is partially supported by the Bulgarian National Science Fund under Project KP-06-N62/1.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ulam, S.M. A Collection of Mathematical Problems; Interscience Tracts in Pure and Applied Mathematics, no. 8; Interscience Publishers: New York, NY, USA, 1960. [Google Scholar]
  2. Hyers, D.H. On the stability of the linear functional equation. Proc. Nat. Acad. Sci. USA 1941, 27, 222–224. [Google Scholar] [CrossRef]
  3. Tamilvanan, K.; Alkhaldi, A.H.; Jakhar, J.; Chugh, R.; Jakhar, J.; Rassias, J.M. Ulam Stability Results of Functional Equations in Modular Spaces and 2-Banach Spaces. Mathematics 2023, 11, 371. [Google Scholar] [CrossRef]
  4. Ciepliński, K. On Ulam Stability of a Functional Equation. Results Math. 2020, 75, 151. [Google Scholar] [CrossRef]
  5. Jung, S.-M. Hyers–Ulam–Rassias Stability of Functional Equations in Nonlinear Analysis; Springer: New York, NY, USA, 2011. [Google Scholar]
  6. Cadariu, L. Stabilitatea Ulam–Hyers–Bourgin Pentru Ecuatii Functionale; Universitatea de Vest din Timisoara: Timisoara, Romania, 2007. [Google Scholar]
  7. Hyers, D.H.; Isac, G.; Rassias, T. Stability of Functional Equations in Several Variables; Birkhauser: Basel, Switzerland, 1998. [Google Scholar]
  8. Agarwal, R.; Hristova, S.; O’Regan, D. Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam Type-Stability. Mathematics 2022, 10, 2315. [Google Scholar] [CrossRef]
  9. Batista, M. On stability of non-inflectional elastica. Comptes Rendus. Mec. 2020, 348, 137–148. [Google Scholar] [CrossRef]
  10. Joshi, H.; Jha, B.K. 2D memory-based mathematical analysis for the combined impact of calcium influx and efflux on nerve cells. Comput. Math. Appl. 2023, 134, 33–44. [Google Scholar] [CrossRef]
  11. Rus, I.A. Ulam stability of ordinary differential equations. Stud. Univ.-Babes-Bolyai Math. 2009, LIV, 125–133. [Google Scholar]
  12. Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpathian J. Math. 2010, 26, 103–107. [Google Scholar]
  13. Tripathy, A.K. Hyers-Ulam Stability of Ordinary Differential Equations; Chapman and Hall/CRC: New York, NY, USA, 2021. [Google Scholar]
  14. Liao, Y.; Wang, J. A note on stability of impulsive differential equations. Bound Value Probl. 2014, 2014, 67. [Google Scholar] [CrossRef]
  15. Zada, A.; Faisal, S.; Li, Y. On the Hyers-Ulam stability of first-order impulsive delay differential equations. J. Funct. Spaces 2006, 2006, 8164978. [Google Scholar] [CrossRef]
  16. Hristova, R.A.S.; O’Regan, D. Ulam type stability results for non-instantaneous impulsive differential equations with finite state dependent delay. Dyn. Syst. Appl. 2019, 28, 47–61. [Google Scholar]
  17. Wang, J.R.; Feckan, M.; Zhou, Y. Ulam’s type stability of impulsive ordinary differential equations. J. Math. Anal. Appl. 2012, 395, 258–264. [Google Scholar] [CrossRef]
  18. Wang, J.; Lv, L.; Zhou, Y. Ulam stability and data dependence for fractional differential equations with Caputo derivative. Electron. J. Qual. Theory Differ. Equ. 2011, 63, 1–10. [Google Scholar] [CrossRef]
  19. Limpanukorn, N.; Ngiamsunthorn, P.S. Existence and Ulam stability of solution to fractional order hybrid differential equations of variable order. Thai J. Math. 2020, 18, 453–463. [Google Scholar]
  20. Wang, X.; Luo, D.; Zhu, Q. Ulam-Hyers stability of Caputo type fuzzy fractional differential equations with time-delays. Chaos Solitons Fractals 2022, 156, 111822. [Google Scholar] [CrossRef]
  21. Chen, C.; Bohner, M.; Jia, B. Ulam-Hyers stability of Caputo fractional difference equations. Math. Meth. Appl. Sci. 2019, 42, 7461–7470. [Google Scholar] [CrossRef]
  22. Liu, K.; Feckan, M.; O’Regan, D.; Wang, J.R. Hyers–Ulam Stability and Existence of Solutions for Differential Equations with Caputo–Fabrizio Fractional Derivative. Mathematics 2019, 7, 333. [Google Scholar] [CrossRef]
  23. Liu, H.; Li, Y. Hyers–Ulam stability of linear fractional differential equations with variable coefficients. Adv. Differ. Equ. 2020, 2020, 404. [Google Scholar] [CrossRef]
  24. Wang, J.; Zhou, Y.; Feckan, M. Nonlinear impulsive problems for fractional differential equations and Ulam stability. Comput. Math. Appl. 2012, 64, 3389–3405. [Google Scholar] [CrossRef]
  25. Harikrishnan, S.; Ibrahim, R.W.; Kanagarajan, K. Fractional Ulam-stability of fractional impulsive differential equation involving Hilfer-Katugampola fractional differential operator. Univers. J. Math. Appl. 2018, 1, 106–112. [Google Scholar] [CrossRef]
  26. Agarwal, R.; Hristova, S.; O’Regan, D. Existence and Ulam type stability for nonlinear Riemann–-Liouville fractional differential equations with constant delay. Electron. J. Qual. Theory Differ. Equ. 2020, 67, 1–18. [Google Scholar] [CrossRef]
  27. Vanterler da, C.; Sousa, J.; Capelas de Oliveira, E.; Rodrigues, F.G. Ulam-Hyers stabilities of fractional functional differential equations. AIMS Math. 2020, 5, 1346–1358. [Google Scholar] [CrossRef]
  28. O’Regan, D.; Hristova, S.; Agarwal, R.P. Ulam-Type Stability Results for Variable Order Ψ-Tempered Caputo Fractional Differential Equations. Fractal Fract. 2024, 8, 11. [Google Scholar] [CrossRef]
  29. El-Sayed, A.M.A.; Al Issa, S.M.; Elmiari, M. Ulam-type Stability for a Boundary Value Problem of Implicit Fractional-orders Differential Equation. Adv. Dynam. Syst. Appl. (ADSA) 2021, 16, 75–89. [Google Scholar] [CrossRef]
  30. Rahman, G.; Agarwal, R.P.; Ahmad, D. Existence and stability analysis of n-th order multi term fractional delay differential equation. Chaos Solitons Fractals 2022, 155, 111709. [Google Scholar] [CrossRef]
  31. 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. [Google Scholar] [CrossRef]
  32. Shah, K.; Tunc, S. Existence theory and stability analysis to a system of boundary value problem. J. Taibah Uni. Sci. 2017, 11, 1330–1342. [Google Scholar] [CrossRef]
  33. Zada, A.; Ali, S.; Li, Y. Ulam-type stability for a class of implicit fractional differential equations with non-instantaneous integral impulses and boundary condition. Adv. Differ. Equ. 2017, 2017, 317. [Google Scholar] [CrossRef]
  34. Salim, A.; Benchohra, M.; Lazreg, J.E.; N’Guerekata, G. Boundary Value Problem for Nonlinear Implicit Generalized Hilfer-Type Fractional Differential Equations with Impulses. Abstr. Appl. Anal. 2021, 2021, 5592010. [Google Scholar] [CrossRef]
  35. Benkerrouche, A.; Souid, M.S.; Stamov, G.; Stamova, I. Multiterm Impulsive Caputo–Hadamard Type Differential Equations of Fractional Variable Order. Axioms 2022, 11, 634. [Google Scholar] [CrossRef]
  36. Zamir, M.; Shah, K.; Chohan, M.I. Stability theory and existence of solution to a multi-point boundary value problem of fractional diferential equations. Math. Sci. 2019, 13, 53–60. [Google Scholar] [CrossRef]
  37. Treanbucha, C.; Sudsutad, W. Stability analysis of boundary value problems for Caputo proportional fractional derivative of a function with respect to another function via impulsive Langevin equation. AIMS Math. 2021, 6, 6647–6686. [Google Scholar] [CrossRef]
  38. Agarwal, R.; Hristova, S.; O’Regan, D. Boundary Value Problems for Fractional Differential Equations of Caputo Type and Ulam Type Stability: Basic Concepts and Study. Axioms 2023, 12, 226. [Google Scholar] [CrossRef]
  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. Wahash, H.A.; Panchal, S.K.; Abdo, M.S. Existence and stability of a nonlinear fractional differential equation involving a ψ-Caputo operator. Adv. Theory Nonlinear Anal. Its Appl. 2020, 4, 266–278. [Google Scholar] [CrossRef]
  41. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 460–481. [Google Scholar] [CrossRef]
Figure 1. Graph of | x ( t ) ν ( t ) | for various values of ε .
Figure 1. Graph of | x ( t ) ν ( t ) | for various values of ε .
Mathematics 12 01626 g001
Figure 2. Graph of | t D t 0.3 ( t ) | and the bound 1.
Figure 2. Graph of | t D t 0.3 ( t ) | and the bound 1.
Mathematics 12 01626 g002
Figure 3. Graph of | x ( t ) ν ( t ) | for various values of ε .
Figure 3. Graph of | x ( t ) ν ( t ) | for various values of ε .
Mathematics 12 01626 g003
Figure 4. Graph of | t + 1 E 0.3 ( t 0.3 ) | and the bound | 2 E 0.3 ( 1 ) | .
Figure 4. Graph of | t + 1 E 0.3 ( t 0.3 ) | and the bound | 2 E 0.3 ( 1 ) | .
Mathematics 12 01626 g004
Figure 5. Graph of | D t α t 2 1.5 t 2 | = | 2 t 2 α Γ ( 3 α ) 1.5 t 2 | and the bound 1.
Figure 5. Graph of | D t α t 2 1.5 t 2 | = | 2 t 2 α Γ ( 3 α ) 1.5 t 2 | and the bound 1.
Mathematics 12 01626 g005
Figure 6. Graph of | E α ( 1.5 t α ) ε t 2 | for various ε ’s.
Figure 6. Graph of | E α ( 1.5 t α ) ε t 2 | for various ε ’s.
Mathematics 12 01626 g006
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

Agarwal, R.P.; Hristova, S.; O’Regan, D. Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them. Mathematics 2024, 12, 1626. https://doi.org/10.3390/math12111626

AMA Style

Agarwal RP, Hristova S, O’Regan D. Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them. Mathematics. 2024; 12(11):1626. https://doi.org/10.3390/math12111626

Chicago/Turabian Style

Agarwal, Ravi P., Snezhana Hristova, and Donal O’Regan. 2024. "Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them" Mathematics 12, no. 11: 1626. https://doi.org/10.3390/math12111626

APA Style

Agarwal, R. P., Hristova, S., & O’Regan, D. (2024). Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them. Mathematics, 12(11), 1626. https://doi.org/10.3390/math12111626

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