Next Article in Journal
Computational Model of Noradrenaline Modulation of Astrocyte Responses to Synaptic Activity
Next Article in Special Issue
High-Order Approximation to Generalized Caputo Derivatives and Generalized Fractional Advection–Diffusion Equations
Previous Article in Journal
Delicate Comparison of the Central and Non-Central Lyapunov Ratios with Applications to the Berry–Esseen Inequality for Compound Poisson Distributions
Previous Article in Special Issue
On the Preservation with Respect to Nonlinear Perturbations of the Stability Property for Nonautonomous Linear Neutral Fractional Systems with Distributed Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Stability Criteria Induced by the Resolvent Kernel for a Fractional Neutral Linear System with Distributed Delays

by
Ekaterina Madamlieva
1,*,
Marian Milev
2 and
Tsvetana Stoyanova
3
1
Department of Mathematical Analysis and Differential Equations, Faculty of Applied Mathematics and Informatics, Technical University of Sofia, 1756 Sofia, Bulgaria
2
Department of Mathematics and Physics, University of Food Technology, 4000 Plovdiv, Bulgaria
3
Department of Management and Administration, University of National and World Economy, 1700 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(3), 626; https://doi.org/10.3390/math11030626
Submission received: 28 December 2022 / Revised: 19 January 2023 / Accepted: 20 January 2023 / Published: 26 January 2023
(This article belongs to the Special Issue Stability Analysis of Fractional Systems-II)

Abstract

:
We consider an initial problem (IP) for a linear neutral system with distributed delays and derivatives in Caputo’s sense of incommensurate order, with different kinds of initial functions. In the case when the initial functions are with bounded variation, it is proven that this IP has a unique solution. The Krasnoselskii’s fixed point theorem, a very appropriate tool, is used to prove the existence of solutions in the case of the neutral systems. As a corollary of this result, we obtain the existence and uniqueness of a fundamental matrix for the homogeneous system. In the general case, without additional assumptions of boundedness type, it is established that the existence and uniqueness of a fundamental matrix lead existence and uniqueness of a resolvent kernel and vice versa. Furthermore, an explicit formula describing the relationship between the fundamental matrix and the resolvent kernel is proven in the general case too. On the base of the existence and uniqueness of a resolvent kernel, necessary and sufficient conditions for the stability of the zero solution of the homogeneous system are established. Finally, it is considered a well-known economics model to describe the dynamics of the wealth of nations and comment on the possibilities of application of the obtained results for the considered systems, which include as partial case the considered model.

1. Introduction

Fractional calculus is not a new mathematical tool, but in view of its application for modeling many real-world phenomena, it has attracted considerable attention in recent decades. Comprehensive information about the fractional calculus theory, fractional differential equations and its applications can be found in the monographs [1,2,3]. For a practical oriented exposition of this theme, we refer to [4,5].
As is mentioned in [6], a predictable process can be physically realized only if it is stable in some suitable natural sense. It must be noted that the ascertainment of this fact is obtained in general from practical experience. This fact explains why the investigation of the property stability of the models of processes is so important. The same is true concerning the existence of a fundamental matrix and the existence of an integral representations of their solutions, since they are the main tools for investigation of the stability properties. It is well-known, generally speaking, that for the delayed fractional differential equations, stability properties are more complicated to be studied in comparison with fractional differential equations and systems without delay.
An excellent historical overview of the works related to this theme up until 2011 can be obtained from the survey [7] and the references therein. Concerning the recent works related to this theme devoted to fractional equations and systems without delay we refer to [8,9]. The autonomous case for retarded (delayed) and neutral fractional equations and systems on this theme are considered in [10] for single constant delay, while [11,12] discusses distributed delay. Regarding distributed delay, systems have also been studied using Riemann–Liouville type derivatives. The singular case is treated in [13] and the case of distributed order fractional derivatives is considered in [14]. The integral representation for the neutral case with distributed delay is studied in [15] and some explicit conditions for stability of the same type systems in terms of logarithmic norm (Losinskii measure) are given in [16]. From the works devoted to the integral representation of the solutions for retarded fractional equations and systems in the nonautonomous case, we refer to [17], which treated a system with single delay, and for the case of distributed delay, we refer to [18]. The same problem in the neutral case is considered in [19] (single variable delay), while in [20] discusses the case of distributed delay. The important problem of the existence and uniqueness of a fundamental matrix and its smoothness is studied for the retarded systems in [21] and for neutral systems in [22]. Concerning the investigation of the different kinds of stability of the zero solution for linear retarded fractional systems, we note the remarkable work of [23], devoted to the case of multiple concentrated delays, [24], which studies the asymptotic stability of a system with distributed delays, and [25], devoted to the finite time stability. Different kinds of stability criteria for a linear system with multiple delays are also considered in [26], while [27] considers the neutral case with distributed delays and with the Riemann–Liouville type derivatives. The stability of the zero solution of nonlinear fractional systems were studied under different approaches. In the work [28], the asymptotic stability is studied via linearization, and in [29,30], for retarded and neutral systems with distributed delays, respectively, the preservation of the asymptotic stability of linear systems is studied under nonlinear perturbation. For numerical aspects, we refer to [31,32].
The present work is inspired by [33] and is mainly devoted of the solving of an open problem posed in this work, which investigated a neutral system with fractional derivatives in Caputo sense, of incommensurate orders belonging to the interval ( 0 , 1 ) and distributed delays. Our approach is based on a direct construction of a fundamental matrix under as minimal as possible restrictions, which leads to the solving of an auxiliary matrix IP for the studied system with discontinuous initial matrix-valued piecewise continuous or bounded variation on every compact interval in R initial functions. As an application of the obtained results, we obtain some analytical properties of the fundamental matrix. Note that the conditions and the obtained results are similar to those in the case of delayed systems with integer derivatives.
The paper is organized as follows. In Section 2, we recall the definitions of Riemann–Liouville and Caputo fractional derivatives; some with additional theorems and notations. Section 3 is devoted to the problem of the existence and uniqueness of the fundamental matrix of the homogeneous system. Section 4 is devoted to an application of the established results. We prove the existence and uniqueness of the corresponding resolvent kernel, without assumptions of uniform boundness in t. In addition, we also study the relationship between the resolvent kernel and the fundamental matrix and establish some analytical properties. In Section 5, a well-known economics model is considered, describing the dynamic of the wealth of nations, and some comments are made on the possibilities of the application of the obtained results for the considered systems, which include as partial case the considered model. In this section, some conclusions are also given.

2. Preliminaries and Problem Statement

