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
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
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 and arbitrarily, for each , the left-sided fractional integral operator and the left side Riemann–Liouville of order are defined via and , respectively.
The Caputo fractional derivative of same order by , where , denote the real linear space of all locally Lebesgue integrable functions and for the subspace of all locally bounded functions.
Below the following notations will be used: , , , , , , , , , , which are the identity and zero matrix, respectively, and is the zero vector-column.
For , , we denote , , , which denote the functions , which, for every fixed , has bounded variation in on every compact interval and denotes and
,
denote the functions with bounded variation in every compact interval
and
It is well known that for , the gamma function has a minimum at , where it attains the value (truncated).
With BL, we denote the Banach space of all vector functions , which are bounded and Lebesgue measurable on the interval with norm and denotes the set of all jumps point of arbitrary .
denote the subspace of all right continuous for piecewise continuous functions, , , and all of them are endowed with the same sup norm.
For
consider the inhomogeneous neutral linear delayed system with and distributed delays in the form:
and the corresponding homogeneous neutral linear system:
where the differential orders,
, can be incommensurate,
,
,
,
,
,
,
,
,
,
,
,
for every
, where
denotes the left side Caputo fractional derivative. For clarity, we rewrite the system (
2) in more detail form:
where
,
,
,
and for arbitrary
BL, we introduce the following initial condition:
In our consideration below, we need the following auxiliary system:
where
.
Definition 1. We say that for the kernels , the conditions(S)are fulfilled if for every , the following conditions hold (see [6,33]): - (S1)
The functions , are measurable in and normalized so that for , for for and every The kernels and are continuous from left in θ on and , respectively, for and for every fixed
- (S2)
For , the functionsand the kernel is uniformly nonatomic at zero (i.e. for every there exists such that for each , we have ([34]). - (S3)
For and , the Lebesgue decompositions of the kernels have the form: , where the indexes denoted the jump, the absolutely continuous and the singular part, respectively, in the Lebesgue decompositions and are the continuous part of these decompositions. In addition, where is the Heaviside function: for every
- (S4)
For each , the relations:hold, and there exists such that the kernels and are continuous in t when and - (S5)
The sets , , for every PC, do not have limit points.
Definition 2. The vector function is a solution of the initial problem (IP) (
1), (
3)
in if satisfies the system (
1)
for all ( ) and the initial condition (
3)
for , Definition 3. The vector function is a solution of the initial problem (IP) (
4),
where (
3)
in if , satisfies the system (
4)
for all () and the initial condition (
3)
for In this paper, we assume that the condition
(S) holds. For arbitrary fixed
, we define for
two initial matrix functions:
and introduce the auxiliary matrix initial problem:
Definition 4. For each , the matrix-valued function is called a solution of the IP (
5), (
6)
for if is continuous in t on , satisfying the matrix Equation (
5)
for , as well as the initial condition (6). Definition 5. For each the matrix-valued function is called a solution of the IP (
5), (
7)
for if is continuous in t on , satisfying the matrix Equation (
5)
for as well as the initial condition (
7)
. The solution of the IP (
5), (
6)
is called the fundamental matrix of the system (
2)
, and obviously, 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 is arbitrary fixed number.
- (ii)
The functions
Then, the function is continuous in t for
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 , the following conditions hold: - (i)
Tis contraction with constant ;
- (ii)
Sis continuous and the set is contained in a compact set;
- (iii)
for every
Then, there exist a with
Theorem 3. (Theorem 4 in [16]) Let the following conditions hold: - (i)
The conditions(S)hold;
- (ii)
;
- (iii)
For each solution of IP (
4), (
3)
with arbitrary initial function the following inequality holds:for where and are nondecreasing,
Then, the solution of IP (
4), (
3)
is Mittag–Leffler (ML) bounded of order for , i.e. 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
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 in both cases, i.e. and
Definition 6. [20] The low terminal a will be called a noncritical point (noncritical jump point) for some initial function relative to the delay if the equality implies that there exists a constant (eventually depending from ), such that for Definition 7. [20] The low terminal a for arbitrary function with () will be called a critical point (critical jump point) relative to some delay if the equality implies that there exists a constant (eventually depending from ), such that for It is simple to see that without loss of generality, we can renumber all delays , so that those for which () is a noncritical or critical point (noncritical or critical jump point), to have the numbers 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 (), or a critical jump point (), relative the delays
Then, either or there exists such that for and
Proof. Assume the contrary, that Then, there exists such that for
From condition
(S2), it follows that
such that for
, we have
Since , we have for too.
For definiteness, we assume that the delays with critical points nave numbers from 1 to q and we will consider both cases: , then a critical jump point for the delays , and , which is a critical point for these delays. In both cases, two possibilities exist:
- (a)
There exists such that for and ;
- (b)
There exists a monotone decreasing sequence with such that for some
Consider first the case when is a critical jump point for some Then, in the case (b), the set includes infinitely many points, and hence, we will have at least one accumulation point which contradicts condition (S5).
For case (a), let
be arbitrary with
Then, for
, we have:
From the other side for
, we have:
which contradicts with (
8).
Let 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
In case (b), we have two possibilities: the first is that for all
, there exists
, such that
for
Then, we obtain that
holds for
and
Then, we have:
The second possibility is that there exists a monotone decreasing sequence , with , such that for some and there exists , such that for , we have that where and
Then, for
and
, we have:
From the other side, as in (
8) for
, we obtain:
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 such that the equation has no roots for each and Then, for and , we have that either (noncritical case), or (critical case). The answer is that in the case when from the condition , it follows that there are no other possibilities. Indeed, if we assume that for every , the equation has at least one root for some then the set will have an accumulation point, which contradicts with condition In the case when , it is possible that the case (b) can appear in addition. Thus, from Theorem 4, it follows that the Condition 3 in Theorem 2 in [20] is unnecessary. Let
, with
being an arbitrary finite particle of
and defining the real linear subspaces of
and
as follows:
It is clear that since , then is a subspace of
Following the approach introduced in [
18], we can define the real linear space:
and its linear subspace:
For arbitrary fixed
, the real linear space is as follows:
endowed with the norm
, which is the real Banach space and considers the linear subspace of
:
The subspace
is endowed with the following norm:
which is a Banach space concerning the norm
for every
(see [
25]). The important question which needs an answer is as follows: is the space
endowed with the norm
a closed subspace of
? 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 , the space endowed with the norm is a closed subspace of
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 be an arbitrary fixed number and let be an arbitrary Cauchy sequence under the norm That means that for each , there exists a number , such that for every , we have , and since is a Banach space, there exists , such that
Thus, there exists a number , such that for every , we have where , and hence, has finitely many jumps and
For an arbitrary partition
of the interval
we have:
for every
and arbitrary fixed
where
Then, there exists a number
, such that
for each
and the inequality
holds for
Thus, for
, the inequalities:
holds too.
For each
and
, by adding the inequalities (
11), we obtain:
From the inequalities (
11) and (
13) for each
, the following estimation is obtained:
Therefore, the sequence is bounded from above and lets us denote
We will prove that has bounded the variation on
Since
, there exists a number
, such that for each
, we have that
, and for
, the following inequalities
hold. Then, as above, we can prove that:
hence,
□
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 of ( ) and each initial function there exists such that the IP (
4), (
3)
has at least one solution in the sense of Definition 3 with an interval of existence Proof. Let be the arbitrary finite partition of , with , and be an arbitrary fixed function, introducing the subset The set is nonempty, convex and closed (in virtue of the Lemma 1) subset of for arbitrary concerning the norm
Define the operator
for arbitrary
and
as follows:
For arbitrary
and
, define the operator
via the following equalities:
and for shortness, we rewrite the system (
4) in the form
As first, we will prove that for arbitrary fixed function , and all are small enough.
Let
and
be arbitrary. Then, we have:
Let be the constant existing according to condition (S4), and hence, we have and the continuous functions for
Then, by integrating, by parts, the third addend in the right side of (
16), we have:
Since
is right continuous at
, then from condition
(S4), it follows that there exists
such that the first addend in the right side of (
17) are continuous functions for
For the second and third addends, we have that
is continuous in
while the functions
and
have bounded variation for
and
, respectively. Then, for
, and hence, for
, the second and third addends in the right side of (
17), according to Theorem 1, are continuous functions in
t for
Thus, the third addend in the right side of (
16) is also a continuous function for
. Below, for definiteness, we will assume that
where the critical points are numbered from 1 to
the noncritical points are numbered from
to
Then, for the second addend in the right side of (
16) we have:
where
for
From Definition 7, it follows that there exists a constant
such that
for
and
Since
for
and taking into account Lemma 3 in [
25] and Definition 7, we conclude that there exists a constant
such that
for
and
is a continuous function in
and
Thus, the left side of (
18) is a continuous function for
(only right continuous at
), and hence,
for all
is also continuous function. From (
11), it follows that
and
Let
be an arbitrary particle of
, and, then we have:
where
is a constant not depending from
Thus,
, and hence,
.
To verify that condition 1 of Theorem 2 holds, we will prove that for arbitrary fixed function and all are small enough so that the operator is a contraction in
Let
be arbitrary, and then for all
, we have:
Applying Lemma 1, we have
, and there exists
such that for
, we have
Condition
(S2) implies that there exists
, such that for
and
we have
and
when
Hence, from (
18) for the first addend in the right side of (
19) when
, it follows that:
For the second addend in the right side of (
19), taking into account that
and
(
23) hold, we obtain:
Thus, from (
24) and (
25), it follows that the operator
is a contraction in
, since
To verify that condition 2 of Theorem 2 holds, we must prove that is continuous and the set is a relative compact set.
For arbitrary
from condition
(S2), it follows that:
hence, the integral on the left side is at least locally bounded by the Lebesgue integrable function in
t for
Estimating the following integral:
Taking into account (
23) for every
and
, we have that:
where
Then, from (
24), it follows that
and we can conclude that
Let
be arbitrary and
Then, using (
24), we obtain:
hence, from (
25), it follows that the map
is continuous.
Let
be an arbitrary fixed number, denoted by:
From (
24) for arbitrary
, we have the following estimation:
hence, the set
is uniformly bounded, and thus,
maps every bounded subset of
in a uniformly bounded subset of
To apply Theorem 2, we must prove that the set is relatively compact, and according to Arzela–Ascoli’s theorem, it is enough to prove that the set is equicontinuous.
Let
be arbitrary, and for definiteness, assume that
Then, for every
, when
, we have the following estimation:
Thus, we proved that the set is equicontinuous, and hence, the map is compact.
Let
and
be arbitrary. Then,
, where the function
is continuous,
, and hence,
, i.e. the condition 3 of Theorem 2 holds. Thus, the system (
4) has at least one fixed point in
□
Theorem 6. Let the conditions(S)hold.
Then, for each , there exists , such that the IP (
4)
and (
3)
has a unique solution in the interval Proof. Let
be arbitrary. Then, according to Theorem 5, the IP (
4) and (
3) has at least one solution
in the interval
, and we can assume that there exist two different solutions, i.e.
and
, of the IP (
4) and (
3) in the same interval.
Then, the function
for
, which is a continuous solution of the IP (
4) and (
3), with
and
Then, from (
4) and (
3), it follows that the following inequality holds:
hence:
Since in (
26), we have that
; then, from Theorem 3, it follows that
which contradicts our assumption. Thus, the IP (
4) and (
3) has a unique solution in the interval
□
Corollary 1. Let the conditions(S)hold.
Then, for each , the IP (
4)
and (
3)
has a unique solution with the interval of existence and Proof. From Theorems 3 and 4 in [
19], it follows that the IP (
4) and (
3) has a unique solution
in the sense of Definition 3 with the interval of existence
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 , the matrix IP (
5), (
6)
has a unique solution with interval of existence and ; - (ii)
For each , the matrix IP (
5), (
7)
has a unique solution with the interval of existence and
Proof. (i) Let
and
be arbitrary fixed numbers, define the initial functions
where
is the
j-th column of the matrix function
, and consider the IP (
2) and (
3).
Then, according to Theorems 5 and 6, the IP (
4) and (
3) has unique solution
, where (
2) and (
3) have a unique solution
, with Definition 2, with the interval of existence
and from Corollary 1, it follows that
This matrix
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
corresponding to the kernel
defined via (
28), is equivalent to the problem of the existence of a unique fundamental matrix
of the system (
2). Some stability results are also established.
Definition 8. ( [33]) The function is called Stieltjes-Volterra type kernel on ( ), if the following conditions(K)hold: - (K1)
The function is measurable in t for each fixed s, right continuous in s on and for ;
- (K2)
is bounded, and the total variation in s of for every fixed t is uniformly bounded in s on J too.
With , we denote the set of kernels that restrictions to an arbitrary compact subset belong to
Definition 9. A kernel () is called a Stieltjes-Volterra resolvent of type , corresponding to a kernel () if for (), it satisfies the following system:where the integrals in (
27)
are understood in the sense of Lebesgue-Stieltjes and is an arbitrary compact subset. According to Lemma 1 in [
27] for arbitrary
, the solution
of the IP (
4), (
3) satisfies for
the Volterra-Stieltjes equation
where the function
is given (i.e. it depends only on the kernels
U and
V, as well as on the initial functions
and
) and the kernel
has the following form:
The condition
(S) implies that
(see [
27]) when the kernel
, defined with (
28), and for arbitrary fixed
we introduce the matrix functions
via the following relation:
where
for
and
when
The next theorem is a generalization of Theorem 2 in [
27] for the case when kernel
is defined via (
28) and satisfies only the conditions
, but possibly does not satisfy the conditions
and solves the open problem stated in the same work. Practically, we establish that the problem of the existence of a unique resolvent kernel
corresponding to the kernel
defined via (
28) is equivalent to the problem of the existence of a unique fundamental matrix
of the system (
2).
Theorem 8. Let the following conditions hold:
- (i)
The conditions are fulfilled;
- (ii)
The kernel have the form (
28).
Then, the relation (
29)
holds if and only if when is the unique fundamental matrix of (
2)
and the function defined via (
29)
is the unique solution of the resolvent Equation (
27)
corresponding to the kernel defined with (
28).
Proof. Sufficiency: Let
is the unique fundamental matrix of (
2), existing according to Theorem 7 and
, with the form (
28). Define the matrix functions
via (
29), and then, for
(
), we have that:
hence, we obtain:
Thus, (
30) for
(
) implies that
is the unique solution of the resolvent Equation (
27). Obviously, since
and
for
, then
for
From the conditions
(S) and since
and (
29), it follows that
, and hence,
is the unique resolvent kernel corresponding to the kernel
defined via (
28). The necessity can be proved in a reverse way. □
Corollary 2. Let the following conditions hold:
- (i)
The conditions are fulfilled;
- (ii)
The kernel have the form (
28).
Then, we have the resolvent kernel
Proof. Since according Theorem 8, there exists a unique fundamental matrix
then the statement follows from (
29) and Theorem 8. □
Definition 10. [6] The zero solution of (
2)
is called: - (a)
Stable for a given if for any there is a such that for any initial function with and In the opposite case, the solution is called unstable.
- (b)
Uniformly stable if for any , there is a , such that for any initial function Φ with and
- (c)
Locally asymptotically stable (LAS) if for a given if it is stable and there is a , such that for any initial function Φ with The set of all initial functions Φ for which is called the attraction domain of the zero solution for initial time The zero solution is said to be uniformly LAS, if is independent from the initial time.
- (d)
Globally asymptotically stable (GAS) if it is uniformly stable and 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 are fulfilled;
- (ii)
The kernel have the form (
28);
- (iii)
The zero solution of system (
2)
is stable.
Then, the corresponding resolvent kernel , for each fixed , is bounded in t, i.e.
Proof. Let
be the unique fundamental matrix of (
2) existing according Theorem 7 and the kernel
is defined via (
28). Then, according to Theorem 8, there exists the corresponding resolvent kernel
, which satisfies (
29) for
, and hence,
Condition 2 of the theorem implies that for each
, there exists
such that for arbitrary initial function
with
, the corresponding solution
of the IP (
2), (
3) satisfies the inequality
for
Let
be an arbitrary fixed number. Then, obviously, if
, or
and
, the statement of the theorem holds. Consider
and define the matrix
. Then, we obtain that
for
, and hence, from (
29), it follows that:
which completes the proof. □
Theorem 10. Let the following conditions hold:
- (i)
The conditions are fulfilled;
- (ii)
The kernel has the form (
28);
- (iii)
The corresponding resolvent kernel is uniformly bounded, i.e. and
The zero solution of system (
2)
is stable. Proof. Condition 3 implies that there exists a constant
such that:
and for arbitrary initial function
, the solution
of the IP (
2) and (
3) has the following integral representation (see [
15,
17]):
where
and
From condition 3 and (
31), it follows that:
Let
be arbitrary and let
be an arbitrary fixed number. Then, since
from (
32), it follows that for arbitrary initial function
with
, we have that
, 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:
under the initial condition
, where
and
The function
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
describes the past (historical data) of
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
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.