Next Article in Journal
On q-Hermite–Hadamard Type Inequalities via s-Convexity and (α,m)-Convexity
Next Article in Special Issue
Variable-Order Fractional Linear Systems with Distributed Delays—Existence, Uniqueness and Integral Representation of the Solutions
Previous Article in Journal
Existence Result for Coupled Random First-Order Impulsive Differential Equations with Infinite Delay
Previous Article in Special Issue
Lyapunov Functions and Stability Properties of Fractional Cohen–Grossberg Neural Networks Models with Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ulam-Type Stability Results for Variable Order Ψ-Tempered Caputo Fractional Differential Equations

1
School of Mathematical and Statistical Sciences, University of Galway, H91 TK33 Galway, Ireland
2
Faculty of Mathematics and Informatics, Plovdiv University “P. Hilendarski”, 4000 Plovdiv, Bulgaria
3
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(1), 11; https://doi.org/10.3390/fractalfract8010011
Submission received: 4 November 2023 / Revised: 9 December 2023 / Accepted: 14 December 2023 / Published: 22 December 2023
(This article belongs to the Special Issue Advances in Variable-Order Fractional Calculus and Its Applications)

Abstract

:
An initial value problem for nonlinear fractional differential equations with a tempered Caputo fractional derivative of variable order with respect to another function is studied. The absence of semigroup properties of the considered variable order fractional derivative leads to difficulties in the study of the existence of corresponding differential equations. In this paper, we introduce approximate piecewise constant approximation of the variable order of the considered fractional derivative and approximate solutions of the given initial value problem. Then, we investigate the existence and the Ulam-type stability of the approximate solution of the variable order Ψ -tempered Caputo fractional differential equation. As a partial case of our results, we obtain results for Ulam-type stability for differential equations with a piecewise constant order of the Ψ -tempered Caputo fractional derivative.

1. Introduction

Fractional derivatives and integrals of variable order are widely studied in the literature, and a comprehensive review and a good mathematical framework of fractional calculus of variable order, based on Scarpi’s approach, is provided in [1]. Surveys of the types of derivatives/integrals of variable order with some physical discussions are given in [2,3]. There are several types of definitions of variable order fractional derivatives/integrals (see, for example, [4,5,6]). These definitions keep the basic property of fractional derivatives—the nonlocality. This property unfortunately does not have the semigroup property when variable order is applied (see, for example, [5,7]). This does not allow us directly to convert the differential equation with variable order fractional derivative to an equivalent fractional integral equation as is possible for fractionals of constant order (see, for example, [8,9,10,11,12,13]). In [14], the author defines a derivative in a new way, again called a fractional derivative. This derivative keeps the locality property of ordinary derivatives and it enables one to apply directly a variable order with a semigroup property.
One of the main problems, existence, is studied by many authors. For the Riemann–Liouville fractional derivative of variable order for differential equations, existence is studied in [7], and the concept of approximate solution of the given problem is introduced. This concept is based on the application of piecewise constant orders of fractional integrals and derivatives used for various problems in [15,16].
To be more general, we consider a tempered Caputo fractional derivative with respect to another function and variable order. This type of derivative has an exponential kernel, and it generalizes the classical Caputo fractional derivative. Furthermore, this derivative is called the variable order ψ -tempered Caputo fractional derivative. In this paper, we study the initial value problem of a nonlinear differential equation with the abovementioned fractional derivative. Based on the tempered Caputo fractional derivative, its application to differential equations, and some known results for this type of derivative, we introduce approximate solutions of the given problem and study existence and Ulam-type stability. Several examples are provided to illustrate the main results. Furthermore, as a partial case, we obtain existence and Ulam-type stability results for differential equations with a piecewise constant order ψ -tempered Caputo fractional derivative.

2. Ψ -Tempered Fractional Calculus of Variable Order

There are several definitions for fractional integrals and derivatives of variable order depending on the time variable (see, for example, [5]).
In this paper, we combine the ideas of variable order fractional integrals and derivatives (see [5]) with Ψ -tempered fractional derivatives/integrals to generalize fractional calculus in variable order.
Let 0 < T < be a fixed number; ψ : [ 0 , T ] ( 0 , ) be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , T ] ; λ > 0 be a given constant; the function δ : [ 0 , T ] ( 0 , 1 ) be locally integrable; and υ A C ψ ( [ 0 , T ] , R ) where
A C ψ ( [ 0 , T ] , R ) = { h : [ 0 , T ] R : h ( t ) ψ ( t ) A C ( [ 0 , T ] , R ) } .
Definition 1.
The tempered Riemann–Liouville fractional integral of variable order with respect to the function ψ (TFIVO) is defined by (here t ( 0 , T ] )
0 I ψ ( t ) δ ( t ) , λ υ ( t ) = 1 Γ ( δ ( t ) ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) δ ( t ) 1 ψ ( s ) υ ( s ) d s .
We will also generalize the Caputo-type fractional derivative of variable order combining the ideas of [17] and the type III Caputo fractional derivative of variable order (see Definition 2 [5]).
Definition 2.
The tempered Caputo fractional derivative of variable order with respect to the function ψ ( t ) (GFDVO) is defined by
0 C D ψ ( t ) δ ( t ) , λ υ ( t ) = 0 I ψ ( t ) 1 δ ( t ) , λ 1 ψ ( t ) d d t + λ υ ( t ) = 1 Γ ( 1 δ ( t ) ) ( 0 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) δ ( t ) ν ( s ) d s + λ 0 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) δ ( t ) ψ ( s ) ν ( s ) d s ) , t ( 0 , T ] .
Remark 1.
The tempered derivatives/integrals with respect to another function, given in Definitions 1 and 2, are called Ψ-tempered derivatives/integrals.
Remark 2.
TFIVO and GFDVO are generalizations of the well-known and studied fractional integrals and derivatives (see, for example, [17,18,19]).

2.1. Some Results on Caputo-Type Fractional Derivatives of Constant Order

