Next Article in Journal
Multiple Critical Points for Symmetric Functionals without upper Growth Condition on the Principal Part
Next Article in Special Issue
Univalence Conditions for Gaussian Hypergeometric Function Involving Differential Inequalities
Previous Article in Journal
On the Application of a Design of Experiments along with an ANFIS and a Desirability Function to Model Response Variables
Previous Article in Special Issue
Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Value Problems of Hadamard Fractional Differential Equations of Variable Order

1
Faculty of Mathematics and Informatics, Plovdiv University, Plovdiv 4000, Bulgaria
2
Laboratory ACEDP, Djillali Liabes University, Sidi Bel Abbès 22000, Algeria
3
Department of Economic Sciences, University of Tiaret, Tiaret 14000, Algeria
*
Author to whom correspondence should be addressed.
Submission received: 22 April 2021 / Revised: 5 May 2021 / Accepted: 12 May 2021 / Published: 18 May 2021
(This article belongs to the Special Issue Functional Equations and Analytic Inequalities)

Abstract

:
A boundary value problem for Hadamard fractional differential equations of variable order is studied. Note the symmetry of a transformation of a system of differential equations is connected with the locally solvability which is the same as the existence of solutions. It leads to the necessity of obtaining existence criteria for a boundary value problem for Hadamard fractional differential equations of variable order. Also, the stability in the sense of Ulam–Hyers–Rassias is investigated. The results are obtained based on the Kuratowski measure of noncompactness. An example illustrates the validity of the observed results.

1. Introduction

The idea of fractional calculus is to take real number powers or complex number powers of the differentiation and integration operators ([1,2]). Recently, fractional integration and differentiation of variable orders have also been studied (see, for example, [3,4,5,6,7,8,9,10,11]). A good review of variable-order (VO) fractional calculus and its practical applications in the context of scientific modeling is given in [12]; some physical discussion on VO fractional integral and derivative models is provided in [13].
One of the main problems in differential equations is the solvability. There are many different methods use to study the existence such as the Lie group symmetry ([14,15]). In this paper we use the integral equivalence to prove the existence result for the boundary value problem for the Hadamard fractional differential equation of variable order. The boundary value problems for various types of fractional differential equations have been set up and studied by many authors, such as for a nonlinear fractional differential equations of order α ( 2 , 3 ) ([16]), for ψ -Hilfer fractional derivatives on b-metric spaces ([17]), for two-point boundary value problems with singular differential equations of variable order ([9]) and for impulsive integro-differential equations with Riemann–Liouville boundary conditions ([18]).
While several research studies have been performed on investigating the existence of solutions to fractional constant-order problems, the existence of solutions to variable-order problems is rarely discussed in the literature and there have been very few studies on the stability of solutions; we refer to [9,11,19,20,21,22,23]. Therefore, investigating this interesting special research topic makes our results novel and worthy.
In addition, all results in this work show a great potential to be applied in various applications of multidisciplinary sciences. Further investigations on this open research problem could be also possible with the help of our original results in this research paper. In other words, in the future one could extend the proposed BVP to other complicated real mathematical fractional models by terms of newly-introduced operators with non-singular kernels.
Recently, the Hadamard-type operators originally introduced in [24] and later generalized to variable fractional order have been investigated in [25,26]. Some existence and Ulam stability properties were studied in [27]. Bai et al. [28] studied the following initial value problem for a Caputo–Hadamard differential equation of constant order 0 < u 1
c D 0 + u x ( t ) = f ( t , x ( t ) , I 0 + u x ( t ) ) , t [ a , b ] , u ] 0 , 1 ] , x ( a ) = x a ,
where 0 < a < b < and c D 0 + u , I 0 + u stand for the Caputo–Hadamard derivative and Hadamard integral operators of constant order u, respectively, and f is a given continuous function, x a R .
In this paper we will study the boundary value problem (BVP) for the Hadamard fractional differential equation of variable order (VOHFDE)
H D 1 + u ( t ) x ( t ) = f ( t , x ( t ) , H I 1 + u ( t ) x ( t ) ) , t J , x ( 1 ) = x ( T ) = 0 ,
where J = [ 1 , T ] , 1 < T < , u ( t ) : J ( 1 , 2 ] is the variable order of the fractional derivatives, f : J × R × R R is a given function and the left Hadamard fractional integral (HFI) of variable-order u ( t ) for function x ( t ) is (see, for example, [25,26])
H I 1 + u ( t ) x ( t ) = 1 Γ ( u ( t ) ) 1 t ( ln t s ) u ( t ) 1 x ( s ) s d s , t J ,
where Γ ( z ) = 0 x z 1 e x d x is the Gamma function and the left Hadamard fractional derivative (HFD) of variable-order u ( t ) for function x ( t ) is (see, for example, [25,26])
H D 1 + u ( t ) x ( t ) = t 2 Γ ( 2 u ( t ) ) d 2 d t 2 1 t ( ln t s ) 1 u ( t ) x ( s ) s d s , t J .
The main purpose of this paper is the study of some qualitative properties of the solutions of BVPs for the VOHFDE (1). Initially some preliminary results about HFIs of variable order are provided. Further, based on the Kuratowski measure of noncompactness technique, an appropriate partition of the considered interval J is applied. The studied BVP for VOHFDE (1) is divided into appropriate BVPs for Hadamard fractional differential equations of constant orders. The existence result for (1) is proved.
The relevance between the concepts of stability and symmetry in differential equation is an old question studied by many authors for various types of differential equations such as for irreversible thermodynamics in [29], time-reversal symmetry in nonequilibrium statistical mechanics ([30]). The interplay of symmetries and stability in the analysis and control of nonlinear dynamical systems and networks is studied in [31]. Consequently, this question arises in differential equations with fractional derivatives as a topical problem. It is known that the stability of solutions is one of the most important properties of functional differential equations. In this paper we study the stability in the sense of Ulam–Hyers–Rassias of (1). An example is provided to illustrate the theoretical results.

2. Preliminaries

