Next Article in Journal
A Quasilinearization Approach for Identification Control Vectors in Fractional-Order Nonlinear Systems
Next Article in Special Issue
Conformal and Non-Minimal Couplings in Fractional Cosmology
Previous Article in Journal
The Existence Results of Solutions to the Nonlinear Coupled System of Hilfer Fractional Differential Equations and Inclusions
Previous Article in Special Issue
Some Results on Fractional Boundary Value Problem for Caputo-Hadamard Fractional Impulsive Integro Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fundamental Matrix, Integral Representation and Stability Analysis of the Solutions of Neutral Fractional Systems with Derivatives in the Riemann—Liouville Sense

1
Faculty of Mathematics and Informatics, University of Plovdiv Paisii Hilendarski, 24 Tzar Asen, 4000 Plovdiv, Bulgaria
2
Faculty of Economics and Business Administration, Sofia University “St. Kliment Ohridski”, 15 Tsar Osvoboditel Blvd., 1504 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(4), 195; https://doi.org/10.3390/fractalfract8040195
Submission received: 31 December 2023 / Revised: 21 March 2024 / Accepted: 27 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Advances in Fractional Modeling and Computation)

Abstract

:
The paper studies a class of nonlinear disturbed neutral linear fractional systems with derivatives in the the Riemann–Liouville sense and distributed delays. First, it is proved that the initial problem for these systems with discontinuous initial functions under some natural assumptions possesses a unique solution. The assumptions used for the proof are similar to those used in the case of systems with first-order derivatives. Then, with the obtained result, we derive the existence and uniqueness of a fundamental matrix and a generalized fundamental matrix for the homogeneous system. In the linear case, via these fundamental matrices we obtain integral representations of the solutions of the homogeneous system and the corresponding inhomogeneous system. Furthermore, for the fractional systems with Riemann–Liouville derivatives we introduce a new concept for weighted stabilities in the Lyapunov, Ulam–Hyers, and Ulam–Hyers–Rassias senses, which coincides with the classical stability concepts for the cases of integer-order or Caputo-type derivatives. It is proved that the zero solution of the homogeneous system is weighted stable if and only if all its solutions are weighted bounded. In addition, for the homogeneous system it is established that the weighted stability in the Lyapunov and Ulam–Hyers senses are equivalent if and only if the inequality appearing in the Ulam–Hyers definition possess only bounded solutions. Finally, we derive natural sufficient conditions under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under nonlinear disturbances.

1. Introduction

Practically, it has been shown that many real-world phenomena in various fields of science can be represented more accurately through mathematical models, including fractional differential equations [1,2,3]. More detailed information on fractional calculus theory and fractional differential equations can be seen in the monographs of Kilbas et al. [4] and Podlubny [5]. Compared to fractional equations with Caputo-type derivatives, the fractional equations with Riemann–Liouville-type derivatives have been studied significantly less. The first obstacle was the lack of meaningful geometric and physical interpretations of the Riemann–Liouville type of integration and differentiation, which was overcome with the appearance of Ref. [6]. The other problem was to derive an appropriate formulation of the initial conditions for the initial problem. Some considerations about them are given in [7,8], but as was mentioned in [9] “The initial value problem is a subject that remains quite up-to-date”.
It is well known that the existence of explicit solutions or an integral representation (variation of constants formula) of the solutions of linear fractional differential equations and/or systems (ordinary or delayed) is a main tool in executing their qualitative analysis. That is why establishing either explicit solutions (see [10]) or integral representation, for which existence of a fundamental matrix is needed (see [11]), are important tasks for stability analysis, especially in the case of equations with Riemann–Liouville derivatives. But it is surprising that there are not many articles devoted to this problem. As far as we know, a survey concerning stability results for retarded and neutral fractional differential equations with Riemann–Liouville-type derivatives does not exist. In general, partially this gap can be fulfilled by reading the overview [12] and the references therein. From the recent works devoted to the discussed themes concerning fractional differential equations with Riemann–Liouville-type derivatives we refer to [13] for equations, to [14,15,16] where retarded fractional differential equations are considered, and to [17,18] where the neutral case is considered. We suppose that the mentioned works with their references give a good enough picture of the studies in this area.
In the present article, we study a class of nonlinear disturbed neutral linear fractional systems with derivatives in the the Riemann–Liouville sense and distributed delay. We first study the important problem of the existence and uniqueness of the solutions of an initial problem (IP) for these systems in the case of discontinuous initial functions. As far as we know, there are no results devoted to the initial problem with discontinuous initial functions for neutral differential equations with derivatives in the Riemann–Liouville sense. Since the classically stability concepts are not directly applicable to systems with derivatives with the property that the derivative of a constant is not equal to zero (like the Riemann–Liouville fractional derivative, for example), we introduce a new concept for weighted stabilities in the Lyapunov, Ulam–Hyers, and Ulam–Hyers–Rassias senses, which coincides with the classical stability concepts for the cases of integer-order or Caputo-type derivatives.
The following abbreviations will be used in this manuscript: BV—bounded variation; GAS—globally asymptotically stable; UH—Ulam–Hyers; LAS—locally asymptotically stable; LT—Laplace transform; UHR—Ulam–Hyers–Rassias; IP—initial problem; PC—piecewise continuous; RL—Riemann–Liouville; WML—weighted Mittag–Leffler.
The paper is organized as follows: Section 2 presents the necessary definitions and properties concerning the RL and Caputo fractional derivatives, the problem statement, and the needed auxiliary definitions and facts for our exposition. Section 3 is devoted to the existence and the uniqueness of the solutions of the initial problem (IP) for the studied class of nonlinear neutral systems in the case when the initial function is discontinuous. In Section 4, as a consequence of the derived result we prove the existence and uniqueness of a fundamental matrix and a generalized fundamental matrix for the linear homogeneous system, as well as establishing an integral representation of the solutions of the IP for the homogeneous system and the corresponding inhomogeneous system. Section 5 is devoted to a new concept for weighted stabilities in the Lyapunov, UH, and UHR senses, the definition of which coincides with the classical stability concepts for the cases of integer-order or Caputo-type derivatives. It is proved that the zero solution of the homogeneous system is weighted stable if and only if all its solutions are weighted bounded. In addition, it is established that the weighted stabilities in the Lyapunov and UH senses are equivalent if and only if the inequality appearing in the UH definition possesses only bounded solutions. In Section 6, natural sufficient conditions are obtained under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under nonlinear disturbances. Finally, in Section 7 an illustrative example is presented. In the last Section 8, we summarize some conclusions concerning the obtained results and propose some open problems.

2. Preliminaries and Problem Statement

As is standard to avoid possible misunderstandings, we recall some properties concerning the RL-type derivatives, as well as the needed definitions, conditions, and auxiliary results necessary for the exposition below. For more comprehensive information on fractional calculus we refer to Refs. [4,5].
For a function g : R R , we will say that some property holds locally if it is fulfilled on every compact subinterval [ b , c ] R . We will use the following notation for the real linear spaces: L 1 l o c ( R , R ) consists of all locally Lebesgue integrable functions, B L 1 l o c ( R , R ) L 1 l o c ( R , R ) and B V l o c ( R , R ) B L 1 l o c ( R , R ) are, respectively, the subspaces of all functions which are locally bounded and of all functions which have locally bounded variation.
We define the left-sided fractional RL integral operator of arbitrary order α ( 0 , 1 ) with lower limit a R for any g L 1 l o c ( R , R ) , via the relation
( I a + α g ) ( t ) = 1 Γ ( α ) a t ( t s ) α 1 g ( s ) d s .
The corresponding left-hand side RL and Caputo fractional derivatives of arbitrary order α ( 0 , 1 ) with lower limit a are defined for any t > a by
D a + α R L g ( t ) = d d t ( I a + 1 α g ( t ) ) ,
D a + α C g ( t ) = D a + α R L g ( t ) D a + α R L g ( a ) .
For shortness, we will use the notation D a + α = D a + α R L .
Consider the neutral fractional nonlinear system with RL-type derivatives of order α ( 0 , 1 ) with lower limit a R and distributed delays in the following general form:
D a + α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) + F ( t , X τ T ( θ ) ) ,
where h > 0 , J = [ a , ) , R + = [ 0 , ) , J 0 = ( a , ) , R + 0 = ( 0 , ) , n = { 1 , 2 , , n } , m 0 = m { 0 } , X ( t ) = c o l ( x 1 ( t ) , , x n ( t ) ) : J 0 R n , F ( t , X t T ) = c o l ( f 1 ( t , X t T ) , , f n ( t , X t T ) ) : J × P C R n , (the notation c o l means column and superscript T denotes the transposed vector), U ( t , θ ) = { u k j ( t , θ ) } k , j = 1 n : J × R R n × n , V ( t , θ ) = { v k j ( t , θ ) } k , j = 1 n : J × R R n × n , D a + α X ( t ) = c o l ( D a + α x 1 ( t ) , , D a + α x n ( t ) ) , and X t ( θ ) = X ( t + θ ) (Krasovskii-type functional) for t J and θ [ h , 0 ] . In addition, I , Θ R n × n denote the identity and the zero matrices, 0 R n denotes the zero vector, and I γ ( ( t a ) ) = ( t a ) γ I , γ [ 1 , 1 ] .
For Y : J × R R n × n , Y ( t , θ ) = { y k j ( t , θ ) } k , j = 1 n , we use the norm
| Y ( t , θ ) | = k , j = 1 n | y k j ( t , θ ) | , V a r [ b , c ] Y ( t , · ) = { V a r [ b , c ] y k j ( t , · ) } k , j = 1 n ,
P C ( [ h , 0 ] , R n ) = { Φ = ( ϕ 1 , , ϕ n ) T : [ h , 0 ] R n | Φ is piecewise continuous } and S Φ denotes the set of all jump points for any Φ ( t ) P C ( [ h , 0 ] , R n ) . We will use the following Banach spaces of initial functions:
PC = { Φ ( t ) = c o l ( φ 1 ( t ) , , φ n ( t ) ) : [ h , 0 ] R n | Φ ( t ) P C ( [ h , 0 ] , R n ) } , PC * = P C B V ( [ h , 0 ] , R n ) and C ( [ h , 0 ] , R n ) ,
endowed with the norm Φ = k n sup s [ h , 0 ] | ϕ k ( s ) | < . In addition, we assume that in all spaces the functions Φ ( t ) are right continuous at t S Φ (for Φ C this assumption ultimately holds).
For any Φ PC we define the initial condition for the system (1) as follows:
X ( t ) = Φ ( t a ) t [ a h , a ] , D a + α 1 X ( a + 0 ) = Φ ( 0 ) , h > 0 .
More details concerning other types of initial conditions are given in [8] for the case when the initial functions are continuous.
Let us introduce the auxiliary integral system (whose system is considered in detail in Lemma 4 below) for any Φ PC and t J 0 :
X ( t ) = Φ ( 0 ) ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] X ( t + θ ) + I 1 ( Γ ( α ) ) a t I α 1 ( t τ ) h 0 [ d θ U ( τ , θ ) ] X ( τ + θ ) + F ( τ , X τ T ) d τ .
Let b > a and γ [ 0 , 1 ] be arbitrary.
Definition 1
([4]). The function G ( t ) = ( g 1 ( t ) , , g n ( t ) ) C ( J 0 , R n ) , γ [ 0 , 1 ] will be called right γ-continuous at a if the function I γ ( t a ) G ( t ) = c o l ( ( t a ) γ g 1 ( t ) , , ( t a ) γ g n ( t ) ) C ( [ a , ) , R n ) , i.e., the function I γ ( t a ) G ( t ) is right continuous at a. By C γ we denote the real linear space of all right γ-continuous-at-a functions and for any b R + by C b γ the real Banach space
C b γ = { R ( t ) C ( ( a , b ] , R n ) | R ( t ) = G ( t ) | ( a , b ] , G ( t ) C γ }
with norm R b γ = k n sup t [ a , b ] ( t a ) γ | r k ( t ) | .
Definition 2
([11]). The vector function c o l X ( t ) = ( x 1 ( t ) , , x n ( t ) ) is a solution of IP (1), (2) in ( a , b ] ( J 0 ) if X ( t ) C b 1 α ( C 1 α ) satisfies the system (1) for all t ( a , b ] ( J 0 ) and the initial condition (2) too.
Definition 3
([11]). The vector function c o l X ( t ) = ( x 1 ( t ) , , x n ( t ) ) is a solution of IP (3), (2) in ( a , b ] ( J 0 ) if X ( t ) C b 1 α ( C 1 α ) satisfies system (3) for all t ( a , b ] ( J 0 ) and initial condition (2) too.
The hypotheses ( S ) stated below in Definition 4, as in the cases of systems with derivatives of integer order or Caputo-type fractional order, will play a major role in the solvability of IP (1), (2) (IP (3), (2)) (see [19,20,21]).
Definition 4
([19,20,21]). We say that for the kernels U , V : J × R R n × n the hypotheses ( S ) hold if the following conditions are fulfilled for any i m 0 and q l , l , m N :
( S 1)
The functions ( t , θ ) U ( t , θ ) , ( t , θ ) V ( t , θ ) are measurable in ( t , θ ) J × R and normalized so that U ( t , θ ) = 0 , V ( t , θ ) = 0 for θ 0 , U ( t , θ ) = U ( t , σ ) for θ σ ,
V ( t , θ ) = V ( t , τ ) for θ τ , σ , τ > 0 , h = max ( σ , τ ) for any t J .
( S 2)
For any fixed t J the kernels U ( t , θ ) and V ( t , θ ) are left continuous in θ on ( σ , 0 ) and ( τ , 0 ) , U ( t , · ) , V ( t , · ) B V 0 l o c ( J × R , R n × n ) and V a r [ h , 0 ] U ( t , · ) , V a r [ h , 0 ] V ( t , · ) B L 1 l o c ( J , R + ) , V a r [ h , 0 ] V ( t , · ) is uniformly bounded in t J and is uniformly nonatomic at zero [21] ( i.e., for every ϵ > 0 , there exists δ ( ϵ ) > 0 such that for any t J , we have that V a r [ δ , 0 ] V ( t , · ) < ϵ ).
( S 3)
For any fixed t J , the Lebesgue decomposition of the kernels U ( t , θ ) and V ( t , θ ) has the form:
U ( t , θ ) = U j ( t , θ ) + U a c ( t , θ ) + U s ( t , θ ) , V ( t , θ ) = V j ( t , θ ) + V a c ( t , θ ) + V s ( t , θ ) , U j ( t , θ ) = i m A i ( t ) H ( θ + σ i ( t ) ) , A i ( t ) = { a k j i ( t ) } k , j = 1 n B L 1 l o c ( J , R n × n ) ,
V j ( t , θ ) = l r A ¯ l ( t ) H ( θ + τ l ( t ) ) , A ¯ l ( t ) = { a ¯ k j l ( t ) } k , j = 1 n C ( J , R n × n ) ,
σ 0 ( t ) 0 , σ i ( t ) , τ l ( t ) C ( J , [ 0 , h ] ) , U a c ( t , · ) , V a c ( t , · ) A C ( [ h , 0 ] , R n × n ) , U s ( t , · ) , V s ( t , · ) C ( [ h , 0 ] , R n × n ) and H ( t ) is the Heaviside function.
( S 4)
The sets S Φ i = { t J | t σ i ( t ) S Φ , i m } , S Φ l = { t J | t τ l ( t ) S Φ , l r } do not have limit points and the relations
lim t t * h 0 | U i ( t , θ ) U i ( t * , θ ) | d θ = 0 , lim t t * h 0 | V l ( t , θ ) V l ( t * , θ ) | d θ = 0
hold for any t * J .
Definition 5
([20,21]). We say that the vector-valued functional F ( t , Y ) : J × P C R n satisfies the modified Caratheodory conditions ( C ) if the following conditions hold:
( C 1)
For almost all fixed t J the functional ( t , Y ) F ( t , Y ) is continuous in arbitrary Y PC ( PC * ) and for each fixed function Y PC ( PC * ) the function F ( t , Y ) B L 1 l o c ( J , R n ) .
( C 2)
(Local Lipschitz-type condition) For any ( t , Y ) J × PC and for some its vicinity O ( t , Y ) J × P C , there exists a function 𝓁 ( t ) B L 1 l o c ( J , R + ) such that the inequalities
| F ( t , Y 1 ) F ( t , Y 2 ) | 𝓁 ( t ) | Y 1 ( t ) Y 2 ( t ) |
hold for every ( t , Y 1 ) , ( t , Y 2 ) O ( t , Y ) .
In our investigations below, we will use the following auxiliary result:
Lemma 1
([4]). Let α ( 0 , 1 ) and let y ( t ) be a Lebesgue measurable function on J.
(a)
If there exists a.e. (almost everywhere) the limit lim t a + 0 [ ( t a ) 1 α y ( t ) ] = c R , then there also exists a.e. the limit ( D a α 1 y ) ( a + 0 ) = ( I a 1 α y ) ( a + 0 ) = lim t a + 0 ( I a 1 α y ) ( t ) = c Γ ( α ) .
(b)
If there exist a.e. the limit lim t a + 0 [ ( t a ) 1 α y ( t ) ] and lim t a + 0 ( I a 1 α y ) ( t ) = c * , then we have that lim t a + 0 [ ( t a ) 1 α y ( t ) ] = c * Γ ( α ) .
Definition 6
([22]). The low terminal a will be called a noncritical point (noncritical jump point) for some initial function Φ PC relative to the delay τ l ( t ) , l r if the equality τ l ( a ) = 0 implies that there exists a constant b 0 > a (eventually depending on τ l ( t ) ) such that t τ l ( t ) < a for t ( a , b 0 ] .
Definition 7
([22]). The low terminal a for a function Φ PC with a S Φ ( a S Φ ) will be called a critical point (critical jump point) relative to some delay τ l ( t ) , l r if the equality τ l ( a ) = 0 implies that there exists a constant b 0 > a (eventually depending on τ l ( t ) ) such that t a τ l ( t ) for t ( a , b 0 ] .
Lemma 2
([22]). Let the hypotheses ( S ) hold and l r τ l ( a ) > 0 . Then, there exists b ( a , a + h ] (eventually depending from τ l ), such that Φ ( t τ l ( t ) ) is continuous for t ( a , b ] .
Lemma 3
([22]). Let the hypotheses ( S ) hold and l r τ l ( a ) = 0 .
Then, for any initial function Φ PC with S Φ = { a } , one of the following statements holds:
(i)
The statement of Lemma 2 holds.
(ii)
The low terminal a is a critical jump point for Φ relative to the kernel V ( t , θ ) for some l 0 r .
Theorem 1
([23], Krasnosel’skii’s fixed point theorem). Let ( E , · E ) be a Banach space, H E be a nonempty, closed, and convex subset of E and the maps T , N : H E satisfy the following conditions:
(i)
The operator T is a contraction with constant γ ( 0 , 1 ) ;
(ii)
The operator K is continuous and the set K ( B ) is contained in a compact set;
(iii)
For any x , y H , we have that T x + K y H .
Then, there exists a z H with T z + K z = z .
Theorem 2
([24], Corollary 2). Suppose that α ( 0 , 1 ) and the following conditions hold:
1.
The functions a ( t ) , g ( t ) , y ( t ) L 1 l o c ( J , R + ) .
2.
The functions a ( t ) , g ( t ) are nondecreasing, g ( t ) is bounded on J, and y ( a ) = 0 .
3.
The inequality y ( t ) a ( t ) + g ( t ) a t ( t s ) β 1 y ( s ) d s holds for t J .
Then, for any t J we have that y ( t ) a ( t ) E α ( g ( t ) Γ ( α ) t α ) .
Remark 1.
The main difference between the systems (even in the linear case) with different types of fractional derivatives from the point of view of their applicability as model tools of real-world phenomena is the answer for a concrete type of derivative, whether the derivative of a constant is identically equal to zero or not. Typical representatives of these two classes are the fractional derivatives of Caputo and RL types. So, for the study of retarded and neutral systems with Caputo derivatives a lot of the ideas and techniques known from the systems with integer-order derivatives can be used, because of the continuity of the corresponding Krasovskii functional. Note that in some important special cases it is possible that the Krasovskii functional can have a jump only of the first kind at the initial point, but this obstacle can be overcome with the help of some variants of Banach’s contraction principle. In contrast, for retarded and neutral systems with RL derivatives we have at least in the technical aspect several complications, mainly based on the availability of a discontinuity (jump of second kind) of the solutions at the low terminal, a fact which greatly complicates the use of the Riesz theorem for the representation of linear continuous functionals (Krasovskii functionals) on C via a Lebesgue–Stieltjes integral. Even for the retarded systems with RL derivatives, more sophisticated techniques such as Weissinger theorem combined with application of the Mittag–Lafleur function must be used. We emphasize that the important neutral case with RL derivatives, the study of which is the aim of this article, is essentially more complicated even in comparison with the retarded systems with the same kind of derivatives, not only in the technical aspects, but in their applications as models.