In our paper, we will use some known results for tempered Caputo fractional derivative with reference to another function and a constant order.
Lemma 1
(Theorem 5.2 [18]). Let α ( 0 , 1 ) be a given constant, λ R , the function ψ : [ 0 , b ] ( 0 , ) , b < be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , b ] , and f C ( [ 0 , b ] × R , R ) . Then, the IVP for the fractional differential equation with the ψ-tempered Caputo fractional derivative of constant order
0 C D ψ ( t ) α , λ y ( t ) = f ( t , y ( t ) ) , t ( 0 , b ] , y ( 0 ) = y 0 ,
is equivalent to the fractional integral equation
y ( t ) = e λ ( ψ ( t ) ψ ( 0 ) ) y 0 + 1 Γ ( α ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) α 1 ψ ( s ) f ( s , y ( s ) ) d s , t [ 0 , b ] ,
both to be solved for functions y A C ψ ( [ 0 , T ] , R ) .
Lemma 2
(Lemma 5.4 [18]). Let α ( 0 , 1 ) , λ R , ψ be a smooth monotonic function on [0, b] with ψ > 0 almost everywhere, the function f C ( [ 0 , b ] × R , R ) with 0 < b < and
| f ( t , u ) f ( t , v ) | L | u v | , u , v R , t [ 0 , b ] .
Then, the IVP (1) has an unique solution y A C ψ ( [ 0 , b ] , R ) , provided that L < Γ ( α + 1 ) ( ψ ( b ) ψ ( 0 ) ) α e λ ( ψ ( b ) ψ ( 0 ) ) .
Lemma 3
(Corollary 2 [20]). Let α > 0 , b ( 0 , ) , ψ : [ 0 , T ] ( 0 , ) ) be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , T ] , C , K > 0 , u C ( [ 0 , b ] , [ 0 , ) ) . If
u ( t ) C + K 0 t ( ψ ( t ) ψ ( s ) ) α 1 ψ ( s ) u ( s ) d s , t ( 0 , b ] ,
then we have
u ( t ) C E α ( K Γ ( α ) ( ψ ( t ) ψ ( 0 ) ) α ) , t ( 0 , b ] .
We will use the following result for the Caputo fractional derivative with respect to another function (a special case of GFDVO with λ = 0 and δ ( t ) α ( 0 , 1 ) being a constant).
Lemma 4
(Lemma 2 [21]). Let μ R , α ( 0 , 1 ) and ψ : [ 0 , T ] ( 0 , ) be a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , T ] . Then, the function E α ( μ ( ψ ( t ) ψ ( 0 ) ) α ) is a solution of the equation 0 C D ψ ( t ) α y ( t ) = μ y ( t ) where 0 C D ψ ( t ) α is the Caputo fractional derivative with reference to another function given in Definition 2 with λ = 0 and δ ( t ) α .
Lemma 5
(Lemma 1 [21]). Given β > 1 , we have
0 C D ψ ( t ) α ( ψ ( t ) ψ ( 0 ) ) β 1 = Γ ( β ) Γ ( β α ) ( ψ ( t ) ψ ( 0 ) ) β α 1 .

2.2. Some Results for Variable Order of Fractional Derivatives

For a general function δ ( t ) , the ψ -tempered Riemann–Liouville variable order fractional integral does not have the semigroup property.
Example 1.
Let T = 2 , γ ( t ) = 0.75 t + 0.15 , δ ( t ) = 0.25 ( t + 1 ) , λ = 0 , ψ ( t ) = t and υ ( t ) 1 .
We will compare 0 I ψ ( t ) γ ( t ) , λ 0 I ψ ( t ) δ ( t ) , λ υ ( t ) | t = 1 and 0 I ψ ( t ) γ ( t ) + δ ( t ) , λ υ ( t ) | t = 1 applied in Definition 1. From Definition 1 we have
0 I ψ ( t ) γ ( t ) + δ ( t ) , λ υ ( t ) = 1 Γ ( γ ( t ) + δ ( t ) ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) γ ( t ) + δ ( t ) 1 ψ ( s ) υ ( s ) d s = 1 Γ ( γ ( t ) + δ ( t ) ) 0 t ψ ( t ) ψ ( s ) γ ( t ) + δ ( t ) 1 ψ ( s ) d s
and
0 I ψ ( t ) γ ( t ) , λ 0 I ψ ( t ) δ ( t ) , λ υ ( t ) = 1 Γ ( γ ( t ) ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) γ ( t ) 1 ψ ( s ) 0 I ψ ( s ) δ ( s ) , λ υ ( s ) d s = 1 Γ ( γ ( t ) ) 0 t 1 Γ ( δ ( s ) ) t s γ ( t ) 1 0 s t s δ ( s ) 1 d σ d s = 1 Γ ( γ ( t ) ) 0 t Γ ( 0.25 + 0.25 s ) Γ ( 1.25 + 0.25 s ) Γ ( 0.25 s + 0.25 ) t s γ ( t ) 1 s 0.25 s + 0.25 d s .
Note
0 I ψ ( t ) γ ( t ) + δ ( t ) , λ υ ( t ) | t = 1 0.805043 ,
0 I ψ ( t ) γ ( t ) , λ 0 I ψ ( t ) δ ( t ) , λ υ ( t ) | t = 1 1.16595 Γ ( 0.9 ) = 1.09107 .
Thus, 0 I ψ ( t ) γ ( t ) , λ 0 I ψ ( t ) δ ( t ) , λ υ ( t ) | t = 1 0 I ψ ( t ) γ ( t ) + δ ( t ) , λ υ ( t ) | t = 1 .
Remark 3.
Without the semigroup property of the ψ-tempered Riemann–Liouville variable order fractional integral (see Example 1), the integral presentation of Lemma 1 is not true for the case when the order is variable (see, for example, Equation (7) [22]).
In our study, we will use the ideas of Lemma 6 [7], which will be slightly modified:
Lemma 6.
Let δ C ( [ 0 , T ] , ( 0 , 1 ) ) . Then, for any ϵ > 0 there exist a natural number m = m ( ϵ ) and points T k = T k ( ϵ ) : 0 = T 0 < T 1 < T 2 < < T m 1 < T m = T such that
| δ ( t ) α ( t ) | ϵ , t [ 0 , T ] ,
where α ( t ) = k = 1 m δ ( T k 1 ) I k ( t ) for t [ 0 , T ] , and I k ( t ) is the indicator of the interval [ T k 1 , T k ] , k = 1 , 2 , , m , i.e.,
I k ( t ) = 1 , i f t ( T k 1 , T k ] , 0 , o t h e r w i s e .
Proof. 
Let ϵ > 0 be a given fixed number. Note that δ is uniformly continuous on [ 0 , T ] , so there exists a θ > 0 with | δ ( t ) δ ( s ) | ϵ for t , s [ 0 , T ] and | t s | θ .
Take p 1 = δ ( 0 ) . If | δ ( t ) p 1 | ϵ for t ( 0 , T ] , then we choose m = 1 and T 1 = T and stop the process. Otherwise, there exists a number T 1 ( 0 , T ] (note T 1 θ ) such that | δ ( t ) p 1 | ϵ , t [ 0 , T 1 ] . Take p 2 = δ ( T 1 ) . If | δ ( t ) p 1 | ϵ for t ( T 1 , T ] , then we choose m = 2 , T 2 = T and stop the process. Otherwise, there exists a number T 2 ( T 1 , T ] (note T 2 T 1 θ ) such that | δ ( t ) p 2 | ϵ , t ( T 1 , T 2 ] . Take p 3 = δ ( T 2 ) and continue this process (the process will stop after a finite number since there exists a positive integer l with l θ T ). Thus, we have constructed (for some m) the function
α ( t ) = k = 1 m p k I k ( t ) = δ ( T 0 ) , i f t [ T 0 , T 1 ] , δ ( T 1 ) , i f t ( T 1 , T 2 ] , δ ( T m ) , i f t ( T m 1 , T m ] .
Definition 3.
The piecewise constant function α defined by (4) will be called an ϵ-approximation of the variable order δ ( t ) .
Remark 4.
Note that the partition as well as the piecewise function α, defined by (4), depend on the number ϵ.
Example 2.
Let T = 2 , δ ( t ) = 0.5 sin ( t ) + 0.4 ( 0 , 1 ) .
Let ϵ = 0.2 . We will construct an 0.2 -approximation of the fractional order δ ( t ) . Since | δ ( t ) δ ( 0 ) | = | 0.5 sin ( t ) | 0.2 for t ( 0 , sin 1 ( 0.2 0.5 ) = sin 1 ( 0.4 ) ] we take T 1 = sin 1 ( 0.4 ) ( 0 , 2 ) . Then, from | δ ( t ) δ ( T 1 ) | = | 0.5 sin ( t ) 0.5 ( 0.4 ) | 0.2 for t ( T 1 , T 2 ) with T 2 = sin 1 ( 0.8 ) ( T 1 , 2 ) . Furthermore, | δ ( t ) δ ( T 2 ) | = | 0.5 sin ( t ) 0.5 ( 0.8 ) | 0.2 for t [ T 2 , 2 ] .
Therefore, the partition is 0 < s i n 1 ( 0.4 ) < sin 1 ( 0.8 ) < 2 , and the 0.2 -approximation of δ ( t ) is
α ( t ) = 0.4 , i f t [ 0 , s i n 1 ( 0.4 ) ] , 0.6 , i f t ( s i n 1 ( 0.4 ) , s i n 1 ( 0.8 ) ] , 0.8 , i f t ( s i n 1 ( 0.8 ) , 2 ] .
Let ϵ = 0.1 . Since | δ ( t ) δ ( 0 ) | = | 0.5 sin ( t ) | 0.1 for t ( 0 , sin 1 ( 0.1 0.5 ) = sin 1 ( 0.2 ) ] we take T 1 = sin 1 ( 0.2 ) ( 0 , 3 ) . Then, from | δ ( t ) δ ( T 1 ) | = | 0.5 sin ( t ) 0.5 ( 0.2 ) | 0.1 for t ( T 1 , T 2 ) with T 2 = sin 1 ( 0.2 0.5 ) = sin 1 ( 0.4 ) ( T 1 , 2 ) . Furthermore, | δ ( t ) δ ( T 2 ) | = | 0.5 sin ( t ) 0.5 ( 0.4 ) | 0.1 for t [ T 2 , T 3 ] with T 3 = sin 1 0.3 0.5 = sin 1 ( 0.6 ) . Furthermore, | δ ( t ) δ ( T 3 ) | = | 0.5 sin ( t ) 0.5 ( 0.6 ) | 0.1 for t [ T 3 , T 4 ] with T 4 = s i n 1 0.4 0.5 = sin 1 0.8 .
Therefore, the partition is 0 < sin 1 ( 0.2 ) < sin 1 ( 0.4 ) < sin 1 ( 0.6 ) < sin 1 ( 0.8 ) < 2 and the 0.1 -approximation of δ ( t ) is
α ( t ) = 0.4 , i f t [ 0 , s i n 1 ( 0.2 ) ] , 0.5 , i f t ( s i n 1 ( 0.2 ) , s i n 1 ( 0.4 ) ] , 0.6 , i f t ( s i n 1 ( 0.4 ) , s i n 1 ( 0.6 ) ] , 0.7 , i f t ( s i n 1 ( 0.6 ) , s i n 1 ( 0.8 ) ] , 0.8 , i f t ( s i n 1 ( 0.8 ) , 2 ] .
The example illustrates the ϵ-approximation of the fractional order depends on ϵ.
Note that the claim of Lemma 6 could be proved on a half real line.
Lemma 7.
Let δ C ( [ 0 , ) , ( 0 , 1 ) ) be such that lim t δ ( t ) = ξ ( 0 , 1 ) . Then, for any ϵ > 0 , there exist a natural number m = m ( ϵ ) and points T k = T k ( ϵ ) : 0 < T 1 < T 2 < < T m + 1 < such that
| δ ( t ) α ( t ) | < ϵ , t [ 0 , ) ,
where α ( t ) = k = 1 m + 1 δ ( T k 1 ) I k ( t ) + ξ I T ( t ) for t [ 0 , T ] , I k ( t ) is the indicator of the interval [ T k 1 , T k ] , k = 1 , 2 , , m + 1 , and I T ( t ) is the indicator of the interval ( T m + 1 , ) (here T 0 = 0 ).
Proof. 
Let ϵ > 0 be a given fixed number. From lim t δ ( t ) = ξ ( 0 , 1 ) it follows that there exists a number T > 0 such that
| δ ( t ) ξ | < ϵ , t > T .
Consider the finite interval [ 0 , T ] and apply Lemma 6 to it. Then, the piecewise constant function
α ( t ) = k = 1 m p k I k ( t ) = δ ( T 0 ) , i f t [ T 0 , T 1 ] , δ ( T 1 ) , i f t ( T 1 , T 2 ] , δ ( T m ) , i f t ( T m , T m + 1 ] , ξ , i f t > T m + 1 ,
is the function needed. □

3. Differential Equations with Variable Order of the Fractional Derivative

Consider the initial value problem (IVP) for the nonlinear scalar differential equation with the ψ -tempered Caputo fractional derivative of variable order
0 C D ψ ( t ) δ ( t ) , ϱ ω ( t ) = F t , ω ( t ) , t ( 0 , T ] , ω ( 0 ) = V 0 ,
where V 0 R , ϱ > 0 , 0 < T < , δ ( . ) : [ 0 , T ] ( 0 , 1 ) , 0 C D ψ ( t ) δ ( t ) , ϱ denotes the ψ -tempered Caputo fractional derivatives with variable order δ ( · ) , and F : [ 0 , T ] × R R is a given function.

Definition of Approximate Solutions of the Initial Value Problem (8)

We will assume ϵ > 0 is a fixed number and δ C ( [ 0 , T ] , ( 0 , 1 ) ) .
According to Lemma 6 for the given fixed number ϵ , there exists a partition P = { ( T k 1 , T k ] , k = 1 , 2 , , m } of the given interval [ 0 , T ] and an ϵ -approximation of δ , defined by (4) with | δ ( t ) α ( t ) | ϵ , t [ 0 , T ] . Note that for the given number ϵ , the partition P is not unique. In this section, we will consider a fixed partition P, i.e., fixed points 0 = T 0 < T 1 < T 2 < T m 1 < T m = T .
Let t ( T k 1 , T k ] , where k is an arbitrary integer 1 k m . Then, α ( t ) = p k with p k = δ ( T k 1 ) , and the fractional derivative can be written
0 C D ψ ( t ) α ( t ) , λ ω ( t ) = 1 Γ ( 1 α ( t ) ) ( 0 t e λ ( ψ ( t ) ψ ( s ) ) ψ ( t ) ψ ( s ) α ( t ) ω ( s ) d s + λ 0 t e λ ( ψ ( t ) ψ ( s ) ) ψ ( t ) ψ ( s ) α ( t ) ψ ( s ) ω ( s ) d s ) = 1 Γ ( 1 p k ) 0 t e λ ( ψ ( t ) ψ ( s ) ) ψ ( t ) ψ ( s ) p k 1 ψ ( s ) d d t + λ ω ( s ) d s = 0 C D ψ ( t ) p k , λ ω ( t ) .
Applying (9), we could define for any integer k : 1 k m , the following IVP for the differential equation with tempered Caputo fractional derivative of a constant order p k , which is deeply connected to the studied IVP (8) and the ϵ -approximation of the variable order δ ( t ) :
0 C D ψ ( t ) p k , λ ν ( t ) = F t , ν ( t ) f o r t ( 0 , T k ] , ν ( 0 ) = V 0 .
With the help of the IVP (10), we will define an approximate solution of (8).
Definition 4.
For a given ϵ > 0 the function
ω ( t ) = ν 1 ( t ) , i f t [ T 0 , T 1 ] , ν 2 ( t ) , i f t ( T 1 , T 2 ] , ν m ( t ) , i f t ( T m , T m + 1 ] ,
is called an ϵ-approximate solution of IVP (8), where ν k ( t ) C ( [ 0 , T k ] , R ) is a solution of IVP (10), k = 1 , 2 , , m .
Example 3.
Consider the following linear IVP:
0 C D ψ ( t ) δ ( t ) ω ( t ) = 0.2 ω ( t ) , t [ 0 , 2 ] , ω ( 0 ) = V 0 ,
where δ ( t ) = 0.5 sin ( t ) + 0.4 and 0 C D ψ ( t ) δ ( t ) is the variable order Caputo fractional derivative with reference to another function given in Definition 2 with λ = 0 .
Let ϵ = 0.2 . Then, according to Example 2, the 0.2 -partition of the fractional order δ ( t ) is given by (5), and the 0.2 -approximate solution of (12) according to Lemma 4 and Definition 4 is
ω 0.2 ( t ) = V 0 E 0.4 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.4 ) , i f t [ 0 , sin 1 ( 0.4 ) ] , V 0 E 0.6 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.6 ) , i f t ( sin 1 ( 0.4 ) , sin 1 ( 0.8 ) ] , V 0 E 0.8 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.8 ) , i f t ( sin 1 ( 0.8 ) , 2 ] .
Let ϵ = 0.1 . Then, according to Example 2 the 0.1 -partition of the fractional order δ ( t ) is given by (6) and the 0.1 -approximate solution of (12) according to Lemma 4 and Definition 4 is
ω 0.1 ( t ) = V 0 E 0.4 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.4 ) , i f t [ 0 , sin 1 ( 0.2 ) ] , V 0 E 0.5 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.5 ) , i f t ( sin 1 ( 0.2 ) , sin 1 ( 0.4 ) ] , V 0 E 0.6 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.6 ) , i f t ( sin 1 ( 0.4 ) , sin 1 ( 0.6 ) ] , V 0 E 0.7 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.7 ) , i f t ( sin 1 ( 0.6 ) , sin 1 ( 0.8 ) ] , V 0 E 0.8 ( 0.2 ( ψ ( t ) ψ ( 0 ) ) 0.8 ) , i f t ( sin 1 ( 0.8 ) , 2 ] .
Remark 5.
If the piecewise constant function α ( t ) is an ϵ-approximation of the variable order δ ( t ) ; ν k ( t ) C ( [ 0 , T k ] , R ) , k = 1 , 2 , , m are solutions of the initial value problems; (10) and ω ( t ) is an ϵ-approximate solution of (8), then, according to (9), we have 0 C D ψ ( t ) α ( t ) , λ ω ( t ) = 0 C D ψ ( t ) p k , λ ν k ( t ) = F ( t , ν k ( t ) ) = F ( t , ω ( t ) ) for any t ( T k 1 , T k ] with p k = δ ( T k 1 ) , k = 1 , 2 , , m .
We now introduce condition (H):
Hypothesis 1
(H1). The function ψ : [ 0 , T ] ( 0 , ) is a smooth increasing function with ψ ( t ) > 0 almost everywhere in [ 0 , T ] , λ > 0 , and the function δ C ( [ 0 , T ] , ( 0 , 1 ) ) .
Hypothesis 2
(H2). The function F C ( [ 0 , T ] × R , R ) , and it is globally Lipschitz with a constant L such that
L e λ ( ψ ( T ) ψ ( 0 ) ) max 1 , ( ψ ( T ) ψ ( 0 ) ) p < A ,
where p = max t [ 0 , T ] δ ( t ) ( 0 , 1 ) , A = min p [ 0 , 1 ] Γ ( 1 + p ) > 0 .
Theorem 1.
Let ϵ > 0 be a given number, condition (H1) be satisfied; P = { ( T k 1 , T k ] , k = 1 , , m } be the the partition of the interval [ 0 , T ] , defined in Lemma 6 for ϵ; the function F C ( [ 0 , T ] × R , R ) ; and there exist constants L k > 0 , k = 1 , 2 , , m , such that
| F ( t , y ) F ( t , z ) | L k | y z | , t [ 0 , T k ] , y , z R
with
L k < Γ ( 1 + p k ) ( ψ ( T k ) ψ ( 0 ) ) p k e λ ( ψ ( T k ) ψ ( 0 ) ) ) , k = 1 , 2 , , m .
Then, IVP (8) has an ϵ-approximate solution.
Proof. 
According to Lemma 2 with L = L k , b = T k , f = F , the initial value problem with the fractional derivative of constant order (10) has a solution ν k A C ψ ( [ 0 , T k ] , R ) for all k = 1 , 2 , , m . Then, the function ω ( t ) defined by (11) is an ϵ -approximate solution of the IVP (8). □
We now define an approximate solution of (8).
Definition 5.
The initial value problem (8) has an approximate solution if for any ϵ > 0 it has an ϵ-approximate solution.
Remark 6.
The approximate solution of (8) depends on the initial value.
Theorem 2.
Let condition (H) be satisfied. Then, IVP (8) has an approximate solution.
Proof. 
Let ϵ be an arbitrary given number. According to Lemma 6 for the given number ϵ , there exist a partition P of the interval [ 0 , T ] and a piecewise constant function α ( t ) defined by (4) which is an ϵ -approximation of the variable order δ ( t ) .
Then, for any k = 1 , 2 , , m the inequalities
L ( ψ ( T k ) ψ ( 0 ) ) p k Γ ( 1 + p k ) L ( ψ ( T ) ψ ( 0 ) ) p k A
hold.
Furthermore,
-
for ψ ( T ) ψ ( 0 ) > 1 we have ( ψ ( T ) ψ ( 0 ) ) p k < ( ψ ( T ) ψ ( 0 ) ) p ;
-
for ψ ( T ) ψ ( 0 ) 1 we have ( ψ ( T ) ψ ( 0 ) ) p k 1 .
Thus, e λ ( ψ ( T k ) ψ ( 0 ) ) ( ψ ( T ) ψ ( 0 ) ) p k e λ ( ψ ( T ) ψ ( 0 ) ) max 1 , ( ψ ( T ) ψ ( 0 ) ) p . Therefore, L ( ψ ( T k ) ψ ( 0 ) ) p k Γ ( 1 + p k ) e λ ( ψ ( T k ) ψ ( 0 ) ) < 1 , and, according to Theorem 1, the IVP (8) has an ϵ -approximate solution (the function ω ϵ ( t ) defined by (11)). Since ϵ is an arbitrary number, we are finished. □

4. Ulam-Type Stability of Approximate Solutions

Note that the Ulam-type stability for differential equations with the Caputo-tempered fractional derivative with respect to another function is studied in [5] for the case of constant order. In the case of constant order, the fractional differential equation is equivalent to a fractional integral Equation (see Lemma 1), but it is not true for the variable order. We will now study the stability of approximate solutions of (8). In this case, we change the definitions for Ulam-type stability given in [5] (Definition 5.1).
Definition 6.
Let condition (H1) be satisfied, and F C ( [ 0 , T ] × R , R ) . Now, (8) is said to be approximately Hyers–Ulam-stable if there exists a constant C 0 > 0 such that for any ϵ > 0 , an ϵ-approximation α ( t ) of the fractional order δ ( t ) (defined by (4)) with a partition P = { ( T k 1 , T k ] , k = 1 , , m } and for any function η k A C ψ ( [ 0 , T k ] , R ) satisfying the inequality
| 0 C D ψ ( t ) p k , λ η k ( t ) F ( t , η k ( t ) ) | ϵ , t [ 0 , T k ] , k = 1 , 2 , , m ,
with p k = δ ( T k 1 ) , k = 1 , 2 , , m , there exists an ϵ-approximate solution ω ϵ ( t ) of (8) with
| ω ϵ ( t ) η ( t ) | C 0 ϵ , t [ 0 , T ] ,
where
η ( t ) = η 1 ( t ) , i f t [ T 0 , T 1 ] , η 2 ( t ) , i f t ( T 1 , T 2 ] , η m ( t ) , i f t ( T m , T m + 1 ] .
Theorem 3
(UHS). Let condition (H) be satisfied. Then, (8) is approximately Hyers–Ulam-stable.
Proof. 
Let ε > 0 be an arbitrary number. According to Lemma 6 there exist a partition P = { ( T k 1 , T k ] , k = 1 , , m } and an ϵ -approximation α ( t ) of the fractional order δ ( t ) , defined by (4). Let the function η k A C ψ ( [ 0 , T k ] , R ) be a solution of inequality (15), k = 1 , 2 , , m . Then, for any k = 1 , 2 , , m we consider the functions g k C ( [ 0 , T k ] , R ) : | g k ( t ) | ϵ such that
0 C D ψ ( t ) p k , λ η k ( t ) = F ( t , η k ( t ) ) + g k ( t ) , t [ 0 , T k ] .
According to Lemma 1, the integral equality
η k ( t ) = e λ ( ψ ( t ) ψ ( 0 ) ) η k ( 0 ) + 1 Γ ( p k ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) p k 1 ψ ( s ) ( F ( s , η k ( s ) ) + g k ( s ) ) d s , = e λ ( ψ ( t ) ψ ( 0 ) ) η k ( 0 ) + 0 I ψ ( t ) p k , λ ( F ( s , η k ( s ) ) + g k ( s ) ) , t [ 0 , T k ] ,
holds.
Denote η 0 = min k = 1 , 2 , m η k ( 0 ) and consider the initial value problem (8) with V 0 R : V 0 < η 0 + ϵ . According to Theorem 1 the initial value problem (8) has an ϵ -approximate solution ω ϵ ( t ) .
Let k : 1 k m be an arbitrary fixed integer. Then, for t ( T k 1 , T k ] the equality ω ϵ ( t ) = ν k ( t ) holds where ν k ( t ) , t [ 0 , T k ] , is the solution of (10) with V 0 = η 0 .
Then, for any t [ 0 , T k ] , we obtain
| ν k ( t ) η k ( t ) | | η k ( t ) e λ ( ψ ( t ) ψ ( 0 ) ) η k ( 0 ) 0 I ψ ( t ) p k , λ F ( s , η k ( s ) ) | + | ν k ( t ) e λ ( ψ ( t ) ψ ( 0 ) ) V 0 0 I ψ ( t ) p k , λ F ( s , ν k ( s ) ) | + e λ ( ψ ( t ) ψ ( 0 ) ) | V 0 η k ( 0 ) | + | 0 I ψ ( t ) p k , λ ( F ( s , η ( s ) ) F ( s , ν k ( s ) ) ) | | 0 I ψ ( t ) p k , λ g k ( t ) | + e λ ( ψ ( t ) ψ ( 0 ) ) | V 0 η k ( 0 ) | + L Γ ( p k ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) p k 1 ψ ( s ) | η k ( s ) ) ν k ( s ) ) | d s ϵ 1 Γ ( p k ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) p k 1 ψ ( s ) d s + ϵ + L Γ ( p k ) 0 t e λ ψ ( t ) ψ ( s ) ψ ( t ) ψ ( s ) p k 1 ψ ( s ) | η k ( s ) ) ν k ( s ) ) | d s ϵ Γ ( 1 + p k ) ( ψ ( t ) ψ ( 0 ) ) p k + ϵ + L Γ ( p k ) 0 t ψ ( t ) ψ ( s ) p k 1 ψ ( s ) | η k ( s ) ) ν k ( s ) ) | d s ϵ 1 + ( ψ ( T k ) ψ ( 0 ) ) p k Γ ( 1 + p k ) + L Γ ( p k ) 0 t ψ ( t ) ψ ( s ) p k 1 ψ ( s ) | η k ( s ) ) ν k ( s ) ) | d s .
Denote M = 1 + max 1 , ( ψ ( T ) ψ ( 0 ) ) p A . We have ( ψ ( T k ) ψ ( 0 ) ) p k max 1 , ( ψ ( T ) ψ ( 0 ) ) p (see the proof of Theorem 2).
Therefore, for any t [ 0 , T k ] we have
| ν k ( t ) η k ( t ) | ϵ M + L Γ ( p k ) 0 t ψ ( t ) ψ ( s ) p k 1 ψ ( s ) | η k ( s ) ) ν k ( s ) ) | d s .
According to Lemma 3 with u ( t ) = | ν k ( t ) η k ( t ) | C [ 0 , T k ] , [ 0 , ) ) , C = ϵ M , K = L Γ ( p k ) , α = p k , b = T k and inequalities E γ ( x γ ) < e x γ , x > 0 , γ ( 0 , 1 ) , and L p k ( ψ ( T ) ψ ( 0 ) ) < A p k A p because A ( 0 , 1 ] (see, condition (H)), we obtain
| ν k ( t ) η k ( t ) | ϵ M E p k ( L ( ψ ( t ) ψ ( 0 ) ) p k ) ϵ M E p k ( ( L p k ( ψ ( T ) ψ ( 0 ) ) ) p k ) ϵ M e L p k ( ψ ( T ) ψ ( 0 ) ) ξ ϵ M e A p ξ , t [ 0 , T k ] ,
where ξ = min t [ 0 , T ] δ ( t ) > 0 .
From inequality (21), it follows that
| ω ϵ ( t ) η ( t ) | = | ν k ( t ) η k ( t ) | C 0 ϵ , t ( T k 1 , T k ] , k = 1 , 2 , , m ,
with C 0 = M e A p ξ .
We will illustrate the approximate Hyers–Ulam stability with a simple example.
Example 4.
Consider the linear fractional differential Equation (12) with variable order δ ( t ) = 0.5 sin ( t ) + 0.4 , ψ ( t ) = t 2 , and 0 C D ψ ( t ) δ ( t ) -the variable order Caputo fractional derivative with reference to another function given in Definition 2 with λ = 0 .
The assumptions (H1) and (H2) are satisfied with L = 0.2 because in this case L e λ ( ψ ( T ) ψ ( 0 ) ) max { 1 , ( ψ ( T ) ψ ( 0 ) ) p } = 0.2 ( 4 ) = 0.8 < min p ( 0 , 1 ) Γ ( 1 + p ) = 0.8872 . According to Theorem 3 Equation (12) is approximately Hyers–Ulam-stable. We will illustrate this for a particular value of ϵ. For example, let ϵ = 0.2 .
In Example 3 several ϵ-approximate solutions of (12) are given.
Then, the inequality (15) is reduced to
| 0 C D ψ ( t ) p k η k ( t ) + 0.2 η k ( t ) | 0.2 , t [ 0 , T k ] , k = 1 , 2 , 3 ,
with T 1 = sin 1 ( 0.4 ) , T 2 = sin 1 ( 0.8 ) , T 3 = 2 and p 1 = 0.4 , p 2 = 0.6 , p 3 = 0.8 (see Example 3).
Then, according to Lemma 5 with β = 2 , ψ ( t ) = t 2 , α = 0.4 , 0.6 , 0.8 , we obtain (see Figure 1, Figure 2 and Figure 3)
η 1 ( t ) = 0.3 ( t 2 ) 2 1 , t [ 0 , sin 1 ( 0.4 ) ] , η 2 ( t ) = 0.1 ( t 2 ) 2 1 , t [ 0 , sin 1 ( 0.8 ) ] , η 3 ( t ) = 0.05 ( t 2 ) 2 1 , t [ 0 , 2 ] ,
According to (13), the 0.2 -approximate solution is
ω 0.2 ( t ) = V 0 E 0.4 ( 0.2 ( t 2 ) 0.4 ) , i f t [ 0 , sin 1 ( 0.4 ) ] , V 0 E 0.6 ( 0.2 ( t 2 ) 0.6 ) , i f t ( sin 1 ( 0.4 ) , sin 1 ( 0.8 ) ] , V 0 E 0.8 ( 0.2 ( t 2 ) 0.8 ) , i f t ( sin 1 ( 0.8 ) , 2 ] ,
where V 0 = min ( 0.3 , 0.1 , 0.05 ) + 0.2 .
Then, | ω 0.2 ( t ) η ( t ) | < C 0 0.2 with C 0 = 2 and η ( t ) defined by (16) with m = 2 (see Figure 4).

5. Piecewise Constant Order of the Fractional Derivative

Consider the case when the partition 0 = T 0 < T 1 < < T m 1 < T m = T of the interval [ 0 , T ] is given initially, and the fractional order of the fractional derivative is defined by the equality
δ ( t ) = p 1 , i f t [ T 0 , T 1 ] , p 2 , i f t ( T 1 , T 2 ] , p m , i f t ( T m 1 , T m ] ,
where the constants p k ( 0 , 1 ) , k = 1 , 2 , , m are initially given.
Let t ( T k 1 , T k ] , where k is an arbitrary integer 1 k m . Then, δ ( t ) = p k and from Definition 2 we have
0 C D ψ ( t ) δ ( t ) , λ υ ( t ) = 1 Γ ( 1 p k ) ( 0 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ν ( s ) d s + λ 0 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ψ ( s ) ν ( s ) d s ) = 1 Γ ( 1 p k ) ( i = 1 k 1 T i 1 T i e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ν ( s ) d s T k 1 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ν ( s ) d s + λ i = 1 k 1 T i 1 T i e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ψ ( s ) ν ( s ) d s ) + λ T k 1 t e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ψ ( s ) ν ( s ) d s ) = T k 1 C D ψ ( t ) p 2 , λ υ ( t ) + 1 Γ ( 1 p k ) i = 1 k 1 ( T i 1 T i e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ν ( s ) d s + λ T i 1 T i e λ ( ψ ( t ) ψ ( s ) ) ( ψ ( t ) ψ ( s ) ) p k ψ ( s ) ν ( s ) d s ) , t ( T k 1 , T k ] .
Therefore,
0 C D ψ ( t ) δ ( t ) , λ υ ( t ) i = 1 k 1 T i 1 C D ψ ( t ) p i , λ υ ( t ) | t = T i + T k 1 C D ψ ( t ) p k , λ υ ( t ) , t ( T k 1 , T k ] ,
or its equivalent
0 C D ψ ( t ) δ ( t ) , λ υ ( t ) i = 1 k 1 T i 1 C D ψ ( t ) δ ( t ) , λ υ ( t ) | t = T i + T k 1 C D ψ ( t ) δ ( t ) , λ υ ( t ) , t ( T k 1 , T k ] ,
hold and their corresponding equalities cannot be applied (see, for example, Equation (6) [23], Equation (5) [16], Equation (9) [24], Equation (3.1) [25]).
In connection with the above, we will apply our approach to study IVP (8) when the order of the fractional derivative is a piecewise constant function σ ( t ) . Our approach is based on (9).
We now introduce assumption (A):
Assumption 1
(A1). The function ψ C 1 ( [ 0 , T ] , ( 0 , ) ) is an increasing function with ψ ( t ) > 0 almost everywhere, λ > 0 is a given constant, and the function δ : [ 0 , T ] ( 0 , 1 ) is defined by (26).
Assumption 2
(A2). The function F C ( [ 0 , T ] × R , R ) and there exist constants L k > 0 , k = 1 , 2 , , m , such that
| F ( t , y ) F ( t , z ) | L k | y z | , t [ 0 , T k ] , y , z R , k = 1 , 2 , , m ,
with
L k e λ ( ψ ( T k ) ψ ( 0 ) ) ( ψ ( T k ) ψ ( 0 ) ) p k < L ( 1 + p k ) .
We obtain the following existence result for the initial value problem (8) with piecewise constant order, defined by (26).
Theorem 4.
Let conditions (A1) and (A2) hold. Then, the initial value problem (8) has a solution.
Proof. 
Consider the partition P = { ( T k 1 , T k ] , k = 1 , 2 , } defined by the given point T k , k = 1 , 2 , , m 1 . For any k = 1 , 2 , , m , we consider the initial value problem (10). According to Lemma 2, it has a unique solution ν k ( t ) , t ( 0 , T k ] for any k = 1 , 2 , , m . Construct the function
ω ( t ) = ν 1 ( t ) , i f t [ T 0 , T 1 ] , ν 2 ( t ) , i f t ( T 1 , T 2 ] , ν m ( t ) , i f t ( T m 1 , T m ] .
The function ω ( t ) is a solution of (8) because ω ( 0 ) = V 0 and for any t ( T k 1 , T k ] we have
0 C D ψ ( t ) δ ( t ) , λ ω ( t ) = 0 C D ψ ( t ) p k , λ ν k ( t ) = F ( t , ν k ( t ) = F ( t , ω ( t ) ) .
Remark 7.
The solution of (8) with a fractional derivative of piecewise constant order is not continuous on the whole interval [ 0 , T ] .
The definition of Ulam-type stability for fractional derivatives of a piecewise constant order reduces to:
Definition 7.
Let condition (A1) hold and F C ( [ 0 , T ] × R , R ) . Now, (8) is said to be Hyers–Ulam-stable if there exists a constant C 0 > 0 such that for any ϵ > 0 and for any function η k C ψ ( [ 0 , T k ] , R ) , k = 1 , 2 , , m , satisfying the inequality
| 0 C D ψ ( t ) p k , λ η ( t ) F ( t , η ( t ) ) | ϵ , t [ 0 , T k ] ,
with p k = δ ( T k 1 ) , k = 1 , 2 , , m , there exists an ϵ-approximate solution ω ϵ ( t ) of (8) with
| ω ϵ ( t ) η ( t ) | C 0 ϵ , t [ 0 , T ] ,
where the function η is defined by (16).
Theorem 5
(UHS). Let condition (A1) and (A2) hold. Then, (8) is Hyers–Ulam-stable.
The proof of Theorem 5 is similar to the one of Theorem 3 without applying ϵ .

6. Conclusions

In this paper, we consider differential equations with the variable order Caputo-type fractional derivative with respect to another function and we study the general case of continuous variable order of the fractional derivative. We define in appropriate way an ϵ -solution of the given initial value problem and an approximate solution.
Furthermore, we define and study the Hyers–Ulam stability. As a partial case, we obtain results for stability of the differential equations with fractional derivatives of piecewise constant order. In future work, we hope to appropriately define some other types of Ulam stability, such as Hyers–Ulam–Rassias stability, and obtain sufficient conditions. Furthermore, the existence on the the interval [ 0 , ) could be investigated. In addition, we could consider the general case of the variable order applying a local definition for fractional derivative [14] in the future.

Author Contributions

Conceptualization, D.O., S.H. and R.P.A.; Methodology, R.P.A.; Formal analysis, D.O., S.H. and R.P.A.; Writing—original draft, D.O. and S.H.; Writing—review & editing, D.O., S.H. and R.P.A.; Funding acquisition, 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-PN62/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 conflict of interest.

References

  1. Garrappa, R.; Giusti, A.; Mainardi, F. Variable-order fractional calculus: A change of perspective. Commun. Nonlinear Sci. Numer. Simul. 2021, 102, 105904. [Google Scholar] [CrossRef]
  2. Sun, H.G.; Chang, A.; Chang, Y.; Chen, W. A review of variable-order fractional differential equations: Mathematical foundations, physical models, numebriacl methods and applications. Frac. Calc. Appl. Anal. 2019, 22, 27–59. [Google Scholar] [CrossRef]
  3. Patnaik, S.; Hollkamp, J.P.; Semperlotti, F. Applications ofvariable-order fractional operators: A review. Proc. R. Soc. A 2020, 476, 20190498. [Google Scholar] [CrossRef] [PubMed]
  4. Samko, S.; Ross, B. Integration and differentiation to a variable fractional order. Integral Transform. Spec. Funct. 1993, 4, 277–300. [Google Scholar] [CrossRef]
  5. Tavares, D.; Almeida, R.; Torres, D.F.M. Caputo derivatives of fractional variable order: Numerical approximations. Commun. Nonlinear Sci. Numer. Simul. 2016, 35, 69–87. [Google Scholar] [CrossRef]
  6. Lorenzo, C.F.; Hartley, T.T. Initialization, conceptualization, and ap plication in the generalized fractional calculus. Crit. Rev. Biomed. Eng 2007, 35, 477–553. [Google Scholar] [CrossRef]
  7. Zhang, S.; Hu, L. Unique existence result of approximate solution to initial value problem for fractional differential equation of variable order involving the derivative arguments on the half-axis. Mathematics 2019, 7, 286. [Google Scholar] [CrossRef]
  8. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  9. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  10. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Drivatives; Gordon and Breach Science Publishers: Longhorne, PA, USA, 1993. [Google Scholar]
  11. 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]
  12. 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]
  13. Agarwal, R.P.; Hristova, S. Boundary Value Problem for Multi-Term Nonlinear Delay Generalized Proportional Caputo Fractional Differential Equations. Fractal Fract. 2022, 6, 691. [Google Scholar] [CrossRef]
  14. Fedotov, A.I. Substantiation of a Quadrature-Difference Method for Solving Integro-Differential Equations with Derivatives of Variable Order. Comput. Math. Math. Phys. 2022, 62, 548–563. [Google Scholar] [CrossRef]
  15. Telli, B.; Souid, M.; Alzabut, J.; Khan, H. Existence and uniqueness theorems for a variable-order fractional differential equation with delay. Axioms 2023, 12, 339. [Google Scholar] [CrossRef]
  16. Refice, A.; Souid, M.; Yakar, A. Some qualitative properties of nonlinear fractional integro-differential equations of variable order. Intern. J. Optim. Contr. Theor. Appl. 2021, 11, 68–78. [Google Scholar] [CrossRef]
  17. Fahad, H.M.; Fernandez, A.; ur Rehman, M.; Siddiqi, M. Tempered and Hadamard-type fractional calculus with respect to functions. Mediterr. J. Math. 2021, 18, 143. [Google Scholar] [CrossRef]
  18. Mali, A.D.; Kucche, K.D.; Fernandez, A.; Fahad, H.M. On tempered fractional calculus with respect to functions and the associated fractional differential equations. Math. Meth. Appl. Sci. 2022, 45, 11134–11157. [Google Scholar] [CrossRef]
  19. Fernandez, A.; Ustaoglu, C. On some analytic properties of tempered fractional calculus. J. Comput. Appl. Math. 2020, 366, 112400. [Google Scholar] [CrossRef]
  20. Sousa, J.V.D.C.; Oliveira, E.C.D. A Gronwall inequality and the Cauchy-type problem by means of ψ-Hilfer operator. Differ. Equ. Appl. 2019, 11, 87–106. [Google Scholar]
  21. 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]
  22. Sarwar, S. On the Existence and Stability of Variable Order Caputo Type Fractional Differential Equations. Fractal Fract. 2022, 6, 51. [Google Scholar] [CrossRef]
  23. Refice, A.; Souid, M.S.; Stamova, I. On the boundary value problems of Hadamard fractional differential equations of variable order via Kuratowski MNC technique. Mathematics 2021, 9, 1134. [Google Scholar] [CrossRef]
  24. Hristova, S.; Benkerrouche, A.; Souid, M.S.; Hakem, A. Boundary value problems of Hadamard fractional differential equations of variable order. Symmetry 2021, 13, 896. [Google Scholar] [CrossRef]
  25. Benia, K.; Souid, M.S.; Jarad, F.; Alqudah, M.A.; Abdeljawad, T. Boundary value problem of weighted fractional derivative of a function with a respect to another function of variable order. J. Ineq. Appl. 2023, 2023, 127. [Google Scholar] [CrossRef]
Figure 1. Graphs of | 0 C D t 2 0.4 η 1 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [ 0 , sin 1 ( 0.4 ) ] .
Figure 1. Graphs of | 0 C D t 2 0.4 η 1 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [ 0 , sin 1 ( 0.4 ) ] .
Fractalfract 08 00011 g001
Figure 2. Graphs of | 0 C D t 2 0.6 η 2 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [ 0 , sin 1 ( 0.8 ) ] .
Figure 2. Graphs of | 0 C D t 2 0.6 η 2 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [ 0 , sin 1 ( 0.8 ) ] .
Fractalfract 08 00011 g002
Figure 3. Graphs of | 0 C D t 2 0.8 η 3 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [0,2].
Figure 3. Graphs of | 0 C D t 2 0.8 η 3 ( t ) + 0.5 η 1 ( t ) | and ϵ = 0.2 on [0,2].
Fractalfract 08 00011 g003
Figure 4. Graphs of | ω ( t ) η ( t ) | and C 0 ϵ = 2 ( 0.2 ) on [ 0 , 2 ] .
Figure 4. Graphs of | ω ( t ) η ( t ) | and C 0 ϵ = 2 ( 0.2 ) on [ 0 , 2 ] .
Fractalfract 08 00011 g004
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

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. https://doi.org/10.3390/fractalfract8010011

AMA Style

O’Regan D, Hristova S, Agarwal RP. Ulam-Type Stability Results for Variable Order Ψ-Tempered Caputo Fractional Differential Equations. Fractal and Fractional. 2024; 8(1):11. https://doi.org/10.3390/fractalfract8010011

Chicago/Turabian Style

O’Regan, Donal, Snezhana Hristova, and Ravi P. Agarwal. 2024. "Ulam-Type Stability Results for Variable Order Ψ-Tempered Caputo Fractional Differential Equations" Fractal and Fractional 8, no. 1: 11. https://doi.org/10.3390/fractalfract8010011

APA Style

O’Regan, D., Hristova, S., & Agarwal, R. P. (2024). Ulam-Type Stability Results for Variable Order Ψ-Tempered Caputo Fractional Differential Equations. Fractal and Fractional, 8(1), 11. https://doi.org/10.3390/fractalfract8010011

Article Metrics

Back to TopTop