Denote by C ( J , R ) the Banach space of continuous functions ϰ : J R with the norm ϰ = s u p { | ϰ ( t ) | : t J } .
For any number δ [ 0 , 1 ] we define the set
C δ ( J , R ) = { h : h ( . ) C ( J , R ) , ( ln ( . ) ) δ h ( . ) C ( J , R ) } .
Remark 1.
In (2) the variable order u : J ( 1 , 2 ] but the HFI could be defined for any u : J ( 0 , ) .
Remark 2.
In the case of a constant order u in Equations (2) and (3), the HFI and HFD coincide with the standard Hadamard derivative and standard Hadamard integral, respectively (see [2,25,26]).
Remark 3.
Note that the semigroup property is satisfied for a standard Hadamard integral with constant orders but it is not fulfilled for the general case of variable orders u ( t ) , v ( t ) , i.e., H I 1 + u ( t ) ( H I 1 + v ( t ) ) x ( t ) H I 1 + u ( t ) + v ( t ) x ( t ) .
Example 1.
Let J = [ 1 , 2 ] and the function x ( t ) 1 for t J . Consider the following functions as orders of HFI: v ( t ) 2 and u ( t ) = t for t J .
Then for any t J we obtain for the HFI defined by (2)
H I 1 + v ( t ) h ( t ) = 1 Γ ( 1 ) 1 t ( ln t s ) 1 1 1 s d s = ln t ,
H I 1 + u ( t ) H I 1 + v ( t ) x ( t ) = 1 Γ ( t ) 1 t ( ln t s ) t 1 ln s s d s ,
and
H I 1 + u ( t ) + v ( t ) x ( t ) = 1 Γ ( t + 1 ) 1 t ( ln t s ) t 1 s d s .
For t = 1.5 we obtain
H I 1 + u ( t ) H I 1 + v ( t ) x ( t ) | t = 1.5 0.027916
and
H I 1 + u ( t ) + v ( t ) x ( t ) | t = 1.5 0.0418739
Therefore, the semigroup property is not satisfied for the general case of HFIs of variable orders.
Recall the following pivotal observation.
Lemma 1.
([2]) Let α 1 , α 2 > 0 , a 1 > 1 , h L 1 ( a 1 , a 2 ) , and H D a 1 + α 1 h L 1 ( a 1 , a 2 ) . Then, the differential equation
H D a 1 + α 1 h = 0
has a solution
h ( t ) = ω 1 ( ln t a 1 ) α 1 1 + ω 2 ( ln t a 1 ) α 1 2 + + ω n ( ln t a 1 ) α 1 n ,
and
H I a 1 + α 1 ( H D a 1 + α 1 ) h ( t ) = h ( t ) + ω 1 ( ln t a 1 ) α 1 1 + ω 2 ( ln t a 1 ) α 1 2 + + ω n ( ln t a 1 ) α 1 n ,
with n 1 < α 1 n , ω R , = 1 , 2 , , n .
Furthermore,
H D a 1 + α 1 ( H I a 1 + α 1 ) h ( t ) = h ( t ) ,
and
H I a 1 + α 1 ( H I a 1 + α 2 ) h ( t ) = H I a 1 + α 2 ( H I a 1 + α 1 ) h ( t ) = H I a 1 + α 1 + α 2 h ( t ) .
We will give some results about HFI defined by (2) which will be used later.
Lemma 2.
If u C ( J , ( 1 , 2 ] ) and there exists a number δ [ 0 , 1 ] such that h C δ ( J , R ) , then the variable-order fractional integral H I 1 + u ( t ) h ( t ) exists for t J .
Proof. 
The function Γ ( u ( t ) ) is continuous non-zero function on J. Denote M u = max t J | 1 Γ ( u ( t ) ) | > 0 and u * = max t J | ( u ( t ) ) | . Thus, for 1 s t T , we have
( ln t s ) u ( t ) 1 1 , i f 1 t s e ,
( ln t s ) u ( t ) 1 ( ln t s ) u * 1 , i f t s > e .
Then, for 1 t s < + , we obtain
( ln t s ) u ( t ) 1 max { 1 , ( ln t s ) u * 1 } = M * .
For t J by (2), applying that the function ( ln ( . ) ) δ is an increasing function on J for any δ [ 0 , 1 ] we obtain
| H I 1 + u ( t ) h ( t ) | = 1 Γ ( u ( t ) ) 1 t ( ln t s ) u ( t ) 1 | h ( s ) | s d s M u 1 t ( ln t s ) u ( t ) 1 ( ln s ) δ ( ln s ) δ | h ( s ) | s d s M u M * ( ln T ) δ h 1 t 1 s ( ln s ) δ d s M u M * h ( ln T ) 1 δ < .
where h = max t J | h ( t ) | . It yields that the variable-order fractional integral H I 1 + u ( t ) h ( t ) exists for any t J . □
Lemma 3.
Let u C ( J , ( 1 , 2 ] ) , then H I 1 + u ( t ) h C ( J , R ) for h C ( J , R ) .
Proof. 
For t , t 0 ( 1 , T ] , t 0 t and h C ( J , R ) , using (2) and the substitutions s = r ( t 1 ) + 1 , d s = ( t 1 ) d r and s = r ( t 0 1 ) + 1 , d s = ( t 0 1 ) d r , respectively, we obtain
| H I 1 + u ( t ) h ( t ) H I 1 + u ( t 0 ) h ( t 0 ) | = | 0 1 1 Γ ( u ( t ) ) ( t 1 ) r ( t 1 ) + 1 ( ln t r ( t 1 ) + 1 ) u ( t ) 1 h ( r ( t 1 ) + 1 ) d r 0 1 1 Γ ( u ( t 0 ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 h ( r ( t 0 1 ) + 1 ) d r | .
Denote h = max t J | h ( t ) | and obtain the following upper bounds of the absolute values
| 0 1 [ 1 Γ ( u ( t ) ) ( t 1 ) r ( t 1 ) + 1 ( ln t r ( t 1 ) + 1 ) u ( t ) 1 h ( r ( t 1 ) + 1 ) 1 Γ ( u ( t ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t r ( t 1 ) + 1 ) u ( t ) 1 h ( r ( t 1 ) + 1 ) ] d r | h 0 1 1 Γ ( u ( t ) ) ( ln t r ( t 1 ) + 1 ) u ( t ) 1 | ( t 1 ) r ( t 1 ) + 1 ( t 0 1 ) r ( t 0 1 ) + 1 | d r
| 0 1 [ 1 Γ ( u ( t ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 h ( r ( t 1 ) + 1 ) 1 Γ ( u ( t 0 ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 h ( r ( t 1 ) + 1 ) ] d r | h 0 1 ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 | 1 Γ ( u ( t ) ) 1 Γ ( u ( t 0 ) ) | d r ,
| 0 1 [ 1 Γ ( u ( t 0 ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 h ( r ( t 1 ) + 1 ) 1 Γ ( u ( t 0 ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 h ( r ( t 0 1 ) + 1 ) ] d r | 0 1 1 Γ ( u ( t 0 ) ) ( t 0 1 ) r ( t 0 1 ) + 1 ( ln t 0 r ( t 0 1 ) + 1 ) u ( t 0 ) 1 | h ( r ( t 1 ) + 1 ) h ( r ( t 0 1 ) + 1 ) | d r ,
From (4)–(7) taking into account the continuity of functions ( t 1 ) r ( t 1 ) + 1 , ( ln t r ( t 1 ) + 1 ) u ( t ) 1 , 1 Γ ( u ( t ) ) , h ( t ) on J, we obtain that the integral H I 1 + u ( t ) h ( t ) is continuous at point t 0 , and therefore, H I 1 + u ( t ) h C ( J , R ) for h C ( J , R ) . □
Definition 1.
([5,10,11]) Let the set I R .
-
The set I is called a generalized interval if it is either an interval or a point or the empty set.
-
The finite set P of generalized intervals is called a partition of I if each x in I lies in exactly one of the generalized intervals E in P .
-
The function g : I R is called a piecewise constant with respect to partition P of I if for any E P , g is constant on E.
In addition, we will provide some necessary background information about the Kuratowski measure of noncompactness, defined in [32]). These results will be used later in the proof of the main result.
Proposition 1.
([32,33]). Let X be a Banach space and D , D 1 , D 2 are bounded subsets of X, then
  • ζ ( D ) = 0 D is relatively compact.
  • ζ ( ϕ ) = 0 .
  • ζ ( D ) = ζ ( D ¯ ) = ζ ( c o n v D ) .
  • D 1 D 2 ζ ( D 1 ) ζ ( D 2 ) .
  • ζ ( D 1 + D 2 ) ζ ( D 1 ) + ζ ( D 2 ) .
  • ζ ( λ D ) = | λ | ζ ( D ) , λ R .
  • ζ ( D 1 D 2 ) = M a x { ζ ( D 1 ) , ζ ( D 2 ) } .
  • ζ ( D 1 D 2 ) = M i n { ζ ( D 1 ) , ζ ( D 2 ) } .
  • ζ ( D + x 0 ) = ζ ( D ) for any x 0 X .
Lemma 4.
([34]) If U C ( J , R ) is an equicontinuous and bounded set, then
(i) 
the function ζ ( U ( t ) ) is continuous for t J , and ζ ^ ( U ) = sup t J ζ ( U ( t ) ) .
(ii) 
ζ 0 T x ( θ ) d θ : x U 0 T ζ ( U ( θ ) ) d θ ,
where for any s J it is defined U ( s ) = { x ( s ) : x U } .
Theorem 1.
([32] ) Let Λ be a nonempty, closed, bounded and convex subset of a Banach space X and ϝ : Λ Λ is a continuous operator satisfying
ζ ( ϝ ( S ) ) k ζ ( S ) , f o r a n y ( S ) Λ , k [ 0 , 1 ) ,
i.e., ϝ d e n o t e s k setcontractions .
Then, ϝ has at least one fixed point in Λ.
Following the ideas in [27] we will give a definition for Hyers–Ulam–Rassias stability of the BVP for VOHFDE (1).
Definition 2.
Let the function ϑ C ( J , R + ) . The BVP for VOHFDE (1) is Hyers–Ulam–Rassias stable with respect to ϑ (UHR) if there exists a constant c f > 0 , such that for any ϵ > 0 and for every z C ( J , R ) such that
| H D 1 + u ( t ) z ( t ) f ( t , z ( t ) , I 1 + u ( t ) z ( t ) ) | ϵ ϑ ( t ) , t J ,
there exists a solution x C ( J , R ) of the BVP for VOHFDE (1) with
| z ( t ) x ( t ) | c f ϵ ϑ ( t ) , t J .

3. Existence of Solutions of BVP for VOHFDE

Let us introduce the following assumption:
Hypothesis 1 (H1).
Let n N be an integer and the finite sequence of points { T k } k = 0 n be given such that 1 = T 0 < T k < T n = T , k = 1 , , n 1 . Denote J k : = ( T k 1 , T k ] , k = 1 , 2 , , n . Then P = k = 1 n J k is a partition of the interval J.
For each l = 1 , 2 , , n , the symbol E l = C ( J l , R ) , indicates the Banach space of continuous functions x : J l R equipped with the norm x E l = sup t J l | x ( t ) | .
Let u ( t ) : J ( 1 , 2 ] be a piecewise constant function with respect to P , i.e., u ( t ) = l = 1 n u l I l ( t ) , where 1 < u l 2 are constants and I l is the indicator of the interval J l , l = 1 , 2 , , n :
I l ( t ) = 1 , f o r t J l , 0 , e l s e w h e r e .
Then, for any t J l , l = 1 , 2 , , n , the left Hadamard fractional derivative of variable-order u ( t ) for function x ( t ) C ( J , R ) , defined by (3), could be presented as a sum of left Hadamard fractional derivatives of constant orders u k , k = 1 , 2 , , l
H D 1 + u ( t ) x ( t ) = t 2 Γ ( 2 u ( t ) ) d 2 d t 2 1 t ( ln t s ) 1 u ( t ) x ( s ) s d s = t 2 Γ ( 2 u ( t ) ) ( k = 1 l 1 d 2 d t 2 T k 1 T k ( ln t s ) 1 u k x ( s ) s d s + d 2 d t 2 T l 1 t ( ln t s ) 1 u x ( s ) s d s ) .
Thus, the BVP for VOHFDE (1) can be written for any t J l , l = 1 , 2 , , n in the form
t 2 Γ ( 2 u ( t ) ) k = 1 l 1 d 2 d t 2 T k 1 T k ( ln t s ) 1 u k x ( s ) s d s + d 2 d t 2 T l 1 t ( ln t s ) 1 u x ( s ) s d s = f ( t , x ( t ) , H I 1 + u ( t ) x ( t ) ) .
Let the function x ˜ C ( J , R ) be such that x ˜ ( t ) 0 on t [ 1 , T 1 ] and it solves integral Equation (10). Then (10) is reduced to
H D T 1 + u x ˜ ( t ) = f ( t , x ˜ ( t ) , H I T 1 + u x ˜ ( t ) ) , t J .
Taking into account the above for any = 1 , 2 , , n , we consider the following auxiliary BVP for Hadamard fractional differential equations of constant order
H D T 1 + u y ( t ) = f ( t , y ( t ) , H I T 1 + u y ( t ) ) , t J y ( T 1 ) = 0 , y ( T ) = 0 .
Lemma 5.
Let { 1 , 2 , n } be a natural number, f C ( J × R × R , R ) and there exists a number δ ( 0 , 1 ) such that ( ln t ) δ f C ( J × R × R , R ) .
Then the function y E is a solution of (11) if and only if y solves the integral equation
y ( t ) = ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 H I T 1 + u f ( t , y ( t ) , H I T 1 + u y ( t ) ) | t = T + H I T 1 + u f ( t , y ( t ) , H I T 1 + u y ( t ) ) .
Proof. 
Let y E be a solution of BVP (11). Employing the operator H I T 1 + u to both sides of (11), we find (see Lemma 1)
y ( t ) = ω 1 ( ln t T 1 ) u 1 + ω 2 ( ln t T 1 ) u 2 + H I T 1 + u f ( t , y ( t ) , H I T 1 + u y ( t ) ) , t J ,
where ω 1 , ω 2 are constants.
From the conditions of function f and the boundary condition y ( T 1 ) = 0 , we conclude that ω 2 = 0 .
From the boundary condition y ( T ) = 0 we obtain
ω 1 = ( ln T T 1 ) 1 u H I T 1 + u f ( T , x ( T ) , H I T 1 + u x ( T ) )
Then, we find y solves integral Equation (12).
Conversely, let y E be a solution of integral Equation (12). Regarding the continuity of function ( ln t ) δ f , we deduce that y is the solution of BVP (11). □
We will prove the existence result for the BVP for Hadamard fractional differential equations of constant order (11). This result is based on Theorem 1.
Theorem 2.
Let the conditions of Lemma 5 be satisfied and there exist constants K , L > 0 such that
( ln t ) δ | f ( t , x 1 , y 1 ) f ( t , x 2 , y 2 ) | K | x 1 x 2 | + L | y 1 y 2 | , x i , y i , R , i = 1 , 2 , t J
and the inequality
μ < 1 2 ,
holds, where
μ = ( ln T ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u .
Then, the BVP for Hadamard fractional differential equations of constant order (11) possesses at least one solution in E .
Remark 4.
According to the remark on page 20 of [35], we can easily show that the inequality (13) and the following inequality
ζ ( ( ln t ) δ f ( t , B 1 , B 2 ) ) K ζ ( B 1 ) + L ζ ( B 2 )
are equivalent for any bounded sets B 1 , B 2 R and for each t J .
Proof. 
For any function x E we define the operator
W x ( t ) = 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 f ( s , x ( s ) , H I T 1 + u x ( s ) ) s d s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 f ( s , x ( s ) , H I T 1 + u x ( s ) ) s d s .
It follows from the properties of fractional integrals and from the continuity of function ( ln t ) δ f that the operator W : E E defined by (14) is well defined.
Let R = 2 f Γ ( u + 1 ) ( ln T T 1 ) u 1 2 μ with f = sup t J | f ( t , 0 , 0 ) | .
Consider the set
B R = { x E , x E R } .
For all { 1 , 2 , , n } the ball B R is a nonempty, bounded, closed convex subset of E .
Now, we check the assumption of the Theorem 1 for the operator W. We shall prove it in four steps.
Step 1: Claim: W ( B R ) B R .
For x B R and by the conditions of function f, we obtain
| W x ( t ) | 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s 2 Γ ( u ) T 1 T ( ln T s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s 2 Γ ( u ) T 1 T ( ln T s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) f ( s , 0 , 0 ) | d s s + 2 Γ ( u ) T 1 T ( ln T s ) u 1 | f ( s , 0 , 0 ) | d s s 2 Γ ( u ) T 1 T ( ln T s ) u 1 ( ln s ) δ K | x ( s ) | + L | H I T 1 + u x ( s ) | d s s + 2 f Γ ( u ) T 1 T ( ln T s ) u 1 d s s 2 Γ ( u ) ( ln T T 1 ) u 1 K x E + L H I T 1 + u x E T 1 T ( ln s ) δ d s s + 2 f Γ ( u + 1 ) ( ln T T 1 ) u 2 Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u x E T 1 T ( ln s ) δ d s s + 2 f Γ ( u + 1 ) ( ln T T 1 ) u 2 μ R + 2 f Γ ( u + 1 ) ( ln T T 1 ) u = R ,
which means that W ( B R ) B R .
Step 2:Claim: W is continuous.
Let x k E , k = 1 , 2 , . Assume the sequence { x k } k = 1 is convergent to x E . Then for any k = 1 , 2 , we have
| W x k ( t ) W x ( t ) | 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 | f ( s , x k ( s ) , H I T 1 + u x k ( s ) ) f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 | f ( s , x k ( s ) , H I T 1 + u x k ( s ) ) f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s 2 Γ ( u ) ( ln t T 1 ) u 1 T 1 T | f ( s , x k ( s ) , H I T 1 + u x k ( s ) ) f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s 2 Γ ( u ) ( ln t T 1 ) u 1 T 1 T ( ln s ) δ K | x k ( s ) x ( s ) | + L H I T 1 + u | x k ( s ) x ( s ) | d s s 2 Γ ( u ) K x k x E + L H I T 1 + u ( x k x ) E ( ln t T 1 ) u 1 T 1 T ( ln s ) δ d s s 2 [ ( ln T ) 1 δ ( ln T 1 ) 1 δ ] ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u x k x E ,
i.e., we obtain
W x k W x E 0 a s k .
Ergo, the operator W is continuous on E .
Step 3: Claim: W is bounded and equicontinuous.
By Step 1 for any x B R we have W ( x ) E R , i.e., W ( B R ) is bounded. We will prove W ( B R ) is equicontinuous.
Let t 1 , t 2 J , t 1 < t 2 and x B R . Then we have
| W x ( t 2 ) W x ( t 1 ) | 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 × T 1 T ( ln T s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s + 1 Γ ( u ) T 1 t 1 ( ln t 2 s ) u 1 ( ln t 1 s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s + 1 Γ ( u ) t 1 t 2 ( ln t 2 s ) u 1 | f ( s , x ( s ) , H I T 1 + u x ( s ) ) | d s s
or
| W x ( t 2 ) W x ( t 1 ) | 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 × T 1 T ( ln T s ) u 1 ( ln s ) δ ( K | x ( s ) | + L H I T 1 + u | x ( s ) | ) d s s + f Γ ( u ) ( ln T T 1 ) 1 u ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 T 1 T ( ln T s ) u 1 d s s + 1 Γ ( u ) T 1 t 1 ( ln t 2 t 1 ) u 1 ( ln s ) δ ( K | x ( s ) | + L H I T 1 + u | x ( s ) | ) d s s + f Γ ( u ) T 1 t 1 ( ln t 2 t 1 ) u 1 d s s + f Γ ( u ) t 1 t 2 ( ln t 2 s ) u 1 d s s + 1 Γ ( u ) t 1 t 2 ( ln t 2 s ) u 1 ( ln s ) δ ( K | x ( s ) | + L H I T 1 + u | x ( s ) | ) d s s 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln T T 1 ) u 1 ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 × K x E + L H I T 1 + u x E × T 1 T ( ln s ) δ d s s + f Γ ( u + 1 ) ( ln T T 1 ) 1 u ( ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 ) ( ln T T 1 ) u + ( ln t 2 ) 1 δ ( ln t 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln t 2 t 1 ) u 1 K x E + L H I T 1 + u x E + f Γ ( u ) ( ln t 2 t 1 ) u 1 ( ln t 1 ln T 1 ) + ( ln t 2 ) 1 δ ( ln t 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln t 2 t 1 ) u 1 K x E + L H I T 1 + u x E + f Γ ( u + 1 ) ( ln t 2 t 1 ) u .
Therefore, we obtain
| W x ( t 2 ) W x ( t 1 ) | ( ln T ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) K + L Γ ( u + 1 ) ( ln T T 1 ) u x E + f Γ ( u + 1 ) ( ln T T 1 ) × ( ln t 2 T 1 ) u 1 ( ln t 1 T 1 ) u 1 + [ 2 ( ( ln t 2 ) 1 δ ( ln t 1 ) 1 δ ) ( 1 δ ) Γ ( u ) K + L Γ ( u + 1 ) ( ln T T 1 ) u x E + f Γ ( u ) ( ln t 1 ln T 1 ) ] ( ln t 2 t 1 ) u 1 + f Γ ( u + 1 ) ( ln t 2 t 1 ) u .
Hence | W x ( t 2 ) W x ( t 1 ) | 0 as | t 2 t 1 | 0 . This implies that W ( B R ) is equicontinuous.
Step 4:Claim: W is k-set contraction.
For U B R , t J , we obtain,
ζ ( W ( U ) ( t ) ) = ζ ( ( W y ) ( t ) , y U ) { 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 ζ f ( s , x ( s ) , H I T 1 + u x ( s ) ) d s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 ζ f ( s , x ( s ) , H I T 1 + u x ( s ) ) d s , y U } .
Then Remark 4 implies that, for each s J ,
ζ ( W ( U ) ( t ) ) { 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 × K ζ ^ ( U ) s δ + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) s δ d s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 K ζ ^ ( U ) s δ + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) s δ d s , y U } . { 1 Γ ( u ) ( ln T T 1 ) u 1 T 1 T K ζ ^ ( U ) s δ + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) s δ d s + 1 Γ ( u ) ( ln T T 1 ) u 1 T 1 t K ζ ^ ( U ) s δ + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) s δ d s , y U } . ( ln T ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) + ( ln t ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) 2 [ ( ln T ) 1 δ ( ln T 1 ) 1 δ ] ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U )
Thus
ζ ^ ( W U ) 2 [ ( ln T ) 1 δ ( ln T 1 ) 1 δ ] ( 1 δ ) Γ ( u ) ( ln T T 1 ) u 1 K + L Γ ( u + 1 ) ( ln T T 1 ) u ζ ^ ( U ) .
From inequality (13) it follows that W is a k-set contraction.
Therefore, all conditions of Theorem 1 are fulfilled and thus there exists x ˜ B R , such that W ( x ˜ ) = x ˜ , which is a solution of the BVP for Hadamard fractional differential equations of constant order (11). Since B R E the claim of Theorem 2 is proved. □
Now, we will prove the existence result for the BVP for VOHFDE (1).
Introduce the following assumption:
Hypothesis 2 (H2).
Let f C ( J × R × R , R ) and there exists a number δ ( 0 , 1 ) such that ( ln t ) δ f C ( J × R × R , R ) and there exist constants K , L > 0 such that
( ln t ) δ | f ( t , x 1 , y 1 ) f ( t , x 2 , y 2 ) | K | x 1 x 2 | + L | y 1 y 2 |
for any x 1 , x 2 , y 1 , y 2 , R and t J .
Theorem 3.
Let the conditions (H1), (H2) and inequality (13) be satisfied for all { 1 , 2 , , n } .
Then, the BVP for VOHFDE (1) possesses at least one solution in C ( J , R ) .
Proof. 
For any { 1 , 2 , , n } , according to Theorem 2 the BVP for Hadamard fractional differential equations of constant order (11) possesses at least one solution x ˜ E .
For any { 1 , 2 , , n } we define the function
x = 0 , t [ 1 , T 1 ] , x ˜ , t J .
Thus, the function x C ( [ 1 , T ] , R ) solves the integral Equation (10) for t J , which means that x ( 1 ) = 0 , x ( T ) = x ˜ ( T ) = 0 and solves (10) for t J , { 1 , 2 , , n } .
Then the function
x ( t ) = x 1 ( t ) , t J 1 , x 2 ( t ) , t J 2 , x n ( t ) , t J n = [ 1 , T ]
is a solution of the BVP for VOHFDE (1) in C ( J , R ) . □

4. Ulam–Hyers–Rassias Stability of VOHFDE

We introduce the following assumption:
Hypothesis 3 (H3).
The function ϑ C ( J , R + ) is increasing and there exists λ ϑ > 0 such that
H I T 1 + u ϑ ( t ) λ ϑ ϑ ( t ) , f o r t J , = 1 , 2 , , n .
Theorem 4.
Let the conditions (H1), (H2), (H3) and inequality (13) be satisfied.
Then, the BVP for VOHFDE (1) is UHR stable with respect to ϑ.
Proof. 
Let ϵ > 0 be an arbitrary number and the function z ( t ) from C ( J , R ) satisfy inequality (8).
For any { 1 , 2 , , n } we define the functions z 1 ( t ) z ( t ) , t [ 1 , T 1 ] and for = 2 , 3 , , n :
z ( t ) = 0 , t [ 1 , T 1 ] , z ( t ) , t J .
For any { 1 , 2 , , n } according to equality (3) for t J we obtain
H D 1 + u ( t ) z ( t ) = t 2 Γ ( 2 u ( t ) ) d 2 d t 2 T 1 t ( ln t s ) 1 u z ( s ) s d s .
We take the HFI H I T 1 + u of both sides of the inequality (8), apply (H3) and obtain
| z ( t ) + 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 T 1 T ( ln T s ) u 1 f ( s , z ( s ) , H I T 1 + u z ( s ) ) s d s 1 Γ ( u ) T 1 t ( ln t s ) u 1 f ( s , z ( s ) , H I T 1 + u z ( s ) ) s d s | ϵ H I T 1 + u ϑ ( t ) ϵ λ ϑ ϑ ( t ) .
According to Theorem 3, the BVP for VOHFDE (1) has a solution x C ( J , R ) defined by x ( t ) = x ( t ) for t J , = 1 , 2 , , n , where
x = 0 , t [ 1 , T 1 ] , x ˜ , t J .
and x ˜ E is a solution of (11). According to Lemma 5 the integral equation
x ˜ ( t ) = ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 H I T 1 + u f ( t , x ˜ ( t ) , H I T 1 + u x ˜ ( t ) ) | t = T + H I T 1 + u f ( t , x ˜ ( t ) , H I T 1 + u x ˜ ( t ) ) , t J .
holds.
Let t J where { 1 , 2 , , n } . Then by Equations (13) and (14) we obtain
| z ( t ) x ( t ) | = | z ( t ) x ( t ) | = | z ( t ) x ˜ ( t ) | | z ( t ) + 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 × T 1 T ( ln T s ) u 1 f ( s , z ( s ) , H I T 1 + u z ( s ) ) s d s 1 Γ ( u ) T 1 t ( ln t s ) u 1 f ( s , z ( s ) , H I T 1 + u z ( s ) ) s d s | + 1 Γ ( u ) ( ln T T 1 ) 1 u ( ln t T 1 ) u 1 × T 1 T 1 s ( ln T s ) u 1 | f ( s , z ( s ) , H I T 1 + u z ( s ) ) f ( s , x ˜ ( s ) , H I T 1 + u x ˜ ( s ) ) | d s + 1 Γ ( u ) T 1 t 1 s ( ln t s ) u 1 | f ( s , z ( s ) , H I T 1 + u z ( s ) ) f ( s , x ˜ ( s ) , H I T 1 + u x ˜ ( s ) ) | d s × T 1 T ( ln T s ) u 1 1 s ( ln s ) δ ( K | z ( s ) x ˜ ( s ) | + L H I T 1 + u | z ( s ) x ˜ ( s ) | ) d s + 1 Γ ( u ) T 1 t ( ln t s ) u 1 1 s ( ln s ) δ ( K | z ( s ) x ˜ ( s ) | + L H I T 1 + u | z ( s ) x ˜ ( s ) | ) d s λ ϑ ϵ ϑ ( t ) + 1 Γ ( u ) ( ln t T 1 ) u 1 T 1 T 1 s ( ln s ) δ ( K | z ( s ) x ˜ ( s ) | + L H I T 1 + u | z ( s ) x ˜ ( s ) | ) d s + 1 Γ ( u ) ( ln t T 1 ) u 1 T 1 t 1 s ( ln s ) δ ( K | z ( s ) x ˜ ( s ) | + L H I T 1 + u | z ( s ) x ˜ ( s ) | ) d s λ ϑ ϵ ϑ ( t ) + ( ln T ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln t T 1 ) u 1 ( K z x ˜ E + L H I T 1 + u z x ˜ E ) + ( ln t ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u ) ( ln t T 1 ) u 1 ( K z x ˜ E + L H I T 1 + u z x ˜ E ) λ ϑ ϵ ϑ ( t ) + μ z x J .
Then,
z x J 1 μ λ ϑ ϵ ϑ ( t )
or for any t J
| z ( t ) x ( t ) z x J λ ϑ 1 μ ϵ ϑ ( t ) .
Therefore, the BVP for VOHFDE (1) is UHR stable with respect to ϑ . □

5. Example

Let J : = [ 1 , e ] , T 0 = 1 , T 1 = 2 , T 2 = e . Consider the scalar BVP for VOHFDE
D 1 + u ( t ) x ( t ) = 7 5 π ( ln t ) u ( t ) + ( ln t ) 1 3 t + 3 x ( t ) + π ln t t 2 + 1 H I 1 + u ( t ) x ( t ) , t J , x ( 1 ) = 0 , x ( e ) = 0 .
where
u ( t ) = 1.3 , t J 1 : = [ 1 , 2 ] , 1.7 , t J 2 : = ] 2 , e ] .
Denote
f ( t , x , y ) = 7 5 π ( ln t ) u ( t ) + ( ln t ) 1 3 t + 3 x + π ln t t 2 + 1 y , ( t , x , y ) [ 1 , e ] × R × R .
Then, we have
( ln t ) 1 3 | f ( t , x 1 , y 1 ) f ( t , x 2 , y 2 ) | = | 1 t + 3 x 1 + π ( ln t ) 4 3 t 2 + 1 y 1 1 t + 3 x 2 π ( ln t ) 4 3 t 2 + 1 y 2 | 1 t + 3 | x 1 x 2 | + π ( ln t ) 4 3 t 2 + 1 | y 1 y 2 | 1 4 | x 1 x 2 | + π 2 | y 1 y 2 | .
Thus, assumption (H2) holds with δ = 1 3 , K = 1 4 , L = π 2 .
By (20), according to (11) we consider two auxiliary BVPs for Hadamard fractional differential equations of constant order
D 1 + 1.3 x ( t ) = 7 5 π ( ln t ) 1.3 + ( ln t ) 1 3 t + 3 x ( t ) + π ln t t 2 + 1 H I 1 + 1.3 x ( t ) , t J 1 , x ( 1 ) = 0 , x ( 2 ) = 0 .
and
D 2 + 1.7 x ( t ) = 7 5 π ( ln t ) 1.7 + ( ln t ) 1 3 t + 3 x ( t ) + π ln t t 2 + 1 H I 1 + 1.7 x ( t ) , t J 2 , x ( 2 ) = 0 , x ( e ) = 0 .
Next, we prove that the condition (13) is fulfilled for = 1 . Indeed,
( ln T 1 ) 1 δ ( ln T 0 ) 1 δ ( 1 δ ) Γ ( u 1 ) ( ln T 1 T 0 ) u 1 1 K + L Γ ( u 1 + 1 ) ( ln T 1 T 0 ) u 1 = ( ln 2 ) 2 3 ( 2 3 ) Γ ( 1.3 ) ( ln 2 ) 0.3 1 4 + π 2 Γ ( 2.3 ) ( ln 2 ) 1.3 0.3809 < 1 2 .
Let ϑ ( t ) = ( ln t ) 1 2 . Then, we obtain
H I 1 + u 1 ϑ ( t ) = 1 Γ ( 1.3 ) 1 t ( ln t s ) 1.3 1 ( ln s ) 1 2 s d s 1 Γ ( 1.3 ) 1 t ( ln t s ) 0.3 1 s d s 0.75 Γ ( 2.3 ) ( ln t ) 1 2 = λ ϑ ϑ ( t ) ,
where λ ϑ = 0.75 Γ ( 2.3 ) . Thus, condition (H3) is satisfied.
By Theorem 2, the BVP (21) has a solution x ˜ 1 E 1 .
We prove that the condition (13) is fulfilled for = 2 . Indeed,
( ln T 2 ) 1 δ ( ln T 1 ) 1 δ ( 1 δ ) Γ ( u 2 ) ( ln T 2 T 1 ) u 2 1 K + L Γ ( u 2 + 1 ) ( ln T 2 T 1 ) u 2 = 1 ( ln 2 ) 2 3 ( 2 3 ) Γ ( 1.7 ) ( ln e 2 ) 0.7 1 4 + π 2 Γ ( 2.7 ) ( ln e 2 ) 1.7 0.0242 < 1 2 .
Accordingly the condition (13) is achieved. We obtain
H I 2 + u 2 ϑ ( t ) = 1 Γ ( 1.7 ) 2 t ( ln t s ) 1.7 1 ( ln s ) 1 2 s d s 1 Γ ( 1.7 ) 2 t ( ln t s ) 0.7 1 s d s 1 Γ ( 2.7 ) ( ln t ) 1 2 = λ ϑ ϑ ( t ) .
where λ ϑ = 1 Γ ( 2.7 ) . Thus condition (H3) is fulfilled.
According to Theorem 2, the BVP (22) possesses a solution x ˜ 2 E 2 .
Thus, according to Theorem 3 the BVP (19) has a solution
x ( t ) = x ˜ 1 ( t ) , t J 1 , x 2 ( t ) , t J 2 ,
where
x 2 ( t ) = 0 , t J 1 x ˜ 2 ( t ) , t J 2 .
According to Theorem (4), the BVP for VOHFDE (19) is UHR stable with respect to ϑ .

6. Conclusions

Our proposed multiterm BVP has been successfully investigated in this work via the Kuratowski measure of noncompactness technique to prove the existence of solutions to VOHFDEs. A numerical example was given at the end to support and validate the potentiality of all our obtained results. Therefore, all results in this work show a great potential to be applied in various applications of multidisciplinary sciences. Further investigations on this open research problem could be also possible with the help of our original results in this research paper. In other words, in the future one could extend the proposed BVP to other complicated real mathematical fractional models by terms of newly-introduced operators with non-singular kernels.

Author Contributions

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

Funding

S.H. was funded by the Science Fund of Plovdiv University under Project FP21FMI002.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdeljawad, A.; Agarwal, R.P.; Karapinar, E.; Kumari, P.S. Solutions of the nonlinear integral equation and fractional differential equation using the technique of a fixed point with a numerical experiment in extended b-Metric Space. Symmetry 2019, 11, 686. [Google Scholar] [CrossRef] [Green Version]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differenatial Equations; North-Holland Mathematics Studies, 204; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. de Oliveira, E.C.; Sousa, J.; Sousa, J.V.D. Ulam-Hyers-Rassias stability for a class of fractional integro-differential equations. Results Math. 2018, 73, 111. [Google Scholar] [CrossRef] [Green Version]
  4. Chikh, S.B.; Amara, A.; Etemad, S.; Rezapour, S. On Ulam-Hyers-Rassias stability of a generalized Caputo type multi-order boundary value problem with four-point mixed integro-derivative conditions. Adv. Differ. Eq. 2020, 2020, 680. [Google Scholar] [CrossRef]
  5. Jiahui, A.; Pengyu, C. Uniqueness of solutions to initial value problem of fractional differential equations of variable-order. Dyn. Syst. Appl. 2019, 28, 607–623. [Google Scholar]
  6. Ragusa, M.A. Embeddings for Morrey-Lorentz spaces. J. Optim. Theory Appl. 2012, 154, 491–499. [Google Scholar] [CrossRef]
  7. Samko, S.G. Fractional integration and differentiation of variable order. Anal. Math. 1995, 21, 213–236. [Google Scholar] [CrossRef]
  8. Samko, S.G.; Boss, B. Integration and differentiation to a variable fractional order. Integr. Transform. Spec. Funct. 1993, 1, 277–300. [Google Scholar] [CrossRef]
  9. Zhang, S. Existence of solutions for two point boundary value problems with singular differential equations of variable order. Elect. J. Differ. Eq. 2013, 245, 1–16. [Google Scholar]
  10. Zhang, S. The uniqueness result of solutions to initial value problems of differential equations of variable-order. Rev. Real Acad. Cienc. Exactas Físicas Nat. Ser. A Mat. 2018, 112, 407–423. [Google Scholar] [CrossRef]
  11. Zhang, S.; Hu, L. The existence of solutions and generalized Lyapunov-type inequalities to boundary value problems of differential equations of variable order. AIMS Math. 2020, 5, 2923–2943. [Google Scholar] [CrossRef]
  12. Patnaik, S.; Hollkamp, J.P.; Semperlotti, F. Applications of variable-order fractional operators: A review. Proc. R. Soc. 2020, A476, 20190498. [Google Scholar] [CrossRef] [Green Version]
  13. Sun, H.G.; Chang, A.; Zhang, Y.; Chen, W. A review on variable order fractional differential equations: Mathematical foundations, Physical models, numerical methods and applications. Frac. Calc. Appl. Anal. 2019, 22, 27–59. [Google Scholar] [CrossRef] [Green Version]
  14. Starrett, J. Solving differential equations by symmetry groups. Mar. Math. Mon. 2007, 114, 778–792. [Google Scholar] [CrossRef]
  15. Chatibi, Y.; Kinani, E.H.E.; Ouhadan, A. Lie symmetry analysis of conformable differential equations. AIMS Math. 2019, 4, 1133–1144. [Google Scholar] [CrossRef]
  16. Adiguzel, R.S.; Aksoy, U.; Karapınar, E.; Erhan, I.M. On the solution of a boundary value problem associated with a fractional differential equation. Math. Meth. Appl. Sci. 2020. [Google Scholar] [CrossRef]
  17. Afshari, H.; Karapınar, E. A discussion on the existence of positive solutions of the boundary value problems via ψ-Hilfer fractional derivative on b-metric spaces. Adv. Differ. Eq. 2020, 2020, 616. [Google Scholar] [CrossRef]
  18. Zada, A.; Alzabut, J.; Waheed, H.; Popa, I.-L. Ulam—Hyers stability of impulsive integro-differential equations with Riemann–Liouville boundary conditions. Adv. Differ. Eq. 2020, 2020, 64. [Google Scholar] [CrossRef] [Green Version]
  19. Aguilar, J.F.G. Analytical and numerical solutions of nonlinear alcoholism model via variable-order fractional differential equations. Phys. A 2018, 494, 52–57. [Google Scholar] [CrossRef]
  20. Sousa, J.V.D.C.; de Oliverira, E.C. Two new fractional derivatives of variable order with non-singular kernel and fractional differential equation. Comput. Appl. Math. 2018, 37, 5375–5394. [Google Scholar] [CrossRef]
  21. Tavares, D.; Almeida, R.; Torres, D.F. Caputo derivatives of fractional variable order Numerical approximations. Commun. Nonlinear Sci. Numer. Simul. 2016, 35, 6987. [Google Scholar] [CrossRef] [Green Version]
  22. 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] [Green Version]
  23. Zhang, S.; Sun, S.; Hu, L. The existeness and uniqueness result of solutions to initial value problems of nonlinear diffusion equations involving with the conformable variable derivative. Rev. Real Acad. Cienc. Exactas FíSicas Nat. Ser. Mat. 2019, 113, 1601–1623. [Google Scholar] [CrossRef]
  24. Hadamard, J. Essai sur l’étude des fonctions, données par leur développement de Taylor. J. Mat. Pure Appl. Ser. 1892, 4, 101–186. [Google Scholar]
  25. Almeida, R.; Torres, D.F.M. Computing Hadamard type operators of variable fractional order. Appl. Math. Comput. 2015, 257, 74–88. [Google Scholar] [CrossRef] [Green Version]
  26. Almeida, R.; Tavares, D.; Torres, D.F.M. The Variable-Order Fractional Calculus of Variations; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  27. Benchohra, M.; Lazreg, J.E. Existence and Ulam stability for nonlinear implicit fractional differential equations with Hadamard derivative. Stud. Univ. Babes-Bolyai Math. 2017, 62, 27–38. [Google Scholar] [CrossRef] [Green Version]
  28. Bai, Y.; Kong, H. Existence of solutions for nonlinear Caputo-Hadamard fractional differential equations via the method of upper and lower solutions. J. Nonlinear Sci. Appl. 2017, 10, 5744–5752. [Google Scholar] [CrossRef] [Green Version]
  29. Lavenda, B.H. Concepts of Stability and Symmetry in Irreversible Thermodynamics. I. Found. Phys. 1972, 2, 161–179. [Google Scholar] [CrossRef]
  30. Gallavotti, G. Breakdown and regeneration of time reversal symmetry in nonequilibrium statistical mechanics. Phys. D 1998, 112, 250–257. [Google Scholar] [CrossRef]
  31. Russo, G.; Slotine, J.-J.E. Symmetries, stability, and control in nonlinear systems and networks. Phys. Rev. E 2011, 84, 041929. [Google Scholar] [CrossRef] [Green Version]
  32. Banás, J.; Goebel, K. Measures of Noncompactness in Banach Spaces; Marcel Dekker: New York, NY, USA, 1980. [Google Scholar]
  33. Banaś, J.; Olszowy, L. Measures of noncompactness related to monotonicity. Comment. Math. (Prace Mat.) 2001, 41, 13–23. [Google Scholar]
  34. Guo, D.J.; Lakshmikantham, V.; Liu, X. Nonlinear Integral Equations in Abstract Spaces; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996. [Google Scholar]
  35. Benchohra, M.; Bouriah, S.; Lazreg, J.E.; Nieto, J.J. Nonlinear implicit Hadamard’s fractional differential equations with delay in Banach Space. Acta Uni. Palacki. Olomuc. Facult. Rerum Natur. Math. 2016, 55, 15–26. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hristova, S.; Benkerrouche, A.; Souid, M.S.; Hakem, A. Boundary Value Problems of Hadamard Fractional Differential Equations of Variable Order. Symmetry 2021, 13, 896. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13050896

AMA Style

Hristova S, Benkerrouche A, Souid MS, Hakem A. Boundary Value Problems of Hadamard Fractional Differential Equations of Variable Order. Symmetry. 2021; 13(5):896. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13050896

Chicago/Turabian Style

Hristova, Snezhana, Amar Benkerrouche, Mohammed Said Souid, and Ali Hakem. 2021. "Boundary Value Problems of Hadamard Fractional Differential Equations of Variable Order" Symmetry 13, no. 5: 896. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13050896

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