Next Article in Journal
Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework
Next Article in Special Issue
Mann-Type Inertial Subgradient Extragradient Rules for Variational Inequalities and Common Fixed Points of Nonexpansive and Quasi-Nonexpansive Mappings
Previous Article in Journal
Banach Contraction Principle and Meir–Keeler Type of Fixed Point Theorems for Pre-Metric Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method

1
Department of Computer Programing, Baskale Vocational School, Van Yuzuncu Yil University, 65080 Van, Turkey
2
Department of Electric-Electronic Engineering, Faculty of Engineering, Van Yuzuncu Yil University, 65080 Van, Turkey
3
Department of Mathematics, Faculty of Sciences, Van Yuzuncu Yil University, 65080 Van, Turkey
4
Research Center for Interneural Computing, China Medical University Hospital, China Medical University, Taichung 406, Taiwan
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(2), 58; https://doi.org/10.3390/axioms10020058
Submission received: 19 March 2021 / Revised: 4 April 2021 / Accepted: 6 April 2021 / Published: 9 April 2021
(This article belongs to the Special Issue Special Issue in Honor of the 60th Birthday of Professor Hong-Kun Xu)

Abstract

:
The purpose of this paper is to investigate some qualitative properties of solutions of nonlinear fractional retarded Volterra integro-differential equations (FrRIDEs) with Caputo fractional derivatives. These properties include uniform stability, asymptotic stability, Mittag–Leffer stability and boundedness. The presented results are proved by defining an appropriate Lyapunov function and applying the Lyapunov–Razumikhin method (LRM). Hence, some results that are available in the literature are improved for the FrRIDEs and obtained under weaker conditions via the advantage of the LRM. In order to illustrate the results, two examples are provided.

1. Introduction

In recent years, a large number of books [1,2,3] and papers [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24] have been devoted to the study of various qualitative properties of solutions of scalars and systems of linear and nonlinear Volterra integro-differential equations (IDEs) both without and with delay, and that of some other kinds of differential equations due to their important applications in population growth models, mathematical models of biological species living together, mathematical models in physics, control engineering and signal processing, mathematical models of heat transfer and radiation, standard closed electric RLC circuits, and so on.
In the relevant literature three methods, which are called the second Lyapunov method, Lyapunov–Krasovskiĭ method and Lyapunov–Razumikhin method, come to the forefront to investigate qualitative properties of solutions of linear and nonlinear integro-differential equations both without and with retardation. Among these methods, the second Lyapunov method and Lyapunov–Krasovskiĭ method are extensively used to study various qualitative behaviors of solutions of integro-differential equations of integer order (see, [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]). To the best of our knowledge, the Lyapunov–Razumikhin method is less used during that kind of investigation [23,25,26]. However, when it is used for the appropriate problems, it is more effective than the other two methods mentioned, the second Lyapunov method and Lyapunov–Krasovskiĭ method. To the best of our knowledge from the relevant literature, the disadvantages of the Lyapunov second method and Lyapunov–Krasovskiĭ method are that both of these methods require the construction or definition of suitable Lyapunov function(s) and Lyapunov–Krasovskiĭ functional(s), which can include double integrals and additional terms. The construction of suitable Lyapunov function(s) and Lyapunov–Krasovskiĭ functional(s) for nonlinear functional differential equations remains an open problem in the literature at this time. This case is known as a disadvantage. In addition, the time derivatives of double integrals leads to stronger conditions for the negative or negative -semi definite of the time derivative(s) of function(s) or functional(s) used as basic tool(s) in the proof(s).
From this point of view, we would like to present the related work of Du [27]. Indeed, in 1995, the author investigated the uniformly asymptotic stability of trivial solutions of the system of nonlinear RIDEs of the form:
x ˙ ( t ) = f ( t , x ( t ) ) + g ( t , x ( t τ ) ) + t τ t h ( t , s , x ( s ) ) d s
or its equivalent system
x ˙ i ( t ) = f i ( t , x ( t ) ) + g i ( t , x ( t τ ) ) + t τ t h i ( t , s , x ( s ) ) d s , ( i = 1 , 2 , . . . , n ) .
In this paper, we consider the following initial value problem (IVP) for the system of nonlinear fractional retarded Volterra integro-differential equations (FrRIDEs) with Caputo derivative:
t o C D t q x ( t ) = f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + t τ t h ( t , s , x ( s ) ) d s + t ρ t p ( t , s , x ( s ) ) d s + q ( t , x ( t ) , x ( t τ ) , x ( t ρ ) ) , 0 < q < 1 ,
x ( t 0 + θ ) = ϕ ( θ ) , x ( t 0 ) = ϕ ( 0 ) = x 0 , θ [ τ , 0 ] [ ρ , 0 ] ,
where x = ( x 1 , . . . , x n ) T R n , t R , s [ τ , ) [ ρ , ) , τ and ρ are positive constants, i.e., they are constant retardations, f = ( f 1 , , f n ) T C ( R × R n , R n ) , f i ( t , x ( t ) ) = f i ( t , x 1 ( t ) , , x n ( t ) ) , g = ( g 1 , , g n ) T C ( R × R n × C H , R n ) , C H = ϕ : ϕ C and ϕ t 0 H < , g i ( t , x ( t ) , x ( t τ ) ) = g i ( t , x 1 ( t ) , , x n ( t ) , x 1 ( t τ ) , , x n ( t τ ) ) , h C ( R × [ τ , ) × C H , R n ) , p C ( R × [ ρ , ) × C H , R n ) and q C ( R × R n × C H × C H , R n ) . It is supposed that f ( t , 0 ) = 0 , g ( t , 0 , 0 ) = 0 , p ( t , s , 0 ) = 0 , and h ( t , s , 0 ) = 0 . Then, the system of Volterra FrRIDEs (2) with a Caputo derivative includes the zero solution, when q ( . ) 0 .
In this article, motivated by the system of nonlinear RIDEs at Equation (1), i.e., the result of Du [27] (Theorem 4), and those in the bibliography of this paper, we consider the system of nonlinear FrRIDEs at Equation (2) with a Caputo derivative. As indicated above, we plan to investigate the uniformly stability, asymptotic stability, and Mittag–Leffler stability of the zero solution of Equation (2) with q 0 , and the boundedness of all solutions of Equation (2) with q 0 , by using the Razumikhin method (see [25,26,28,29,30]). It should be noted that the Caputo derivative is applicable to continuously differentiable quadratic Lyapunov functions to study qualitative properties of solutions of fractional differential equations and fractional delay differential equations, etc. (see, for example, [25,26,31,32,33,34,35,36,37]).
It is known that the presence of the fractional derivatives in the system requires that we use appropriately defined fractional derivatives of Lyapunov functions. In the literature, four types of fractional derivatives are commonly applied to calculate the derivatives of Lyapunov functions; these are the Caputo fractional derivative, the Caputo fractional Dini derivative, the Riemann–Liouville fractional derivative, and the Grünwald–Letnikov fractional derivative [32,36]. Not all of these will be employed here. The results presented below are new contributions to the literature on delay fractional integro-differential equations with Caputo derivatives.

2. Preliminaries

We begin by considering a system of fractional retarded differential equations (FrRDEs) with a Caputo derivative of order q ( 0 , 1 ) :
t o C D t q x ( t ) = F ( t , x t ) , t J , J = [ t 0 τ , T ) , T + , 0 t 0 t ,
where x R n , F ( t , ϕ ) J × C ( [ r , 0 ] , R n ) , F ( t , 0 ) = 0 , x ( t 0 + s ) = ϕ ( s ) for s [ r , 0 ] , x ( t 0 + ) = ϕ ( 0 ) , ϕ C ( [ r , 0 ] , R n ) , r > 0 is the constant retardation. For ϕ C ( [ r , 0 ] , R n ) , we use the usual Euclidean norms . and . t 0 defined by
x t = sup r s 0 x ( t + s ) and ϕ t 0 = sup t 0 r s t 0 ϕ ( s ) ,
respectively.
Since the function F is continuous, for any initial data ( t 0 , ϕ ) R + × C ( [ r , 0 ] , R n ) , the initial value problem for the system of FrRDEs in Equation (4) has at the least one solution x ( t ) = x ( t , t 0 , ϕ ) C 1 ( [ t 0 , ) , R n ) . If the function F satisfies a Lipschitz condition in x, then the solution is unique.
The following lemmas and other concepts are needed in the remainder of this paper.
Firstly, we give Lemma 1, which is a consequence of (Theorem 2) [31].
Lemma 1.
Assume that for any initial data x ( t 0 , ϕ 0 ) R + × C ( [ τ , 0 ] , R n ) , the system of FrRIDEs in Equation (2) has a solution. If there exists a Lyapunov function V and strictly increasing u , v C ( R + , R + ) with u ( 0 ) = v ( 0 ) = 0 and
u ( x ) V ( t , x ) v ( x ) for all t t 0 τ and all x R n
and such that for any initial data ( t 0 , ϕ 0 ) R + × C ( [ τ , 0 ] , R n ) and any point s > t 0 with
V ( s + ξ , x ( s + ξ ) ) < V ( s , x ( s ) ) for all ξ [ τ , 0 ) ,
the inequality
t 0 C D t q V ( t , x ( t ) ) 0 for all t ( t 0 , s ]
holds, then the zero function of Equation (2) with a zero initial condition is uniformly stable.
Lemma 2
([31]). The zero solution of the FrRDEs in Equation (4) is asymptotically stable if there exist a continuous function V ( t , x ) , continuous increasing and positive definite functions u , v , ω and a continuous non-decreasing function p ( s ) > s for s > 0 such that the following conditions hold for all t J :
V ( t , 0 ) 0 , u ( x ) V ( t , x ) v ( x ) for all t J and all x R n ; t o C D t q V ( t , x ( t ) ) ω ( | x ( t ) | ) for all t ( t 0 , s ] ;
and
V ( t + s , x ( t + s ) ) < p V ( t , x ( t ) ) for all s [ τ , 0 ] .
Lemma 3
([38] Lemma 1). Let x ( t ) R n be a vector of differentiable functions. Then for any t t 0 ,
1 2 t o C D t q ( x T x ) x T ( t ) t o C D t q x ( t ) for all q ( 0 , 1 ] ,
Lemma 4
([38] Lemma 4). Let x ( t ) R n be a vector of differentiable functions. Then, for any t t 0 ,
1 2 t o C D t q ( x T P x ) x T ( t ) P t o C D t q x ( t ) for all q ( 0 , 1 ] ,
where P R n × n is a constant, symmetric and positive definite matrix.
Definition 1
([35] Definiton 3.1). The trivial solution of the system of FrRIDEs in Equation (2) is said to be Mittag–Leffler stable provided the solution x ( . , ϕ ) of (2) satisfies
x ( t , ϕ ) [ m ( ϕ ) E q ( λ ( t t 0 ) q ) ] b ,
where q ( 0 , 1 ) , λ 0 , b > 0 ,
ϕ = max θ [ τ , 0 ] ϕ ( θ ) ,
m ( 0 ) = 0 , m is a locally Lipschitz function and is non-negative, and
E q ( z ) = k = 0 z k Γ ( q k + 1 )
is the one-parameter Mittag–Leffler function, and Γ denotes the Gamma function.
Lemma 5
([35] Lemma 2.1). Let x R n be a vector of differentiable functions. If a continuous function V : [ t 0 , ) × R n R + satisfies
t o C D t q V ( t , x ( t ) ) α V ( t , x ( t ) ) ,
then
V ( t , x ( t ) ) V ( t 0 , x ( t 0 ) ) E q ( α ( t t 0 ) q ) ,
where α > 0 and 0 < q < 1 .
Lemma 6
([33] Property 1).
t o C D t q ( a x ( t ) + b y ( t ) ) = a t o C D t q x ( t ) + b t o C D t q y ( t ) ,
where q ( 0 , 1 ] .
The contents of the next lemma are well known.
Lemma 7.
Let x R n , n N , n 1 , and M R n × n be a positive definite symmetric n × n - matrix such that
λ M λ i ( M ) λ m , ( i = 1 , 2 , . . . , n ) ,
where λ i ( M ) denotes the eigenvalues of M. Then
λ M x 2 M x , x λ m x 2 ,
where λ M and λ m are the greatest and least eigenvalues of the matrix M, respectively.
We know that λ M and λ m are real and positive since M is a positive definite symmetric matrix.

3. Razumikhin Analyses of Solutions

In the system of Volterra FrRIDEs in Equation (2), let q 0 , i.e., we consider the system in Equation (2) with Equation (3) replaced by
t o C D t q x ( t ) = f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + t τ t h ( t , s , x ( s ) ) d s + t ρ t p ( t , s , x ( s ) ) d s .
We will use the following hypotheses in our main results.
Hypothesis 1.
f ( t , 0 ) = g ( s , 0 , 0 ) = h ( t , s , 0 ) 0 , x i f i ( t , x ) > 0 as x i 0 , for t R + , all x R n ;
Hypothesis 2.
The functions H and P satisfy the local Lipschitz condition in x, with
H ( t , s , x ) t τ t h ( t , s , x ( s ) ) d s , h ( t , s , 0 ) = 0 , h ( t , s , x ( s ) ) h 0 x for s t
and
P ( t , s , x ) t ρ t p ( t , s , x ( s ) ) d s , p ( t , s , 0 ) = 0 , p ( t , s , x ( s ) ) p 0 x for s t ,
where h 0 > 0 , p 0 > 0 , h 0 , p 0 R ;
Hypothesis 3.
f ( t , x ) g ( t , x , x ( t τ ) ) ( τ h 0 + ρ p 0 ) x 0 f o r t R +
and all x R n , x ( t τ ) C H ;
Hypothesis 4.
f ( t , x ) g ( t , x , x ( t τ ) ) ( τ h 0 + ρ p 0 ) x ρ 1 x for t R +
and all x R n , x ( t τ ) C H , where ρ 1 > 0 , ρ 1 R ;
Hypothesis 5.
There exists q 0 C ( R + , R ) such that
q ( t , x , x ( t τ ) , x ( t ρ ) ) q 0 ( t ) x for all t t 0 , x R n , x ( t τ ) , x ( t ρ ) C H ,
and
f ( t , x ) g ( t , x , x ( t τ ) ) ( q 0 ( t ) + τ h 0 + ρ p 0 ) x 0 for t R +
and all x R n , x ( t τ ) C H .
Theorem 1.
The zero solution of the system of FrRIDEs in Equation (5) with Caputo derivative is uniformly stable if the conditions of Hypotheses 1–3 hold.
Proof. 
We define a Lyapunov function W : = W ( t , x ) = W ( t , x ( t ) ) by
W ( t , x ) : = x = i = 1 n x i = x 1 + + x n .
For arbitrary initial data ( t 0 , ϕ ) R + × C ( [ τ , 0 ] [ ρ , 0 ] , R n ) and a point t > t 0 , it follows that W ( t , x ) satisfies the Razumikhin condition (see [28,29,30])
W ( t , x ( t ) ) > W ( t + s , x ( t + s ) )
on the initial set [ τ , 0 ] [ ρ , 0 ] , i.e.,
x ( t ) > x ( t + s ) for all s [ τ , 0 ] [ ρ , 0 ] .
Let x ( t ) = x ( t , t 0 , ϕ ) denote the solution of the IVP of Equation (5) such that x ( t 0 + + s ) = ϕ ( s ) for s [ τ , 0 ] [ ρ , 0 ] . From this point, it is clear that W ( t , x ) in Equation (6) satisfies the relations
W ( t , 0 ) = 0 , 1 2 x 1 + . . . + 1 2 x n = 1 2 x W ( t , x ) ,
and
W ( t , x ) 5 4 x 1 + . . . + 5 4 x n = 5 4 x .
Taking the Caputo fractional derivative of the Lyapunov function W ( t , x ) in Equation (6) along the system of FrRIDEs in Equation (5), making use the conditions of Hypotheses 1 and 2 and some elementary calculations, we obtain
t o C D t q W ( t , x ( t ) ) = t o C D t q ( x 1 ( t ) + x 2 ( t ) + . . . + x n ( t ) ) = t o C D t q x 1 ( t ) + t o C D t q x 2 ( t ) + + t o C D t q x n ( t ) = ( s i g n x 1 ( t ) ) t o C D t q x 1 ( t ) + ( s i g n x 2 ( t ) ) t o C D t q x 2 ( t ) + . . . + ( s i g n x n ( t ) ) t o C D t q x n ( t ) = i = 1 n x i ( t ) t o C D t q x i ( t ) = i = 1 n x i ( t ) [ f i ( t , x ( t ) ) + g i ( t , x ( t ) , x ( t τ ) ) ] + i = 1 n x i ( t ) [ t τ t h i ( t , s , x ( s ) ) d s + t ρ t p i ( t , s , x ( s ) ) d s ] i = 1 n [ f i ( t , x ( t ) ) + g i ( t , x ( t ) , x ( t τ ) ) ] + i = 1 n [ t τ t h i ( t , s , x ( s ) ) d s + t ρ t p i ( t , s , x ( s ) ) d s ] = f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + t τ t h ( t , s , x ( s ) ) d s + t ρ t p ( t , s , x ( s ) ) d s f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + h 0 t τ t x ( s ) d s + p 0 t ρ t x ( s ) d s .
Consider the integral terms such that
h 0 t τ t x ( s ) d s + p 0 t ρ t x ( s ) d s ,
which are included in the inequality of Equation (7).
Letting s t = ξ gives d s = d ξ . Hence, for s = t τ and s = t , it follows that ξ = τ and ξ = 0 , respectively. Similarly, by the same transformation and way, for s = t ρ and s = t , we have ξ = ρ and ξ = 0 , respectively. In view of these estimates, using the Razumikhin condition on the set s [ τ , 0 ] [ ρ , 0 ] , we get
h 0 t τ t x ( s ) d s + p 0 t ρ t x ( s ) d s = h 0 τ 0 x ( t + ξ ) d ξ + p 0 ρ 0 x ( t + ξ ) d ξ < h 0 τ 0 x ( t ) d ξ + p 0 ρ 0 x ( t ) d ξ = h 0 x ( t ) τ 0 d ξ + p 0 x ( t ) ρ 0 d ξ = h 0 τ x ( t ) + p 0 ρ x ( t ) .
Then, from Equations (7) and (8), it follows that
t o C D t q W ( t , x ( t ) ) f ( t , x ( t ) ) g ( t , x ( t ) , x ( t τ ) ) ( τ h 0 + ρ p 0 ) x ( t ) 0 ,
that is, using the condition (H3), we have
t o C D t q W ( t , x ( t ) ) 0 .
Thus, from Lemma 1, the zero solution of the system of FrRIDEs in Equation (5) is uniformly stable. □
Our next result deals with the asymptotic stability of the system in Equation (5).
Theorem 2.
The zero solution of the system of FrRIDEs in Equation (5) is asymptotically stable if the conditions of Hypotheses 1, 2 and 4 hold.
Proof. 
With W ( t , x ) defined as in Equation (6), from the conditions (H1), (H2), and (H4) we easily conclude that
1 2 x W ( t , x ) 5 4 x
and
t o C D t q W ( t , x ( t ) ) ρ 1 x ( t ) .
Hence, the zero solution of the system of FrRIDEs in Equation (5) is asymptotically stable by Lemma 2. □
The following theorem shows the Mittag–Leffler stability of the system FrRIDEs in Equation (5).
Theorem 3.
The zero solution of the system of FrRIDEs in Equation (5) is Mittag–Leffler stable if the conditions of Hypotheses 1, 2 and 4 hold
Proof. 
Again with the Lyapunov function W ( t , x ) defined as in Equation (6), from the conditions of Hypotheses 1, 2 and 4, it is clear that
t o C D t q W ( t , x ( t ) ) ρ 1 x ( t ) = ρ 1 W ( t , x ( t ) )
holds.
Using Lemma 5, we obtain
x ( t ) = W ( t , x ( t ) ) W ( t 0 , x ( t 0 ) ) E q ( ρ 1 ( t t 0 ) q ) = x ( t 0 ) E q ( ρ 1 ( t t 0 ) q ) = m ( x ( t 0 ) ) E q ( ρ 1 ( t t 0 ) q ) m ( ϕ ) E q ( ρ 1 ( t t 0 ) q )
with m ( x ) = x ( t ) , which is locally Lipschitz. Thus, the proof of Theorem 3 is completed by using Definition 1. □

4. Boundedness of Solutions of System in Equation (2)

We now turn our attention to the perturbed system in Equation (2).
Theorem 4.
The solutions of the system of FrRIDEs in Equation (2) are bounded if the conditions of Hypotheses 1, 2 and 5 hold.
Proof. 
We again consider the Lyapunov function defined in Equation (6). Calculating the time derivative of the Lyapunov function W ( t , x ) along the system of FrRIDEs in Equation (2) and using the conditions in Hypotheses 1, 2 and 5, we obtain
t o C D t q W ( t , x ( t ) ) [ f ( t , x ( t ) ) g ( t , x ( t ) , x ( t τ ) ) ( τ h 0 + ρ p 0 ) ] x + q ( t , x , x ( t τ ) , x ( t ρ ) ) f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + ( q ( t ) + τ h 0 + ρ p 0 ) x 0 .
Hence, we have
W ( t , x ( t ) ) W ( t 0 , ϕ ( t 0 ) ) .
As a result of this inequality, it follows that
W ( t , x ( t ) ) = x ( t ) = x 1 ( t ) + . . . + x n ( t ) x ( t 0 ) = x 1 ( t 0 ) + . . . + x n ( t 0 ) = W ( t 0 , ϕ ( t 0 ) ) .
Let
K 0 = x ( t 0 ) = x 1 ( t 0 ) + + x n ( t 0 ) .
Hence, we obtain
x ( t ) = x 1 ( t ) + + x n ( t ) K 0 for t R + .
Hence, it is clear that if t , then x ( t ) K 0 . This inequality completes the proof of Theorem 4. □
Remark 1.
Here, if q = 1 , the boundedness of solutions as t was proved without using the Gronwall inequality, see Theorem 4. By this fact, we have removed some unnecessary conditions, and we can obtain some boundedness results in the literature under less restrictive conditions (see, for example, [12,13] and the bibliography therein). Here, we will not state the details of the discussions.
We now give the following example and solve the given system using MATLAB software. In fact, the problem was solved using the 4th order Runge–Kutta method in MATLAB. Here, the graphs of Figure 1, Figure 2, Figure 3 and Figure 4 show the behaviors of paths of the solutions of Example 1 for different values of fractional order q.
Example 1.
Consider the system of nonlinear Volterra IFrRDEs with Caputo derivative of order q ( 0 , 1 ) :
t o C D t q x 1 ( t ) t o C D t q x 2 ( t ) = 12 x 1 ( t ) + x 1 ( t ) 1 + t 2 + x 1 2 ( t ) 12 x 2 ( t ) + x 2 ( t ) 1 + t 2 + x 2 2 ( t ) + x 1 ( t ) 1 + t 2 + x 1 2 ( t 1 10 ) x 2 ( t ) 1 + t 2 + x 2 2 ( t 1 10 ) + t 1 10 t sin x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) sin x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s + t 1 5 t x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s ,
where t 1 10 , τ = 1 10 and ρ = 1 5 are the constant retardations and x ( t ) = x R 2 .
Comparing the system of IFrRDEs with Caputo derivative in Equation (11) with that given by Equation (5), we have the following formulas:
f ( t , x ) = f ( t , x 1 , x 2 ) = 12 x 1 + x 1 1 + t 2 + x 1 2 12 x 2 + x 2 1 + t 2 + x 2 2 ,
f ( t , 0 ) = f ( t , 0 , 0 ) = 0 ,
f 1 ( t , x 1 , x 2 ) = 12 x 1 + x 1 1 + t 2 + x 1 2 ,
x 1 f 1 ( t , x 1 , x 2 ) = 12 x 1 2 + x 1 2 1 + t 2 + x 1 2 > 0 ,
x 1 0 ,
f 2 ( t , x 1 , x 2 ) = 12 x 2 + x 2 1 + t 2 + x 2 2 ,
x 2 f 2 ( t , x 1 , x 2 ) = 12 x 2 2 + x 2 2 1 + t 2 + x 2 2 > 0 ,
x 2 0 ,
x 1 ( t + 0 ) f 1 ( t , x 1 , x 2 ) x 2 ( t + 0 ) f 2 ( t , x 1 , x 2 )
= 12 x 1 s g n x 1 ( t + 0 ) x 1 s g n x 1 ( t + 0 ) 1 + t 2 + x 1 2
12 x 2 s g n x 2 ( t + 0 ) x 2 s g n x 2 ( t + 0 ) 1 + t 2 + x 2 2
11 x 1 11 x 2 = 11 x f ( t , x ) ,
g ( t , x , x ( t τ ) ) = g ( t , x 1 , x 2 , x 1 ( t 1 10 ) , x 2 ( t 1 10 ) ) = x 1 1 + t 2 + x 1 2 ( t 1 10 ) x 2 1 + t 2 + x 2 2 ( t 1 10 ) ,
g ( t , x , x ( t τ ) ) = g ( t , x 1 , x 2 , x 1 ( t 1 10 ) , x 2 ( t 1 10 ) ) = x 1 1 + t 2 + x 1 2 ( t 1 10 ) x 2 1 + t 2 + x 2 2 ( t 1 10 )
x 1 1 + t 2 + x 1 2 ( t 1 10 ) + x 2 1 + t 2 + x 2 2 ( t 1 10 )
x 1 + x 2 = x ,
h ( t , s , x ) = h ( t , s , x 1 , x 2 ) = sin x 1 1 + t 2 + s 2 + x 1 2 sin x 2 1 + t 2 + s 2 + x 2 2 ,
h ( t , s , 0 ) = h ( t , s , 0 , 0 ) = 0 ,
h ( t , s , x ) = h ( t , s , x 1 , x 2 ) = sin x 1 1 + t 2 + s 2 + x 1 2 sin x 2 1 + t 2 + s 2 + x 2 2
= sin x 1 1 + t 2 + s 2 + x 1 2 + sin x 2 1 + t 2 + s 2 + x 2 2
sin x 1 + sin x 2 x 1 + x 2 = x ,
where h 0 = 1 ,
H ( t , s , x ) t τ t h ( t , s , x ( s ) ) d s = t 1 10 t sin x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) sin x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s
= t 1 10 t sin x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) d s
+ t 1 10 t sin x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s
t 1 10 t sin x 1 ( s ) d s + t 1 10 t sin x 2 ( s ) d s
t 1 10 t x 1 ( s ) d s + t 1 10 t x 2 ( s ) d s = t 1 10 t x ( s ) d s .
Let s t = ξ , which implies d s = d ξ . Then, for s = t 1 10 , we derive ξ = 1 10 , and similarly for s = t , we have ξ = 0 .
In view of these findings, using the given Razumikhin condition [28,29,30] on the initial segment [ 1 10 , 0 ] , it follows that
t 1 10 t x ( s ) d s = 1 10 0 x ( t + ξ ) d ξ < 1 10 0 x ( t ) d ξ = 1 10 x .
For this step, we consider the term:
P ( t , s , x ) t ρ t p ( t , s , x ( s ) ) d s
with
p ( t , s , x ) = p ( t , s , x 1 , x 2 ) = x 1 1 + t 2 + s 2 + x 1 2 x 2 1 + t 2 + s 2 + x 2 2 , ρ = 1 5 .
In that case, we derive
p ( t , s , 0 , 0 ) = 0
and
p ( t , s , x ) = p ( t , s , x 1 , x 2 ) = x 1 1 + t 2 + s 2 + x 1 2 x 2 1 + t 2 + s 2 + x 2 2
= x 1 1 + t 2 + s 2 + x 1 2 + x 2 1 + t 2 + s 2 + x 2 2 x 1 + x 2 = x , where p 0 = 1 .
For the next step, it follows that
P ( t , s , x ) t ρ t p ( t , s , x ( s ) ) d s = t 1 5 t x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s
t 1 5 t x 1 ( s ) d s + t 1 5 t x 2 ( s ) d s = t 1 5 t x ( s ) d s .
Let s t = ξ , which implies d s = d ξ . Then, for s = t 1 5 , we derive ξ = 1 5 . Similarly, for s = t , we have ξ = 0 . Then,
t 1 5 t x ( s ) d s = 1 5 0 x ( t + ξ ) d ξ < 1 5 0 x ( t ) d ξ = 1 5 x .
Hence, bringing together the above results, we derive
f ( t , x ) g ( t , x , x ( t τ ) ) ( τ h 0 + ρ p 0 ) x f ( t , x ) g ( t , x , x ( t 1 2 ) ) ( 1 10 + 1 5 ) x 11 x x 3 10 x = ( 9.7 ) x , where ρ 1 = 97 10 .
In the light of the above discussion, the conditions of Hypotheses 1–3 of Theorem 1, and the conditions of Hypotheses 1, 2 and 4 of Theorems 2 and 3 hold. For this reason, the zero solution of the system of FrRIDEs in Equation (11) with Caputo derivative is uniformly stable, asymptotically stable and Mittag–Leffler stable.
In Figure 1, Figure 2, Figure 3 and Figure 4, the system of FrRIDEs (11) was solved and the orbits of the solutions x 1 ( t ) , x 2 ( t ) were drawn for τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
For the case q ( . ) 0 , we now give the second example and solve it using MATLAB software.
Example 2.
Consider the system of nonlinear Volterra FrRIDEs with Caputo derivative of order q ( 0 , 1 ) :
t o C D t q x 1 ( t ) t o C D t q x 2 ( t ) = 12 x 1 ( t ) + x 1 ( t ) 1 + t 2 + x 1 2 ( t ) 12 x 2 ( t ) + x 2 ( t ) 1 + t 2 + x 2 2 ( t ) + x 1 ( t ) 1 + t 2 + x 1 2 ( t 1 10 ) x 2 ( t ) 1 + t 2 + x 2 2 ( t 1 10 ) + t 1 10 t sin x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) sin x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s + t 1 5 t x 1 ( s ) 1 + t 2 + s 2 + x 1 2 ( s ) x 2 ( s ) 1 + t 2 + s 2 + x 2 2 ( s ) d s + exp ( t ) x 1 1 + exp ( 2 t ) + x 1 ( t 1 10 ) + x 1 ( t 1 5 ) exp ( t ) x 2 1 + exp ( 2 t ) + x 2 ( t 1 10 ) + x 2 ( t 1 5 ) ,
where t 1 10 , τ = 1 10 and ρ = 1 5 are the constant delay terms and x ( t ) = x R 2 .
Comparing the systems of FrRIDEs in Equation (13) with Caputo derivative and Equation (2), we note that the functions f ( t , x ) , g ( t , x , x ( t 1 10 ) ) , h ( t , s , x ) and p ( t , s , x ) are the same as those in Example 1. Then, the satisfaction of the conditions of Hypotheses 1 and 2 have been shown in Example 1. For the verification of the condition of Hypothesis 5, we consider the last term of Equation (13):
q ( t , x , x ( t 1 10 ) , x ( t 1 5 ) ) = exp ( t ) x 1 1 + exp ( 2 t ) + x 1 ( t 1 10 ) + x 1 ( t 1 5 ) exp ( t ) x 2 1 + exp ( 2 t ) + x 2 ( t 1 10 ) + x 2 ( t 1 5 )
Clearly, it follows that
| q ( t , x , x ( t 1 10 ) , x ( t 1 5 ) ) | = exp ( t ) x 1 1 + exp ( 2 t ) + x 1 ( t 1 10 ) + x 1 ( t 1 5 ) exp ( t ) x 2 1 + exp ( 2 t ) + x 2 ( t 1 10 ) + x 2 ( t 1 5 ) = exp ( t ) x 1 1 + exp ( 2 t ) + x 1 ( t 1 10 ) + x 1 ( t 1 5 ) + exp ( t ) x 2 1 + exp ( 2 t ) + x 2 ( t 1 10 ) + x 2 ( t 1 5 ) exp ( t ) x 1 1 + exp ( 2 t ) + exp ( t ) x 2 1 + exp ( 2 t ) = exp ( t ) 1 + exp ( 2 t ) x 1 + x 2 = q 0 ( t ) x ,
where
q 0 ( t ) = exp ( t ) 1 + exp ( 2 t ) 1 2 , x 1 + x 2 = x .
In view of Equations (12) and (14), it is clear that
f ( t , x ) g ( t , x , x ( t τ ) ) ( q 0 ( t ) + τ h 0 + ρ p 0 ) x ( 9.7 ) x exp ( t ) 1 + exp ( 2 t ) x ( 9.7 ) x 1 2 x = ( 9.2 ) x .
As a consequence of this inequality, the condition of Hypothesis 5 holds. Thus, the solutions of the system of FrRIDEs in Equation (13) with Caputo derivative are bounded as t .
The following graphs of Figure 5, Figure 6, Figure 7 and Figure 8 show the behaviors of paths of the solutions of Example 2 for different values of fractional order q.

5. Discussions

We would like to explain the contributions of this paper to the relevant literature as the following.
(1)
To the best of our knowledge, in the literature, there are numerous papers on the uniform stability, asymptotic stability, Mittag–Leffer stability and boundedness of fractional differential equations of integer order both with and without delay. However, there are no papers in the literature on the asymptotic stability, Mittag–Leffer stability and boundedness of the FrRIDEs in Equation (2) with Caputo fractional derivative, except the two papers of Hristova and Tunç [25,26], which include some results on the uniform stability. Next, qualitative behaviors of the FrRIDEs in Equation (2) have not been discussed in the literature yet. Therefore, the results of this paper are new, original and they have scientific novelty.
(2)
If q = 1 in the FrRIDEs in Equation (2), then we have the system of RIDEs
x ˙ ( t ) = f ( t , x ( t ) ) + g ( t , x ( t ) , x ( t τ ) ) + t τ t h ( t , s , x ( s ) ) d s + t ρ t p ( t , s , x ( s ) ) d s + q ( t , x ( t ) , x ( t τ ) , x ( t ρ ) ) .
It is clear that the system of RIDEs in Equation (15) includes, extends and improves the system of RIDEs in Equation (1). This is a contribution to the topic and the relevant literature.
(3)
In Du [27] (Theorem 4), the uniform asymptotic stability of the zero solution of the system of RIDEs in Equation (1) was proved using the Lyapunov–Krasovskiĭ functional:
V ( t , x ( . ) ) : = x + t τ t g ( s , x ( s ) ) d s + t τ t t h ( u , s , x ( s ) ) d u d s .
We can prove the same result, [27] (Theorem 4) using the Lyapunov-Razumikhin method and the Lyapunov function
W ( t , x ) : = x = i = 1 n x i = x 1 + . . . + x n .
Clearly, this Lyapunov function does not include the term t τ t g ( s , x ( s ) ) d s . The time derivative of this term gives
d d t t τ t g ( s , x ( s ) ) d s = g ( t , x ( t ) ) g ( t τ , x ( t τ ) ) .
Based on this approach, we can obtain the result of Du [27] (Theorem 4) under weaker conditions. Namely, we remove the following hypothesis from Du [27] (Theorem 4):
g ( t τ , x ( t τ ) ) g ( t , x ( t τ ) ) 0 , t R + .
To the best of our information, this is a stronger condition and the satisfaction of this hypothesis can be difficult. Removing this condition from that of Du [27] (Theorem 4) leads to an important and strong advantage during the applications of that kind of equation.
(4)
Du [27] (Theorem 4) proved the related theorem without giving an example in a particular case, which verifies the hypotheses of [27] (Theorem 4). In this paper, we provided two examples and solved them with MATLAB software, which verifies the applicability of the results of this paper.

6. Conclusions

This paper has proposed an effective way to discuss some qualitative properties of solutions of nonlinear Volterra integro-differential equations with Caputo fractional derivatives and multiple constant retardations. Here, a new mathematical model consisting of non-linear fractional Volterra integro-differential equations with Caputo fractional derivatives and two constant retardations was considered. New sufficient conditions for the uniform stability, asymptotic stability and Mittag–Leffer stability of the zero solution, as well as the boundedness of the solutions were obtained. The presented results were proved by defining an appropriate Lyapunov function and applying the Lyapunov–Razumikhin method. An advantage of the new function and method used here is that they eliminate using Gronwall’s inequality. Compared to related results in the literature, the conditions here are new, more general, simple and convenient to apply. Examples to show the application of the theorems have been included.

Author Contributions

Investigation, O.T., Ö.A., C.T. and J.-C.Y. Writing—Original Draft: O.T., Ö.A., C.T. and J.-C.Y. All authors have read andagreed to the published version of the manuscript.

Funding

The research of J. C. Yao was partially supported by the Grant MOST 108-2115-M-039-005-MY3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the three anonymous referees and the handling for many useful comments and suggestions, leading to a substantial improvement in the presentation of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burton, T.A. Volterra integral and differential equations. In Mathematics in Science and Engineering, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2005; p. 202. [Google Scholar]
  2. Rahman, M. Integral Equations and Their Applications; WIT Press: Southampton, UK, 2007. [Google Scholar]
  3. Wazwaz, A.M. Linear and nonlinear integral equations. In Methods and Applications; Higher Education Press: Beijing, China; Springer: Heidelberg, Germany, 2011. [Google Scholar]
  4. Alahmadi, F.; Raffoul, Y.; Alharbi, S. Boundedness and stability of solutions of nonlinear Volterra integro-differential equations. Adv. Dyn. Syst. Appl. 2018, 13, 19–31. [Google Scholar]
  5. Andreev, A.S.; Peregudova, O.A. On the stability and stabilization problems of Volterra integro-differential equations. Russ. J. Nonlinear Dyn. 2018, 14, 387–407. [Google Scholar] [CrossRef] [Green Version]
  6. Berezansky, L.; Domoshnitsky, A. On stability of a second order integro-differential equation. Nonlinear Dyn. Syst. Theory 2019, 19, 117–123. [Google Scholar]
  7. Chang, Y.-K.; Ponce, R. Uniform exponential stability and applications to bounded solutions of integro-differential equations in Banach spaces. J. Integr. Equ. Appl. 2018, 30, 347–369. [Google Scholar] [CrossRef]
  8. Anh, N.P.H.; Tran, A. The New stability criteria for nonlinear Volterra integro-differential equations. Acta Math. Vietnam. 2018, 43, 485–501. [Google Scholar]
  9. Raffoul, Y.; Rai, H. Uniform stability in nonlinear infinite delay Volterra integro-differential equations using Lyapunov functionals. Nonauton. Dyn. Syst. 2016, 3, 14–23. [Google Scholar] [CrossRef]
  10. Tunç, C. Properties of solutions to Volterra integro-differential equations with delay. Appl. Math. Inf. Sci. 2016, 10, 1775–1780. [Google Scholar] [CrossRef]
  11. Tunç, C. Qualitative properties in nonlinear Volterra integro-differential equations with delay. J. Taibah Univ. Sci. 2017, 11, 309–314. [Google Scholar] [CrossRef] [Green Version]
  12. Tunç, C. Asymptotic stability and boundedness criteria for nonlinear retarded Volterra integro-differential equations. J. King Saud Univ. Sci. 2016, 30, 3531–3536. [Google Scholar] [CrossRef]
  13. Tunç, C. Stability and boundedness in Volterra-integro differential equations with delays. Dyn. Syst. Appl. 2017, 26, 121–130. [Google Scholar]
  14. Tunç, O. On the qualitative analyses of integro-differential equations with constant time lag. Appl. Math. Inf. Sci. 2020, 14, 57–63. [Google Scholar]
  15. Tunç, C.; Tunç, O. New results on the stability, integrability and boundedness in Volterra integro-differential equations. Bull. Comput. Appl. Math. 2018, 6, 41–58. [Google Scholar]
  16. Tunç, C.; Tunç, O. New qualitative criteria for solutions of Volterra integro-differential equations. Arab J. Basic Appl. Sci. 2018, 25, 158–165. [Google Scholar] [CrossRef]
  17. Tunç, C.; Tunç, O. A note on the qualitative analysis of Volterra integro-differential equations. J. Taibah Univ. Sci. 2019, 13, 490–496. [Google Scholar]
  18. Vasundhara, D.J.; Mc Rae, F.A.; Drici, Z. Variational Lyapunov method for fractional differential equations. Comput. Math. Appl. 2012, 64, 2982–2989. [Google Scholar]
  19. Wang, K. Uniform asymptotic stability in functional-differential equations with infinite delay. Ann. Differ. Equ. 1993, 9, 325–335. [Google Scholar]
  20. Wang, Q. The stability of a class of functional differential equations with infinite delays. Ann. Differ. Equ. 2000, 16, 89–97. [Google Scholar]
  21. Wang, Y.; Li, T. Stability analysis of fractional-order nonlinear systems with delay. Math. Probl. Eng. 2014. [Google Scholar] [CrossRef]
  22. Wen, Y.; Zhou, X.F.; Zhang, Z.; Liu, S. Lyapunov method for nonlinear fractional differential systems with delay. Nonlinear Dyn. 2015, 82, 1015–1025. [Google Scholar] [CrossRef]
  23. Zhou, B.; Egorov, A.V. Razumikhin and Krasovskii stability theorems for time-varying time-delay systems. Autom. J. IFAC 2016, 71, 281–291. [Google Scholar] [CrossRef]
  24. Agarwal, R.; Bohner, M.; Domoshnitsky, A.; Goltser, Y. Floquet theory and stability of nonlinear integro-differential equations. Acta Math. Hung. 2005, 109, 305–330. [Google Scholar] [CrossRef]
  25. Hristova, S.; Tunç, C. Stability of nonlinear Volterra integro-differential equations with Caputo fractional derivative and bounded delays. Electron. J. Differ. Equ. 2019, 30, 11. [Google Scholar]
  26. Hristova, S.; Tunç, C. On the stability properties of retarded Volterra integro-fractional differential equations with Caputo derivative. AIP Conf. Proc. 2021, 2321, 030013. [Google Scholar]
  27. Du, X.T. Some kinds of Liapunov functional in stability theory of RFDE. Acta Math. Appl. Sin. 1995, 11, 214–224. [Google Scholar] [CrossRef]
  28. Hale, J.K.; Verduyn Lunel, S.M. Introduction to functional-differential equations. In Applied Mathematical Sciences; Springer: New York, NY, USA, 1993; p. 99. [Google Scholar]
  29. Razumihin, B.S. On stability of systems with retardation. Russ. Prikl. Mat. Meh. 1956, 20, 500–512. [Google Scholar]
  30. Razumihin, B.S. The application of Lyapunov’s method to problems in the stability of systems with delay. Automat. Remote Control 1960, 21, 515–520. [Google Scholar]
  31. Agarwal, R.; Hristova, S.; O’Regan, D. Lyapunov functions and stability of Caputo fractional differential equations with delays. Differ. Equ. Dyn. Syst. 2020, in press. [Google Scholar] [CrossRef]
  32. Agarwal, R.; O’Regan, D.; Hristova, S. Stability of Caputo fractional differential equations by Lyapunov functions. Appl. Math. 2015, 60, 653–676. [Google Scholar] [CrossRef] [Green Version]
  33. Hu, J.-B.; Lu, G.-P.; Zhang, S.-B.; Zhao, L.-D. Lyapunov stability theorem about fractional system without and with delay. Commun. Nonlinear Sci. Numer. Simul. 2015, 20, 905–913. [Google Scholar] [CrossRef]
  34. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, S.; Jiang, W.; Li, X.; Zhou, X.-F. Lyapunov stability analysis of fractional nonlinear systems. Appl. Math. Lett. 2016, 51, 13–19. [Google Scholar] [CrossRef]
  36. Noeiaghdam, S.; Sidorov, D. Caputo-Fabrizio fractional derivative to solve the fractional model of energy supply-demand System. Math. Model. Eng. Probl. 2020, 7, 359–367. [Google Scholar] [CrossRef]
  37. Podlubny, I. Fractional differential equations. An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. In Mathematics in Science and Engineering; Academic Press, Inc.: San Diego, CA, USA, 1999; p. 198. [Google Scholar]
  38. Duarte-Mermoud, M.A.; Aguila-Camacho, N.; Gallegos, J.A.; Castro-Linares, R. Using general quadratic Lyapunov functions to prove Lyapunov uniform stability for fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 650–659. [Google Scholar] [CrossRef]
Figure 1. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 1 ( t ) of the system of of fractional retarded Volterra integro-differential equations (FrRIDEs) in Equation (11) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 1. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 1 ( t ) of the system of of fractional retarded Volterra integro-differential equations (FrRIDEs) in Equation (11) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g001
Figure 2. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 2 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 2. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 2 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g002
Figure 3. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 1 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 3. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 1 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g003
Figure 4. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 2 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values and different initial values when t 1 10 .
Figure 4. The behaviors of uniformly, asymptotically and Mittag–Leffler stable solution x 2 ( t ) of the system of of FrRIDEs in Equation (11) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values and different initial values when t 1 10 .
Axioms 10 00058 g004
Figure 5. The boundedness of solution x 1 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 5. The boundedness of solution x 1 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g005
Figure 6. The boundedness of solution x 2 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 6. The boundedness of solution x 2 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g006
Figure 7. The boundedness of solution x 1 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 7. The boundedness of solution x 1 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.9 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g007
Figure 8. The boundedness of solution x 2 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Figure 8. The boundedness of solution x 2 ( t ) of the system of IFrRDEs in Equation (13) for q = 0.5 , τ = 1 10 , ρ = 1 5 and different initial values when t 1 10 .
Axioms 10 00058 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tunç, O.; Atan, Ö.; Tunç, C.; Yao, J.-C. Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method. Axioms 2021, 10, 58. https://doi.org/10.3390/axioms10020058

AMA Style

Tunç O, Atan Ö, Tunç C, Yao J-C. Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method. Axioms. 2021; 10(2):58. https://doi.org/10.3390/axioms10020058

Chicago/Turabian Style

Tunç, Osman, Özkan Atan, Cemil Tunç, and Jen-Chih Yao. 2021. "Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method" Axioms 10, no. 2: 58. https://doi.org/10.3390/axioms10020058

APA Style

Tunç, O., Atan, Ö., Tunç, C., & Yao, J. -C. (2021). Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method. Axioms, 10(2), 58. https://doi.org/10.3390/axioms10020058

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