3. The Initial Problem (1), (2) with Discontinuous Initial Functions

The next Lemma is well known for the case of systems with Caputo-type derivatives and for delayed systems with RL-type derivatives, and therefore, we will sketch only the differences appearing in the neutral case with RL-type derivatives and distributed delays.
Lemma 4.
Let the following conditions be fulfilled:
1.
The hypotheses S hold.
2.
For any Φ PC , the vector-valued functional F : J × PC R n B L 1 l o c ( J , R n ) .
Then, every solution X ( t ) of IP (1), (2) is a solution of IP (3), (2) and vice versa.
Proof. 
Let Φ PC be arbitrary and X ( t ) be the corresponding solution of IP (1), (2). Then, according to formula 2.1.40 in [5], we have
I a + α D a + α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) lim t a + 0 I a + 1 α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) Φ ( 0 ) ( t a ) α 1 + lim t a + 0 I a + 1 α h 0 [ d θ V ( t , θ ) ] X ( t + θ ) .
Then, obviously X ( t ) will be a solution of IP (3), (2) if we prove
lim t a + 0 I a + 1 α h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = 0 .
Condition ( S 2 ) implies that there exists a constant V h > 0 such that | V a r θ [ h , 0 ] V ( t , θ ) | V h for t [ a , a + h ] , and since V a r θ [ h , 0 ] V ( t , θ ) is uniformly nonatomic at zero in t, then for any ϵ > 0 there exists δ ( 0 , ϵ ) such that | V a r θ [ h , 0 ] V ( t , θ ) | < ϵ Γ ( α ) 2 Φ uniformly for t [ a , a + h ] . Then, for t [ a , a + δ ] we have
lim t a + 0 I a + 1 α h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = 1 Γ ( α ) lim t a + 0 ( t a ) 1 α h a t [ d θ V ( t , θ ) ] Φ ( t + θ ) + a t 0 [ d θ V ( t , θ ) ] X ( t + θ ) 1 Γ ( α ) lim t a + 0 ( t a ) 1 α Φ | V a r θ [ h , 0 ] V ( t , θ ) | + 1 Γ ( α ) lim t a + 0 ( t a ) 1 α sup θ [ δ , 0 ] | X ( t , θ ) | V a r θ [ δ , 0 ] V ( t , θ )
and, hence, for the first addend on the right-hand side of (4) the relation
lim t a + 0 ( t a ) 1 α Φ | V a r θ [ h , 0 ] V ( t , θ ) | = 0
holds. Since X ( t ) C 1 α , then there exists δ ¯ ( 0 , δ ) such that for t [ a , a + δ ] and any θ [ h , 0 ] we have that | X ( t + θ ) | | X ( t + θ ) X ( t ) | + | X ( t ) | 2 | X ( t ) | . Then, for the second addend on the right-hand side of (4) for t [ a , a + δ ] we obtain that
1 Γ ( α ) lim t a + 0 ( t a ) 1 α sup θ [ δ , 0 ] | X ( t , θ ) | V a r θ [ δ , 0 ] V ( t , θ ) 2 Γ ( α ) lim t a + 0 ( t a ) 1 α | X ( t ) | ϵ Γ ( α ) 2 Φ 2 Φ ( 0 ) ϵ 2 Φ < ϵ ,
which implies that
lim t a + 0 ( t a ) 1 α sup θ [ δ , 0 ] | X ( t , θ ) | V a r θ [ δ , 0 ] V ( t , θ ) = 0 .
Thus, X ( t ) is a solution of IP (3), (2).
Let X ( t ) is a solution of IP (3), (2), and then, applying the operator D a + α to both sides of (3) we obtain immediately the opposite statement. □
Modifying the approach in [11], for an arbitrary fixed number α ( 0 , 1 ) we introduce the following real linear space:
E = { G : [ a h , ) | G ( t ) | [ a h , a ] = Φ ( t a ) , Φ PC ; G | J 0 C 1 α ( J 0 , R n ) , lim t a + 0 ( t a ) 1 α G ( t ) = Φ ( 0 ) }
and for any b > a introduce the linear subspaces:
E b = { R : [ a h , b ] | R ( t ) = G | [ a h , b ] , G E }
endowed with the norm R b = Φ + G | ( a , b ] b 1 α .
Let Φ 0 PC with S Φ 0 = { a } be arbitrary fixed and define the nonempty, closed, and convex subset E b Φ 0 E b as follows:
E b Φ 0 = { R E b | R ( t ) | [ a h , a ] = Φ 0 ( t a ) , lim t a + 0 ( t a ) 1 α R ( t ) = Φ 0 ( 0 ) } .
It is simple to see that for any R E b we have that the norms R b = Φ 0 + R b 1 α and R b = R b 1 α are equivalent.
For any b > a , t ( a , b ] and R E b define the operator R : E b Φ 0 E b as follows:
( R R ) ( t ) = Φ ( 0 ) ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] R ( t + θ ) + I 1 ( Γ ( α ) ) a t I α 1 ( t s ) h 0 [ d U ( s , θ ) ] R ( s + θ ) d s + a t I α 1 ( t s ) F ( t , R s ) d s ,
with the additional conditions
( R R ) ( t ) = Φ 0 ( t a ) , t [ a h , a ] ; lim t b 0 ( R R ) ( t ) = ( R R ) ( b ) .
In addition, for arbitrary b > h , t ( a , b ] and R E b we define the operators T , K : E b Φ 0 E b as follows:
( T R ) ( t ) = Φ 0 ( 0 ) ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] R ( t + θ ) ( T R ) ( t ) = Φ 0 ( t a ) , t [ a h , a ] ; lim t b 0 ( T R ) ( t ) = ( T R ) ( b ) ,
( K R ) ( t ) = I 1 ( Γ ( α ) ) a t I α 1 ( t s ) h 0 [ d θ U ( s , θ ) ] R ( s + θ ) d s + a t I α 1 ( t s ) F ( s , R s ) d s ( K R ) ( t ) = 0 , t [ a h , a ] ; lim t b 0 ( K R ) ( t ) = ( K R ) ( b ) .
Then, the operator R , defined via (5), (6), for any R E b Φ 0 has the form
( R R ) ( t ) = ( T R ) ( t ) + ( K R ) ( t ) .
The next technical lemma is an immediate generalization of Theorem 4 in [25]. The lemma below is more appropriate for applications and can be used in the cases of RL derivatives too. For simplicity we will consider only the case when the initial function Φ PC has only one jump point, but it is simple to see that the proof provided below can be used in the case of an arbitrary but finite number of jumps too.
Lemma 5.
Let the following conditions be fulfilled:
1.
The hypotheses ( S ) hold and Φ PC with S Φ = { t j u m p } , t j u m p ( a h , a ] are arbitrary.
2.
The low terminal a is a jump point for Φ relative to the delay τ l 0 ( t ) for some l 0 r .
Then, lim t a + 0 A ¯ l 0 ( t ) = Θ , l 0 r .
Proof. 
We will consider the more complicated case when S Φ = { t j u m p = a } and a is a critical jump point. Then, since τ l ( t ) C ( J , [ 0 , h ] ) for all l r we have that lim t a + 0 t τ l 0 ( t ) = a . Let assume that there exists a sequence { t i } i = 1 ( a , b 0 ) such that t i τ l 0 ( t i ) = a , and hence, the set S Φ has a limit point which contradicts condition ( S 4 ) . Then, there exists b 0 ( a , a + h ) such that t τ l 0 ( t ) > a for any t ( a , b 0 ] . Let us assume that A ¯ l 0 ( a ) > 0 , ϵ 0 , A ¯ l 0 ( a ) 4 and δ ( o , ϵ ) is the number existing according to condition ( S 2 ) such that V a r [ δ , 0 ] V ( t , · ) < ϵ for any t J . Then, for any t J θ ( 0 , δ ] with t + θ a
A ¯ l 0 ( a ) H ( a ) = A ¯ l 0 ( a ) V d ( t , θ ) V d ( t , 0 ) + V a r θ [ δ , 0 ] V d ( t , θ ) = V a r θ [ δ , 0 ] V d l ( t , θ ) < ϵ < A ¯ l 0 ( a ) 4 ,
which is impossible. Thus, we have that A ¯ l 0 ( a ) = 0 , and hence, lim t a + 0 = A ¯ l 0 ( t ) = 0 .
The proofs of the cases when a is a noncritical jump point and when t j u m p a are simpler versions of the proof in the case when t j u m p = a (a is a critical jump point), and therefore, will be omitted. □
It is clear that without loss of generality it is possible to renumber all concentrated delays in the jump part V j ( t , θ ) = l r A ¯ l ( t ) H ( θ + τ l ( t ) ) and give the first l ¯ numbers, 0 l ¯ r of these delays, relative to which a is a critical jump point. Everywhere below we will assume that this renumbering is made.
Theorem 3.
Let the following conditions be fulfilled:
1.
The conditions of Lemma 4 hold.
2.
sup θ [ h , 0 ] | V θ ( t , θ ) | , sup θ [ h , 0 ] | U θ ( t , θ ) | B L 1 l o c ( J , R n × n ) , where V θ , U θ denote the partial derivatives with respect to θ.
3.
The condition ( C ) holds.
Then, for any Φ 0 PC with S Φ = { t j u m p } , t j u m p ( a h , a ] there exists b 0 > a such that IP (3), (2) has at least one solution X ( t ) E b Φ 0 with an interval of existence ( a , b 0 ] .
Proof. 
As above, we will consider the more complicated case when Φ PC with S Φ = { a } . Let b ( a , a + h ) is arbitrary. As a first step we will prove that T ( E b Φ 0 ) E b Φ 0 .
For any R E b Φ 0 from (6) and Lemmas 2 and 3 it follows that
( T R ) ( t ) = Φ 0 ( 0 ) ( t a ) α 1 + h a t [ d θ V ( t , θ ) ] Φ 0 ( t a + θ ) + a t 0 [ d θ V ( t , θ ) ] R ( t + θ ) = Φ 0 ( 0 ) ( t a ) α 1 + h a t [ d θ V ( t , θ ) ] Φ 0 ( t a + θ ) + l = l ¯ + 1 r A ¯ l ( t ) Φ 0 ( t τ l ( t ) a ) + l = 1 l ¯ A ¯ l ( t ) R ( t τ l ( t ) ) + a t 0 [ d θ V c ( t , θ ) ] R ( t + θ ) ,
where V c ( t , θ ) = V a c ( t , θ ) + V s ( t , θ ) , and t ( a , b ] .
The first addend on the right-hand side of (10) is continuous for any t J 0 , and hence, belongs to C 1 α . Since
lim t a + 0 ( t a ) 1 α Φ 0 ( 0 ) ( t a ) α 1 = Φ 0 ( 0 ) ,
we can prolong Φ 0 ( 0 ) ( t a ) α 1 as a function from E b Φ 0 . According to Lemma 1 in [11] and Lemma 4 in [22], since the initial function Φ 0 ( t a ) is a continuous function for t [ a h , a ) , the second and the third addends on the right-hand side of (10) are continuous functions at t ( a , b ] . Taking into account that a is a critical jump point relative to the delays with numbers l l ¯ , we conclude that the fourth addend is a continuous function at t ( a , b ] too. For the fifth addend, via the substitution t + θ = s we obtain
a t 0 [ d θ V c ( t , θ ) ] R ( t + θ ) = a t [ d s V c ( t , s t ) ] R ( s ) .
Because V c ( t , s t ) and R ( s ) for any t ( a , b ] are continuous at s ( a , t ] , then the fifth addend is a continuous function at t ( a , b ] by virtue of lemma 1 in [19]. Since the right-hand side of (10) is a continuous function at t = b , then the second additional relation in (7) holds too, and hence, from (7) it follows that T ( E b Φ 0 ) E b Φ 0 . Thus, T ( E b Φ 0 ) E b Φ 0 .
As a second step we will prove that the operator T is a contraction. Denote  L b = sup t [ a , b ] 𝓁 ( t ) ,
V b = max sup θ [ h , 0 ] t [ a , b ] , | θ V ( t , θ ) | , sup t [ a , b ] V a r θ [ h , 0 ] V ( t , θ )
and
U b = max sup t [ a , a + b ] V a r θ [ h , 0 ] U ( t , θ ) , sup θ [ h , 0 ] t [ a , b ] , | θ U ( t , θ ) | .
Then, for arbitrary R , R ¯ E b Φ 0 from (7) and Lemmas 2 and 3 it follows that
( T ( R ) ) ( t ) ( T ( R ¯ ) ) ( t ) = h 0 [ d θ V ( t , θ ) ] ( R ( t + θ ) R ¯ ( t + θ ) )
h 0 [ d θ V j ( t , θ ) ] ( R ( t + θ ) R ¯ ( t + θ ) ) + h 0 [ d θ V c ( t , θ ) ] ( R ( t + θ ) R ¯ ( t + θ ) )
l = 1 l ¯ A ¯ l ( t ) ( R ( t τ l ( t ) ) R ¯ ( t τ l ( t ) ) ) + l = l ¯ l = r A ¯ l ( t ) ( Φ 0 ( t a τ l ( t ) ) Φ 0 ¯ ( t a τ l ( t ) ) ) + a t 0 [ d θ V c ( t , θ ) ] ( R ( t + θ ) R ¯ ( t + θ ) ) l = 1 l ¯ A ¯ l ( t ) ( R ( t τ l ( t ) ) R ¯ ( t τ l ( t ) ) ) + t + θ = s a t [ d s V ( t , s t ) ] ( R ( s ) R ¯ ( s ) ) .
For the first addend on the right-hand side of (12), since τ l ( a ) = 0 and the low terminal a is a critical jump point for Φ relative to the delays τ l ( t ) , l l ¯ , by virtue of Lemma 5 there exists b 1 ( a , b ] such that for any t ( a , b 1 ] we have the estimation
l = 1 l ¯ A ¯ l ( t ) ( R ( t τ l ( t ) ) R ¯ ( t τ l ( t ) ) ) ( t a ) α 1 l = 1 l ¯ A ¯ l ( t ) ( t a ) 1 α ( R ( t τ l ( t ) ) R ¯ ( t τ l ( t ) ) ) ( t a ) 1 α 6 sup s ( a , t ] ( s a ) 1 α R ( s ) R ¯ ( s ) .
Hypothesis ( S 2 ) implies that for ϵ = 1 6 there exists δ ( 0 , ϵ ) such that
V a r θ [ δ , 0 ] V c ( t , θ ) < 1 6
for any t ( a , b 2 ] , where b 2 = min ( b 1 , a + δ ) , we obtain
| V ( t , s t ) | | V ( t , 0 ) | + V a r s t [ δ , 0 ] V c ( t , s t ) 1 6 .
Then, the second addend in (12) for any t ( a , b 2 ] has the estimation
a t [ d s V ( t , s t ) ] ( R ( s ) R ¯ ( s ) ) = a t [ s V ( t , s t ) ] ( s a ) α 1 ( s a ) 1 α ( R ( s ) R ¯ ( s ) ) d s ( t a ) α V b α sup s ( a , t ] ( s a ) 1 α R ( s ) R ¯ ( s )
and, hence, for t ( a , b 3 ] , where b 3 a , min 1 , b 2 , a + α 6 V b 1 α we obtain that
a t [ d s V ( t , s t ) ] ( R ( s ) R ¯ ( s ) ) 1 6 sup s [ a , t ] ( s a ) 1 α R ( s ) R ¯ ( s ) .
Multiplying both sides of (12) by ( t a ) 1 α and then from (12), (13), and (14) it follows that for any t ( a , b 3 ] we obtain that T ( R ) T ( R ¯ ) b 3 Φ 0 1 3 R R ¯ b 3 Φ 0 and, hence, T is a contraction operator in E b 3 Φ 0 .
For arbitrary R , R ¯ E b 3 Φ 0 , from (7) we obtain that
( K R ) ( t ) ( K R ¯ ) ( t ) n Γ ( α ) a t I α 1 ( t s ) h 0 [ d θ U ( s , θ ) ] ( R ( s + θ ) R ¯ ( s + θ ) ) d s + n Γ ( α ) a t I α 1 ( t s ) F ( s , R s ) F ( s , R ¯ s ) d s n 2 Γ ( α ) a t ( t s ) α 1 a s 0 [ d θ U ( s , θ ) ] ( R ( s + θ ) R ¯ ( s + θ ) ) d s + n 2 Γ ( α ) a t ( t s ) α 1 F ( s , R s ) F ( s , R ¯ s ) d s n 2 Γ ( α ) a t ( t s ) α 1 a s 0 | θ U ( s , θ ) | | R ( s + θ ) R ¯ ( s + θ ) | d s + n 2 Γ ( α ) a t ( t s ) α 1 𝓁 ( s ) R ( s ) R ¯ ( s ) d s n 2 U b Γ ( α ) a t ( t s ) α 1 a s | R ( η ) R ¯ ( η ) | d η d s + L b a t ( t s ) α 1 | R ( s ) R ¯ ( s ) | d s .
For the first addend on the right-hand side of (15), for any t ( a , b 3 ] the following estimation holds:
n 2 U b Γ ( α ) a t ( t s ) α 1 a s | R ( η ) R ¯ ( η ) | d η d s = n 2 U b Γ ( α ) a t a s | R ( η ) R ¯ ( η ) | d η d ( t s ) α n 2 U b Γ ( α ) a t ( t s ) α | R ( s ) R ¯ ( s ) | d s = t s = y a n 2 U b Γ ( α ) a t ( y a ) α ( y a ) α 1 ( y a ) 1 α | R ( t y + a ) R ¯ ( t y + a ) | d y ( 2 α ) 1 n 2 U b Γ ( α ) ( t a ) 2 α sup s [ a , t ] ( s a ) 1 α | R ( s ) R ¯ ( s ) | n 2 U b 2 Γ ( 1 + α ) sup s [ a , t ] ( s a ) 1 α | R ( s ) R ¯ ( s ) |
and analogously, for the second one we obtain
n 2 L b Γ ( α ) a t ( t s ) α 1 | R ( s ) R ¯ ( s ) | d s = n 2 L b Γ ( 1 + α ) a t | R ( s ) R ¯ ( s ) | d ( t s ) α n 2 L b Γ ( 1 + α ) a t ( t s ) α d | R ( s ) R ¯ ( s ) | n 2 L b ( t a ) 2 α 1 Γ ( 1 + α ) sup s [ a , t ] ( s a ) 1 α | R ( s ) R ¯ ( s ) | .
Multiplying both sides of (15) by ( t a ) 1 α , and then, from (15), (16), and (17), it follows that for any t ( a , b 3 ] we obtain that
K ( R ) K ( R ¯ ) b 3 Φ 0 n 2 2 Γ ( 1 + α ) ( U b + 2 L b ( b 3 a ) α ) R R ¯ b 3 Φ 0
and, hence, K : E b 3 Φ 0 E b 3 is a continuous operator.
Let r 0 R + 0 and R 0 E b 3 Φ 0 be arbitrary fixed and consider the ball
B ( R 0 , r 0 ) = { R E b 3 Φ 0 | R 0 R b 3 Φ 0 r 0 } E b 3 Φ 0 .
For any R B ( R 0 , r 0 ) , from (18) it follows that
K ( R ) b 3 Φ 0 K ( R 0 ) b 3 Φ 0 + K ( R ) K ( R 0 ) b 3 Φ 0 n 2 2 Γ ( 1 + α ) ( U b + 2 L b ( b 3 a ) α ) R R 0 b 3 Φ 0 + N ( R 0 ) b 3 Φ 0 n 2 r 0 2 Γ ( 1 + α ) ( U b + 2 L b ( b 3 a ) α ) + N ( R 0 ) b 3 Φ 0 = C ( r 0 , b 3 )
and, hence, the set K ( B ( R 0 , r 0 ) ) is uniformly bounded.
To apply Theorem 1 we must prove that the set N ( B ( R 0 , r 0 ) ) is equicontinuous, and hence, it is at least relative compact. For any R B ( R 0 , r 0 ) and arbitrary t 1 , t 2 ( a , b 3 ] we have the estimation
| ( K R ) ( t 2 ) ( K R ) ( t 1 ) | n 2 Γ ( α ) a t 1 ( ( t 2 s ) α 1 ( t 1 s ) α 1 ) h 0 [ d θ U ( s , θ ) ] R ( s + θ ) + F ( s , R s ) d s + n 2 Γ ( α ) t 1 t 2 ( t 2 s ) α 1 h 0 [ d θ U ( s , θ ) ] R ( s + θ ) + F ( s , R s ) d s n 2 Γ ( α ) a t 1 ( t 2 s ) α 1 ( t 1 s ) α 1 h 0 [ d θ U ( s , θ ) ] R ( s + θ ) + F ( s , R s ) d s + n 2 Γ ( α ) t 1 t 2 ( t 2 s ) α 1 h 0 [ d θ U ( s , θ ) ] R ( s + θ ) + F ( s , R s ) d s .
Since for t ( a , b 3 ] from conditions ( C ) and (19) we have that
h 0 [ d θ U ( s , θ ) ] R ( s + θ ) + F ( s , R s ) ( U b + L b ) C ( r 0 , b 3 ) ,
then from (20) and (21) for any t ( a , b 3 ] we obtain
| ( K R ) ( t 2 ) ( K R ) ( t 1 ) | n 2 C ( r 0 , b 3 ) ( U b + L b ) Γ ( α ) a t 1 ( t 2 s ) α 1 ( t 1 s ) α 1 d s + t 1 t 2 ( t 2 s ) α 1 d s n 2 C ( r 0 , b 3 ) ( U b + L b ) α Γ ( α ) ( 2 ( t 2 t 1 ) α + ( t 1 a ) α ( t 2 a ) α ) n 2 C ( r 0 , b 3 ) ( U b + L b ) Γ ( 1 + α ) ( 2 ( t 2 t 1 ) α + | ( t 1 a ) α ( t 2 a ) α | ) .
Then, since ( t a ) α is uniformly continuous at t [ a , a + h ] , for any ϵ > 0 there exists δ ( 0 , ϵ ) such that if | t 2 t 1 | < δ we have that
| ( t 1 a ) α ( t 2 a ) α | < ϵ Γ ( 1 + α ) n 2 C ( r 0 , b 3 ) ( U b + L b )
and, hence, from (22) it follows that the set K ( B ( R 0 , r 0 ) ) is equicontinuous. Thus, by virtue of the Arzella–Ascoli theorem, K ( B ( R 0 , r 0 ) ) is at least relative compact. For arbitrary R , R ¯ E b 3 Φ 0 we consider the sum ( T R ) ( t ) + ( K R ¯ ) ( t ) . Taking into account that lim t a + 0 ( K R ¯ ) ( t ) = 0 and can be prolonged as 0 for any t [ a h , a ] , from (7) it follows that ( T R ) ( t ) + ( K R ¯ ) ( t ) = Φ 0 ( t a ) . Since lim t a + 0 ( K R ¯ ) ( t ) = 0 , then lim t a + 0 ( t a ) 1 α ( K R ¯ ) ( t ) = 0 too. Considering the fact that
lim t a + 0 ( t a ) 1 α ( T R ) ( t ) = Φ 0 ( 0 ) + lim t a + 0 ( t a ) 1 α h 0 [ d θ V ( t , θ ) ] R ( t + θ ) = Φ 0 ( 0 ) ,
then, we obtain that
lim t a + 0 ( t a ) 1 α ( ( T R ) ( t ) + ( K R ¯ ) ( t ) ) = Φ 0 ( t a )
and, hence, ( T R ) ( t ) + ( K R ¯ ) ( t ) E b 3 Φ 0 . Applying Theorem 1, we obtain that there exists at least one fixed point R * E b 0 Φ 0 , i.e., R * ( t ) = ( T R * ) ( t ) + ( K R * ) ( t ) , t ( a , b 0 ] , b 0 = b 3 . □
Theorem 4.
Let the conditions of Theorem 3 be fulfilled.
Then, for any Φ 0 PC with S Φ = { t j u m p } , t j u m p ( a h , a ] the solution X ( t ) E b 0 Φ 0 of IP (3), (2) is unique in the interval of its existence ( a , b 0 ] .
Proof. 
Assume the contrary, that IP (3), (2) possess two different solutions X 1 ( t ) , X 2 ( t ) E b Φ 0 . Then, from (3) for any t ( a , b 0 ] we obtain that
X 1 ( t ) X 2 ( t ) = Y ( t ) ( T Y ) ( t ) + ( K Y ) ( t ) ( T Y ) ( t ) + ( K Y ) ( t ) q | Y ( t ) | + n 2 Γ ( α ) a t ( t τ ) α 1 h 0 θ U ( τ , θ ) Y ( τ + θ ) d θ + F ( τ , X τ 1 , T ) F ( τ , X τ 2 , T ) d τ q | Y ( t ) | + n 2 ( b 0 a ) 1 α U b 0 Γ ( α ) a t ( t τ ) α 1 ( sup s [ a , τ ] | Y ( s ) | + 𝓁 ( τ ) | Y ( s ) | ) d τ .
Since (2) implies that | Y ( a ) | = 0 , then from (23) it follows that
sup s [ a , t ] | Y ( s ) | n 2 ( b 0 a ) 1 α U b 0 Γ ( α ) ( 1 q ) a t ( t τ ) α 1 ( sup s [ a , τ ] | Y ( s ) | + 𝓁 ( τ ) | Y ( s ) | ) d τ n 2 ( b 0 a ) 1 α U b 0 ( 1 + L b 0 ) Γ ( α ) ( 1 q ) a t ( t τ ) α 1 sup s [ a , τ ] | Y ( s ) | d τ .
Denoting g ( t ) = g 0 = n 2 ( b 0 a ) 1 α U b 0 ( 1 + L b 0 ) Γ ( α ) ( 1 q ) , a ( t ) 0 and applying Theorem 2, we obtain for t ( a , b 0 ] the estimation sup s [ a , t ] | Y ( s ) | a ( t ) E α ( g 0 Γ ( α ) t α ) 0 which contradicts our assumption. The proof of the case when t j u m p a is almost the same as the proof in the case when t j u m p = a , and therefore, will be omitted. □
Definition 8.
We say that the solution of IP (3), (2) X 2 ( t ) E b 2 Φ 0 is a continuation of a solution X 1 ( t ) E b 1 Φ 0 of IP (3), (2) if ( a , b 1 ] ( a , b 2 ] and X 2 ( t ) X 1 ( t ) on ( a , b 1 ] . The solution X M a x ( t ) E b M a x Φ 0 of IP (3), (2) which coincides with all of its continuations will be called the maximal solution.
Theorem 5.
Let the conditions of Theorem 3 be fulfilled.
Then, for any Φ 0 PC with S Φ = { t j u m p } , t j u m p ( a h , a ] the unique solution X ( t ) E b 0 Φ 0 with interval of existence ( a , b 0 ] of IP (3), (2) can be continued as a unique solution X ( t ) E of IP (3), (2) with interval of existence J 0 .
Proof. 
Assume the contrary, that there exists a maximal solution X M a x ( t ) E b M a x Φ 0 of IP (3), (2) with the interval of its existence ( a , b M a x ] (which is closed from the right), b M a x < , and define the operator R : E b Φ 0 E b with b > b M a x via (5) with the additional conditions ( R R ) ( t ) = X M a x ( t ) , t [ a h , a ] ; lim t b 0 ( R R ) ( t ) = ( R R ) ( b ) . Then, we can prove, fully analogous to the proof of Theorem 3, that for any Φ 0 PC with S Φ = { t j u m p } , t j u m p ( a h , a ] there exists b ¯ > b M a x such that IP (3), (2) has at least one solution X b ¯ ( t ) E b ¯ Φ 0 with interval of existence ( a , b ¯ ] , and obviously, X b ¯ ( t ) X M a x ( t ) for any t ( a h , b M a x ] . Thus, we obtain that X b ¯ ( t ) is a continuation of X M a x ( t ) , which is impossible.
Let assume that that there exists a maximal solution X M a x ( t ) of IP (3), (2) and its interval of existence is open and finite, i.e., ( a , b M a x ) , b M a x < . Then, since X M a x ( t ) satisfies (3) for any t ( a , b M a x ) and the additional condition lim t b M a x 0 ( R R ) ( t ) = ( R R ) ( b ) too, we conclude that X M a x ( t ) satisfies (3) for any t ( a , b M a x ] , which contradicts our assumption. Thus, we have that b M a x = , which completes the proof. □
Remark 2.
We note that from Theorems 3–5 it follows that the requirement in Theorem 2 [11], that in the Lebesgue decomposition of U ( t , θ ) a singular part does not exist, is unnecessary. So, the results proved in Theorems 3–5 generalize the statement of Theorem 2 [11] even in the retarded case.

4. The Linear Case F ( t , X t T ) = F ( t ) – Fundamental Matrices and Integral Representation of the Solutions

In this section, we establish the existence of two different types of fundamental matrices and of their bases; we obtain two types integral representations for the solutions of IP (24), (2) for different kinds of initial functions. In the case when F ( t , X t T ) = F ( t ) , system (1) becomes an inhomogeneous linear system and has the following form:
D a + α Y ( t ) h 0 [ d θ V ( t , θ ) ] Y ( t + θ ) = h 0 [ d θ U ( t , θ ) ] Y ( t + θ ) + F ( t ) ,
where Y ( t ) = ( y 1 ( t ) , , y n ( t ) ) T : J R n and F ( t ) = ( f 1 ( t ) , , f n ( t ) ) T B L 1 l o c ( J , R n ) . Consider the corresponding homogeneous linear system of (1) and (24):
D a + α Y ( t ) h 0 [ d θ V ( t , θ ) ] Y ( t + θ ) = h 0 [ d θ U ( t , θ ) ] Y ( t + θ ) ,
and following [11] introduce for any j n and s ¯ J the initial function
Y j ( t , s ¯ ) = 0 , t [ s ¯ h , s ¯ ) ; Y ( t , s ¯ ) = I j , t = s ¯ ,
where I j denotes the j-th column of the identity matrix I .
From Theorem 5 it follows that for any j n and s ¯ J IP (25), (2) has a unique solution H j ( t , s ) , which satisfies (25) for any t ( s ¯ , ) and the initial condition (2) with the initial function (26) too. Then, the matrix H ( t , s ) = ( H 1 ( t , s ) , , H n ( t , s ) ) will be called fundamental matrix of the system (25).
Additionally, we introduce the following initial function for any j n and s ¯ [ a h , a ] :
Y j ( t , s ¯ ) = I j , a h s ¯ t a ; Y ( t , s ¯ ) = 0 , t < s ¯ ; Y ( t , s ¯ ) = 0 , s ¯ < a h .
As above, for any j n and s ¯ J we can conclude that IP (24), (2) has a unique solution Q j ( t , s ) which satisfies (24) for any t J 0 and the initial condition (2) with the initial function (27) too. The matrix Q ( t , s ) = ( Q 1 ( t , s ) , , Q n ( t , s ) ) we will call the generalized fundamental matrix.
Both matrices play crucial roles in the construction of different kinds of integral representations of the solutions of IP (24), (2) and IP (25), (2), having a lot of applications in the qualitative theory of the fractional systems.
It is clear that H ( t , a ) Q ( t , a ) if and only if s = a . First, we will establish some analytical properties of the matrices Q ( t , a ) and H ( t , a ) , needed to obtain appropriate integral representations of the solutions of (24) and (25).
As in the case of Caputo-type derivatives [26], for any Φ PC * we define the function
X ˜ ( t ) = a h a [ d s Φ ˜ ( s a ) ] Q ( t , s ) ,
where Φ ˜ ( s a ) = Φ ( s a ) , s ( a h , a ] , and Φ ˜ ( h ) = 0 , s = a h . It is necessary to establish for any j n that the solutions Q j ( t , s ) to IP (25), (2) are Lebesgue integrable in s on s [ a h , a ] for any fixed t J 0 .
Remark 3.
Note that to prove the statement that Q j ( t , s ) , j n are Lebesgue integrable in s on s [ a h , a ] for any fixed t J 0 is not trivial, in contrast with the Caputo case where this fact can be proved via elementary application of a generalized Bellman–Gronwall inequality (see Theorem 2). The main tool in all these proofs of theorems devoted to integral representations of the solutions of fractional systems with Caputo-type derivatives is the Fubini theorem. Thus, from our point of view it is important to justify its correct application. The greatest obstacle in the neutral case with RL derivatives is the problem with the integrability of the integrands relative to the product measures, which are defined as the products of two Lebesgue–Stieltjes measures. Note that this problem arises in our case, since the integrand is with a singularity of order 1 α at the point a in contrast to the case of neutral systems with Caputo-type fractional derivatives, in which case the integrand is a Lebesgue measurable and locally bounded function. This difference leads at least to some additional technical complications, but not only that.
To overcome this obstacle, we need the statements of the next two lemmas.
Lemma 6.
Let the following conditions be fulfilled:
1.
Conditions 1 and 2 of Theorem 3 hold.
2.
sup t J 0 V a r θ [ h , 0 ] V ( t , θ ) < 1
Then, for any compact interval [ a ¯ , b ] J 0 , t [ a ¯ , b ] , and j n we have that
Q j ( t , · ) B L 1 l o c ( [ a h , a ] , R n ) .
Proof. 
Let t J 0 and j n be arbitrary fixed. Then, from (25), since F ( τ , X τ T ) 0 we have
| Q j ( t , s ) | | I j 1 [ s , a ] ( t a ) | ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] Q j ( t , s ) + n 2 Γ ( α ) a t ( t τ ) α 1 h 0 [ d θ U ( τ , θ ) ] Q j ( τ + θ , s ) d τ | I j 1 [ s , a ] ( t a ) | ( t a ) α 1 + V sup η ( a , t ] Q j ( t , s ) = n 2 U b Γ ( α ) ( 1 V ) a t ( t τ ) α 1 sup η ( a , τ ] Q j ( η , s ) d τ
and, hence, by virtue of Theorem 2 we obtain
sup s [ h , 0 ] sup η ( a , t ] Q j ( η , s ) | I j 1 [ s , a ] ( t a ) | ( t a ) α 1 1 V + n 2 U b Γ ( α ) ( 1 V ) a t ( t τ ) α 1 sup s [ h , 0 ] sup η ( a , τ ] Q j ( η , s ) d τ | I j 1 [ s , a ] ( t a ) | ( t a ) α 1 1 V E α n 2 U b ( t a ) α ( 1 V ) ,
where in (29) E α ( z ) = k = 1 z k Γ ( α k + 1 ) is the one-parameter Mittag–Laffler function. This completes the proof. □
By virtue of Lemma 6 we can state that the integral on the right-hand side of (28) is correctly defined.
Corollary 1.
Let the conditions of Lemma 6 hold.
Then, for arbitrary fixed s ¯ [ h , 0 ) , Q ( t , s ) is continuous in t for t a [ h , s ¯ ) ( s ¯ , 0 ) or when s ¯ = 0 for t a [ h , 0 ) . At t a = s ¯ , Q ( t , s ) has a first-kind jump.
Lemma 7.
Let the conditions of Lemma 6 hold.
Then, for any t J 0 we have that Q ( t , · ) B V ( [ a h , a ] , R n × n ) .
Proof. 
Let Φ * ( s ) = ( φ 1 * ( s ) , , φ n * ( s ) ) T PC * with inf s [ a h , a ] φ j * ( s ) c 0 > 0 , j n . Then, integrating by parts formally the integral in (28) we have that
a h a [ d s Φ ˜ * ( s a ) ] Q ( t , s ) = Φ ˜ * ( 0 ) Q ( t , a ) Φ ˜ ( h ) ] Q ( t , a h ) a h a [ d s Q ( t , s ) ] Φ ˜ * ( s a ) ,
and since the integral on the left-hand side of (30) exists, then the integral on the right-hand side exists too (see [27], point 5, page 229). Let t J 0 be arbitrary, ϵ > 0 and q = { s 0 = a h , s 1 , , s q = a } [ a h , a ] be an arbitrary partition. Then, there exists a number q ϵ N such that for any q q ϵ we have that
a h a [ d s Q ( t , s ) ] Φ ˜ * ( s a ) k = 1 q Φ ˜ * ( s k a ) Q j ( t , s k ) Q j ( t , s k 1 ) < ϵ ,
and hence,
0 c 0 k = 1 q Q j ( t , s k ) Q j ( t , s k 1 ) inf s [ a h , a ] Φ ˜ * ( s a ) k = 1 q Q j ( t , s k ) Q j ( t , s k 1 ) k = 1 q Φ ˜ * ( s k ) Q j ( t , s k ) Q j ( t , s k 1 ) a h a [ d s Q ( t , s ) ] Φ ˜ * ( s a ) + ϵ .
Thus, Q ( t , · ) B V ( [ a h , a ] , R n × n ) . □
Theorem 6.
Let the conditions of Lemma 6 hold.
Then, for arbitrary initial function Φ PC * the function X ˜ ( t ) , defined via (28), is the unique solution of IP (25), (2) with interval of existence t J 0 .
Proof. 
The proof is based on ideas used in [22]. Thereupon, we will emphasize in detail those differences that arise from the influence on the neutral system of the Riemann–Liouville-type derivatives.
Theorem 5 implies that Q ( · , s ) C ( J 0 , R n × n ) for any s [ a h , a ] . Then, for the function X ˜ ( t ) , defined via (28), by virtue of Lemma 1 in [22] we have that X ˜ ( t ) C ( J 0 , R n ) too.
More concretely, for arbitrary fixed t J 0 the kernels V ( t , θ ) of U ( t , θ ) define two Lebesgue–Stieltjes measures μ V θ = μ V ( ( θ 1 , θ 2 ] ) = V ( t , θ 2 ) V ( t , θ 1 ) and μ U θ = μ U ( ( θ 1 , θ 2 ] ) = U ( t , θ 2 ) U ( t , θ 1 ) for any ( θ 1 , θ 2 ] [ h , 0 ] , as well as the function Φ ¯ ( s ) defining a measure μ Φ ¯ s = μ Φ ¯ ( ( s 1 , s 2 ] ) = Φ ¯ ( s 1 ) Φ ¯ ( s 2 ) for any ( s 1 , s 2 ] [ a h , a ] . We introduce the product measures
( μ V θ × μ Φ ¯ s ) ( ( θ 1 , θ 2 ] × ( s 1 , s 2 ] ) = μ V ( ( θ 1 , θ 2 ] μ Φ ¯ ( ( s 1 , s 2 ] )
and
( μ U θ × μ Φ ¯ s ) ( ( θ 1 , θ 2 ] × ( s 1 , s 2 ] ) = μ U ( ( θ 1 , θ 2 ] μ Φ ¯ ( ( s 1 , s 2 ] )
of the rectangles in P = [ h , 0 ] × [ a h , a ] . Then, to use the proposition 5.15 in [28] we need to prove that the relations P Q ( t + θ , s ) μ U θ × μ Φ ¯ s < and P Q ( t + θ , s ) μ V θ × μ Φ ¯ s < hold. Since according to Lemmas 5 and 6 we have Q ( t , · ) = B V ( [ a h , a ] , R n × n ) and Q ( t + θ , · ) = L 1 l o c ( [ a h , ) , R n × n ) for any t J 0 and θ [ h , 0 ] concerning the first argument, then taking into account that μ V θ × μ Φ ¯ s = μ Φ ¯ s × μ V θ and μ U θ × μ Φ ¯ s = μ Φ ¯ s × μ U θ we obtain
P Q ( t + θ , s ) μ U θ × μ Φ ¯ s = h 0 [ d θ U ( t , θ ) ] a h a [ d s Φ ( s a ) ] Q ( t + θ , s ) < , P Q ( t + θ , s ) μ V θ × μ Φ ¯ s = h 0 [ d θ V ( t , θ ) ] a h a [ d s Φ ( s a ) ] Q ( t + θ , s ) < .
Then, applying proposition 5.15 in [28] we have for t J 0 that
h 0 [ d θ V ( t , θ ) ] X ˜ ( t + θ ) = h 0 [ d θ V ( t , θ ) ] a h a [ d s Φ ( s a ) ] Q ( t + θ , s ) = a h a [ d s Φ ( s a ) ] h 0 [ d θ V ( t , θ ) ] Q ( t + θ , s )
and in the same way we obtain
h 0 [ d θ U ( t , θ ) ] X ˜ ( t + θ ) = a h a [ d s Φ ( s a ) ] h 0 [ d θ U ( t , θ ) ] Q ( t + θ , s ) .
From (25) and (31) via the Fubini theorem it follows that
D a + α X ˜ j ( t ) h 0 [ d θ V ( t , θ ) ] X ˜ j ( t + θ ) = 1 Γ ( 1 α ) d d t a t ( t η ) α a h a [ d s Φ ˜ ( s a ) ] Q ( η , s ) d η 1 Γ ( 1 α ) d d t a t ( t η ) α h 0 [ d θ V ( η , θ ) ] a h a [ d s Φ ˜ ( s a ) ] Q ( η + θ , s ) d η = a h a [ d s Φ ˜ ( s a ) ] 1 Γ ( 1 α ) d d t a t ( t η ) α Q ( η , s ) d η a h a [ d s Φ ˜ ( s a ) ] 1 Γ ( 1 α ) d d t a t ( t η ) α h 0 [ d θ V ( η , θ ) ] Q ( η + θ , s ) d η = a h a [ d s Φ ˜ ( s a ) ] D α Q ( t , s ) h 0 [ d θ V ( t , θ ) ] Q ( t + θ , s ) .
Hence, from (25), (32), and (33) it follows that X ˜ ( t ) , defined via (28), satisfies (25) for any t J 0 .
Let s ˜ [ a h , a ] be an arbitrary number, and then, from (28) when t = s ˜ , it follows that
X ˜ ( s ˜ ) = a h a [ d s Φ ˜ ( s a ) ] Q ( t , s ) = a h s ˜ [ d s Φ ˜ ( s a ) ] Q ( t , s ) + s ˜ a [ d s Φ ˜ ( s a ) ] Q ( t , s ) = a h s ˜ [ d s Φ ˜ ( s a ) ] I = Φ ˜ ( s ˜ ) Φ ˜ ( a h ) = Φ ˜ ( s ˜ ) ,
i.e., X ˜ ( t ) satisfies the initial condition (2), and hence, it is the unique solution of the IP (25), (2) with interval of existence t J 0 . □
Following the idea in [22], we introduce
X ( t ) = a t H ( t , s ) D a + 1 α F 0 ( s ) d s ,
where F 0 ( t ) F ( t ) , t J 0 , and F 0 ( 0 ) = 0 .
Theorem 7.
Let the conditions of Lemma 6 hold.
Then, the function X ( t ) introduced via (34) is the unique solution of IP (24), (2) for the initial function Φ ( t a ) 0 , t [ a h , a ] and interval of existence t J 0 .
Proof. 
As in the theorem above, the proof is based on some ideas used in [15]. Therefore, we will only sketch the similar parts and we will emphasize in detail those differences that arise from the influence on the neutral system of the Riemann–Liouville-type derivatives.
For any t J 0 and θ [ h , 0 ] , since H ( t , s ) = 0 when t < s , we have that
t + θ t H ( t + θ , s ) D a + 1 α F 0 ( s ) d s = 0 .
Then, substituting X ( t ) on the left-hand side of (24) for any t J 0 via the Fubini theorem and using formula (2.211) in [5] and (35) we obtain that
D a + α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = 1 Γ ( 1 α ) d d t a t ( t η ) α a η H ( η , s ) D a + 1 α F 0 ( s ) d s d η 1 Γ ( 1 α ) d d t a t ( t η ) α h 0 [ d θ V ( η , θ ) ] a η + θ H ( η + θ , s ) D a + 1 α F 0 ( s ) d s d η = 1 Γ ( 1 α ) d d t a t s t ( t η ) α H ( η , s ) d η D a + 1 α F 0 ( s ) d s 1 Γ ( 1 α ) d d t a t ( t η ) α h 0 [ d θ V ( η , θ ) ] a η H ( η + θ , s ) D a + 1 α F 0 ( s ) d s d η = 1 Γ ( 1 α ) d d t a t a t ( t η ) α H ( η , s ) d η D a + 1 α F 0 ( s ) d s 1 Γ ( 1 α ) d d t a t a t ( t η ) α h 0 [ d θ V ( η , θ ) ] H ( η + θ , s ) d η D a + 1 α F 0 ( s ) d s = a t D a + 1 α F 0 ( s ) D a + α ( H ( t , s ) ) d s a t D a + 1 α F 0 ( s ) D a + α h 0 [ d θ V ( t , θ ) ] H ( t + θ , s ) d s + lim s t 0 D a + α 1 R L D a + 1 α F 0 ( s ) H ( t , s ) lim s t 0 D a + α 1 R L D a + 1 α F 0 ( s ) h 0 [ d θ V ( t , θ ) ] H ( t + θ , s ) .
For the third addend on the right-hand side of (36), we have that
lim s t 0 D a + α 1 R L D a + 1 α F 0 ( s ) ( H ( t , s ) ) = lim s t 0 1 Γ ( α ) a t ( t τ ) α H ( τ , s ) D a + 1 α F 0 ( s ) d τ = 1 Γ ( α ) a t ( t τ ) α lim s t 0 H ( τ , s ) D a + 1 α F 0 ( s ) d τ = 1 Γ ( α ) a t ( t τ ) α H ( τ , τ ) D a + 1 α F 0 ( τ ) d τ = 1 Γ ( α ) a t ( t τ ) α I D a + 1 α F 0 ( τ ) d τ = D a + α 1 R L D a + 1 α R L F 0 ( t ) = F 0 ( t ) F 0 ( 0 ) = F 0 ( t ) .
Since V ( τ , 0 ) = Θ and H ( t + θ , t ) = 0 , for any t J 0 and θ [ h , 0 ) we obtain
lim s t 0 D a + α 1 R L D a + 1 α F 0 ( s ) h 0 [ d θ V ( τ , θ ) ] H ( τ + θ , s ) = lim s t 0 1 Γ ( α ) a t ( t τ ) α D a + 1 α F 0 ( s ) h 0 [ d θ V ( τ , θ ) ] H ( τ + θ , s ) d τ = 1 Γ ( α ) a t ( t τ ) α D a + 1 α F 0 ( τ ) h 0 [ d θ V ( τ , θ ) ] H ( τ + θ , τ ) d τ = 0 ,
and hence, from (36), (37), and (38) it follows that
D a + α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) = a t D a + 1 α F 0 ( s ) D a + α H ( t , s ) h 0 [ d θ V ( t , θ ) ] H ( t + θ , s ) d s + F 0 ( t ) .
Similarly as above, for the right-hand side of (24) we obtain
h 0 [ d θ U ( t , θ ) ] X ( t + θ ) + F 0 ( t ) = h 0 [ d θ U ( t , θ ) ] a t + θ H ( t + θ , s ) D a + 1 α F 0 ( s ) d s + F 0 ( t ) = a t D a + 1 α F 0 ( s ) h 0 [ d θ U ( t , θ ) ] H ( t , s ) d s + F 0 ( t ) ,
and since H ( t , s ) is a fundamental matrix, then the statement of the theorem follows from (39) and (40). □
Corollary 2.
Let the following conditions be fulfilled:
1.
The conditions of Lemma 6 hold.
2.
For any initial function Φ P C * and F ( t ) B L 1 l o c ( J , R n ) , the functions X ˜ ( t ) and X ( t ) are defined via (28) and (34), respectively.
Then, the function X ( t ) = X ˜ ( t ) + X ( t ) is the unique solution of IP (24), (2) and has the following representation for any t J 0 :
X ( t ) = a t H ( t , s ) D a + 1 α F ( s ) d s + a h a Q ( t , s ) d s Φ ˜ ( s a ) ,
where Φ ˜ ( s a ) Φ ( s a ) for s [ a h , a ) and Φ ˜ ( a h ) = 0 .
Proof. 
It follows from the superposition principle and Theorems 5 and 6. □

5. Weighted Stabilities

It is worth noting that the standard definitions for stability in the Lyapunov or Ulam–Hyers senses introduced for the systems with integer-order derivatives (without or with delays) are applicable directly for systems with fractional derivatives only when the fractional derivative of a constant is equal to zero (as Caputo-type derivatives, etc.). As we have mentioned previously [11], in general when the fractional derivative of a constant is not equal to zero (as for RL-type derivatives), then the standard definitions are not applicable since the solutions of the systems with RL-type derivatives have singularity at a + 0 (the low terminal) of power order α 1 . That is why new types of definitions for the different kinds of stabilities applicable to systems with RL-type derivatives are needed. The aim of this section is to introduce definitions of weighted stability in the Lyapunov sense, as well as UH and UHR weighted stabilities for fractional systems (equations) with RL-type derivatives and to study the relations with the classical definitions and with the concept “stability in time” in the Lyapunov sense, introduced in [9].
Let assume that F ( t , 0 ) 0 for any t J . Below, we recall the classical definitions for Lyapunov stability, applicable only for the case when the fractional derivatives are, in the Caputo sense, of order α ( 0 , 1 ] (and not applicable for RL derivatives):
Definition 9
([29]). The zero solution of IP (1), (2) (with derivatives in the Caputo sense!) is said to be:
(a)
Stable if for a given low terminal a R and any ϵ > 0 there is a δ ( ϵ , a ) > 0 such that for each initial function Φ PC ( PC * ) with Φ < δ the corresponding solution X ( t ) satisfies for each t J the inequality | X ( t ) | ϵ .
(b)
Stable (uniformly) if for a given low terminal a R and any ϵ > 0 there is a δ ( ϵ ) > 0 such that for every initial function Φ PC ( PC * ) , with Φ < δ the corresponding solution X ( t ) satisfying for each t J the inequality | X ( t ) | ϵ .
(c)
Locally asymptotically stable (LAS) if for a given low terminal a R there is a Δ ( a ) > 0 such that for every initial function Φ PC ( PC * ) with Φ < Δ the corresponding solution X ( t ) is stable and lim t | X ( t ) | = 0 .
(d)
Globally asymptotically stable (GAS) if for every for every initial function Φ PC ( PC * ) the corresponding solution X ( t ) is stable and lim t | X ( t ) | = 0 .
Both definitions below are applicable not only for initial problems with Caputo-type derivatives, but also for initial problems with fractional derivatives in the RL sense.
Define for any Φ PC the weighted function w ( t ) : J R as follows:
w ( t ) = 1 , t = a , ( t a ) 1 α , t ( a , a + 1 ] , 1 , t > a + 1 .
Definition 10.
The zero solution of IP (1), (2) is said to be weighted stable (W-stable), weighted uniformly stable (W-uniformly stable), weighted locally asymptotically stable (W-LAS stable), or weighted globally asymptotically stable (W-GAS stable) if for the corresponding solution X ( t ) the product w ( t ) X ( t ) satisfies the conditions (a), (b), (c), or (d) from Definition 9.
Definition 11
([9]). The zero solution of IP (1), (2) is said to be:
(i)
Stable in time (Lyapunov in time stable) if for an arbitrary ϵ > 0 there exists a point t ϵ J 0 and number δ ( ϵ , t ϵ ) > 0 such that for any initial function Φ PC ( PC * ) with Φ < δ , the corresponding solution X ( t ) satisfies for each t t ϵ the inequality | X ( t ) | ϵ .
(ii)
Asymptotically stable in time if for any initial function Φ PC ( PC * ) the corresponding solution X ( t ) is stable in time and additionally lim t | X ( t ) | = 0 .
Remark 4.
It is simple to see that if α = 1 (the order of differentiation) then Definition 10 for the weighted stability coincides with the classical one. In addition, when a solution is weighted stable then we can take t ϵ = a + 1 for any ϵ > 0 in Definition 11, and hence, the weighted stability implies stability in time, as well as the weighted global asymptotically stability implying asymptotic stability in time.
Let ϵ > 0 be an arbitrary number and consider the inequality
D a + α X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) h 0 [ d θ U ( t , θ ) ] X ( t + θ ) + F ( t , X τ T ( θ ) ) ϵ .
Definition 12.
The function Y ( t ) : [ a h , ) R n is a solution of IP (42), (2) in J 0 if Y ( t ) | t J 0 C 1 α satisfies the inequality (42) for t J 0 and the condition (2) for t [ a h , a ] , with initial function Y ( t a ) = Φ Y ( t a ) PC ( PC * ) .
Definition 13.
The system (1) is said to be Ulam–Hyers weighted stable in J 0 if for any ϵ > 0 there exists a number δ ( 0 , ϵ ) and constant C ϵ > 0 such that for any solution Y ( t ) of (42) there exists an initial function Φ ( t ) PC ( PC * ) , with | Φ Y ( t a ) Φ ( t a ) | δ for t [ a h , a ] and for the corresponding unique solution X ( t ) of IP (1), (2), with initial function Φ ( t ) , the inequality
w ( t ) Y ( t ) | J X ( t ) C ϵ ϵ
holds for any t J 0 .
The first result in this section clarifies the relations between the boundedness of all solutions and the weighted stabilities in the Lyapunov sense.
Theorem 8.
Let the conditions of Lemma 6 hold.
Then, the zero solution of the system (25) is W-stable in the Lyapunov sense if and only if for all solutions X ( t ) of IP (25), (2) with initial function Φ ( t ) PC , the product w ( t ) X ( t ) is bounded in J 0 .
Proof. Necessity: 
Let us assume that the zero solution is W-stable in the Lyapunov sense. Assume the contrary, that there exists Φ ¯ ( t ) PC such that for the existing, according to Theorem 3, unique solution X ¯ ( t ) of IP (25), (2) we have that w ( t ) X ¯ ( t ) is unbounded in J 0 . Since the zero solution is W-stable, then for any ϵ > 0 there exists a number δ ( 0 , ϵ ) such that for any Φ ( t ) PC with Φ < δ we have that the corresponding solution X ( t ) satisfies the inequality | w ( t ) X ( t ) | ϵ for any t J 0 . Consider the following initial function Φ * = δ Φ ¯ 2 Φ ¯ with norm Φ * = δ Φ ¯ 2 Φ ¯ = δ 2 < δ . Then, since the system (25) is a linear system the function X * ( t ) = δ 2 Φ ¯ X ¯ ( t ) is the corresponding solution of IP (25) (2) for the initial function Φ * ( t ) PC , and hence, for t J 0 we have that | w ( t ) X * ( t ) | = δ 2 Φ ¯ w ( t ) X ( t ) ϵ , which implies that | w ( t ) X ¯ ( t ) | 2 Φ ¯ ϵ δ for t J 0 , which contradicts our assumption.
Sufficiency: Let for any Φ ( t ) PC the product w ( t ) X ( t ) be bounded in J 0 , where X ( t ) is the corresponding solution of IP (25), (2) with this initial function, and hence, according to Lemma 4, we have that X ( t ) satisfies IP (3), (2) with F ( t , X t T ) 0 for t J and vice versa. Let t * J 0 be arbitrary and let b > a * > a be fixed numbers with t * [ a * , b ] . Then, for arbitrary initial functions Φ ¯ , Φ PC , their corresponding unique solutions of IP (3), (2) are X ¯ ( t ) , X ( t ) , respectively, and for any t * [ a * , b ] we have
| X ¯ ( t ) X ( t ) | Φ ¯ ( 0 ) Φ ( 0 ) ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] X ¯ ( t + θ ) X ( t + θ ) + n Γ ( α ) a t ( t τ ) 1 α h 0 [ d θ U ( τ , θ ) ] X ¯ ( τ + θ ) X ( τ + θ ) d τ Φ ¯ Φ ( t a ) α 1 + V sup s [ a , t ] | X ¯ ( s ) X ( s ) | + V Φ ¯ Φ + n Γ ( α ) a t ( t τ ) α 1 h 0 [ θ U ( τ , θ ) ] X ¯ ( τ + θ ) X ( τ + θ ) d θ d τ V sup s [ a , t ] | X ¯ ( s ) X ( s ) | + Φ ¯ Φ ( t a ) α 1 + V + n U b Γ ( 1 + α ) a t ( t τ ) α 1 sup η [ a , τ ] | X ¯ ( η ) X ( η ) | d τ + n U b ( t a ) α Γ ( 1 + α ) Φ ¯ Φ V sup s [ a , t ] | X ¯ ( s ) X ( s ) | + Φ ¯ Φ ( t a ) α 1 + n U b ( t a ) α Γ ( 1 + α ) + V + n U b Γ ( 1 + α ) a t ( t τ ) α 1 sup η [ a , τ ] | X ¯ ( η ) X ( η ) | d τ .
From (44), for any t [ a * , b ] it follows that
sup s [ a , t ] w ( t ) | X ¯ ( s ) X ( s ) | Φ ¯ Φ sup s [ a , t ] w ( s ) ( s a ) α 1 1 V + n U b b α ( 1 V ) Γ ( 1 + α ) + V ( 1 V ) + n U b Γ ( 1 + α ) a t ( t τ ) α 1 sup η [ a h , τ ] | X ¯ ( η ) X ( η ) | d τ .
Then, applying corollary 2 in [24] to (45), we obtain for any t [ a * , b ] the estimation
sup s [ a , t ] | X ¯ ( s ) X ( s ) | Φ ¯ Φ sup s [ a , t ] w ( s ) ( s a ) α 1 1 V + n U b ( b a ) α ( 1 V ) Γ ( 1 + α ) + V ( 1 V ) ( E α ( n U b t α ) ) .
Let us define for any t J 0 a family of functionals t : PC R n via the equality t ( Φ ) = w ( t ) X ( t ) , where X ( t ) denotes the corresponding solution of IP (25),(2) for initial function Φ ( t ) PC . From (46), it follows that t ( Φ ) is a continuous functional at any Φ ( t ) PC for any t J 0 i.e., for every Φ ( t ) PC we have that | t ( Φ ) | t Φ , where t denotes the norm of the functional. Then applying the Banach–Steinhaus theorem we obtain that the norms of all functionals are uniformly bounded, i.e., there exists a constant > 0 such that t for any t J . Let ϵ > 0 be arbitrary and choose δ = ϵ 2 > 0 . Then, for any Φ ( t ) PC with Φ < δ for the corresponding solution X ( t ) we have that | w ( t ) X ( t ) | = | t ( Φ ) | Φ = ϵ 2 < ϵ for any t J 0 , which implies that the zero solution is W-stable in the Lyapunov sense. □
Remark 5.
It is worth emphasizing that the result from the application of any kind of stability definitions concerning the differential equation (system) with all types of derivatives (integer or fractional order) essentially depends on the functional type of the set of all solutions of the studied object. As an example, for equations (systems) with first-order derivatives the solutions are usually either continuous differentiable or absolutely continuous functions. So, the conclusions could be true for all solutions which are absolutely continuous or only for these solutions which are continuous differentiable. For retarded or neutral equations (systems) the situation is more complicated. The space of the initial functions plays a leading role in its functional type, because the analytical type of the solutions essentially depends on them. This is very important for the stability in the Ulam–Hyers sense for the following reason: the functional type of the initial function and the solutions of system (1) must be the same for the corresponding inequality (42), since all solutions of (1) are solutions of inequality (42) too.
As a consequence, from the proved results, we obtain two necessary conditions for W-stability of the zero solution of the system (25).
Theorem 9.
Let the following conditions be fulfilled:
1.
The conditions of Lemma 6 hold.
2.
The zero solution of the system (25) is W-stable in the Lyapunov sense.
Then, the following relations hold:
Q = sup t J 0 | w ( t ) Q ( t ) | < , H = sup t J 0 | H ( t ) | < ,
where Q ( t ) = V a r s [ a h , a ] Q ( t , s ) and H ( t ) = sup s [ a , t ] | w ( t ) H ( t , s ) | .
Proof. 
Let us, as above for Φ ( t ) PC * and any t J 0 , define a family of functionals t : PC R n via the equality t ( Φ ) = w ( t ) X ( t ) , where X ( t ) denotes the corresponding solution of IP (25), (2) for this initial function. Then, according Theorem 6 and Lemma 6, for any t J 0 we have that
| t ( Φ ) | = | w ( t ) X ( t ) | = w ( t ) a h a [ d s Φ ˜ ( s a ) ] Q ( t , s ) sup s [ h , 0 ] Q ( t , s ) Φ Q ( t ) Φ V a r s [ h , 0 ] | Q ( t , s ) | Φ ,
and hence, t ( Φ ) is a continuous functional at any Φ ( t ) PC * for any t J 0 with norm t = V a r s [ h , 0 ] | Q ( t , s ) | . Then, by virtue of the Banach-Steinhaus theorem we have that there exists a constant > 0 such that t for any t J 0 . Thus, we have that the W-boundedness ( sup t J 0 ( V a r s [ h , 0 ] | Q ( t , s ) | ) < ) of the matrix w ( t ) Q ( t ) is a necessary condition for the W-stability of the zero solution of system (25).
By virtue of Theorem 8 we have that for any s J 0 and j n , the relation sup s [ a , t ] | w ( t ) H j ( t , s ) | < holds, and taking into account that | w ( t ) H j ( t , s ) | = 0 when t < s , we can conclude that H = sup t J 0 | H ( t ) | < . □
Corollary 3.
Let the following conditions be fulfilled:
1.
The conditions of Lemma 6 hold.
2.
The zero solution of the system (25) is W-GAS stable in the Lyapunov sense.
Then, the following relations
lim t | Q ( t ) | = 0 , lim t | H ( t ) | = 0
hold, where H ( t ) = sup s [ a , t ] | w ( t ) H ( t , s ) | .
Proof. 
The statement follows immediately from Theorem 9, since according to condition 2 of Corollary 3, for any j n we have lim t | Q j ( t , s ) | = 0 for s [ a h , a ] and lim t | H j ( t , s ) | = 0 for any s [ a , t ] . □
Remark 6.
Note that the differences in the proofs of Theorems 8 and 9 are caused from the differences in the spaces of the initial functions Φ ( t ) PC and Φ ( t ) PC * , respectively.
The next theorem makes clear the relation between the W-stability in the Lyapunov sense and Ulam–Hyers W-stability.
Theorem 10.
Let the following conditions be fulfilled:
1.
The conditions of Lemma 6 hold.
2.
The system (25) is Ulam–Hyers W-stable for any initial functions Φ ( t ) PC .
Then, the zero solution of the system (25) is stable in the Lyapunov sense if and only if for any Φ ( t ) PC and arbitrary ϵ > 0 IP (42), (2) has only W-bounded solutions.
Proof. Necessity: 
Let the zero solution of the system (25) be stable in the Lyapunov sense and assume that IP (42), (2) has a W-unbounded solution Y ( t ) | t J 0 C 1 α for some Φ Y PC and ϵ > 0 . Then, since (25) is Ulam–Hyers W-stable for any ϵ > 0 , there exists δ = δ ( ϵ ) > 0 and Φ * ( t ) C with | Φ * ( t ) Φ Y ( t ) | < δ for t [ a h , a ] such that for t J 0 inequality (43) holds for Y ( t ) | t J 0 and X * ( t ) , where X * ( t ) is the corresponding solution of IP (25), (3) with initial function Φ * ( t ) . On the other hand, since according to Theorem 8 the function X * ( t ) is a W-bounded solution of IP (25), (3), and since Y ( t ) | t J 0 is a W-unbounded solution of (42) according to our assumption, we conclude that inequality (43) is impossible to be fulfilled for Y ( t ) | t J 0 and X * ( t ) , which contradicts our assumption. Thus, all solutions of the inequality (43) are W-bounded.
Sufficiency: Let IP (42), (2) have only bounded solutions. Then, IP (25), (2) has, for any Φ PC , only W-bounded solutions. Then, by virtue of Theorem 8, the zero solution of system (25) is W-stable in the Lyapunov sense. □

6. On the Preservation of the Stability Properties

The aim of this section is, based on the obtained results in the previous sections, to make clear the preservation of the W-stability properties of the system (25) under nonlinear disturbances.
Definition 14.
The function X C 1 α will be called weighted Mittag–Leffler bounded (WML-bounded) of order ω R + if there exist t 0 J 0 , a constant C X > 0 , and a function a ( t ) P C ( J 0 , R + ) with a ( t ) = O ( t α ) such that | w ( t ) X ( t ) | a ( t ) E ω ( C X Γ ( ω ) t ω ) for t t 0 . By W M L ω , we denote the subset of all WML-bounded functions of order ω in C 1 α .
Definition 15.
We say that the vector-valued functional F : J × P C R n is a damper of order ω for system (1) if for any X ( t ) W M L ω for the function F X ( t ) = ( F 1 X ( t ) , , F n X ( t ) ) T = F ( t , X t T ) there exist a point t 0 J 0 and a constant C 0 = C 0 ( t 0 , X ) > 0 such that the following estimation | F X ( t ) | C 0 t ( ω + 1 ) holds for t t 0 .
The next theorem establishes an a priori estimate of the solutions of IP (1), (2) with initial functions Φ ( t ) PC .
Theorem 11.
Let the following conditions be fulfilled.
1.
The conditions of Lemma 6 hold.
2.
The conditions C hold.
3.
F ( t , 0 ) = 0 and sup t J | 𝓁 ( t ) | = L < .
Then, any solution X C 1 α of IP (3), (2) with initial function Φ ( t ) PC is WML-bounded of order α.
Proof. 
Let Φ ( t ) PC be an arbitrary initial function and X C 1 α be the corresponding unique solution of IP (3), (2). Then, for t J 0 from (3) we obtain that
| X ( t ) | Φ ( t a ) α 1 + h 0 [ d θ V ( t , θ ) ] X ( t + θ ) + n Γ ( α ) a t ( t τ ) 1 α h 0 [ d θ U ( τ , θ ) ] X ( t + θ ) + F ( τ , X τ T ) d τ Φ ( t a ) α 1 + V sup s [ a , t ] | X ( s ) | + V Φ + n Γ ( α ) a t ( t τ ) α 1 h 0 [ θ U ( τ , θ ) ] X ( t + θ ) d θ + F ( τ , X τ T ) d τ V sup s [ a , t ] | X ( s ) | + Φ ( ( t a ) α 1 + V ) + n ( U b + L ) Γ ( α ) a t ( t τ ) α 1 sup s [ a h , τ ] | X ( s ) | d τ V sup s [ a , t ] | X ( s ) | + Φ ( ( t a ) α 1 + V ) + n ( U b + L ) Γ ( α ) a t ( t τ ) α 1 sup s [ a , τ ] | X ( s ) | d τ ( t a ) α Φ V sup s [ a , t ] | X ( s ) | + Φ ( t a ) α 1 + n ( U b + L ) ( t a ) α Γ ( α ) + V + n ( U b + L ) Γ ( α ) a t ( t τ ) α 1 sup s [ a , τ ] | X ( s ) | d τ .
Without loss of generality, we can assume that sup s [ a , a + 1 ] X ( s ) α Φ = C * 1 when Φ > 0 . Then, from (47) we have
| X ( t ) | Φ ( t a ) α 1 ( 1 V ) + V Γ ( 1 + α ) + n ( U b + L ) ( t a ) α ( 1 V ) Γ ( α ) = n ( U b + L ) ( 1 V ) Γ ( α ) a a + 1 ( t τ ) α 1 sup s [ a , τ ] | X ( s ) | d τ + a + 1 t ( t τ ) α 1 sup s [ a , a + 1 ] | X ( s ) | d τ + a + 1 t ( t τ ) α 1 sup s [ a + 1 , τ ] | X ( s ) | d τ | Φ ( t a ) α 1 ( 1 V ) + V Γ ( 1 + α ) + n ( U b + L ) ( t a ) α ( 1 V ) Γ ( α ) + n ( U b + L ) ( 1 V ) Γ ( α ) α 1 sup s [ a , a + 1 ] X ( s ) ( t a ) α + a + 1 t ( t τ ) α 1 sup s [ a + 1 , τ ] | X ( s ) | d τ | Φ ( t a ) α 1 ( 1 V ) + n ( U b + L ) ( 1 + C * ) ( t a ) α + V Γ ( 1 + α ) ( 1 V ) Γ ( α ) + n ( U b + L ) ( 1 V ) Γ ( α ) a + 1 t ( t τ ) α 1 sup s [ a + 1 , τ ] | X ( s ) | d τ .
Since ( t a ) α 1 is monotone decreasing for t ( a , a + 1 ] , then sup s [ a + 1 , t ] | X * ( s ) | sup s [ a , t ] w ( s ) | X * ( s ) | , and hence,
sup s [ a + 1 , t ] | X ( s ) | Φ sup t a + 1 ( t a ) α 1 ( 1 V ) + n ( U b + L ) ( 1 + C * ) ( t a ) α + V Γ ( α ) ( 1 V ) Γ ( α ) + n ( U b + L ) ( 1 V ) Γ ( α ) a + 1 t ( t τ ) α 1 sup s [ a + 1 , τ ] | X ( s ) | d τ Φ 1 + V ( 1 V ) + n ( U b + L ) ( 1 + C * ) ( t a ) α ( 1 V ) Γ ( α ) + n ( U b + L ) ( 1 V ) Γ ( α ) a + 1 t ( t τ ) α 1 sup s [ a + 1 , τ ] | X ( s ) | d τ sup s [ a + 1 , t ] | X ( s ) | Φ 1 + V 1 V + n ( U b + L ) ( 1 + C * ) ( t a ) α ( 1 V ) Γ ( α ) E α ( n ( U b + L ) ( 1 V ) t α ) .
Applying to (48) Theorem 2, we obtain that
sup s [ a + 1 , t ] | X ( s ) | Φ 1 + V ( 1 V ) + n ( U b + L ) ( 1 + C * ) ( t a ) α ( 1 V ) Γ ( α ) E α ( n ( U b + L ) ( 1 V ) t α ) ,
which completes the proof. □
Remark 7.
We note that Theorem 11 establishes that for Lipschitz-type nonlinear disturbances, all solutions of IP (3), (2) with initial functions Φ ( t ) PC are WML-bounded of order α.
The proof of the next theorem is based on the following practical variant of the well-known Barbalat’s lemma.
Lemma 8
([30]). Let f L 1 l o c ( R + , R + ) be a bounded function and lim t 0 t f ( s ) d s < . Then, the relation lim t f ( t ) = 0 holds.
Theorem 12.
Let the following conditions be fulfilled:
1.
The conditions of Theorem 11 hold.
2.
The zero solution of system (25) is W-GAS.
3.
There exists ω > α such that the vector-valued functional F : J × P C R n is a damper of order ω.
Then, the zero solution of IP (1), (2) (or IP (3), (2)) with initial function Φ ( t ) PC * is W-GAS.
Proof. 
Let Φ * ( t ) PC * be arbitrary and X * ( t ) be the corresponding unique solution of IP (1), (2) (or (3), (2)). By virtue of Theorem 11 we have that X * ( t ) is WML-bounded of order α . We introduce the function F * ( t ) F ( t , X t * T ) , t J 0 , and using (28) and (34) define the functions X ˜ ( t ) and X ( t ) as follows:
X ˜ ( t ) = a h a [ d s Φ ˜ * ( s a ) ] Q ( t , s ) ,
X ( t ) = a t H ( t , s ) D a + 1 α F * ( s ) d s .
Then, from Theorem 6 it follows that X ˜ ( t ) is the unique solution of IP (25), (2) with initial function Φ * ( t ) , and by virtue of Theorem 7 the function X ( t ) is the unique solution of IP (24), (2) with F ( t ) F * ( t ) and initial function Z ( t a ) 0 , t [ a h , a ] . Thus, X ˜ ( t ) + X ( t ) is a solution of IP (2), (1) (IP (3), (2)), and hence, we have that
X * ( t ) = X ˜ ( t ) + X ( t ) = a h a [ d s Φ ˜ ( s a ) ] Q ( t , s ) + a t H ( t , s ) D a + 1 α F * ( s ) d s .
From (49) and Corollary 3 it follows that
lim t | X ˜ ( t ) | = lim t a h a [ d s Φ ˜ ( s a ) ] Q ( t , s ) lim t ( Φ ˜ ) | Q ( t ) | Φ lim t | Q ( t ) | = 0 .
Consideringthe function
G ( t ) = a t a τ H ( τ , s ) D a + 1 α F * ( s ) d s d τ ,
obviously the integrand X ( τ ) = a τ H ( τ , s ) D a + 1 α F * ( s ) d s is a continuous function at any τ J 0 , and hence, d d t G ( t ) = X ( t ) for t J 0 too. From (53) it follows that
| X ( τ ) | b τ H ( τ , s ) D a + 1 α F * ( s ) d s Γ 1 ( α ) | H ( t ) | a τ d d s a s ( s η ) α 1 F * ( η ) d η d s = Γ 1 ( α ) | H ( τ ) | a τ ( τ η ) α 1 F * ( η ) d η ( α Γ ( α ) 1 ) | H ( τ ) | a τ | F * ( η ) | d ( τ η ) α Γ 1 ( 1 + α ) | H ( τ ) | sup s [ a , τ ] ( | F * ( s ) | s α ) Γ 1 ( 1 + α ) | H ( τ ) | sup s [ a , b ] ( | F * ( s ) | s α ) + sup s [ b , τ ] ( | F * ( s ) | s α ) Γ 1 ( 1 + α ) | H ( τ ) | sup s [ a , b ] ( | F * ( s ) | s α ) + sup s [ b , τ ] ( C 0 s ( ω + 1 ) s α ) .
Since ω > α and the zero solution of the system (25) is W-GAS, then by virtue of Theorem 9 we have that H = sup t J 0 | H ( t ) | < , and hence, | g ( τ ) | is bounded on J 0 . Then, from (53) and (54) it follows that
Γ 1 ( 1 + α ) | H ( τ ) | sup s [ a , b ] ( | F * ( s ) | s α ) + sup s [ b , τ ] ( C 0 s ( ω + 1 ) s α )
| G ( t ) | = a t b τ H ( τ , s ) D a + 1 α F * ( s ) d s d τ Γ 1 ( 1 + α ) a t | H ( τ ) | sup s [ a , τ ] ( | s α F * ( s ) | ) d τ
Γ 1 ( 1 + α ) b α sup s [ a , b ] | H ( τ ) | | F * ( s ) | ( b a ) + Γ 1 ( 1 + α ) b t | H ( τ ) | sup s [ b , τ ] ( C 0 s ( ω + 1 ) s α ) d τ Γ 1 ( 1 + α ) b α sup s [ a , b ] | H ( τ ) | | F * ( s ) | ( b a ) + Γ 1 ( 1 + α ) H C 0 t ( ω + 1 ) + α + 1
and, hence, lim t | G ( t ) | < . Then, by virtue of Lemma 8 we have that lim t | X ( τ ) | = 0 , and thus, from (51) and (52) it follows that lim t | X * ( t ) | = 0 , which completes the proof. □

7. An Example

Example 1.
Let n = 1 , the lower terminal a = 0 , α ( 0 , 1 ) , h > 0 , p , b , c R , U ( t , θ ) = a H ( t ) + b H ( t h ) , V ( t , θ ) = p H ( t h ) , and ϕ P C ( [ h , 0 ] ) . Then, IP (25), (2) obtains the following form:
D a + α ( x ( t ) + p x ( t h ) ) = c x ( t ) + b x ( t h ) , t ( 0 , ) ,
x ( t ) = ϕ ( t ) , t [ h , 0 ] , D a + α 1 x ( 0 + ) = ϕ ( 0 ) .
According to the results obtained by us in Section 3, IP (55), (56) has a unique solution, which we will find an explicit representation of in this example.
To obtain an explicit representation of the solution we will use the approach based on the Laplace transform (LT) introduced in [31] for the case of Caputo-type derivatives. Below, we assume that | p | < 1 , and then, by virtue of theorem 3 in [17] the LT can be applied correctly to Equation (55). Denote the LT of x ( t ) by
x ^ ( s ) = L ( x ( t ) ; s ) = 0 e s t x ( t ) d t , ϕ ^ ( s ) = h 0 e s t ϕ ( t ) d t .
Applying the LT to both sides of (55) and taking into account (56) we have
s α x ^ ( s ) + s α e s h p ( x ^ ( s ) + ϕ ^ ( s ) ) ( ϕ ( 0 ) + p ϕ ( h ) ) = c x ^ ( s ) + b e s h ( x ^ ( s ) + ϕ ^ ( s ) ) ,
and hence,
x ^ ( s ) 1 b p s α s α c e s h = ϕ ( 0 ) + p ϕ ( h ) s α c + b p s α s α c e s h ϕ ^ ( s ) .
For any s C with | s | α > | c | , we have the series expansion
1 b p s α s α c e s h 1 = k = 0 ( b p s α ) k ( s α c ) k e k s h
and then, from (57) it follows that
x ^ ( s ) = ( ϕ ( 0 ) + p ϕ ( h ) ) k = 0 ( b p s α ) k ( s α c ) k + 1 e k s h + ϕ ^ ( s ) k = 0 ( b p s α ) k ( s α c ) k e k s h .
Substituting into (58) the expression ( b p s α ) k = i = 0 k k i ( 1 ) i p i s i α b k i we obtain that the LT of x ^ ( s ) has the following form:
x ^ ( s ) = ( ϕ ( 0 ) + p ϕ ( h ) ) k = 0 i = 0 k k i ( p ) i b k i s i α ( s α c ) k + 1 + ϕ ^ ( s ) k = 0 i = 0 k k i ( p ) i b k i s i α ( s α c ) k .
For simplicity we will assume that ϕ ( t ) ϕ 0 R , t [ h , 0 ] , and then, we have that
ϕ ^ ( s ) = ϕ 0 s 1 ( e s h 1 ) .
Consider the three-parameter Mittag–Leffler (Prabhakar) function
E α , β γ ( z ) = Γ 1 ( γ ) j = 0 Γ ( j + γ ) z j j ! Γ ( α j + β )
introduced in [32], where β , γ R . For z = c t α , t 0 , R e s > 0 with | s | α > | c | the LT of the function e α , β γ ( t ; c ) = t β 1 E α , β γ ( c t α ) is L ( e α , β γ ( t ; c ) ; s ) = s α γ β ( s α c ) γ and when h 0 (see theorem 1.31 in [33]) we have that the following relation holds:
L 1 s α γ β ( s α c ) γ e s h , t = e α , β γ ( t h ; c ) , t h , 0 , t < h .
Denoting as usual by [ t ] the greatest integer number less or equal to t for any t R , from (59), (60), and (61), it follows that the solution x ( t ) of IP (55), (56) possess the representation
x ( t ) = ϕ 0 ( 1 + c ) k = 0 t h i = 0 k k i ( p ) i b k i e α , α ( k i ) + 1 k + 1 ( t h k ; c ) ϕ 0 k = 0 t h i = 0 k k i ( p ) i b k i e α , α ( k i ) + 1 k ( t h k ; c ) + ϕ 0 k = 0 t h + 1 i = 0 k k i ( p ) i b k i e α , α ( k i ) + 1 k ( t h k + h ; c ) .
Remark 8.
Note that the approach used by the presentation of the solution x ( t ) of IP (55), (56) can be extended to obtain a presentation of the solutions with step functions (with one jump of the first kind in the interval [ h , 0 ] ) as initial functions, and hence, as a consequence, a presentation of the generalized fundamental matrix Q ( t , s ) can be obtained. But this idea, as well as a numerical simulation via the methods developed in [34,35,36,37], is beyond the area of this article and can be a theme for future research.

8. Conclusions and Comments

In this paper, we have studied a general class of nonlinear disturbed neutral linear fractional systems with derivatives in the Riemann–Liouville sense and distributed delays. As motivation for this study we can refer to the meaningful physical interpretations of models with these Riemann–Liouville fractional derivatives presented in [7]. Mainly, we have studied the most important, from the point of view of applications, and the technically more complicated case, when the lower terminal of the Riemann–Liouville derivatives coincides with the end point of the initial interval. First, it was proved that the initial problem for these systems with discontinuous initial functions possesses a unique solution under some natural assumptions. It is worth mentioning that the assumptions used to derive this result are similar to those used in the case of systems with first-order derivatives, and as far as we know this is the first result on this theme. Then, as a consequence of the obtained result, we have proved the existence and uniqueness of a fundamental matrix and a generalized fundamental matrix for the studied neutral linear homogeneous system. The existence of the fundamental matrices have allowed us to establish an integral representation for the solutions of the initial problem for the homogeneous system and also for the corresponding inhomogeneous system. Furthermore, for fractional systems with Riemann–Liouville derivatives we have introduced a new concept for weighted stabilities in the Lyapunov, Ulam–Hyers and Ulam–Hyers–Rassias senses. Note that the introduced concept coincides with the classical stability concept for the cases of integer-order or Caputo-type derivatives. To prove the applicability of the introduced concept it was proved that the zero solution of the homogeneous system is weighted stable if and only if all its solutions are weighted bounded; this result is well known for systems with first-order derivatives and is also established for the case of systems with Caputo-type derivatives in our former works. In addition, for the homogeneous system it was established that the weighted stability in the Lyapunov sense and weighted stability in the Ulam–Hyers sense are equivalent if and only if the inequality appearing in the Ulam–Hyers definition possesses only bounded solutions. Finally, we have derived natural sufficient conditions under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under appropriate nonlinear disturbances.
In our point of view the main contributions in the article can be highlighted as follows:
  • A new concept was introduced for weighted stabilities in the Lyapunov, Ulam–Hyers and Ulam–Hyers–Rassias senses; it coincides with the classical stability concept for the cases of integer-order or Caputo-type fractional derivatives.
  • Sufficient conditions have been obtained which guarantee that the weighted stability in the Lyapunov sense and the weighted stability in the Ulam–Hyers sense are equivalent.
  • Sufficient conditions have been given under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under appropriate nonlinear disturbances.
As a future perspective for research we think that it will be important from the point of view of applications to establish explicit-type sufficient conditions which guarantee weighted stability in the Lyapunov or Ulam–Hyers senses.

Author Contributions

Conceptualization, H.K., M.M., S.I.C. and A.Z.; methodology, H.K., M.M., S.I.C. and A.Z.; validation, H.K., M.M., S.I.C. and A.Z.; formal analysis, H.K., M.M., S.I.C. and A.Z.; investigation, H.K., M.M., S.I.C. and A.Z.; resources, H.K., M.M., S.I.C. and A.Z.; data curation, H.K., M.M., S.I.C. and A.Z.; writing—original draft preparation, H.K., M.M., S.I.C. and A.Z.; writing—review and editing; supervision and project administration, H.K., M.M., S.I.C. and A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The authors of this research have been partially supported as follows: Hristo Kiskinov by the Bulgarian National Science Fund, Grant KP-06-N52/9, Mariyan Milev by EU-NextGeneration, through the NRRPlan of Bulgaria, SUMMIT BG-RRP-2.004-0008-C01 and Andrey Zahariev by the Bulgarian National Science Fund under Grant KP-06-N52/4, 2021.

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. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  2. Saxena, R.K.; Kalla, S.L. On the solutions of certain fractional kinetic equations. Appl. Math. Comput. 2008, 199, 504–511. [Google Scholar] [CrossRef]
  3. Li, C.; Qian, D.; Chen, Y.Q. On Riemann-Liouville and Caputo derivatives. Discret. Dyn. Nat. Soc. 2011, 2011, 562494. [Google Scholar] [CrossRef]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equation; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  6. Podlubny, I. Geometric and physical interpretation of fractional integration and fractional differentiation. Fract. Calc. Appl. Anal. 2002, 5, 367–386. [Google Scholar]
  7. Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol. Acta 2006, 46, 765–771. [Google Scholar]
  8. Kiskinov, H.; Zahariev, A. On fractional systems with Riemann-Liouville derivatives and distributed delays-Choice of initial conditions, existence and uniqueness of the solutions—Choice of initial conditions, existence and uniqueness of the solutions. Eur. Phys. J. Spec. Top. 2017, 9, 3473–3487. [Google Scholar] [CrossRef]
  9. Agarwal, R.; Hristova, S.; O’Regan, D. Stability concepts of Riemann-Liouville fractional-order delay nonlinear systems. Mathematics 2021, 9, 435. [Google Scholar] [CrossRef]
  10. Agarwal, R.; Hristova, S.; O’Regan, D. Explicit solutions of initial value problems for linear scalar Riemann-Liouville fractional differential equations with a constant delay. Mathematics 2020, 8, 32. [Google Scholar] [CrossRef]
  11. Kiskinov, H.; Madamlieva, E.; Zahariev, A. Hyers–Ulam and Hyers–Ulam–Rassias stability for linear fractional systems with Riemann–Liouville derivatives and distributed delays. Axioms 2023, 12, 637. [Google Scholar] [CrossRef]
  12. Li, C.P.; Zhang, F.R. A survey on the stability of fractional differential equations. Eur. Phys. J. Spec. Top. 2011, 193, 27–47. [Google Scholar] [CrossRef]
  13. Liu, S.; Wu, X.; Zhou, X.F.; Jiang, W. Asymptotical stability of Riemann-Liouville fractional nonlinear systems. Nonlinear Dyn. 2016, 86, 65–71. [Google Scholar] [CrossRef]
  14. Li, M.; Wang, J. Representation of solution of a Riemann–Liouville fractional differential equation with pure delay. Appl. Math. Lett. 2018, 85, 118–124. [Google Scholar] [CrossRef]
  15. Liang, C.; Wang, J.; O’Regan, D. Representation of a solution for a fractional linear system with pure delay. Appl. Math. Lett. 2018, 77, 72–78. [Google Scholar] [CrossRef]
  16. 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]
  17. Veselinova, M.; Kiskinov, H.; Zahariev, A. Explicit conditions for stability of neutral linear fractional system with distributed delays. AIP Conf. Proc. 2016, 1789, 040005. [Google Scholar] [CrossRef]
  18. Chartbupapan, W.; Bagdasar, O.; Mukdasai, K. A novel delay-dependent asymptotic stability conditions for differential and Riemann-Liouville fractional differential neutral systems with constant delays and nonlinear perturbation. Mathematics 2020, 8, 82. [Google Scholar] [CrossRef]
  19. Myshkis, A. Linear Differential Equations with Retarded Argument; Nauka: Moscow, Russia, 1972. (In Russian) [Google Scholar]
  20. Hale, J.; Lunel, S. Introduction to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  21. Kolmanovskii, V.; Myshkis, A. Introduction to the Theory and Applications of Functional Differential Equations; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999. [Google Scholar]
  22. Kiskinov, H.; Madamlieva, E.; Veselinova, M.; Zahariev, A. Integral representation of the solutions for neutral linear fractional system with distributed delays. Fractal Fract. 2021, 5, 222. [Google Scholar] [CrossRef]
  23. Krasnosel’skii, M.A. Some problems of nonlinear analysis. Am. Math. Soc. Transl. Ser. 1958, 2, 10. [Google Scholar]
  24. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equations. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef]
  25. Madamlieva, E.; Milev, M.; Stoyanova, T. On stability criteria induced by the resolvent kernel for a fractional neutral linear system with distributed delays. Mathematics 2023, 11, 626. [Google Scholar] [CrossRef]
  26. Kiskinov, H.; Milev, M.; Veselinova, M.; Zahariev, A. Continuous dependence on the initial functions and stability properties in Hyers–Ulam–Rassias sense for neutral fractional systems with distributed delays. Fractal Fract. 2023, 7, 742. [Google Scholar] [CrossRef]
  27. Natanson, I.P. Theory of Functions of a Real Variable, 5th ed.; Frederick Ungar Publishing Co.: New York, NY, USA, 1983. [Google Scholar]
  28. Reitano, R.R. Foundations of Quantitative Finance: 5. General Measure and Integration Theory; Brandeis International Business School: Waltham, MA, USA, 2018. [Google Scholar]
  29. Madamlieva, E.; Kiskinov, H.; Petkova, M.; Zahariev, A. On the preservation with respect to nonlinear perturbations of the stability property for nonautonomous linear neutral fractional systems with distributed delays. Mathematics 2022, 10, 2642. [Google Scholar] [CrossRef]
  30. Sun, Z. A gathering of Barbalat’s lemmas and their (unsung) cousins. arXiv 2023. [Google Scholar] [CrossRef]
  31. Tuan, H.T.; Thai, H.D.; Garrappa, R. An analysis of solutions to fractional neutral differential equations with delay. Commun. Nonlinear Sci. Numer. Simul. 2021, 100, 105854. [Google Scholar] [CrossRef]
  32. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  33. Schiff, J.L. The Laplace Transform: Theory and Applications; Springer: New York, NY, USA, 1999. [Google Scholar]
  34. Garrappa, R.; Popolizio, M. Computing the matrix Mittag-Leffler function with applications to fractional calculus. J. Sci. Comput. 2018, 77, 129–153. [Google Scholar] [CrossRef]
  35. Micula, S. An iterative numerical method for fractional integral equations of the second kind. J. Comput. Appl. Math. 2018, 339, 124–133. [Google Scholar] [CrossRef]
  36. Jafari, H.; Tuan, N.A.; Ganji, R.M. A new numerical scheme for solving pantograph type nonlinear fractional integro-differential equations. J. King Saud Univ.-Sci. 2021, 33, 101185. [Google Scholar] [CrossRef]
  37. Moumen, A.; Mennouni, A.; Bouye, M. A novel Vieta–Fibonacci projection method for solving a system of fractional integrodifferential equations. Mathematics 2023, 11, 3985. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kiskinov, H.; Milev, M.; Cholakov, S.I.; Zahariev, A. Fundamental Matrix, Integral Representation and Stability Analysis of the Solutions of Neutral Fractional Systems with Derivatives in the Riemann—Liouville Sense. Fractal Fract. 2024, 8, 195. https://doi.org/10.3390/fractalfract8040195

AMA Style

Kiskinov H, Milev M, Cholakov SI, Zahariev A. Fundamental Matrix, Integral Representation and Stability Analysis of the Solutions of Neutral Fractional Systems with Derivatives in the Riemann—Liouville Sense. Fractal and Fractional. 2024; 8(4):195. https://doi.org/10.3390/fractalfract8040195

Chicago/Turabian Style

Kiskinov, Hristo, Mariyan Milev, Slav Ivanov Cholakov, and Andrey Zahariev. 2024. "Fundamental Matrix, Integral Representation and Stability Analysis of the Solutions of Neutral Fractional Systems with Derivatives in the Riemann—Liouville Sense" Fractal and Fractional 8, no. 4: 195. https://doi.org/10.3390/fractalfract8040195

APA Style

Kiskinov, H., Milev, M., Cholakov, S. I., & Zahariev, A. (2024). Fundamental Matrix, Integral Representation and Stability Analysis of the Solutions of Neutral Fractional Systems with Derivatives in the Riemann—Liouville Sense. Fractal and Fractional, 8(4), 195. https://doi.org/10.3390/fractalfract8040195

Article Metrics

Back to TopTop