To avoid misunderstandings, below are given the definitions of the Riemann–Liouville and Caputo fractional derivatives as well as their properties. Detailed information can be found in the monographs [1,2].
For a R and α ( 0 , 1 ) arbitrarily, for each t > a , the left-sided fractional integral operator and the left side Riemann–Liouville of order α are defined via ( D a + α f ) ( t ) = Γ 1 ( α ) a t ( t s ) α 1 f ( s ) d s and R L D a + α f ( t ) = 1 Γ ( 1 α ) d d t D a + α 1 f ( t ) , respectively.
The Caputo fractional derivative of same order by C D a + α f ( t ) = R L D a + α [ f ( s ) f ( a ) ] ( t ) , where f L 1 l o c ( R , R ) , L 1 l o c ( R , R ) denote the real linear space of all locally Lebesgue integrable functions f : R R and B L 1 l o c ( R , R ) L 1 l o c ( R , R ) for the subspace of all locally bounded functions.
Below the following notations will be used: R + = ( 0 , ) , R ¯ + = [ 0 , ) , N 0 = { 0 } , k n = { 1 , 2 , . . . , n } , n N , i m 0 { 0 } , J a = [ a , ) , a R , J a + b = [ a , a + b ] , b R + , I , Θ R n × n , which are the identity and zero matrix, respectively, and 0 R n is the zero vector-column.
For W : J a × R R n × n , W ( t , θ ) = { w k j ( t , θ ) } k , j = 1 n , we denote | W ( t , θ ) | = k , j = 1 n | w k j ( t , θ ) | , ( t , θ ) J a × R , W ( t , θ ) B V l o c ( J a × R , R n × n ) , which denote the functions W ( t , θ ) , which, for every fixed t J a , has bounded variation in θ on every compact interval K R and denotes Var θ K W ( t , θ ) = { Var θ K w k j ( t , θ ) } k , j = 1 n and Var θ K W ( t , · ) = k , j = 1 n Var θ K w k j ( t , · ) .
Y ( t ) = ( y 1 ( t ) , . . . , y n ( t ) ) T : J a R n , Y ( t ) B V l o c ( J a , R n ) denote the functions with bounded variation in every compact interval K R and
I β ( Y ( t ) ) = diag ( ( y 1 ( t ) ) β 1 , , ( y n ( t ) ) β n ) , where β = ( β 1 , . . . , β n ) , β k [ 1 , 1 ] , k n .
It is well known that for z R + , the gamma function has a minimum at z m i n + 1.46 , where it attains the value Γ z m i n + 0.8856 (truncated).
With BL, we denote the Banach space of all vector functions Φ = ϕ 1 , . . . , ϕ n T : [ h , 0 ] R n , which are bounded and Lebesgue measurable on the interval [ h , 0 ] with norm | | Φ | | = k n sup s [ h , 0 ] | ϕ k ( s ) | < and S Φ denotes the set of all jumps point of arbitrary Φ BL .
PC = P C ( [ h , 0 ] , R n ) denote the subspace of all right continuous for t S Φ piecewise continuous functions, P C * = PC B V ( [ h , 0 ] , R n ) , C = C ( [ h , 0 ] , R n ) , and all of them are endowed with the same sup norm.
For t > a consider the inhomogeneous neutral linear delayed system with and distributed delays in the form:
D a + α ( X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) + F ( t ) ,
and the corresponding homogeneous neutral linear system:
D a + α ( X ( t ) h 0 [ d θ V ( t , θ ) ] X ( t + θ ) ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) ,
where the differential orders, α = ( α 1 , . . . , α n ) , α k ( 0 , 1 ) , can be incommensurate, U i , V l : J a × R R n × n , V ( t , θ ) = { v k j ( t , θ ) } k , j = 1 n = l r V l ( t , θ ) , r N , U ( t , θ ) = { u k j ( t , θ ) } k , j = 1 n = i m 0 U i ( t , θ ) , F : J a R n , X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T , X t ( θ ) = ( x t 1 ( θ ) , . . . , x t n ( θ ) ) T , X t ( θ ) = X ( t + θ ) , h θ 0 , h R + , F B L 1 l o c ( J a , R ) , D a + α X ( t ) = ( D a + α 1 x 1 ( t ) , . . . , D a + α n x n ( t ) ) T , for every t J a , where D a + α k = C D a + α k denotes the left side Caputo fractional derivative. For clarity, we rewrite the system (2) in more detail form:
D a + α k ( x k ( t ) l = 1 r ( j = 1 n τ l 0 x j ( t + θ ) d θ v k j l ( t , θ ) ) ) = i = 0 m ( j = 1 n σ i 0 x j ( t + θ ) d θ u k j i ( t , θ ) ) ,
where σ i , τ l R + , τ = max l r τ l , σ = max i m 0 σ i , h = max ( τ , σ ) and for arbitrary Φ  BL, we introduce the following initial condition:
X t ( θ ) = X ( t + θ ) = Φ ( t a + θ ) for t + θ a , θ [ h , 0 ] .
In our consideration below, we need the following auxiliary system:
X ( t ) = C Φ ( 0 ) + h 0 [ d θ V ( t , θ ) ] X ( t + θ ) + I 1 ( Γ ( α ) ) a t I α 1 ( t τ ) ( h 0 [ d θ U ( τ , θ ) ] X ( τ + θ ) + F ( τ ) ) d τ ,
where C Φ ( 0 ) = Φ ( 0 ) h 0 [ d θ V l ( a , θ ) ] Φ ( θ ) .
Definition 1. 
We say that for the kernels U i , V l : J a × R R n × n , the conditions(S)are fulfilled if for every l r , i m 0 , the following conditions hold (see [6,33]):
(S1) 
The functions ( t , θ ) U i ( t , θ ) , ( t , θ ) V l ( t , θ ) are measurable in ( t , θ ) J a × R and normalized so that U ( t , θ ) = V ( t , θ ) = 0 for θ 0 , U i ( t , θ ) = U i ( t , σ i ) for θ σ i , V l ( t , θ ) = V l ( t , τ l ) for θ τ l and every t J a . The kernels U i ( t , θ ) and V l ( t , θ ) are continuous from left in θ on ( σ i , 0 ) and ( τ l , 0 ) , respectively, for t J a , and U i ( t , θ ) , V l ( t , θ ) B V l o c ( J a × R , R n × n ) for every fixed t J a .
(S2) 
For t J a , the functions
V * ( t ) = Var [ h , 0 ] V ( t , · ) , U * ( t ) = Var [ h , 0 ] U ( t , · ) B L 1 l o c ( J a , R n × n )
and the kernel V ( t , θ ) is uniformly nonatomic at zero (i.e. for every ε > 0 , there exists δ ( ε ) > 0 such that for each t J a , we have Var [ δ , 0 ] V ( t , · ) < ε ([34]).
(S3) 
For t J a and θ [ h , 0 ] , the Lebesgue decompositions of the kernels have the form: U i ( t , θ ) = U d i ( t , θ ) + U c i ( t , θ ) , V l ( t , θ ) = V d l ( t , θ ) + V c l ( t , θ ) , U c i ( t , θ ) = U a c i ( t , θ ) + U s i ( t , θ ) , V c l ( t , θ ) = V c l ( t , θ ) + V s l ( t , θ ) , where the indexes d , a c , s denoted the jump, the absolutely continuous and the singular part, respectively, in the Lebesgue decompositions and U c i ( t , θ ) , V c l ( t , θ ) are the continuous part of these decompositions. In addition, U d i ( t , θ ) = { a k j i ( t ) H ( θ + σ i ( t ) ) } k , j = 1 n , V d l ( t , θ ) = { a ¯ k j l ( t ) H ( θ + τ l ( t ) ) } k , j = 1 n , where H ( t ) is the Heaviside function:
A i ( t ) = { a k j i ( t ) } k , j = 1 n B L 1 l o c ( J a , R n × n ) , A ¯ l ( t ) = { a ¯ k j l ( t ) } k , j = 1 n C ( J a , R n × n ) ,
σ 0 ( t ) 0 , σ i ( t ) C ( J a , [ 0 , σ i ] ) , τ l ( t ) C ( J a , [ 0 , τ l ] ) for every t J a .
(S4) 
For each t * J a , the relations:
lim t t * σ 0 | U i ( t , θ ) U i ( t * , θ ) | d θ = lim t t * τ 0 | V l ( t , θ ) V l ( t * , θ ) | d θ = 0
hold, and there exists γ ( 0 , ) such that the kernels V a c l ( t , θ ) and V s l ( t , θ ) are continuous in t when t [ a , a + γ ] , θ [ h , 0 ] and l r .
(S5) 
The sets S Φ l = { t [ a , a + h ] | t τ l ( t ) a S Φ } , S Φ i = { t [ a , a + h ] | t σ i ( t ) a S Φ } , for every Φ PC, do not have limit points.
Definition 2. 
The vector function X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T is a solution of the initial problem (IP) (1), (3) in J a + b ( J a ) , b R + , if X | [ a , a + b ] C ( [ a , a + b ] , R n ) , ( X | J a C ( J a , R n ) ) , X ( t ) satisfies the system (1) for all t ( a , b ] ( t ( a , ) ) and the initial condition (3) for t + θ a , θ [ h , 0 ] .
Definition 3. 
The vector function X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T is a solution of the initial problem (IP) (4), where (3) in J a + b ( J a ) if X | [ a , a + b ] C ( [ a , a + b ] , R n ) , ( X | J a C ( J a , R n ) ) , X ( t ) satisfies the system (4) for all t ( a , b ] ( t ( a , ) ) and the initial condition (3) for t + θ a , θ [ h , 0 ] .
In this paper, we assume that the condition (S) holds. For arbitrary fixed s J a , we define for ( t , s ) J a × R two initial matrix functions:
Φ C ( t , s ) = I , t = s Θ , t < s , ( t , s ) J a × R Θ , s < a , Φ Q ( t , s ) = I , h s t a Θ , t < s Θ , s < h
and introduce the auxiliary matrix initial problem:
D a + α ( W ( t , s ) h 0 [ d θ V ( t , θ ) ] W ( t + θ , s ) ) = h 0 [ d θ U ( t , θ ) ] W ( t + θ , s ) , t ( s , ) ;
W ( t , s ) = Φ C ( t , s ) , t ( , s ] ;
W ( t , s ) = Φ Q ( t , s ) , t ( , s ] .
Definition 4. 
For each s J a , the matrix-valued function C ( · , s ) : R R n × n , C ( t , s ) = { c k j ( t , s ) } k , j n is called a solution of the IP (5), (6) for t ( s , ) , if C ( · , s ) is continuous in t on J s , satisfying the matrix Equation (5) for t ( s , ) , as well as the initial condition (6).
Definition 5. 
For each s J a the matrix-valued function Q ( · , s ) : R R n × n , Q ( t , s ) = { q k j ( t , s ) } k , j n is called a solution of the IP (5), (7) for t ( a , ) , if Q ( · , s ) is continuous in t on J a , satisfying the matrix Equation (5) for t ( a , ) , as well as the initial condition (7).
The solution of the IP (5), (6) C ( t , s ) is called the fundamental matrix of the system (2), and obviously, C ( t , a ) = Q ( t , a ) .
In our exposition below, we will need three theorems presented below, one of them in a slightly modified version.
Theorem 1. 
([35], page 17) Let the following conditions be fulfilled:
(i)
The conditions(S)hold and h max ( τ , σ ) is arbitrary fixed number.
(ii)
The functions G C ( [ a h , ) , R n ) .
Then, the function Y ( t ) = h 0 [ d θ V ( t , θ ] G ( t + θ ) is continuous in t for ( t , θ ) J a × [ h , 0 ] .
Theorem 2. 
([36] Krasnosel’skii’s fixed point theorem) Let E be a Banach space with norm · , M be a nonempty, closed and convex subset of E and for the maps T , S : M E , the following conditions hold:
(i)
Tis contraction with constant β ( 0 , 1 ) ;
(ii)
Sis continuous and the set S ( M ) is contained in a compact set;
(iii)
T x + S y M for every x , y M .
Then, there exist a y M with T y + S y = y .
Theorem 3. 
(Theorem 4 in [16]) Let the following conditions hold:
(i)
The conditions(S)hold;
(ii)
sup t J a V a r θ [ h , 0 ] | V ( t , θ ) | = V < 1 ;
(iii)
For each solution of IP (4), (3) with arbitrary initial function Φ C the following inequality holds:
Y ( t ) a ( t ) + g ( t ) a t ( I 1 ( Γ ( α ) ) I α 1 t η ) Y ( η ) d τ
for t [ a , T ) , where a ( t ) L 1 l o c ( [ a , T ) , R + ) and g ( t ) C ( [ a , T ) , R + ) are nondecreasing, a < T .
Then, the solution of IP (4), (3) is Mittag–Leffler (ML) bounded of order α M for t [ a , T ) , i.e.
Y ( t ) a ( t ) E α M ( g ( t ) Γ ( α M ) t α M ) .

3. Main Results

The aim of this section is two-fold: First, with Theorem 4, we provide a positive answer of the open problem stated in [33]. Second, we study the existence and uniqueness of the solutions of the IP (4) and (3) in the cases when Φ P C * .
The next definitions clear the possible interaction between the concentrated delays in the neutral part of the system, i.e. the low terminal a of the fractional derivatives, in the case when Φ PC in both cases, i.e. a S Φ and a S Φ .
Definition 6. 
[20] The low terminal a will be called a noncritical point (noncritical jump point) for some initial function Φ BL relative to the delay τ l ( t ) , l r if the equality τ l ( a ) = 0 implies that there exists a constant ε ( 0 , h ] (eventually depending from τ l ), such that t a < τ l ( t ) t for t ( a , a + ε ] .
Definition 7. 
[20] The low terminal a for arbitrary function Φ BL 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 ε ( 0 , h ] (eventually depending from τ l ), such that t a τ l ( t ) for t [ a , a + ε ] .
It is simple to see that without loss of generality, we can renumber all delays τ l ( t ) , so that those for which a S Φ ( a S Φ ) is a noncritical or critical point (noncritical or critical jump point), to have the numbers τ 1 ( t ) , . . . , τ q ( t ) , q r . In the next exposition, for convenience, we assume that this renumbering is made.
Theorem 4. 
Let the following conditions hold:
(i)
The conditions(S)hold;
(ii)
The low terminal a is either a critical point ( a S Φ ), or a critical jump point ( a S Φ ), relative the delays τ 1 ( t ) , . . . , τ q ( t ) , 1 q r .
Then, either l q A ¯ l ( a ) < 1 or there exists ε ( 0 , h ] such that V d l ( t , θ ) Θ for t [ a , a + ε ] , θ [ h , 0 ] and l q .
Proof. 
Assume the contrary, that l q A ¯ l ( a ) 1 . Then, there exists ε 1 > 0 such that l q A ¯ l ( t ) 1 2 for t [ a , a + ε 1 ] .
From condition (S2), it follows that δ > 0 , such that for t [ a , a + ε 1 ] , we have
l q Var θ [ δ , 0 ] V d l ( t , θ ) 1 4 .
Since V d l ( t , θ ) V d l ( t , 0 ) + V a r θ [ δ , 0 ] V d l ( t , θ ) = V a r θ [ δ , 0 ] V d l ( t , θ ) , we have l q V d l ( t , θ ) l q V a r θ [ δ , 0 ] V d l ( t , θ ) 1 4 for t [ a , a + ε 1 ] too.
For definiteness, we assume that the delays with critical points nave numbers from 1 to q and we will consider both cases: ( i ) a S Φ , then a critical jump point for the delays τ 1 ( t ) , , τ q ( t ) , and ( ii ) a S Φ , which is a critical point for these delays. In both cases, two possibilities exist:
(a)
There exists ε 2 ( 0 , ε 1 ] such that τ l ( t ) > 0 , for t ( a , a + ε 2 ] and l q ;
(b)
There exists a monotone decreasing sequence { t j } j N J a with lim j t j = a such that τ l ( t j ) = τ l ( a ) = 0 , j N for some l q .
( i ) Consider first the case when a S Φ is a critical jump point for some l q . Then, in the case (b), the set S Φ includes infinitely many points, and hence, we will have at least one accumulation point which contradicts condition (S5).
For case (a), let t * ( a , a + ε 2 ] be arbitrary with τ l ( t * ) ( 0 , δ ) . Then, for θ [ τ l ( t * ) , 0 ] , we have:
l q V d l ( t * , θ ) l q V d l ( t * , 0 ) + V a r θ [ τ l ( t * ) , 0 ] V d l ( t * , θ ) = l q V a r θ [ τ l ( t * ) , 0 ] V d l ( t * , θ ) l q V a r θ [ δ , 0 ] V d l ( t * , θ ) 1 4 .
From the other side for θ [ τ l ( t * ) , 0 ] , we have:
l q V d l ( t * , θ ) = l q A ¯ l ( t * ) H ( θ + τ l ( t * ) ) = l q A ¯ l ( t * ) 1 2 ,
which contradicts with (8).
( ii ) Let a S Φ be a critical point for the delays. It is simple to see that case (a) can be treated as case (a) in the former point ( i ) .
In case (b), we have two possibilities: the first is that for all l q , there exists ε > 0 , such that τ l ( t ) 0 for t [ a , a + ε ] . Then, we obtain that H ( θ + τ l ( t ) ) = H ( θ ) = 0 holds for t [ a , a + ε 2 ] and l q . Then, we have:
l q V d l ( t , θ ) = l q A ¯ l ( t ) H ( θ + τ l ( t ) ) = l q A ¯ l ( t ) H ( θ ) = Θ .
The second possibility is that there exists a monotone decreasing sequence { t j } j N J a , with lim j t j = a , such that τ l ( t j ) = τ ( a ) = 0 , j N for some l q and there exists j n 0 N , such that for j j n 0 , we have that δ > τ l ( t j * ) > 0 , where j j n 0 and t j * = max t [ t j + 1 , t j ] τ l ( t ) .
Then, for j j n 0 and θ [ τ l ( t j * ) , 0 ] , we have:
l q V d l ( t j * , θ ) = l q A ¯ l ( t j * ) H ( θ + τ l ( t j * ) ) = l q A ¯ l ( t * ) 1 2 .
From the other side, as in (8) for θ [ τ l ( t j * ) , 0 ] , we obtain:
l q V d l ( t j * , θ ) l q V a r θ [ τ l ( t * ) , 0 ] V d l ( t j * , θ ) l q V a r θ [ δ , 0 ] V d l ( t j * , θ ) 1 4 ,
which contradicts with (10). □
Remark 1. 
Both Definitions 6 and 7 lead to the following essentially question: are there only these two possibilities or do others exist? In fact, in both definitions, we suppose that either there exists ε ( 0 , h ] such that the equation t τ l ( t ) = a has no roots for each t ( a , a + ε ] and l r . Then, for t ( a , a + ε ] and l r , we have that either t τ l ( t ) < a (noncritical case), or t τ l ( t ) a (critical case). The answer is that in the case when a S Φ , from the condition ( S 5 ) , it follows that there are no other possibilities. Indeed, if we assume that for every ε ( 0 , h ] , the equation t τ l ( t ) = a has at least one root for some l r , then the set l r S Φ l will have an accumulation point, which contradicts with condition ( S 5 ) . In the case when a S Φ , it is possible that the case ( ii ) (b) can appear in addition. Thus, from Theorem 4, it follows that the Condition 3 in Theorem 2 in [20] is unnecessary.
Let J L = { t 0 , . . . , t L } , L N , with J L ( a h , a ] being an arbitrary finite particle of [ a h , a ] and defining the real linear subspaces of PC and P C * as follows:
E ( J L ) = { Φ : [ a h , a ] R n | Φ PC , S Φ J L } , E * ( J L ) = { Φ : [ a h , a ] R n | Φ P C * , S Φ J L } .
It is clear that since P C * PC , then E * ( J L ) is a subspace of E ( J L ) .
Following the approach introduced in [18], we can define the real linear space:
E = { G : [ a h , ) R n | G | J a C ( J a , R n ) , G ( t ) = Φ ( t ) , t [ a h , a ] , Φ E ( J L ) } ,
and its linear subspace:
E * = { G ¯ : [ a h , ) , R n | G ¯ | J a C ( J a , R n ) B V l o c ( J a , R n ) , G ¯ ( t ) = Φ ( t ) , t [ a h , a ] , Φ E * ( J L ) } .
For arbitrary fixed M ( 0 , h ] , the real linear space is as follows:
E M = { G : [ a h , a + M ] R n | G = G | [ a h , a + M ] , G E } ,
endowed with the norm G s = sup t [ a h , a + M ] G ( t ) , which is the real Banach space and considers the linear subspace of E M :
E M * = { G : [ a h , a + M ] R n | G = G ¯ | [ a h , a + M ] , G ¯ E * } .
The subspace E M * is endowed with the following norm:
G V a r = G ( a h ) + V a r s [ a h , a + M ] G ( s ) = G ( a h ) + V a r s [ a h , a ] Φ ( s ) + V a r s [ a , a + M ] G ( s ) = G ( a h ) + k = 1 n V a r s [ a h , a ] | ϕ k ( s ) | + k = 1 n V a r s [ a , a + M | g k ( s ) | ,
which is a Banach space concerning the norm · V a r for every M ( 0 , h ] (see [25]). The important question which needs an answer is as follows: is the space E M * endowed with the norm · s a closed subspace of E M ? It must be noted that the following lemma that answers the question above is a generalization of Lemma 1 in [26] in the case of Banach spaces.
Lemma 1. 
For arbitrary fixed M ( 0 , h ] , the space E M * endowed with the norm · s is a closed subspace of E M .
Proof. 
The proof of the statement uses the idea of the proof of Lemma 1 in [26], and in that way, we will mainly be sketching the differences.
Let M ( 0 , h ] be an arbitrary fixed number and let { G q ( t ) = ( g 1 q ( t ) , . . . , g n q ( t , s ) ) T } l = 1 E M * be an arbitrary Cauchy sequence under the norm · s . That means that for each ε > 0 , there exists a number q 0 ( ε ) N , such that for every p N , we have G q 0 G q 0 + p s < ε , and since E M is a Banach space, there exists G 0 ( t , s ) = ( g 1 0 ( t ) , . . . , g n 0 ( t ) ) T E M , such that lim p G 0 G q 0 + p s = 0 .
Thus, there exists a number q 1 ( ε ) q 0 ( ε ) , such that for every p N , we have Φ 0 Φ q 1 + p s < ε , where Φ 0 = G 0 | [ a h , a ] , Φ q = G q | [ a h , a ] , and hence, Φ 0 ( t ) has finitely many jumps and S Φ 0 J N .
For an arbitrary partition P N = { t 0 , . . . , t N } , N N , of the interval [ a h , a + M ] , t 0 = a h , t N = a + M , we have:
i = 0 N 1 | g k q ( t i + 1 ) g k q ( t i ) | V k q ,
for every q N and arbitrary fixed k n , where V q = V a r t [ a h , s + M ] G q ( t ) = k = 1 n V a r t [ a h , s + M ] g k q ( · , s ) = k = 1 n V k q . Then, there exists a number q 2 ( ε ) q 1 ( ε ) , such that G q 2 G q 2 + p s < ε for each p N and the inequality | g k q 2 ( t i ) g k q 2 + p ( t i ) | < ε 2 N holds for t j [ a h , a + M ] . Thus, for j N 1 0 , the inequalities:
g k q 2 ( t j + 1 ) ε 2 N < g k q 2 + p ( t j + 1 , s ) < g k q 2 ( t j + 1 ) + ε 2 N , g k q 2 ( t j ) ε 2 N < g k q 2 + p ( t j ) < g k q 2 ( t j ) + ε 2 N
holds too.
For each p N and j N 1 0 , by adding the inequalities (11), we obtain:
| g k q 2 + p ( t j + 1 ) g k q 2 + p ( t j ) | < | g k q 2 ( t j + 1 ) g k q 2 ( t j ) | + ε N .
From the inequalities (11) and (13) for each p N , the following estimation is obtained:
j = 0 N 1 | g k q 2 + p ( t j + 1 ) g k q 2 + p ( t j ) | < j = 0 N 1 | g k q 2 ( t j + 1 ) g k q 2 ( t j ) | + ε < V k q 2 + ε .
Therefore, the sequence { V q } q N is bounded from above and lets us denote V = sup q N { V q } q N .
We will prove that G 0 ( t ) has bounded the variation on [ a h , a + M ] .
Since lim p G 0 G q 0 + p s = 0 , there exists a number q 3 ( ε ) q 2 ( ε ) , such that for each p N , we have that G 0 G q 2 + p s < ε , and for j N 1 0 , the following inequalities | g k q 3 ( ε ) + p ( t j ) g k 0 ( t j ) | < ε 2 N hold. Then, as above, we can prove that:
j = 0 N 1 | g k 0 ( t j + 1 ) g k 0 ( t j ) | < j = 0 N 1 | g k q 3 ( ε ) + p ( t j + 1 ) g k q 3 ( ε ) + p ( t j ) | + ε V + ε
hence, G 0 E M * .
The statement of the next theorem is essentially based off of the statement of Theorem 4, and in this sense, it is a generalization of Theorem 2 in [20], i.e. the statement is still true under weaker conditions.
Theorem 5. 
Let the conditions(S)hold.
Then, for arbitrary finite partition J L = { t 0 , . . . , t L } , L N of [ a h , a ] ( J L ( a h , a ] ) and each initial function Φ E * ( J L ) , there exists M * R + such that the IP (4), (3) has at least one solution X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T in the sense of Definition 3 with an interval of existence [ a , a + M * ] .
Proof. 
Let J L = { t 0 , . . . , t L } , L N be the arbitrary finite partition of [ a h , a ] , with J L ( a h , a ] , and Φ * E * ( J L ) be an arbitrary fixed function, introducing the subset E M Φ * = { G Φ * E M * | G ( t ) = Φ * ( t ) , t [ a h , a ] } , E M Φ * E M * . The set E M Φ * is nonempty, convex and closed (in virtue of the Lemma 1) subset of E M * for arbitrary M ( 0 , h ] , concerning the norm G s = sup t [ a h , a + M ] G ( t ) , G E M * .
Define the operator T for arbitrary t ( a , a + M ] and G E M * as follows:
( T G ) ( t ) = C Φ ( 0 ) + h 0 [ d θ V ( t , θ ) ] G ( t + θ ) ; ( T G ) ( t ) = Φ ( t a ) , t [ a h , a ] , Φ P C * ; lim t a + M 0 ( T G ) ( t ) = ( T G ) ( a + M ) .
For arbitrary M ( 0 , h ] and G E M * , define the operator S via the following equalities:
( S G ) ( t ) = I 1 ( Γ ( α ) ) a t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] G ( η + θ ) + F ( η ) ) d τ , t ( a , a + M ] ; ( S G ) ( t ) = 0 , t [ a h , a ]
and for shortness, we rewrite the system (4) in the form G ( t ) = ( T G ) ( t ) + ( S G ) ( t ) .
As first, we will prove that T E M Φ * E M Φ * for arbitrary fixed function Φ * E * ( J L ) , and all M ( 0 , h ] are small enough.
Let G E M Φ * and M ( 0 , τ ] be arbitrary. Then, we have:
( T G ) ( t ) = C Φ * ( 0 ) + h 0 [ d θ V d ( t , θ ) ] G Φ * ( t + θ ) + h 0 [ d θ V c ( t , θ ) ] G Φ * ( t + θ ) .
Let γ R + be the constant existing according to condition (S4), and hence, we have V c ( t , θ ) and the continuous functions for ( t , θ ) [ a , a + γ ] × [ h , 0 ] .
Then, by integrating, by parts, the third addend in the right side of (16), we have:
h 0 [ d θ V c ( t , θ ) ] G ( t + θ ) = V c ( t , h ) G ( t h ) + h 0 [ d θ G ( t + θ ) ] V c ( t , θ ) = V c ( t , h ) Φ * ( t a h ) + h 0 [ d θ G ( t + θ ) ] V c ( t , θ ) = V c ( t , h ) Φ * ( t a h ) + h a t [ d θ Φ * ( t a + θ ) ] V c ( t , θ ) + a t 0 [ d θ G ( t + θ ) ] V c ( t , θ ) .
Since Φ * ( t a + θ ) is right continuous at θ = h , then from condition (S4), it follows that there exists M 1 ( 0 , min ( M , γ ) ) such that the first addend in the right side of (17) are continuous functions for t [ a , a + M 1 ] . For the second and third addends, we have that V c ( t , θ ) is continuous in [ a , a + M 1 ] × [ h , 0 ] , while the functions Φ * ( t a + θ ) and G ( t + θ ) have bounded variation for t + θ [ a h , a ] and t + θ [ a , a + M 1 ] , respectively. Then, for ( t , θ ) [ a , a + M 1 ] × [ h , 0 ] , and hence, for t + θ [ a h , a + M 1 ] , the second and third addends in the right side of (17), according to Theorem 1, are continuous functions in t for t [ a , a + M 1 ] . Thus, the third addend in the right side of (16) is also a continuous function for t [ a , a + M 1 ] . Below, for definiteness, we will assume that 1 q 1 q r , where the critical points are numbered from 1 to q 1 the noncritical points are numbered from q 1 + 1 to q . Then, for the second addend in the right side of (16) we have:
l q h 0 [ d θ V d l ( t , θ ) ] G Φ * ( t + θ ) ) = l q 1 A l ( t ) G Φ * ( t τ l ( t ) ) + l q q 1 A l ( t ) Φ * ( t a τ l ( t ) ) + l r q A l ( t ) Φ * ( t a τ l ( t ) ) ,
where τ l ( a ) > 0 for l r q . From Definition 7, it follows that there exists a constant ε 1 ( 0 , h ] , such that t τ l ( t ) a for t [ a , a + ε 1 ] and l q 1 . Since τ l ( a ) > 0 for l r q and taking into account Lemma 3 in [25] and Definition 7, we conclude that there exists a constant ε 2 ( 0 , ε 1 ] such that t τ l ( t ) < a for t ( a , a + ε 2 ] and Φ * ( t a τ l ( t ) ) is a continuous function in t ( a , a + ε 2 ] and l r q 1 . Thus, the left side of (18) is a continuous function for t ( a , a + ε 2 ] (only right continuous at t = a ), and hence, ( T G Φ * ) ( t ) for all M ( 0 , M 2 ] , M 2 = min ( M 1 , ε 2 ) is also continuous function. From (11), it follows that lim t a + 0 ( T G Φ * ) ( t ) = Φ * ( 0 ) and lim t a + M 0 ( T G Φ * ) ( t ) = ( T G Φ * ) ( a + M ) .
Let J N = { t 0 , . . . , t N } , N N be an arbitrary particle of [ a , a + M ] , M ( 0 , M 2 ] , and, then we have:
z N 1 0 h 0 G ( t z + 1 + θ ) d θ V ( t z + 1 , θ ) h 0 G ( t z + θ ) d θ V ( t z , θ ) V a r s [ a h , a + M ] G ( s ) sup t [ a , a + M ] V a r θ [ h , 0 ] V ( t , θ ) ,
where C > 0 is a constant not depending from N . Thus, T E M Φ * E M Φ * , and hence, T E M * E M * .
To verify that condition 1 of Theorem 2 holds, we will prove that for arbitrary fixed function Φ * P C * and all M ( 0 , h ] are small enough so that the operator T is a contraction in E M Φ * .
Let G , G ˜ E M Φ * be arbitrary, and then for all t [ a , a + M 2 ] , M ( 0 , M 2 ] , we have:
T G T G ¯ s = sup t [ a h , a + M ] T G ( t ) T G ˜ ( t ) sup t [ a h , a + M ] h 0 [ d θ V d ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) + sup t [ a h , a + M ] h 0 [ d θ V c ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) .
Applying Lemma 1, we have l q 1 A ¯ l ( a ) = β < 1 , and there exists M 3 ( 0 , M 2 ] such that for t [ a , a + M 3 ] , we have l q 1 A ¯ l ( t ) β + ( 1 β ) 8 . Condition (S2) implies that there exists δ > 0 , such that for θ [ δ , 0 ] and t [ a , a + M 4 ] , M 4 = min ( M 3 , δ ) we have sup t [ a , a + M ] V a r θ [ δ , 0 ] V ( t , θ ) = ( 1 β ) 8 and
Φ * ( t a τ l ( t ) ) = G ( t + τ l ( t ) ) = G ˜ ( t + τ l ( t ) )
when l r q 1 .
Hence, from (18) for the first addend in the right side of (19) when t [ a , a + M 4 ] , it follows that:
sup t [ a h , a + M ] h 0 [ d θ V d ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) ( β + 1 β 8 ) G G ˜ s .
For the second addend in the right side of (19), taking into account that t + θ a and l r (23) hold, we obtain:
sup t [ a h , a + M ] h 0 [ d θ V c ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) sup t [ a h , a + M ] δ 0 [ d θ V c ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) + sup t [ a h , a + M ] h δ [ d θ V c ( t , θ ) ] ( G ( t + θ ) G ˜ ( t + θ ) ) sup t [ a , a + M ] V a r θ [ δ , 0 ] V ( t , θ ) G G ˜ s = ( 1 β ) 8 G G ˜ s .
Thus, from (24) and (25), it follows that the operator T is a contraction in E M Φ * , since ( 1 + 3 β ) 4 < 1 .
To verify that condition 2 of Theorem 2 holds, we must prove that S is continuous and the set S ( E M Φ * ) is a relative compact set.
For arbitrary G E M * from condition (S2), it follows that:
h 0 [ d θ U ( τ , θ ) ] G ( τ + θ ) U * G M , U * = sup t [ a , a + M ] U * ( t )
hence, the integral on the left side is at least locally bounded by the Lebesgue integrable function in t for t [ a , a + M ] , M ( 0 , h ] .
Estimating the following integral:
a t I 1 ( Γ ( α ) ) I α 1 ( t η ) d η = k n Γ 1 ( α k ) a t ( t η ) α k 1 d η = k n ( t a ) α k Γ ( 1 + α k ) ( t a ) α m i n n Γ 1 ( z min ) .
Taking into account (23) for every M ( 0 , M 5 ] and t [ a , a + M ] , M 5 = min ( 1 , M 4 ) , we have that:
( S G ) ( t ) M = sup t [ a h , a + M ] a t I 1 ( Γ ( α ) ) I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] G ( η + θ ) + F ( η ) ) d τ ( U * G M + F M ) sup t [ a h , a + M ] a t I 1 ( Γ ( α ) ) I α 1 ( t η ) d τ ( t a ) α m i n n Γ 1 ( z min ) ( U * G M + F M ) ,
where F M = sup t [ a , a + M ] F ( t ) . Then, from (24), it follows that lim t a + 0 ( S G ) ( t ) = 0 , and we can conclude that S E M * E M * .
Let G 1 , G 2 E M * be arbitrary and t [ a , a + M ] , M ( 0 , h ] . Then, using (24), we obtain:
( S G 1 ) ( t ) ( S G 2 ) ( t ) M sup t [ a , a + M I 1 ( Γ ( α ) ) a t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] G 1 ( η + θ ) G 2 ( η + θ ) ) d τ sup t [ a , a + M a t I 1 ( Γ ( α ) ( I α 1 ( t η ) h 0 [ d θ U ( η , θ ) ] G 1 ( η + θ ) G 2 ( η + θ ) ) d τ d τ U * G 1 G 2 M sup t [ a , a + M a t I 1 ( Γ ( α ) ( I α 1 ( t η ) d τ M α m i n U * n Γ 1 ( z min ) G 1 G 2 M
hence, from (25), it follows that the map S : E M * E M * is continuous.
Let R Φ * be an arbitrary fixed number, denoted by:
B ( G * , R ) = { G E M * | G G * M R } , where G * ( t ) = Φ * ( t a ) , t [ a h , a ] ; Φ * ( 0 ) , t [ a , a + M ] .
From (24) for arbitrary G M B ( G * , R ) , we have the following estimation:
( S G ) ( t ) M M α m i n n Γ 1 ( z min ) ( U * R + F M )
hence, the set S ( B ( G * , R ) ) is uniformly bounded, and thus, S maps every bounded subset of E M * in a uniformly bounded subset of E M * .
To apply Theorem 2, we must prove that the set S ( B ( G * , R ) ) is relatively compact, and according to Arzela–Ascoli’s theorem, it is enough to prove that the set S ( B ( G * , R ) ) is equicontinuous.
Let ε > 0 , t 1 , t 2 [ a , a + M ] , M ( 0 , M 5 ] be arbitrary, and for definiteness, assume that t 1 < t 2 . Then, for every G B ( G * , R ) , when t 2 t 1 < δ , δ < ε C * ( n + 1 ) 1 α m i n , C * = n Γ 1 ( z min ) ( U * R + F M ) , we have the following estimation:
( S G ) ( t 2 ) ( S G ) ( t 1 ) I 1 ( Γ ( α ) ) t 1 t 2 I α 1 ( t 2 η ) ( h 0 [ d θ U ( η , θ ) ] G ( η + θ ) + F ( η ) ) d τ + I 1 ( Γ ( α ) ) a t 1 ( I α 1 ( t 2 η ) I α 1 ( t 1 η ) ) ( h 0 [ d θ U ( η , θ ) ] G ( η + θ ) ) d τ ( t 2 t 1 ) α m i n n Γ 1 ( z min ) ( U * R + F M ) + U * R a t 1 I 1 ( Γ ( α ) ( I α 1 ( t 2 η ) I α 1 ( t 1 η ) ) d τ ( t 2 t 1 ) α m i n n Γ 1 ( z min ) ( U * R + F M ) + U * R k n Γ 1 ( α k ) a t 1 ( t 2 η ) α k 1 ( t 1 η ) α k 1 d η C * ( t 2 t 1 ) α m i n + U * R k n Γ 1 ( α k ) a t 1 ( ( t 1 η ) α k 1 ( t 2 η ) α k 1 ) d η C * ( t 2 t 1 ) α m i n + U * R Γ 1 ( z min ) k n ( ( t 1 a ) α k ( t 2 a ) α k + ( t 2 t 1 ) α k ) C * ( t 2 t 1 ) α m i n + C * k n ( α k ( t 1 t 2 ) ( ξ a ) α k 1 + ( t 2 t 1 ) α k ) C * ( t 2 t 1 ) α m i n + C * k n ( t 2 t 1 ) α k ( 1 α k ( t 2 t 1 ) 1 α k ( t 2 a ) 1 ) C * ( t 2 t 1 ) α m i n ( 1 + k n ( 1 α k ( t 2 t 1 ) 1 α k ( t 2 a ) 1 ) C * ( t 2 t 1 ) α m i n ( n + 1 ) < ε .
Thus, we proved that the set S ( B ( G * , R ) ) is equicontinuous, and hence, the map S : E M * E M * is compact.
Let G E M * and G ˜ E M Φ * be arbitrary. Then, t [ a , a + M ] , M ( 0 , M * ] , M * = M 5 , where the function T G ˜ ( t ) + S G ( t ) is continuous, lim t a + 0 ( T G ˜ ( t ) + S G ( t ) ) = Φ * ( 0 ) , T G ˜ ( t ) + S G ( t ) = Φ * ( t a ) , t [ a h , a ] , lim t a + M 0 ( T G ˜ ( t ) + S G ( t ) ) = ( T G ˜ ) ( a + M ) + S G ( a + M ) , and hence, T G ˜ ( t ) + S G ( t ) E M Φ * , i.e. the condition 3 of Theorem 2 holds. Thus, the system (4) has at least one fixed point in E M Φ * .
Theorem 6. 
Let the conditions(S)hold.
Then, for each Φ * P C * , there exists M * R + , such that the IP (4) and (3) has a unique solution X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T in the interval [ a , a + M * ] .
Proof. 
Let Φ * P C * be arbitrary. Then, according to Theorem 5, the IP (4) and (3) has at least one solution X ( t ) in the interval M * R + , and we can assume that there exist two different solutions, i.e. X 1 and X 2 , of the IP (4) and (3) in the same interval.
Then, the function Y ( t ) = sup s [ a , a + t ] X 2 ( s ) X 1 ( s ) > 0 for t [ a , a + M * ] , which is a continuous solution of the IP (4) and (3), with F ( t ) 0 and Φ ( t ) 0 .
Then, from (4) and (3), it follows that the following inequality holds:
Y ( t ) ( T Y ) ( t ) + S Y ( t ) β Y ( t ) + S Y ( t ) , β ( 0 , 1 )
hence:
Y ( t ) S Y ( t ) U * 1 β a t ( I 1 ( Γ ( α ) ) I α 1 t η ) Y ( η ) d τ .
Since in (26), we have that a ( t ) 0 , ; then, from Theorem 3, it follows that Y ( t ) 0 , which contradicts our assumption. Thus, the IP (4) and (3) has a unique solution in the interval [ a , a + M * ] .
Corollary 1. 
Let the conditions(S)hold.
Then, for each Φ * P C * , the IP (4) and (3) has a unique solution X ( t ) = ( x 1 ( t ) , . . . , x n ( t ) ) T with the interval of existence J a and X ( t ) B V l o c ( J a , R n ) .
Proof. 
From Theorems 3 and 4 in [19], it follows that the IP (4) and (3) has a unique solution X ( t ) in the sense of Definition 3 with the interval of existence J a . Then, the statement of the theorem follows from Lemma 1. □
Theorem 7. 
Let the conditions(S)hold.
Then, the following statements hold:
(i)
For each fixed s J a , the matrix IP (5), (6) has a unique solution C ( · , s ) : R R n × n with interval of existence J s , and C ( t , s ) = { c k j ( t , s ) } k , j n B V l o c ( R , R n ) ;
(ii)
For each s [ a h , a ] , the matrix IP (5), (7) has a unique solution Q ( · , s ) : R R n × n with the interval of existence J a , and Q ( t , s ) = { q k j ( t , s ) } k , j n B V l o c ( R , R n ) .
Proof. 
(i) Let j n and s J a be arbitrary fixed numbers, define the initial functions Φ * ( t , s ) = Φ C j ( t , s ) , where Φ C j ( t , s ) is the j-th column of the matrix function Φ C ( t , s ) , and consider the IP (2) and (3).
Then, according to Theorems 5 and 6, the IP (4) and (3) has unique solution C j ( t , s ) = ( c 1 j ( t , s ) , . . . , c n j ( t , s ) ) T , where (2) and (3) have a unique solution C j ( t , s ) = ( c 1 j ( t , s ) , , c n j ( t , s ) ) T , with Definition 2, with the interval of existence J s and from Corollary 1, it follows that C j ( · , s ) B V l o c ( R , R n ) .
This matrix C ( t , s ) = ( C 1 ( t , s ) , . . . , C n ( t , s ) ) is the unique fundamental matrix for the system (2).
Case (ii) can be treated in an analogical way. □

4. Applications

In this section, for the application of the obtained results concerning the fundamental matrix of the system (2), we establish that the problem of the existence of a unique resolvent kernel R ( t , s ) SVB l o c J × R , R n × n corresponding to the kernel K ( t , s ) SBV l o c ( J a × R , R n × n ) defined via (28), is equivalent to the problem of the existence of a unique fundamental matrix C ( t , s ) B V l o c ( R , R n ) of the system (2). Some stability results are also established.
Definition 8. 
( [33]) The function K ( t , s ) : J a × J a R n × n is called Stieltjes-Volterra type B kernel on J a × J a ( K SV B J a × J a , R n × n ), if the following conditions(K)hold:
(K1) 
The function t , s K ( t , s ) is measurable in t for each fixed s, right continuous in s on ( a , t ) and K ( t , s ) = Θ for s > t ;
(K2) 
K ( t , s ) is bounded, and the total variation in s of K ( t , s ) for every fixed t is uniformly bounded in s on J too.
With K SVB l o c J × R , R n × n , , we denote the set of kernels that restrictions to an arbitrary compact subset J , J J a belong to B .
Definition 9. 
A kernel R SVB l o c J a × R , R n × n ( R SVB J a × R , R n × n ) is called a Stieltjes-Volterra resolvent of type B , corresponding to a kernel K SVB l o c J a × R , R n × n ( K SVB J a × R , R n × n ) if for s R , t J { t s } , ( t J a { t s } ), it satisfies the following system:
R ( t , s ) = K ( t , s ) + s t d η [ K ( t , η ) ] R ( η , s ) = K ( t , s ) + s t d η [ R ( t , η ) ] K ( η , s ) ,
where the integrals in (27) are understood in the sense of Lebesgue-Stieltjes and J J a is an arbitrary compact subset.
According to Lemma 1 in [27] for arbitrary Φ C , the solution X ( t ) of the IP (4), (3) satisfies for t J a the Volterra-Stieltjes equation X ( t ) = s t [ d s K ( t , s ) ] X ( s ) + f ( t ) , where the function f ( t ) is given (i.e. it depends only on the kernels U and V, as well as on the initial functions Φ ( t ) and F ( t ) ) and the kernel K ( t , s ) has the following form:
K ( t , s ) = V ( t , s t ) + I 1 ( Γ ( α ) s t ( U ( τ , s τ ) I α 1 ( t τ ) d τ .
The condition (S) implies that K ( t , s ) SBV l o c J a × R , R n × n (see [27]) when the kernel K ( t , s ) , defined with (28), and for arbitrary fixed s J a , we introduce the matrix functions R ( t , s ) via the following relation:
R ( t , s ) = C ( t , s ) H ( t , s ) ,
where H ( t , s ) = I for t s and H ( t , s ) = Θ when t < s .
The next theorem is a generalization of Theorem 2 in [27] for the case when kernel K ( t , s ) is defined via (28) and satisfies only the conditions ( S ) , but possibly does not satisfy the conditions ( K ) and solves the open problem stated in the same work. Practically, we establish that the problem of the existence of a unique resolvent kernel R ( t , s ) SVB l o c J × R , R n × n corresponding to the kernel K ( t , s ) SBV l o c ( J a × R , R n × n ) defined via (28) is equivalent to the problem of the existence of a unique fundamental matrix C ( t , s ) B V l o c ( R , R n ) of the system (2).
Theorem 8. 
Let the following conditions hold:
(i)
The conditions ( S ) are fulfilled;
(ii)
The kernel K ( t , s ) have the form (28).
Then, the relation (29) holds if and only if when C ( t , s ) B V l o c J a × R , R n × n is the unique fundamental matrix of (2) and the function R ( t , s ) SVB l o c J a × R , R n × n defined via (29) is the unique solution of the resolvent Equation (27) corresponding to the kernel K ( t , s ) SBV l o c ( J a × R , R n × n ) defined with (28).
Proof. 
Sufficiency: Let C ( t , s ) B V l o c J a × R , R n × n is the unique fundamental matrix of (2), existing according to Theorem 7 and K ( t , s ) , with the form (28). Define the matrix functions R ( t , s ) via (29), and then, for s R , t J { t s } , ( t J a { t s } ), we have that:
C ( t , s ) = I + s t [ d η K ( t , η ) ] C ( η , s ) = H ( t , s ) + s t [ d η K ( t , η ) ] C ( η , s )
hence, we obtain:
R ( t , s ) = C ( t , s ) H ( t , s ) = s t [ d η K ( t , η ) ] C ( η , s ) = s t [ d η K ( t , η ) ] C ( η , s ) + H ( t , s ) H ( t , s ) = s t [ d η K ( t , η ) ] H ( t , s ) + s t [ d η K ( t , η ) ] C ( η , s ) H ( t , s ) = K ( t , s ) + s t [ d η K ( t , η ) ] R ( t , s ) .
Thus, (30) for s R , t J { t s } , ( t J a { t s } ) implies that R ( t , s ) is the unique solution of the resolvent Equation (27). Obviously, since C ( s , s ) = I and H ( t , s ) = Θ for t < s , then R ( t , s ) = Θ for t s . From the conditions (S) and since C ( t , s ) B V l o c J a × R , R n × n and (29), it follows that R ( t , s ) SVB l o c J × R , R n × n , and hence, R ( t , s ) is the unique resolvent kernel corresponding to the kernel K SVB l o c J a × R , R n × n defined via (28). The necessity can be proved in a reverse way. □
Corollary 2. 
Let the following conditions hold:
(i)
The conditions ( S ) are fulfilled;
(ii)
The kernel K ( t , s ) have the form (28).
Then, we have the resolvent kernel R ( t , s ) B V l o c J a × R , R n × n .
Proof. 
Since according Theorem 8, there exists a unique fundamental matrix C ( t , s ) B V l o c J a × R , R n × n , then the statement follows from (29) and Theorem 8. □
Definition 10. 
[6] The zero solution of (2) is called:
(a) 
Stable for a given a R if for any ϵ > 0 there is a δ = δ ( ε , a ) > 0 , such that X ( t ; a ; Φ ) ε for any initial function Φ P C with Φ < δ and t J a . In the opposite case, the solution is called unstable.
(b) 
Uniformly stable if for any ϵ > 0 , there is a δ ( ϵ ) > 0 , such that X ( t ; a ; Φ ) ε for any initial function Φ with Φ < δ and t J a .
(c) 
Locally asymptotically stable (LAS) if for a given a R , if it is stable and there is a Δ = Δ ( a ) > 0 , such that lim t X ( t ; a ; Φ ) = 0 for any initial function Φ with Φ < Δ . The set Ω ( a ) PC of all initial functions Φ for which lim t X ( t ; t ¯ ; Φ ) = 0 is called the attraction domain of the zero solution for initial time a . The zero solution is said to be uniformly LAS, if Δ > 0 is independent from the initial time.
(d) 
Globally asymptotically stable (GAS) if it is uniformly stable and lim t X ( t ; a ; Φ ) = 0 for any initial function Φ .
The next result gives a simple but useful necessary condition for the asymptotic stability of the zero solution of the system (2).
Theorem 9. 
Let the following conditions hold:
(i)
The conditions ( S ) are fulfilled;
(ii)
The kernel K ( t , s ) have the form (28);
(iii)
The zero solution of system (2) is stable.
Then, the corresponding resolvent kernel R ( t , s ) B V l o c J a × R , R n × n , for each fixed s R , is bounded in t, i.e. sup t J a R ( t , s ) < .
Proof. 
Let C ( t , s ) B V l o c J a × R , R n × n be the unique fundamental matrix of (2) existing according Theorem 7 and the kernel K ( t , s ) is defined via (28). Then, according to Theorem 8, there exists the corresponding resolvent kernel R ( t , s ) , which satisfies (29) for s R , t J a , and hence, R ( t , s ) B V l o c J a × R , R n × n . Condition 2 of the theorem implies that for each ε > 0 , there exists δ ( ε ) > 0 , such that for arbitrary initial function Φ P C with Φ < δ , the corresponding solution X ( t ; a ; Φ ) of the IP (2), (3) satisfies the inequality X ( t ; a ; Φ ) ε n for t J a .
Let s R be an arbitrary fixed number. Then, obviously, if s J a , or s J a and s > t , the statement of the theorem holds. Consider t J a { t s } and define the matrix C δ ( t , s ) = δ 2 C ( t , s ) = ( δ 2 C 1 ( t , s ) , , δ 2 C n ( t , s ) ) . Then, we obtain that C δ ( t , s ) ε for t J a { t s } , and hence, from (29), it follows that:
R ( t , s ) C ( t , s ) + 1 2 δ C δ ( t , s ) 1 + 2 ε δ ,
which completes the proof. □
Theorem 10. 
Let the following conditions hold:
(i)
The conditions ( S ) are fulfilled;
(ii)
The kernel K ( t , s ) has the form (28);
(iii)
The corresponding resolvent kernel is uniformly bounded, i.e. sup ( t , s ) J a × R R ( t , s ) < and sup s [ a h , a ] Q ( t , s ) sup s R R ( t , s ) < .
The zero solution of system (2) is stable.
Proof. 
Condition 3 implies that there exists a constant Q > 0 such that:
sup ( t , s ) J a × [ a h , a ] Q ( t , s ) Q
and for arbitrary initial function Φ P C * , the solution X ( t ; a ; Φ ) of the IP (2) and (3) has the following integral representation (see [15,17]):
X ( t ; a ; Φ ) = a h a Q ( t , s ) d Φ ˜ ( s ) ,
where Φ ˜ ( a h ) = 0 and Φ ˜ ( s ) Φ ( s ) , s ( a h , a ] .
From condition 3 and (31), it follows that:
sup t J a X ( t ; a ; Φ ) sup ( t , s ) J a × [ a h , a ] Q ( t , s ) V a r s [ a h , a ] Φ ˜ ( s ) Q V a r s [ a h , a ] Φ ( s ) Q Φ V a r .
Let ε > 0 be arbitrary and let δ ( 0 , ε Q ) be an arbitrary fixed number. Then, since Φ V a r Φ s from (32), it follows that for arbitrary initial function Φ P C * with Φ Var < δ , we have that sup t J a X ( t ; a ; Φ ) ε , and hence, the zero solution of system (2) is stable. □

5. Comments and Conclusions

As a motivation of our mathematical consideration, we note some possibilities of application of the studied systems as economics models.
As far we know, firstly, in the remarkable book [37], it was argued persuasively that delay differential equations are more suitable than ordinary differential equations alone or difference equations alone for an adequate treatment of dynamic economic phenomena. It is well known that there are at least two ways that time delays emerge in the dynamics of economic variables: there is some time lag between when the time economic decisions are made and the time the decisions bear fruit (see Chukwu [38]). There is a second "hidden" way, the way of rational expectation, see Fair [39] and Taylor [40] (inclusive expectation of bankruptcy [41]). In the second way, one assumes that expected future values of a variable are functions of the current and the past values of all relevant variables. In the monograph [42], the following model is introduced, describing the dynamic of the wealth of nations:
d X d t X ( t ) λ 0 t A ¯ ( t s ) X ( s ) d s = A 0 X ( t ) + λ 0 t A ( t s ) X ( s ) d s + F ( t ) ,
under the initial condition X ( t ) = Φ ( t ) , t ( [ h , 0 ] , where A ¯ ( t ) , A ( t ) C 1 ( R ¯ + , R n × n ) , A 0 R n × n , n N , X ( t ) , F ( t ) C ( R ¯ + , R n ) and Φ ( t ) C ( [ h , 0 ] , R n ) . The function F ( t ) summarize the government and the private controls and some structural factors. The proposed model was derived from familiar economic principles and was used to study the dynamics of six important economic factors: national income, interest rate, employment, value of capital stock, prices, and cumulative balance of payment. The function Φ ( t ) describes the past (historical data) of X ( t ) , as well as which data have an impact on the dynamics of the state of the economy. Note that for the considered model, the conditions (S) hold, and then the stability criteria proved in the present work can be used to study the stability properties of the model.
As a first result for the considered IP for a linear neutral system with distributed delays and derivatives in Caputo’s sense of incommensurate order, the existence of a unique solution is proven in the case when the initial functions are with a bounded variation. As corollary, we obtain the existence and uniqueness of a fundamental matrix for the homogeneous system, which has a bounded variation on each compact subinterval of J a . Second, without any additional assumptions of the boundedness type, it is established that the existence and uniqueness of a fundamental matrix lead to the existence and uniqueness of a resolvent kernel and vice versa. The explicit formula describing the relationship between the fundamental matrix and the resolvent kernel is proven in the general case too. Furthermore, on the base of the existence and uniqueness of a resolvent kernel, a necessary condition as well as a sufficient condition for the stability of the zero solution of the homogeneous system are established. Finally, a well-known economics model is considered, describing the dynamics of the wealth of nations, and we comment on the possibilities of the application of the obtained results for the considered systems, which include as a partial case the considered model. Note that the validation of our conclusion follows from the results proved in Section 3 and Section 4.

Author Contributions

Conceptualization, E.M., M.M. and T.S. Writing—Review and Editing, E.M., M.M. and T.S. The author contributions in this article are equal. 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: Ekaterina Madamlieva by the Bulgarian Ministry of Education and Science under the National Program “Young Scientists and Postdoctoral Students—2” (approved with RMS No. 206/ 7.04.2022), Stage I, 2022/2023, at the Faculty of Applied Mathematics and Informatics, Technical University of Sofia, Tsvetana Stoyanova by the Bulgarian National Science Fund under the project No. KP-06-H 45/7 entitled “Digitalisation of the Learning Process in Higher Education—Identification and Management Model”.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the anonymous reviewers for their very helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science BV: Amsterdam, The Netherlands, 2006. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equation; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  3. Jiao, Z.; Chen, Y.Q.; Podlubny, I. Distributed-Order Dynamic Systems: Stability, Simulation, Applications and Perspectives; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  4. Diethelm, K. The Analysis of Fractional Differential Equations, an Application-Oriented Exposition Using Differential Operators of Caputo Type; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010; Volume 2004. [Google Scholar]
  5. Stamova, I.; Stamov, G. Functional and Impulsive Differential Equations of Fractional Order: Qualitative Analysis and Applications; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  6. Kolmanovskii, V.; Myshkis, A. Introduction to the Theory and Applications of Functional Differential Equations; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999. [Google Scholar]
  7. Li, C.; Zhang, F. A survey on the stability of fractional differential equations. Eur. Phys. J. Spec. Top. 2011, 193, 27–47. [Google Scholar] [CrossRef]
  8. Kaczorek, T.; Idczak, D. Cauchy formula for the time-varying linear systems with Caputo derivative. Fract. Calc. Appl. Anal. 2017, 20, 494–505. [Google Scholar] [CrossRef]
  9. Gomoyunov, M.I. On representation formulas for solutions of linear differential equations with Caputo fractional derivatives. Fract. Calc. Appl. Anal. 2020, 23, 1141–1160. [Google Scholar] [CrossRef]
  10. Li, K.; Peng, J. Laplace transform and fractional differential equations. Appl. Math. Lett. 2011, 24, 2019–2023. [Google Scholar] [CrossRef] [Green Version]
  11. Krol, K. Asymptotic properties of fractional delay differential equations. Appl. Math. Comput. 2011, 218, 1515–1532. [Google Scholar] [CrossRef] [Green Version]
  12. Veselinova, M.; Kiskinov, H.; Zahariev, A. Stability analysis of linear fractional differential system with distributed delays. AIP Conf. Proc. 2015, 1690, 040013. [Google Scholar] [CrossRef]
  13. Wei, Y. The constant variation formulae for singular fractional differential systems with delay. Comput. Math. Appl. 2010, 59, 1184–1190. [Google Scholar] [CrossRef] [Green Version]
  14. Boyadzhiev, D.; Kiskinov, H.; Veselinova, M.; Zahariev, A. Stability analysis of linear distributed order fractional systems with distributed delays. Fract. Calc. Appl. Anal. 2017, 20, 914–935. [Google Scholar] [CrossRef]
  15. Madamlieva, E.; Konstantinov, M.; Milev, M.; Petkova, M. Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay. Mathematics 2020, 8, 364. [Google Scholar] [CrossRef] [Green Version]
  16. 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]
  17. Zhang, H.; Wu, D. Variation of constant formulae for time invariant and time varying Caputo fractional delay differential systems. J. Math. Res. Appl. 2014, 34, 549–560. [Google Scholar] [CrossRef]
  18. Boyadzhiev, D.; Kiskinov, H.; Zahariev, A. Integral representation of solutions of fractional system with distributed delays. Integral Transform. Spec. Funct. 2018, 29, 725–744. [Google Scholar] [CrossRef]
  19. Zhang, H.; Cao, J.; Jiang, W. General solution of linear fractional neutral differential difference equations. Discret. Dyn. Nat. Soc. 2013, 2013, 489521. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. Zahariev, A.; Kiskinov, H.; Angelova, E. Smoothness of the fundamental matrix of linear fractional system with variable delays. Neural Parall. Sci. Comput. 2019, 27, 71–83. [Google Scholar] [CrossRef]
  22. Kiskinov, H.; Madamlieva, E.; Veselinova, M.; Zahariev, A. Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays. Mathematics 2021, 9, 150. [Google Scholar] [CrossRef]
  23. Cermak, J.; Hornicek, J.; Kisela, T. Stability regions for fractional differential systems with a time delay. Commun. Nonlinear Sci. Numer. Simul. 2015, 31, 108–123. [Google Scholar] [CrossRef]
  24. Veselinova, M.; Kiskinov, H.; Zahariev, A. About stability conditions for retarded fractional differential systems with distributed delays. Commun. Appl. Anal. 2016, 20, 325–334. [Google Scholar]
  25. Li, M.; Wang, J. Finite time stability of fractional delay differential equations. Appl. Math. Lett. 2017, 64, 170–176. [Google Scholar] [CrossRef]
  26. Milev, M.; Zlatev, S. A note about stability of fractional retarded linear systems with distributed delays. Int. J. Pure Appl. Math. 2017, 115, 873–881. [Google Scholar] [CrossRef] [Green Version]
  27. Veselinova, M.; Kiskinov, H.; Zahariev, A. Stability analysis of neutral linear fractional system with distributed delays. Filomat 2016, 30, 841–851. [Google Scholar] [CrossRef] [Green Version]
  28. Cong, N.D.; Doan, T.S.; Siegmund, S.; Tuan, H.T. Linearized asymptotic stability for fractional differential equations. Electron. J. Qual. Theory Differ. Equ. 2016, 39, 1–13. [Google Scholar] [CrossRef]
  29. Zahariev, A.; Kiskinov, H. Asymptotic stability of delayed fractional system with nonlinear perturbation. AIP Conf. Proc. 2018, 2048, 050014. [Google Scholar] [CrossRef]
  30. 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]
  31. Singh, H. Numerical simulation for fractional delay differential equations. Int. J. Dynam. Control 2021, 9, 463–474. [Google Scholar] [CrossRef]
  32. Ahmed, A.I.; Al-Ahmary, T.A. Fractional–Order Chelyshkov Collocation Method for Solving Systems of Fractional Differential Equations. Math. Probl. Eng. 2022, 2022, 1–17. [Google Scholar] [CrossRef]
  33. Kiskinov, H.; Milev, M.; Zahariev, A. About the Resolvent Kernel of Neutral Linear Fractional System with Distributed Delays. Mathematics 2022, 10, 4573. [Google Scholar] [CrossRef]
  34. Hale, J.; Lunel, S. Introduction to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  35. Myshkis, A. Linear Differential Equations with Retarded Argument; Nauka: Moscow, Russia, 1972. (In Russian) [Google Scholar]
  36. Krasnosel’skii, M.A. Some problems of nonlinear analysis. Amer. Math. Soc. Transl. Ser. 1958, 2, 10. [Google Scholar]
  37. Gandolfo, G. Economic Dynamics, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  38. Chukwu, E.N. Stability and Time Optimal Control of Hereditary Systems; Academic Press: Boston, MA, USA, 1992. [Google Scholar]
  39. Fair, R. Specification, Estimation, and Analysis of Macroeconometric Models; Harvard University Press: Cambridge, MA, USA, 1984. [Google Scholar]
  40. Taylor, J.B. Macroeconomic Policy in a World Economy from Econometric Design to Practical Operation; W. W. Norton and Company: New York, NY, USA, 1993. [Google Scholar]
  41. Gigova, T.; Nikolova-Alexieva, V.; Valeva, K. Application of Business Crisis Forecasting Models in Bulgaria. IOP Conf. Ser. Mater. Sci. Eng. 2019, 618, 012071. [Google Scholar] [CrossRef]
  42. Chukwu, E.N. Differential Models and Neutral Systems for Controlling the Wealth of Nations; Series on Advances in Mathematics for Applied Sciences, 54; World Scientific Publishing Co. Pte. Ltd.: New York, NY, USA, 2001. [Google Scholar]
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

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

AMA Style

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(3):626. https://doi.org/10.3390/math11030626

Chicago/Turabian Style

Madamlieva, Ekaterina, Marian Milev, and Tsvetana Stoyanova. 2023. "On Stability Criteria Induced by the Resolvent Kernel for a Fractional Neutral Linear System with Distributed Delays" Mathematics 11, no. 3: 626. https://doi.org/10.3390/math11030626

APA Style

Madamlieva, E., Milev, M., & Stoyanova, T. (2023). On Stability Criteria Induced by the Resolvent Kernel for a Fractional Neutral Linear System with Distributed Delays. Mathematics, 11(3), 626. https://doi.org/10.3390/math11030626

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop