Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays

: The goal of the present paper is to obtain sufﬁcient conditions that guaranty the existence and uniqueness of an absolutely continuous fundamental matrix for a retarded linear fractional differential system with Caputo type derivatives and distributed delays. Some applications of the obtained result concerning the integral representation of the solutions are given too.


Introduction
It is well known that many mathematical models of real world phenomena can be described more accurately through fractional derivative formulation. For more details on fractional calculus theory and fractional differential equations, we recommend the monographs of Kilbas et al. [1] and Podlubny [2]. For distributed order fractional differential equations see Jiao at all [3] and for an application-oriented exposition Diethelm [4]. Impulsive differential and functional differential equations with fractional derivative and some applications are studied by Stamova and Stamov [5].
The theme of the integral representation (variation of constants formula) of the solutions of linear fractional differential equations and/or systems (ordinary or with delay) is an "evergreen" theme for research. This explains why a lot of papers are devoted to different aspects of this problem. For linear fractional ordinary differential equations and systems, we refer the works [1,2,[6][7][8][9][10] and the references therein. Relatively, as far as we know, there are not many works devoted to the variation of constants formula for delayed linear fractional systems [11][12][13][14][15]. The case of neutral fractional systems is studied in [16][17][18][19][20].
The establishment of a fundamental matrix with appropriate properties (for example in [21] smoothness is obtained) is the basis for obtaining any integral representation and is a key tool in the study of different types of stability of linear and nonlinear disturbed systems (see [20]).
In the present work, we consider linear fractional systems with distributed delays and incommensurate order derivatives in the Caputo sense. The first goal of the work is to establish sufficient conditions for existence and uniqueness of a fundamental matrix C(t, s), which is absolutely continuous in t on every compact subinterval of R. The second one is to clarify the analytic properties in s, which are very similar to these in the integer case. As an application of the obtained results, some results concerning the integral representation of the solutions given in [15,18] are improved.
The paper is organized as follows. In Section 2 we recall the definitions of Riemann-Liouville and Caputo fractional derivatives with some of their properties. In the same section is the statement of the problem, as well as some necessary preliminary results used later. Section 3 is devoted to the existence and the uniqueness of the solutions of the Initial Problem (IP) for linear fractional systems with distributed delays and incommensurate order derivatives in the Caputo sense and special type discontinuous initial function. In Section 4 the existence and uniqueness of an absolutely continuous fundamental matrix is proved und its analytical properties are studied. Using the obtained results, in Section 5 we establish new integral representations for the solutions of the studied systems.

Preliminaries and Problem Statement
For readers convenience, below we recall the definitions of Riemann-Liouville and Caputo fractional derivatives as well as some needed properties. For details and other properties we refer to [1][2][3].
Let α ∈ (0, 1) be an arbitrary number and denote by L loc 1 (R, R) the linear space of all locally Lebesgue integrable functions f : R → R. Then for each a ∈ R, each t > a and f ∈ L loc 1 (R, R) the definitions of the left-sided fractional integral operator, the left side Riemann-Liouville and Caputo fractional derivatives of order α and some properties are given below (see [1]): If f ∈ AC(R, R) then the next formula gives a direct definition of the Caputo left side derivative: Everywhere below, the following notations will be used: R + = (0, ∞),R + = [0, ∞), J a = [a, ∞), a ∈ R, J s+M = [s, s + M], s ∈ J a , M ∈ R + , n = {1, 2, . . . , n}, m 0 = m ∪ {0}, n, m ∈ N, I, Θ ∈ R n×n denote the identity and zero matrix, respectively, I k , k ∈ n denotes the k-th column of the identity matrix and 0 ∈ R n is the zero element.

Definition 1.
A function Z(t) = (z 1 (t), . . . , z n (t)) T : J → R n is called piecewise absolutely continuous on some interval J ⊂ R (denoted Z(t) ∈ PAC(J, R n )), if it is with bounded variation (BV) in t on J, there is no singular term in the Lebesgue decomposition of Z(t) and the set of discontinuity points of the function Z(t) has no limit points in J.

Definition 2.
With C * a , a ∈ R we denote the Banach space of right continuous vector functions the subspace of all absolutely continuous functions by The set of jump points of every initial vector function Φ ∈ C * a we denote by S Φ . We emphasize that the set S Φ ∩ K is finite for every compact interval K ⊂ R and the case S Φ = ∅ is not excluded.
For the system (1) or (2) introduce the following initial conditions: We say that for the kernel U : J a × R → R n×n the conditions (S) hold, if the following conditions are fulfilled: (S1) The functions (t, θ) → U(t, θ) are measurable in (t, θ) ∈ J a × R and normalized so that for t ∈ J a , U(t, θ) = 0 when θ ∈R + and U(t, θ) = U(t, −h) for all θ ∈ (−∞, −h], h ∈ R + and Var θ∈[−h,0] U(t, ·) < ∞ for t ∈ J a . (S2) The Lebesgue decomposition of the kernel U(t, θ) for t ∈ J a and θ ∈ [−h, 0] has the form: bounded on J a and U S (t, θ) ∈ C(J a × R, R n×n ). (S3) There exists a locally bounded function z u ∈ L loc for t ∈ J a and for every t * ∈ J a the following relations hold: lim Consider the following auxiliary system in matrix form where I −1 (Γ(α)) = diag(Γ −1 (α 1 ), . . . , Γ −1 (α n )), or for k ∈ n in more detailed form with the initial condition (3).
We will need a slightly modified version of the Weissinger generalization of the Banach's fixed point theorem for complete metric spaces (see [22], Fixpunktsatz, p. 195).

Theorem 1.
Let Ω be a complete metric space with metric d Ω and let the following conditions hold:

2.
The operator T : Ω → Ω satisfies for each q ∈ N and for arbitrary x, y ∈ Ω the inequality Then T has a uniquely fixed point x * ∈ Ω(x * = Tx * ) and for every x ∈ Ω we have that lim q→∞ T q x = x * . Remark 1. This modification of the Weissinger generalization of the Banachs fixed point is not new. It is used in [23] and in the case when Ω is a Banach space in [24]. It is clear that the original Weissinger proof is correct for the presented variant too, with elementary modifications.
Let B be an arbitrary real Banach space.

Definition 4 ([25]
). The function f (t) : R → B is called a regulated function if it has one-sided (left and right) limits at every point t ∈ R.
Remark 2. If f (t) : J K → B, where J K ⊂ R is an arbitrary compact interval and f (t) is a regulated function, then it is assumed that in the left (right) end on the interval J K the function f (t) has only a right (left) limit. Theorem 2 ([25]). Let f (t) : R → B be an arbitrary function.
Then a necessary and sufficient condition for f (t) to be a regulated function is that f (t) must in every compact interval J K ⊂ R be a limit of n uniformly convergent sequence of step-functions (i.e., with respect to the supremum norm · ∞ = sup t∈J K · B ). Theorem 3 ([26]). Let t 0 ∈ J a and α > 0 be arbitrary fixed numbers and the following conditions hold: 1.
Then for t ∈ J a the following inequality holds

Remark 3.
Note that the statement of Theorem 3 is proved in the partial case t 0 = 0, but with small modifications the proof will be correct for arbitrary t 0 ∈ J a .

Existence and Uniqueness of the Solutions
Let it be that for every Φ ∈ C * a consider the corresponding linear space where Φ ∈ C * a is arbitrary and define a metric function d M : M (a) then for every M ∈ R + , according to a well-known result we conclude that E Φ M is a complete metric space concerning the metric It was a very strange for us that we could not find a result from which the statement of the next lemma directly follows. Proof. Let M ∈ R + be an arbitrary fixed number and consider an arbitrary Cauchy Let ε > 0 and k ∈ N be arbitrary numbers. There exist n 0 = n 0 (ε) ∈ N and δ = δ(ε, k) > 0 such that for every n ≥ n 0 we have that |G n (t) Thus G 0 ∈ AC([a, a + M]) for every M ∈ R + and the statement is proved.

Remark 4.
It must be noted that the statement of the next theorem cannot be obtained directly as a corollary from analogical results for the considered initial problems in [15,18,23]. Conditions (S) hold.

2.
The function Φ ∈ C * a is arbitrary. Then the IP (1) and (3) has a unique absolutely continuous solution in the interval [a, a + M] for every M ≥ max{h, 1}.
Proof. In the proof of this theorem we will use the approach introduced in [23].
Let Φ ∈ C * a and M ≥ max{h, 1} be arbitrary. From (5) in the case f k (t) ≡ 0 for k ∈ n we obtain the system For every G(t) = (g 1 (t), . . . , g n (t)) T ∈ E Φ M define the operator ( G)(t) = ( 1 g 1 (t), . . . , n g n (t)) T via the operators k g k (t) for t ∈ (a, a + M], k ∈ n by and the following additional condition: Since Φ is PAC, the conditions (S) imply that for each k, j ∈ n , the functions (see [27], Lemma 1). Then the right side of (7) is absolutely continuous on the interval [a, a + M], which implies that the function k g k (t) is absolutely continuous on the same interval and since lim M is a full metric space concerning the metric d M , it is enough to check that for the operator the conditions of Theorem 1 hold and then in virtue of Theorem 1 we will obtain that the operator has a unique fixed point in E Φ M . We recall that the Γ(z), z ∈ R + , has a local minimum at z min ≈ 1.46163, where it attains the value Γ(z min ) ≈ 0.885603. There exists q 0 ∈ N such that q 0 + 1 > α −1 m ≥ q 0 and for every q ∈ q 0 + 1 we denote with α q , α q ∈ {α 1 , . . . , α n } that number for which θ)) and let G,Ḡ ∈ E Φ M be arbitrary.
Then from (7) for k ∈ n and every t ∈ [a, a + M] we obtain Let assume that for some q ∈ N, k ∈ n and for every t ∈ [a, a + M] the inequality holds. Obviously from (9) it follows that the inequality (10) holds for each k ∈ n and every t ∈ [a, a + M] at least for q = 1. Denoting for simplicity q G Let us assume that (10) holds for each k ∈ n , every t ∈ [a, a + M] and for some q ∈ N. Then a similar way as in (9) from (10) and (11) we obtain Substituting s − a = z(t − a) in the integral in the right side of (12) and using the well-known relation between the gamma and beta functions we obtain and hence (10) holds for every q ∈ N , for each k ∈ n and every t ∈ [a, a + M].
For q ∈ q 0 + 1 from (10) and (13) it follows that and denote γ q = n(MP) q Γ(1+α q q) for q ∈ q 0 + 1 . For all q > q 0 + 1 from (10) and (13) we obtain that and denote γ q = n(MP) q Γ(1+α m q) for q > q 0 + 1. Consider the one parameter Mittag-Leffler function It is simply to see that the series ) < ∞ and then from Theorem 1 it follows that the IP (1) and (3)  Proof. Let Φ ∈ C * a be arbitrary and assume the contrary, that there exists M max < ∞ such that the solution X max (t) in the interval t ∈ [a, a + M max ] does not possess a prolongation. Let M * > M max be an arbitrary number. In virtue of Theorem 4 the IP (1) and (3) has a unique absolutely continuous solution X * (t) in the interval t ∈ [a, a + M * ]. The solution X * (t) obviously is a prolongation of the solution X max (t), which contradicts of our assumption that M max < ∞.
For arbitrary fixed s > a consider the following auxiliary system with the following condition Corollary 2. Let the following conditions hold: 1. Conditions (S) hold.

2.
The functionΦ js has the form Φ js (t) = I j , t = s, j ∈ n 0, t < s Then for each s > a and arbitrary j ∈ n the problem (16) and (17) has a unique solution, which satisfies Equation (16) for t > s, the condition (17) for t ≤ s and is absolutely continuous in (−∞, s) ∪ (s, ∞) with a first kind jump at t = s. Proof. Let s > a be an arbitrary fixed number, introduce the system and consider the IP (17) and (18). Since Theorem 4 is proved for arbitrary a ∈ R andΦ js (t) for arbitrary j ∈ n is PAC on the interval [s − h, s], then from Theorem 4 it follows that the IP (17) and (18) possess a unique solution X j (·) ∈ AC(J s , R n ) and moreover, from (17) it follows also that X j (·) ∈ AC([a, s) ∪ J s , R n ). Then for every k ∈ n we have and hence X j (t) = (x j1 , . . . , x jn ) T satisfies the Equation (16) for t ∈ (s, ∞) and the condition (17) for t ∈ (−∞, s]. Let consider an IP with Equation (16) for t > a and initial condition X(t) = X j (t) = 0 for t ∈ [a − h, a]. Then obviously X j (t) = 0 is its unique solution in t ∈ [a, s). This completes the proof.

Fundamental Matrix
Let s ∈ J a be an arbitrary fixed number and define the following matrix valued functionΦ(t, s) = (φ kj (t, s)) n k,j=1 : R × J a → R n×n with Φ(t, s) = I, t = s Θ, t < s and denoteΦ j (t, s) = (φ 1j (t, s), . . . ,φ nj (t, s)) T , j ∈ n . For arbitrary fixed number s ∈ J a consider the following matrix IP Definition 5. The matrix valued function t → C(t, s) = (C 1 (t, s), . . . , C n (t, s)) = {c j k (t, s)} n k,j=1 , s ∈ J a , is called a solution of the IP (19) and (20) in J s if C(·, s) : [s, ∞) → R n×n is continuous for t ∈ [s, ∞) and satisfies the matrix Equation (19) on t ∈ [s, ∞) as well as the initial condition (20) too. (20) we need only the values ofΦ(·, s) for t ∈ [s − h, s], but for convenience we define C(t, s) = Θ also for t ∈ (−∞, s − h). Then C(t, s) is prolonged as continuous in t function on (−∞, s).

Theorem 5. Let the conditions (S) hold.
Then for every initial point s ∈ J a , the matrix IP (19) and (20) has a unique absolutely continuous solution t → C(t, s) in the interval J s .
Proof. The statement of the Theorem follows immediately from Corollary 2. (19) and (20), will be called fundamental (or Cauchy) matrix for the homogeneous system (1).

Lemma 2.
Let the conditions (S) hold and the matrix valued function t → C(t, s) is the fundamental matrix of the system (1).
Then for everyt ∈ J a the matrix function C(t, ·) : [a,t] → R n×n is locally bounded in s for s ∈ [a,t] and t ∈ (−∞,t] . Proof. Lett ∈ J a be an arbitrary fixed number, s ∈ [a,t] be arbitrary and consider the fundamental matrix C(t, ·) : [a,t] → R n×n . According to Remark 5 for s ∈ [a,t] we have that C(t, s) = Θ for t ∈ (−∞, a) and C(a, a) = I.
Taking into account Theorem 5 and Corollary 2 it is easy to be seen, that the unique solution C(t, s) of IP (19) and (20) is a solution of the equation Introduce for every k, j ∈ n the notations: From Theorem 3 and (22) for a ≤ t ≤ s ≤t we obtain the estimation and thus C(t, s) is locally bounded in s for s ∈ [a,t] and t ∈ (−∞,t].
Theorem 6. Let the conditions (S) hold and the matrix valued function t → C(t, s) be the fundamental matrix of system (1). Then for every fixedt ∈ J a the matrix function C(t, ·) : [a,t] → R n×n is continuous for s ∈ [a,t] when s = t, for s = t possess first kind jumps and hence is Lebesgue integrable in s on [a,t] for each t ∈ (−∞,t]. Proof. Lett ∈ J a be an arbitrary fixed number, s ∈ [a,t] be arbitrary and consider the fundamental matrix C(t, ·) : [a,t] → R n×n . According to Remark 5 for s ∈ (a,t] we have that C(t, s) = Θ for t ∈ (−∞, a] and C(a, a) = I. In virtue of Theorem 2 and Lemma 2 it is enough to prove that C(t, s) has left and right limits for each s ∈ (a,t) and there exist the limits C(t,t − 0) and C(a, a + 0).
Then from (26) and Theorem 3 it follows that for each j ∈ n , k ∈ n and s ∈ (a, s * ) we obtain the following estimation: From (27) and (28) it follows that lim The case when t > s * and s ∈ (s * , t) can be treated fully analogically to obtain that c j k (t, s * + 0) exists and c j k (t, s * + 0) = c j k (t, s * ).
Let s * ∈ [a − h, a] be an arbitrary fixed number and define the following matrix valued and consider following IP: Theorem 7. Let the conditions (S) hold andt ∈ J a be arbitrary. Then the following statements hold: 1.
The matrix IP (29) and (30) has a unique absolutely continuous solution T s * (t, s) in t for t ∈ J a for every s * ∈ [a − h, a].

2.
The matrix function T a−h (t, ·) : [a,t] → R n×n is continuous in s for each s ∈ [a,t] with s = t.

3.
When s = t with t ≥ a, T a−h possess first kind jumps and hence is Lebesgue integrable in s on [a,t].
Proof. 1. Let s * ∈ [a − h, a] be an arbitrary fixed number. Then since Φ * s * is PAC for each s * ∈ [a − h, a], then the statement of point 1 follows from Theorem 4.

2.
Let s ∈ [a,t] with s = t. Then if t < s we have that T a−h (t, s) = 0. Consider the case t > s. Then the same way as in the proof of point (iii) of Theorem 6 we obtain that T a−h (t, s) is continuous in s when s = t.

3.
Let s = t with t ≥ a. Then obviously lim

Applications
We will demonstrate that the obtained results concerning the fundamental matrix allow to improve the integral representation of the solution of the IP (2) and (3) and simplify the proofs.
As usual according the superposition principle we will seek a solution of IP (2) and (3) with initial condition Φ(t) ≡ 0, t ∈ [a − h, a] for the case when the function F ∈ L loc 1 (J a , R n ) is locally bounded. Let where K(t, s) = C(t, s)T(s), T(s) = RL D 1−α a+ F(s).
Theorem 8. Let the following conditions hold.

2.
The function F ∈ L loc 1 (J a , R n ) is locally bounded and D α a+ F(t) ∈ L loc 1 (J a , R n ). Then the vector function X F (t) defined by equality (31) is a solution of IP (2) and Proof. From (31) and Theorem 5 it follows that X F (t) is an absolutely continuous function in J a . Then we have that Taking into account (31) and (32), Lemma 2.5 in [1], Condition 2 of the Theorem and applying Fubini's theorem we obtain For the right side of (33) we have and then the statement of the corollary follows from (33) and (34).

Corollary 3.
Let the following conditions hold.

2.
The function F ∈ L loc 1 (J a , R n ) is locally bounded and F(a) = 0. Then the vector function X F (t) defined by equality (31) is a solution of IP (2) and (3) with initial condition X F (t) = Φ(t) ≡ 0, t ∈ [a − h, a].

Proof. For
Then the proof is the same as in Theorem 8.
Let T a−h (t, s) be a solution of the IP (29) and (30) for s * ∈ [a − h, a]. For arbitrary function Φ(t) ∈ BV([a − h, a], R n ) define the following function: Theorem 9. Let the following conditions hold.

2.
The initial function Φ(t) ∈ BV([a − h, a], R n ) is not a constant and has finitely many jumps.
Then the function X Φ (t) defined with (35) is the unique solution of IP (1) and (3).
and hence X Φ (t) satisfies the initial condition (3). Theorem 7 implies that the matrix valued function T a−h (t, s) is an absolutely continuous function for t ∈ J a and then the vector valued X Φ (t) is defined by equality (35) too. Then we have and hence applying the Fubini's theorem we have that In the right side of (1) for each i ∈ m applying the unsymmetric Fubini theorem [28] we obtain that Theorem 10. Let the following conditions hold.
The function F ∈ L loc 1 (J a , R n ) is locally bounded.
Proof. The statement of the theorem follows immediately from the superposition principle, Corollary 3 and Theorems 8 and 9.
Corollary 4. Let the following conditions hold.

1.
The conditions of Theorem 10 hold.

2.
The Lebesgue decomposition of the function Φ(t) ∈ BV([a − h, a], R n ) does not include a singular term.
Then the function X Φ F (t) defined by (38) possesses the following integral representation: where Φ(t) = Φ J (t) + Φ A (t) and Φ J (t), Φ A (t) are the jump term and the absolutely continuous term, respectively, in its Lebesgue decomposition and the summation is taken over all jump points s i ∈ S Φ .
Proof. The statement of the corollary follows immediately from Theorem 10.
One of the important questions that arise when we use a fundamental matrix is what kind are its analytical properties concerning the variable s. It is well known that for the integer-order linear differential systems without delay, this problem is fully solved, i.e., both variables are symmetric. Generally speaking, this is not true even in the integer case for the delayed differential systems (the symmetry disappears for the non-autonomous systems and in some cases for the autonomous systems too). As far as we know there are no results in this direction for delayed non-autonomous fractional differential systems. The next result is a first attempt to establish some analytical properties of the fundamental matrix in the mentioned case concerning the variable s. Theorem 11. Let the conditions of Theorem 9 hold and the matrix valued function t → C(t, s) is the fundamental matrix of the system (1).
(i) When t ∈ (−∞, s) then for the fundamental matrix C(t, ·) : [a,ā] → R n×n according to Remark 5 for arbitrary s ∈ (a,ā) we have that C(t, s) = Θ and C(s, s) = I. Thus C(t, s) has jumps of first kind for t = s and C(t, s) is absolutely continuous in s for s ∈ (a,ā).
(ii) Let s ∈ (a,ā) andt ∈ (ā, ∞) be an arbitrary fixed number. For purposes of clarity we assume that s, s * ∈ (a,ā) with s < s * . Then as in the proof of Theorem 6 we obtain that (25) holds for t ∈ [a,ā] and introduce the same notations used there.
Since the function p(η, s, s * ) = bC j (η−s) α k −(η−s * ) α k α k is monotonically decreasing in η for η ∈ [s * ,ā] then in a similar way as in the proof of Theorem 6 for each k, j ∈ n from (25) we obtain that when s * − s ≤ 1. Then we can conclude that c j k (t, s) for each k, j ∈ n is α m -Hoelder continuous in s in every closed subinterval [a 1 , a 2 ] ⊂ (a,ā) with a 2 − a 1 ≤ 1. Thus c j k (t, s) is absolutely continuous in s for s ∈ (a,ā). (iii) Let s ∈ (a,ā),t ∈ (s,ā] and s < s * <t. Then as in the former case we conclude that (40) holds and then c j k (t, s) for each k, j ∈ n is α m -Hoelder continuous in s for every closed subinterval [a 1 , a 2 ] ⊂ (a,ā) for which a 2 − a 1 ≤ 1. Thus c j k (t, s) is absolutely continuous in s for every compact subinterval [a 1 , a 2 ] ⊂ (a,ā). Remark 6. It is not difficult to see that under the conditions of Theorem 11 the statement of the theorem holds for T a−h (t, s) too.
Corollary 5. Let the following conditions hold.
The Lebesgue decomposition of the function Φ(t) ∈ BV([a − h, a], R n ) does not include a singular term. 3.
The delays σ i (t) ∈ C 1 (J a ,R + ), σ i (t) < 1 for t ∈ J a and i ∈ m .
Then the unique solution X Φ (t) of IP (1) and (3) defined with (35) has the following representation: Letā ∈ J a be an arbitrary fixed number and [a 1 , a 2 ] ⊂ (a,ā) be an arbitrary compact subinterval. Condition 3 of the theorem implies that for arbitrary fixed number t ∈ R the set S σ t = {s ∈ [a 1 , a 2 ]|s = t − σ i (t), i ∈ m 0 } is finite.
Then the statement of the corollary follows from (42) and (43).

Conclusions
In this article, first the existence and uniqueness of the solutions of an initial problem for linear differential systems with incommensurate order Caputo fractional derivatives and with piecewise absolutely continuous (PAC) initial function is proved.
Then we prove the existence and uniqueness of an absolutely continuous fundamental matrix C(t, s), which has the following properties: 1. C(t, s) is absolutely continuous in t for PAC initial functions; 2. C(t, s) is absolutely continuous in s (with appropriate additional assumptions).
It must be noted that when the fundamental matrix is absolutely continuous in t and in s, the fundamental matrix has integrable derivatives in t and in s and this allows simpler and more applicable formulas to be obtained in the integral representations, as well as simpler and shorter proofs.
As far we know there are no other articles where such properties of the fundamental matrix C(t, s) concerning the variable s for delayed non-autonomous fractional differential systems are obtained. A brief comparison with similar fundamental matrix studies shows that the same system was studied in [15], but there is proof of existence of a continuous fundamental matrix, which is only continuous in t for initial functions with bounded variation. Our result is more general than that obtained in [21] where the smoothness of the fundamental matrix is proven.
Finally, using the properties of the fundamental matrix thus obtained, integral representations are obtained in the paper for the particular solution of the inhomogeneous system with zero initial conditions and for the general solution of the homogeneous system. A comparison, for example, with the integral representations obtained in [15] shows, that all the proofs are shorter and the obtained formula for the general solution of the homogeneous system is simpler and more applicable.
A general comparison with the analogous results for integer order derivatives shows that those obtained in the article results coincide with them at α = 1, which means that they are a generalization of the classical ones.
We hope that the results obtained will be useful both for future research and generalizations from a mathematical point of view, as well as for modeling of real-world phenomena.

Conflicts of Interest:
The authors declare no conflict of interest.