Next Article in Journal
Center-like Subsets in Semiprime Rings with Multiplicative Derivations
Next Article in Special Issue
On Some New Dynamic Hilbert-Type Inequalities across Time Scales
Previous Article in Journal
Quantization of the Rank Two Heisenberg–Virasoro Algebra
Previous Article in Special Issue
C0–Semigroups Approach to the Reliability Model Based on Robot-Safety System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order Sequential Linear Differential Equations with Nabla Derivatives on Time Scales

1
School of Science, Jiangnan University, Wuxi 214122, China
2
School of Mathematics and Statistics, Jiangsu Normal University, Xuzhou 221116, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(7), 447; https://doi.org/10.3390/axioms13070447
Submission received: 30 May 2024 / Revised: 29 June 2024 / Accepted: 30 June 2024 / Published: 1 July 2024
(This article belongs to the Special Issue Infinite Dynamical System and Differential Equations)

Abstract

:
In this paper, we present a general theory for fractional-order sequential differential equations with Riemann–Liouville nabla derivatives and Caputo nabla derivatives on time scales. The explicit solution, in the case of constant coefficients, for both the homogeneous and the non-homogeneous problems, are given using the ∇-Mittag-Leffler function, Laplace transform method, operational method and operational decomposition method. In addition, we also provide some results about a solution to a new class of fractional-order sequential differential equations with convolutional-type variable coefficients using the Laplace transform method.

1. Introduction

Continuous fractional calculus theory has been studied for more than 300 years; now, it has evolved in a respected mathematical discipline with numerous applications in many fields of science and engineering (see, for example, [1,2,3,4], and the references contained therein). However, the study of discrete problems and discrete models is still a relatively new field. The pioneering works [5,6] about fractional operators were introduced for q-calculus and difference calculus, respectively. It was not until the early 21st century that this topic further developed. The papers on this topic [7,8,9,10], which are required reading in order to deeply understand the background of discrete dynamic behaviors, obtained some interesting results by applying discrete fractional calculus to discrete chaos behaviors. In [11,12,13,14], the delta-type discrete fractional calculus was studied. In [15,16], the nabla-type discrete fractional calculus was studied. In order to unify differential equations and difference equations, Hilger [17] proposed was the first to propose the time scale. P.R. Williams [18] provided a definition of fractional integrals and derivatives on time scales to unify three cases of specific time scales, which improved the results in [19]. N.R.D.O. Bastos provided a definition of fractional ∆-integrals and ∆-derivatives on time scales in [20]. The application of delta fractional calculus and Laplace transform on some specific discrete time scales were also discussed in [21,22,23]. In light of the above work, we further studied the theory of fractional integrals and derivatives on general time scales in [24,25,26], where the ∇-Laplace transform, fractional ∇-power function, ∇-Mittag-Leffler function, Riemann-Liouville-type, and Caputo-type fractional ∇-integrals and fractional ∇-differentials were studied. The Δ-power function and fractional Δ-integrals and fractional ∆-differentials on time scales were, respectively, defined. Some of their properties were discussed in detail. After that, through using the Laplace transform method, the existence of the solution and the dependency of the solution upon the initial value for a Cauchy-type problem with Riemann-Liouville fractional ∇-derivatives and ∆-derivatives were studied. With the rise of interdisciplinary research, an increasing number of researchers are shifting their focus from the application of fractional equations to practical problems in time scales [3,27,28,29,30,31,32]. However, the theory of fractional linear differential equations is the foundation of the study of nonlinear fractional differential equation theory; thus, deeper theoretical research, such as [33,34,35,36], should not stagnate.
Inspired by these results, in this paper, we will continue the work we started in [24,25] by presenting a general theory for sequential differential equations of a fractional order on time scales. The aim of this paper is to present a general and unified theory to study the continuous and discrete sequential differential equations of a fractional order on time scales. Specifically, we will provide the explicit solution, in the case of constant coefficients, for both the homogeneous and the non-homogeneous problem using the ∇-Mittag-Leffler function, Laplace transform method, operational method and operational decomposition method. In addition, we studied a class of new fractional sequential differential equations with convolutional-type variable coefficients; in addition, some results about the solutions for this kind of fractional sequential differential equation are given using the Laplace transform method and convolution method.
The structure of this paper is as follows: In Section 2, we provide some preliminaries about time scales, convolution, etc., and we also detail definitions of nabla generalized power functions, fractional nabla integrals and fractional nabla derivatives on time scales, as well as a definition of the ∇-Mittag-Leffler function, which is an important tool for solving fractional differential equations. Then, in Section 3, the basic conception of sequential linear differential equations of a fractional order on time scales and the structure of a general solution are presented. In Section 4, we introduce a method, one that is independent of the Laplace transform, for obtaining a fundamental system of solutions for the equation L n α y ( t ) = 0 , which yields an explicit expression for the general solution to the non-homogeneous equation L n α y ( t ) = f ( t ) . In Section 5, we first prove a fact that the nabla generalized power function is positive, which is also an open problem (see [37]); then, we introduce the concept of α -analyticity on time scales and detail some results about a solution for fractional sequential differential equations with convolution-type variable coefficients using the Laplace transform method.

2. Preliminaries

First, readers can refer to reference [38] for some well-known definitions and theorems about equations on time scale T , where T is an nonempty, closed subset of real numbers. Next, we will briefly introduce some of the results of our previous research [24,25], which are important for the research in this article.
Definition 1.
(see [24]) On time scales, we define the fractional generalized-power function as
h ^ α ( t , t 0 ) = L , t 0 1 { 1 z α + 1 } ( t ) ( α > 1 ) ,
where z C { 0 } , t t 0 . If α > 0 , then
h ^ α ( t 0 , t 0 ) = 0 .
Definition 2.
(see [24]) Let f : [ t 0 , ) T C , then the solution of the shifting problem is
u t ( t , ρ ( τ ) ) = u τ ( t , τ ) , t , τ T , t τ t 0 , u ( t , t 0 ) = f ( t ) , t T , t t 0
which is called the shift of f and is denoted by f ˜ .
Definition 3.
(see [24]) Let f , g : T R , then f g is defined by
f g ( t ) = t 0 t f ˜ ( t , ρ ( τ ) ) g ( τ ) τ , t T ,
where f ˜ is the shift that was introduced in Definition 9.
Definition 4.
(see [24]) On time scales, the fractional generalized-power function h ^ α ( t , s ) is defined as the shift of h ^ α ( t , t 0 ) , which is
h ^ α ( t , s ) = h ^ α ( · , t 0 ) ˜ ( t , s ) ( t s t 0 ) .
Lemma 1.
(see [37]) Let T be an isolated time scale and r > 1 be a rational number. Then, there exists a unique power function h ^ r ( t , s ) on a time scale such that h ^ r ( t , s ) > 0 for t > s and h ^ r ( t , t 0 ) = L , t 0 1 { 1 z r + 1 } .
Further, in accordance with the results in [24], we have the following lemma.
Lemma 2.
Let T be an isolated time scale, and let α > 1 be a real number. Then, on a time scale, there exists a unique power function h ^ α ( t , s ) such that h ^ α ( t , s ) 0 for t > s and h ^ α ( t , t 0 ) = L , t 0 1 { 1 z α + 1 } .
In this paper, we always denote Ω : = [ t 0 , t 1 ] T as a finite interval on a time scale T ( sup T = ) .
Definition 5.
(see [24]) Assume that t , t 0 Ω . Define the Riemann–Liouville fractional-integral I , t 0 α f of order α > 0 as
I , t 0 α f ( t ) : = h ^ α 1 ( t , t 0 ) f ( t ) = t 0 t h ^ α 1 ( · , t 0 ) ˜ ( t , ρ ( τ ) ) f ( τ ) τ = t 0 t h ^ α 1 ( t , ρ ( τ ) ) f ( τ ) τ t > t 0 .
Definition 6.
(see [24]) Assume that t , t 0 Ω . Define the Riemann–Liouville fractional -derivative D , t 0 α f of order α 0 as
D , t 0 α f ( t ) = D m I , t 0 m α f ( t ) m = [ α ] + 1 ; t > t 0 .
Definition 7.
(see [25]) Assume that t , t 0 Ω . The C a p u t o f r a c t i o n a l   d e r i v a t i v e o f o r d e r α 0 is defined via the Riemann–Liouville fractional derivative by
D , t 0 α C f ( t ) : = D , t 0 α [ f ( t ) k = 0 m 1 h ^ k ( t , t 0 ) f k ( t 0 ) ] ( t > t 0 ) ,
where
m = [ α ] + 1 for α N ; m = α for α N .
Definition 8.
(see [18]) A subset I T is called a time scale interval if it is of the form I = A T for some real interval A R . On a time scale interval I, f : I R is said to be left dense and absolutely continuous if for all ε > 0 there exist δ > 0 such that k = 1 n | f ( b k ) f ( a k ) | < ε whenever a disjoint finite collection of sub-time scale intervals ( a k , b k ] T I for 1 k n satisfies k = 1 n | b k a k | < δ . This was then denoted by f A C . If f m 1 A C , then we denoted it as f A C m .
The following properties can be deduced from the above definitions.
Property 1.
(see [25]) Let α 0 and m be given by (3). If f ( t ) A C m ( Ω ) , then the Caputo fractional derivative D , t 0 α C f ( t ) exists almost everywhere on Ω k m .
(a)
If α N , D , t 0 α C f ( t ) is represented by
D , t 0 α C f ( t ) = h ^ m α 1 ( t , t 0 ) f m ( t ) = : I , t 0 m α D m f ( t ) ,
where m = [ α ] + 1 , then α N ,   D , t 0 α C f ( t 0 ) = 0 (where the notation D , t 0 α C f ( t 0 ) denotes the limit of D , t 0 α C f ( t ) as t t 0 + ). In particular, when 0 < α < 1 and f ( t ) A C ( Ω ) , we have
D , t 0 α C f ( t ) = h ^ α ( t , t 0 ) f ( t ) = : I , t 0 1 α f ( t ) .
(b)
If α = m N , then D , t 0 α C f ( t ) is represented by D , t 0 m C f ( t ) = f m ( t )   m N . In particular, we have
D , t 0 0 C f ( t ) = f ( t ) .
Throughout this paper, we denote f n = D n f = D , t 0 n f , n N .
Lemma 3.
(see [24]) Assume that α > 0 , m 1 < α m ( m N ) and f : Ω R for t 0 , t Ω k m with t 0 < t . Then, we have the following:
(1) 
If f L , p ( Ω ) , then
L , t 0 { I , t 0 α f ( t ) } ( z ) = 1 z α L , t 0 { f ( t ) } ( z ) ,
(2) 
and if f A C m ( Ω ) , then
L , t 0 D , t 0 α f ( t ) ( z ) = z α L , t 0 f ( t ) ( z ) j = 1 m z j 1 D , t 0 α j f ( t 0 )
for those regressive z C satisfying lim t { D j I , t 0 m α f ( t ) e ^ ν z ( t , t 0 ) } = 0 , j = 0 , 1 , · · · , m 1 .
Definition 9.
(see [24]) The - M i t t a g - L e f f l e r f u n c t i o n is defined as
F α , β ( λ ; t , t 0 ) = j = 0 λ j h ^ α j + β 1 ( t , t 0 )
provided the right hand series is convergent, where α , β > 0 , λ R .

3. Fundamental Conceptions and the Structure of General Solutions on Sequential Linear Differential Equations of a Fractional Order

In this section, let us first introduce a few necessary definitions.
Definition 10.
For n N , we defined a linear sequential fractional differential equation of order n α with equations of the form
k = 0 n b k ( t ) y ( k α ( t ) = f ( t ) ,
where b k ( t ) are given real functions, y ( 0 ( t ) : = y ( t ) , and y ( k α ( t ) : = D , t 0 k α y ( t ) ( k = 1 , 2 , · · · , n ) (which represents a fractional sequential derivative). For example, for the Riemann–Liouville derivative, D , t 0 α y is
D , t 0 α y : = D , t 0 α y , D , t 0 k α y : = D , t 0 α D , t 0 ( k 1 ) α y ( k = 2 , 3 , · · · , n ) ;
for the Caputo derivative, D , t 0 α y is
D , t 0 α y : = D , t 0 α C y , D , t 0 k α y : = D , t 0 α D , t 0 ( k 1 ) α y ( k = 2 , 3 , · · · , n ) .
In the following section, we mainly discuss Riemann–Liouville derivatives since the study of them is similar to the study of Caputo derivatives.
If t [ t 0 , t 1 ] T , b n ( t ) 0 , then Equation (10) may be expressed in its normal form for the Riemann-Liouville derivative as follows:
L n α y ( t ) : = D , t 0 n α y ( t ) + k = 0 n 1 a k ( t ) D , t 0 k α y ( t ) y ( n α ( t ) + k = 0 n 1 a k ( t ) y ( k α ( t ) = f ( t ) .
Note that Equation (11) reduces to
D , t 0 α Y ( t ) = A ( t ) Y ( t ) + B ( t ) ,
where
A ( t ) = 0 1 0 · · · 0 0 0 1 · · · 0 · · · · · · · · · · · · · · · · · · · · 1 a 0 a 1 a 2 · · · a n 1
and
B ( t ) = 0 0 · · · f ( t ) ; Y ( t ) = y 1 ( t ) y 2 ( t ) · · · y n ( t )
apply by just changing the variables
y 1 ( t ) = y ( t ) ; D , t 0 α y j ( t ) = y j + 1 ( t ) ( j = 1 , 2 , · · · , n 1 ) .
Definition 11.
For α > 0 , we call the α-Wronskian of n functions u j ( t ) ( j = 1 , · · · , n ) those that have fractional sequential derivatives up to order ( n 1 ) α in interval V ( t 0 , t 1 ] T with the following determinant:
W α ( u 1 , · · · , u n ) ( t ) = u 1 ( t ) u 2 ( t ) · · · u n ( t ) u 1 ( α ( t ) u 2 ( α ( t ) · · · u n ( α ( t ) u 1 ( 2 α ( t ) u 2 ( 2 α ( t ) · · · u n ( 2 α ( t ) · · · · · · · · · · · · · · · · · · u 1 ( ( n 1 ) α ( t ) u 2 ( ( n 1 ) α ( t ) · · · u n ( ( n 1 ) α ( t ) .
To simplify the notation, the above will be presented as W α ( t ) .
We will use W α ( t ) for the corresponding Wronskian matrix.
Definition 12.
We defined a fundamental system of solutions of Equation L n α y ( t ) = 0 in V [ t 0 , t 1 ] T as a family of n functions, which are the solutions of this equation and are linearly independent in V.
Definition 13.
The term general solution of the fractional differential equation L n α y ( t ) = f ( t ) refers to any solution to this equation (which only depend on n-independent constants).
Through using the similar iterative approximation method used in [24], we can prove the existence and uniqueness of the global solutions of (12). From these, we can obtain the following two theorems on the existence and uniqueness of global solutions to Equation L n α y ( t ) = f ( t ) with certain initial conditions.
Theorem 1.
Assume that y 0 k R ( k = 0 , 1 , , n 1 ) has been given real numbers, and let a j ( t ) C l d [ t 0 , t 1 ] T ( j = 0 , 1 , , n 1 ) and f ( t ) C l d [ t 0 , t 1 ] T be given left dense and continuous functions on [ t 0 , t 1 ] T . Then, there exists a unique left dense and continuous solution y ( t ) C l d [ t 0 , t 1 ] T of the Cauchy problem
L n α y ( t ) = f ( t ) , D , t 0 k α y ( t ) : = y ( k α ( t 0 ) = y 0 k ( k = 0 , 1 , · · · , n 1 ) .
Theorem 2.
The initial value problem
L n α y ( t ) = 0 , D , t 0 k α y ( t ) : = y ( k α ( t 0 ) = 0 ( k = 0 , 1 , · · · , n 1 )
has only the trivial solution y ( t ) = 0 .
Proposition 1.
Any linear combination of solutions of the homogeneous equation
L n α y ( t ) = 0
is also a solution to this equation.
Proof. 
The result is immediately evident when keeping in mind the linearity of L n α . □
Proposition 2.
If u j ( t ) ( j = 1 , 2 , , n ) are linearly dependent in [ t 0 , t 1 ] T , then | W α ( t ) | = 0 .
Proof. 
Since u j ( t ) ( j = 1 , 2 , , n ) are linearly dependent in [ t 0 , t 1 ] T , then there exists n constants { c j } j = 1 n , of which not all are zero, such that
W α ( t ) c 1 c 2 · · · c n = 0
has a nontrivial solution; thus, | W α ( t ) | = 0 . □
It follows from Proposition 2 that we can obtain the following Proposition 3.
Proposition 3.
If there exists a t such that | W α ( t ) | 0 , then u j ( t ) ( j = 1 , 2 , , n ) are linearly independent in [ t 0 , t 1 ] T .
Assume u j ( t ) ( j = 1 , 2 , , n ) are n solutions to Equation (13) and that the linear dependency of them are closely connected with the corresponding Wronskian matrix.
Proposition 4.
u j ( t ) ( j = 1 , 2 , , n ) are linearly dependent if and only if there exists a t such that | W α ( t ) | = 0 .
Proof. 
As the necessity was proven in Proposition 2, we will now prove the sufficiency. Suppose t t 0 , as well as consider the equation
W α ( t ) C = 0 ,
where
C = c 1 c 2 · · · c n 0 .
Then, as W α ( t ) = 0 , the equation has a nontrivial solution c j ( j = 1 , 2 , · · · , n ) . As such, we will construct a function using the following constants:
u ( t ) = j = 1 n c j u j ( t ) ,
and we obtain that u ( t ) is a solution to Equation (13). From (15), we obtain that u ( t ) satisfies the initial value condition
u ( t ) = u ( α ( t ) = u ( ( n 1 ) α ( t ) = 0 .
However, u ( t ) 0 is also a solution to the equation, satisfying the initial value condition. Via the uniqueness of the solution, we have
j = 1 n c j u j ( t ) = 0 ;
thus, u j ( t ) ( j = 1 , 2 , , n ) are linearly dependent. □
Proposition 5.
If { u j ( t ) } j = 1 n is a fundamental system of solutions to the equation L n α y ( t ) = 0 in a certain interval V ( t 0 , t 1 ] T , then the general solution to this differential equation in V is given by
y g ( t ) = k = 1 n c k u k ( t ) ,
where { c k } k = 1 n are arbitrary constants.
Proof. 
Since y g ( t ) is a solution to (13), we need only show that any other solution is a special case of y g ( t ) . Let u ( t ) be a solution to Equation (13), which satisfies, for a given θ ( t 0 , t 1 ] T , the following conditions:
u ( k α ( θ ) = u 0 k R ( k = 0 , 1 , · · · , n 1 ) .
Since u j ( t ) j = 1 n is a fundamental system of solutions, | W α ( θ ) | 0 , as well as the equation
c 1 u 1 ( θ ) + c 2 u 2 ( θ ) + · · · + c n u n ( θ ) = u 0 0 c 1 u 1 ( α ( θ ) + c 2 u 2 ( α ( θ ) + · · · + c n u n ( α ( θ ) = u 0 1 · · · c 1 u 1 ( ( n 1 ) α ( θ ) + c 2 u 2 ( ( n 1 ) α ( θ ) + · · · + c n u n ( ( n 1 ) α ( θ ) = u 0 n 1 ,
has a unique solution of c 1 , c 2 , · · · , c n . It is evident that k = 1 n c k u k ( t ) satisfies the initial condition u ( k α ( θ ) = u 0 k ( k = 0 , 1 , · · · , n 1 ) . Then, through using Theorem 1, we have
u ( t ) = k = 1 n c k u k ( t ) ,
which proves the theorem. □
Proposition 6.
The set of all solutions of the differential equation L n α y ( t ) = 0 , in a certain interval V ( t 0 , t 1 ] T , is a vector space of n dimensions.
Proof. 
The result follows directly from Proposition 5. □
Proposition 7.
If y p ( t ) is a particular solution to the equation L n α y ( t ) = f ( t ) , then the general solution to this equation is given by
y g ( t ) = y h ( t ) + y p ( t ) ,
where y h ( t ) is the general solution to the associated homogeneous equation L n α y ( t ) = 0 .
Proof. 
The result is evident. □

4. The Solution of Sequential Linear Differential Equations with Constant Coefficients

This section introduces a general theory to solve Riemann–Liouville sequential linear fractional differential equations for fractional cases. Some general methods were developed for the purpose of obtaining the general solution to linear sequential fractional differential equations with constant coefficients, and this was achieved using the roots of the characteristic polynomial of the corresponding homogeneous equation, or in finding a certain Green function in the non-homogeneous case.

4.1. The Solution of Integral-Order Differential Equations in the Homogeneous Case

To provide a method or inspiration to solve linear sequential fractional differential equations, we firstly and predominantly undertook a discussion on integral-order differential equations with constant coefficients in the homogeneous case. First, we studied the following second-order differential equation:
D 2 y ( t ) + a D y ( t ) + b y ( t ) = 0 ,
which is denoted as L y ( t ) = 0 .
We then sought a solution of the form
y ( t ) = e ^ λ ( t , t 0 ) ,
where we can derive
D e ^ λ ( t , t 0 ) = λ e ^ λ ( t , t 0 )
and
D 2 e ^ λ ( t , t 0 ) = λ 2 e ^ λ ( t , t 0 ) .
By substituting them into (16), we can obtain the characteristic equation:
P ( λ ) = λ 2 + a λ + b = 0 .
If we suppose the characteristic roots are λ 1 , λ 2 , then there the following two cases present themselves:
(1)
When λ 1 λ 2 , the linearly independent solutions of (16) are
y 1 ( t ) = e ^ λ 1 ( t , t 0 ) , y 2 ( t ) = e ^ λ 2 ( t , t 0 ) .
(2)
When λ 1 = λ 2 = λ , it is evident that y 1 ( t ) = e ^ λ ( t , t 0 ) is one solution of (16). But, how does one derive another solution y 2 ( t ) that is independent of y 1 ( t ) ?
If we take the derivative of y 1 ( t ) = e ^ λ ( t , t 0 ) for λ , then there is
L e ^ λ ( t , t 0 ) λ = λ L e ^ λ ( t , t 0 ) = λ P ( λ ) e ^ λ ( t , t 0 ) = P ( λ ) e ^ λ ( t , t 0 ) + P ( λ ) e ^ λ ( t , t 0 ) λ .
Since λ = λ 1 = λ 2 is a double root of P ( λ ) = 0 , P ( λ ) = P ( λ ) = 0 , and thus
y 2 ( t ) = e ^ λ ( t , t 0 ) λ
is the solution of (16), then apparently y 1 ( t ) and y 2 ( t ) are linearly independent and therefore we have the following calculations:
Theorem 3.
For differential Equation (16), if λ 1 , λ 2 are characteristic roots of
P ( λ ) = λ 2 + a λ + b = 0 ,
then we have the following:
(1) 
When λ 1 λ 2 , the linearly independent solutions of (16) are
y 1 ( t ) = e ^ λ 1 ( t , t 0 ) , y 2 ( t ) = e ^ λ 2 ( t , t 0 ) .
(2) 
When λ 1 = λ 2 = λ , the linearly independent solutions of (16) are
y 1 ( t ) = e ^ λ ( t , t 0 ) , y 2 ( t ) = e ^ λ ( t , t 0 ) λ .
For a general linear homogeneous differential equation, we can utilize the following:
L y ( t ) = D m y ( t ) + a 1 D m 1 y ( t ) + · · · + a m 1 D y ( t ) + a m y ( t ) = 0 ,
Theorem 4.
For the differential Equation (17), its characteristic equation is
P ( λ ) = λ m + a 1 λ m 1 + · · · + a m 1 λ + a m = 0 .
(1) 
If there are m distinct roots λ 1 , λ 2 , · · · , λ m for the characteristic Equation (18), then
y 1 ( t ) = e ^ λ 1 ( t , t 0 ) , y 2 ( t ) = e ^ λ 2 ( t , t 0 ) , · · · , y m ( t ) = e ^ λ m ( t , t 0 )
are the fundamental solution system of Equation (17).
(2) 
If λ is a root of multiplicity k of the characteristic Equation (18), then the functions
y 1 ( t ) = e ^ λ ( t , t 0 ) , y 2 ( t ) = e ^ λ ( t , t 0 ) λ , y 2 ( t ) = 2 e ^ λ ( t , t 0 ) λ 2 , · · · , y k ( t ) = k e ^ λ ( t , t 0 ) λ k
are linearly independent solutions of Equation (17).
Proof. 
The proof of (1) is evident.
(2) We only need to prove (for l = 1 , 2 , · · · , k , L y l ( t ) = 0 , where y l ( t ) = l 1 e ^ λ ( t , t 0 ) λ l 1 ), as λ is a root of multiplicity k of the characteristic Equation (18),
P ( λ ) = P ( λ ) = · · · = P ( k 1 ) ( λ ) = 0 , P ( k ) ( λ ) 0 ,
as well as
L y l ( t ) = L l 1 e ^ λ ( t , t 0 ) λ l 1 = l 1 λ l 1 L e ^ λ ( t , t 0 ) = l 1 λ l 1 P ( λ ) e ^ λ ( t , t 0 ) = j = 0 l 1 l 1 j P ( j ) ( λ ) l 1 j e ^ λ ( t , t 0 ) λ l 1 j .
From (19), we have L y l ( t ) = 0 , l = 1 , 2 , · · · , k , that is, y l ( t ) , l = 1 , 2 , · · · , k , which are linearly independent solutions of Equation (17). □
In general, when the characteristic root is a multiple, we have the following theorem:
Theorem 5.
Let λ 1 , λ 2 , · · · , λ r be r distinct roots of multiplicity k 1 , k 2 , · · · , k r ( j = 1 r k j = n ) of the characteristic Equation (18). Then, the functions
e ^ λ 1 ( t , t 0 ) , e ^ λ 1 ( t , t 0 ) λ 1 , 2 e ^ λ 1 ( t , t 0 ) λ 1 2 , · · · , k e ^ λ 1 ( t , t 0 ) λ 1 k ; e ^ λ 2 ( t , t 0 ) , e ^ λ 2 ( t , t 0 ) λ 2 , 2 e ^ λ 2 ( t , t 0 ) λ 2 2 , · · · , k e ^ λ 2 ( t , t 0 ) λ 2 k ; · · · e ^ λ r ( t , t 0 ) , e ^ λ r ( t , t 0 ) λ r , 2 e ^ λ r ( t , t 0 ) λ r 2 , · · · , k e ^ λ r ( t , t 0 ) λ r k
are linearly independent solutions of Equation (17).

4.2. General Solution of Sequential Differential Equations of a Fractional Order in the Homogeneous Case

As we were inspired by the thinking of solving integral-order differential equations, we will discuss sequential linear differential equations of a fractional order in this section.
Let
L n α y ( t ) : = D , t 0 n α y ( t ) + k = 0 n 1 a k D , t 0 k α y ( t ) ,
where the coefficients { a j } j = 1 n 1 are real constants.
As in the usual case, we shall seek the solution of (20) in the form y ( t ) = F α ( λ ; t , t 0 ) , where F α ( λ ; t , t 0 ) is defined as
F α ( λ ; t , t 0 ) : = F α , α ( λ ; t , t 0 ) = k = 0 λ k h ^ α k + α 1 ( t , t 0 ) .
In the following, we first give a lemma for F α ( λ ; t , t 0 ) .
Lemma 4.
For F α ( λ ; t , t 0 ) , it is valid that
D , t 0 α F α ( λ ; t , t 0 ) = λ F α ( λ ; t , t 0 ) , λ F α ( λ ; t , t 0 ) = h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) .
Proof. 
According to the definition of F α ( λ ; t , t 0 ) , as well as by taking D , t 0 α h ^ α 1 ( t , t 0 ) = 0 (see Formula (85) of Property 1 in [24]) into account, we have
D , t 0 α F α ( λ ; t , t 0 ) = D , t 0 α k = 0 λ k h ^ α k + α 1 ( t , t 0 ) = D , t 0 α h ^ α 1 ( t , t 0 ) + k = 1 λ k h ^ α k + α 1 ( t , t 0 ) = k = 1 λ k h ^ α k 1 ( t , t 0 ) = k = 0 λ k + 1 h ^ α k + α 1 ( t , t 0 ) = λ k = 0 λ k h ^ α k + α 1 ( t , t 0 ) = λ F α ( λ ; t , t 0 ) .
For the proof of the second formula, through using the convolution property of the generalized polynomial h ^ α ( t , t 0 ) , we have
λ F α ( λ ; t , t 0 ) = λ k = 0 λ k h ^ α k + α 1 ( t , t 0 ) = k = 1 λ k 1 h ^ α k + α 1 ( t , t 0 ) = k = 1 λ k 1 h ^ α 1 ( t , t 0 ) h ^ α k 1 ( t , t 0 ) = h ^ α 1 ( t , t 0 ) k = 1 λ k 1 h ^ α k 1 ( t , t 0 ) = h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) .
The proof is thus completed. □
From Lemma 4, we have
D , t 0 α F α ( λ ; t , t 0 ) = λ F α ( λ ; t , t 0 ) , D , t 0 2 α F α ( λ ; t , t 0 ) = λ 2 F α ( λ ; t , t 0 ) , · · · D , t 0 n α F α ( λ ; t , t 0 ) = λ n F α ( λ ; t , t 0 ) . l λ l F α ( λ ; t , t 0 ) = i = 1 l h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) ,
where i = 1 l h ^ α 1 ( t , t 0 ) = h ^ α 1 ( t , t 0 ) · · · h ^ α 1 ( t , t 0 ) , that is, the m time convolution of h ^ α 1 ( t , t 0 ) . It follows from (20) and the equalities detailed above that
L n α [ F α ( λ ; t , t 0 ) ] = P n ( λ ) F α ( λ ; t , t 0 ) ,
where
P n ( λ ) = λ n + k = 0 n 1 a k λ k
is the characteristic polynomial associated with the equation L n α [ y ( t ) ] = 0 , and the root of P n ( λ ) = 0 is called a characteristic root. Thus, y ( t ) = F α ( λ ; t , t 0 ) is the solution to Equation (20) if and only if λ is the root of (21).
Theorem 6.
For the differential Equation (20), there are the following two cases for a characteristic root that corresponds to the characteristic equation:
(1) 
When all of the characteristic roots are simple, that is, if there are n-distinct roots λ 1 , λ 2 , . . . , λ n for the characteristic polynomial (21), then
y j ( t ) = F α ( λ j ; t , t 0 ) ( j = 1 , 2 , . . . , n )
are the fundamental solution system of Equation (20).
(2) 
When the characteristic root is a multiple, that is, if λ is a root of multiplicity k ( 1 < k n ) of the characteristic polynomial (21), then the functions
i = 1 j h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) , ( j = 0 , 1 , 2 , . . . , k 1 )
are linearly independent solutions of Equation (20).
Proof. 
The result of (1) is evident, and so we only provided a proof for (2). Since λ is a root of multiplicity k of P n ( λ ) = 0 , we have
P n ( λ ) = P n ( λ ) = · · · = P n ( k 1 ) ( λ ) = 0 , P n ( k ) ( λ ) 0 .
In addition, it follows from L n α [ F α ( λ ; t , t 0 ) ] = P n ( λ ) F α ( λ ; t , t 0 ) and the classical Leibniz rule that, for j = 0 , 1 , · · · , k 1 ,
L n α [ j λ j F α ( λ ; t , t 0 ) ] = j λ j L n α [ F α ( λ ; t , t 0 ) ] = j λ j [ P n ( λ ) F α ( λ ; t , t 0 ) ] = m = 0 j j m P n ( m ) ( λ ) [ F α ( λ ; t , t 0 ) ] ( j m ) ,
and, using (22), we have
L n α [ j λ j F α ( λ ; t , t 0 ) ] 0 , j = 0 , 1 , 2 , . . . , k 1 .
Since
j λ j F α ( λ ; t , t 0 ) = i = 1 j h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) ,
we have
L n α [ i = 1 j h ^ α 1 ( t , t 0 ) F α ( λ ; t , t 0 ) ] 0 , j = 0 , 1 , 2 , . . . , k 1 .
The result then follows. □
When there are several multiple roots for the characteristic equation, we have the following results.
Theorem 7.
Let { λ j } j = 1 k be k-distinct roots of multiplicity { μ j } j = 1 k 1 μ j n , j = 1 k μ j = n of the characteristic polynomial (21). Then, the functions
F α ( λ 1 ; t , t 0 ) , h ^ α 1 ( t , t 0 ) F α ( λ 1 ; t , t 0 ) , · · · , i = 1 μ 1 1 h ^ α 1 ( t , t 0 ) F α ( λ 1 ; t , t 0 ) ; F α ( λ 2 ; t , t 0 ) , h ^ α 1 ( t , t 0 ) F α ( λ 2 ; t , t 0 ) , · · · , i = 1 μ 2 1 h ^ α 1 ( t , t 0 ) F α ( λ 2 ; t , t 0 ) ; · · · · · · · · · · · · F α ( λ k ; t , t 0 ) , h ^ α 1 ( t , t 0 ) F α ( λ k ; t , t 0 ) , · · · , i = 1 μ k 1 h ^ α 1 ( t , t 0 ) F α ( λ k ; t , t 0 )
are linearly independent solutions of Equation (20).
Proof. 
The result follows from Theorem 6. □
Example 1.
1. The fractional differential equation y ( 2 α ( t ) y ( t ) = 0 has the following fundamental system of solutions:
F α ( 1 ; t , t 0 ) , F α ( 1 ; t , t 0 ) .
This result follows from Theorem 6 if we take into account that λ 1 = 1 and λ 2 = 1 are the roots of the characteristic polynomial P 2 ( λ ) : = λ 2 1 = 0 .
2. The fractional differential equation y ( 2 α ( t ) 2 y ( α ( t ) + y ( t ) = 0 has the following fundamental system of solutions:
F α ( 1 ; t , t 0 ) , h ^ α 1 ( t , t 0 ) F α ( 1 ; t , t 0 ) .
Since λ 1 = 1 is the root of multiplicity two of the characteristic polynomial P 2 ( λ ) : = λ 2 2 λ + 1 = 0 , then Theorem 6 yields the above result.

4.3. General Solution of the Sequential Differential Equations of a Fractional Order in the Non-Homogeneous Case

Consider the following non-homogeneous linear differential equation of a fractional order:
L n α [ y ( t ) ] = f ( t ) ,
which is defined by (11). As in the case of ordinary differential equations, the general solution to this equation is a sum of a particular solution to Equation (23), as well as a general solution to the corresponding homogeneous Equation (20).
First, we applied the operational method to derive a particular solution y p ( t ) of (23).
Theorem 8.
Let { λ j } j = 1 k be k-different roots of multiplicity { μ j } j = 1 k ( j = 1 k μ j = n ) of the characteristic polynomial P n ( λ ) = 0 , which is associated with the equation L n α [ y ( t ) ] = 0 (where y ( k α ( t ) = D , t 0 k α y ( t ) ) . Then, the particular solution to Equation (23), with a given function f ( t ) on ( t 0 , t 1 ] T , is given by
y p ( t ) = [ j = 1 k ( n = 0 λ j n I , t 0 . ( n + 1 ) α ) μ j f ] ( t )
provided that the series in the right-hand side of (24) are uniformly convergent.
Proof. 
We can consider the inverse operator to the left for L n α in the form
L n α 1 = j = 1 k [ ( D , t 0 α λ j ) 1 ] μ j = j = 1 k D , t 0 α ( 1 λ j I , t 0 . α ) 1 μ j = j = 1 k ( n = 0 λ j n I , t 0 . ( n + 1 ) α ) μ j ,
which yields the formal solution to Equation (23) in (24). Applying the operator
L n α = j = 1 k ( D , t 0 α λ j ) μ j
to the function y p ( t ) in (24), as well as when using the uniform convergence of the series considered, we derive that y p ( t ) is the solution to the equation L n α [ y ( t ) ] = f ( t ) .
Second, we applied the Laplace transform method to explicitly derive a particular solution for the non-homogenous equation
L n α [ y ( t ) ] = f ( t ) .
Via Lemma 1, the Laplace transform L , t 0 of D , t 0 α f ( t ) , for a suitable f ( t ) , is given by
L , t 0 { D , t 0 α f ( t ) } ( z ) = z α L { f } ( z ) ;
therefore,
L , t 0 { D , t 0 n α f ( t ) } ( z ) = z n α L { f } ( z ) .
Theorem 9.
A particular solution to the non-homogenous linear fractional differential Equation (23) with y ( k α ( t ) = D , t 0 k α y ( t ) is given by
y p ( t ) = j = 1 k m = 1 μ j c j , m ( g j , m f ) ( t ) ,
where
g j , m ( t ) = L 1 { 1 ( z α λ j ) m } ( t ) = 1 m 1 ! i = 1 m 1 h ^ α 1 ( t , t 0 ) F α , α ( λ j ; t , t 0 )
and c j , m ( j = 1 , 2 , . . . , k ; m = 1 , 2 , . . . , μ j ) are constants that are defined by the following decomposition into simple fractions:
1 P n ( λ ) = 1 j = 1 k ( λ λ j ) μ j = j = 1 k m = 1 μ j c j , m ( λ λ j ) m ,
where ( g f ) ( t ) represents the convolution of g ( t ) and f ( t ) .
Proof. 
In applying the Laplace transform to Equation (23), we have
L L n α [ y ( t ) ] = P n ( z α ) L y ( t ) = L f ( t ) .
Using (27), we can obtain
L y ( t ) = j = 1 k m = 1 μ j L L 1 c j , m ( z α λ j ) m L f ( t ) = L j = 1 k m = 1 μ j L 1 c j , m ( z α λ j ) m f ( t ) = L j = 1 k m = 1 μ j c j , m L 1 1 ( z α λ j ) m f ( t ) .
By taking the inverse Laplace transform, we find the particular solution to (23) in the form
y p ( t ) = j = 1 k m = 1 μ j c j , m L 1 1 ( z α λ j ) m f ( t ) = j = 1 k m = 1 μ j c j , m ( g j , m f ) ( t ) ,
where g j , m ( t ) = L 1 { 1 ( z α λ j ) m } ( t ) . Through using Theorem 40 in [24], we obtain
g j , m ( t ) = L 1 { 1 ( z α λ j ) m } ( t ) = L 1 { 1 ( m 1 ) ! m 1 λ m 1 1 ( z α λ ) | λ = λ j } ( t ) = 1 ( m 1 ) ! m 1 λ m 1 L 1 { 1 ( z α λ ) } | λ = λ j ( t ) = 1 ( m 1 ) ! m 1 λ m 1 F α , α ( λ ; t , t 0 ) | λ = λ j ( t ) = 1 m 1 ! i = 1 m 1 h ^ α 1 ( t , t 0 ) F α , α ( λ j ; t , t 0 ) .
This yields the result in (26). The proof is thus completed. □
Lastly, we applied the operational decomposition method to derive a particular solution for some non-homogenous equations.
Theorem 10.
Let f ( t ) be a ld-continuous function. Then, the linear differential equation
y ( α ( t ) λ y t = f t
has the general solution
y g t = F α ( λ ; t , t 0 ) + y p ( t ) ,
where y ( α ( t ) = D , t 0 α y ( t ) , y p ( t ) = F α ( λ ; t , t 0 ) f t , which is a particular solution to (28).
Proof. 
Via Theorem 6, y ( t ) = F α ( λ ; t , t 0 ) is the general solution to the corresponding homogeneous equation y ( α ( t ) λ y t = 0 . Therefore, it is sufficient to verify that y p ( t ) is a particular solution to (28).
Let y ( t ) be a solution of (28), we then have
D , t 0 α 1 λ I , t 0 α y t = f t ;
therefore,
y ( t ) = 1 λ I , t 0 α 1 I , t 0 α f t = j = 0 I , t 0 ( j + 1 ) α λ j f t = j = 0 λ j h ^ ( j + 1 ) α 1 ( t , t 0 ) f ( t ) = F α ( λ ; t , t 0 ) f t = y p ( t ) .
Thus, y p ( t ) is a particular solution to (28). The proof is then completed. □
Theorem 11.
Let { λ j } j = 1 k be k-distinct roots of multiplicity { σ j } j = 1 k of the characteristic polynomial (21) for the following non-homogeneous differential equation
L n α y t = j = 1 k ( D , t 0 α λ j ) σ j y t = f t .
Then, the particular solution is given by
y p ( t ) = G α f t ,
where G α t is given by
G α t = j = 1 k l = 1 σ j F α , α ( λ j ; t , t 0 ) .
Proof. 
Let
y 1 ( t ) = j = 2 k ( D , t 0 α λ j ) σ j y t ,
then Equation (29) can be rewritten as
( D , t 0 α λ 1 ) σ 1 y 1 t = f ( t ) .
By successively applying the result of Theorem 10, we can obtain
y 1 t = l = 1 σ 1 F α , α ( λ 1 ; t , t 0 ) f ( t ) .
Let
y l ( t ) = j = l + 1 k ( D , t 0 α λ j ) σ j y t , l = 2 , 3 , · · · , k 1 , y k ( t ) = y ( t ) ,
then
( D , t 0 α λ l ) σ l y l t = y l 1 t , l = 2 , 3 , · · · , k .
As it is similar to the proof of (32), we can prove that
y l ( t ) = l = 1 σ l F α , α ( λ l ; t , t 0 ) y l 1 t , l = 2 , 3 , · · · , k .
Thus, we have
y ( t ) = j = 1 k l = 1 σ j F α , α ( λ j ; t , t 0 ) f ( t ) = ( G α f ) t .
The proof is thus completed. □

5. Solution of Fractional Sequential Differential Equations with Convolution

In this section, we analyze the solutions of the homogeneous equation L n α y ( t ) = 0 , with variable coefficients expressed by a convolution, where
L n α y ( t ) : = D , t 0 n α y ( t ) + k = 0 n 1 a k ( t ) D , t 0 k α y ( t ) .
First, we needed to prove that the power functions are non-negative on general time scales, as well as introduce the concept of α -analyticity.
Let P n = t 0 , t 1 , · · · , t n [ t 0 , t ] T be a partition of [ t 0 , t ] T , where t 0 < t 1 < · · · < t n = t . We denoted P n = max t i t i 1 : t i 1 < ρ ( t i ) for any i 1 , 2 , · · · , n , as well as either t i t i 1 P n , or t i t i 1 > P n , and t i 1 = ρ ( t i ) . P n can be viewed as an isolated time scale; thus, we denoted the graininess function on P n by υ P n . The Laplace transform on P n is denoted by L , t 0 , P n , and the power function on P n is denoted by h ^ α , P n ( t , t 0 ) .
Lemma 5.
Let e ^ z ( t , t 0 ) be the nabla exponential function on time scale T and e ^ z , P n ( t , t 0 ) be the nabla exponential function on time scale P n . It then holds that
(i) 
lim P n 0 e ^ z , P n ( t , t 0 ) = e ^ z ( t , t 0 ) ,
(ii) 
for any α > 1 , h ^ α ( t , t 0 ) 0 .
Proof. 
Via the definition of the nabla exponential function on a time scale (see [38]), we can obtain
e ^ z ( t , t 0 ) = e t 0 t ξ ^ υ ( τ ) ( τ ) τ = e i = 1 n t i 1 t i ξ ^ υ ( τ ) ( τ ) τ
and
e ^ z , P n ( t , t 0 ) = e t 0 t ξ ^ υ P n ( τ ) ( τ ) τ = e i = 1 n t i 1 t i ξ ^ υ P n ( τ ) ( τ ) τ .
If t i is left-scattered in T , and when P n is small enough, we have t i 1 = ρ ( t i ) . This implies that
t i 1 t i ξ ^ υ P n ( τ ) ( τ ) τ = 1 υ P n ( t i ) L o g ( 1 z υ P n ( t i ) ) · υ P n ( t i ) = 1 υ ( t i ) L o g ( 1 z υ ( t i ) ) · υ ( t i ) = t i 1 t i ξ ^ υ ( τ ) ( τ ) τ .
If t i is left-dense in T , then υ ( t i ) = 0 and ξ ^ υ ( τ ) ( t i ) = z . Since lim x 0 + 1 x L o g ( 1 x z ) = z , for any ε > 0 , there exists a δ > 0 , such that
1 x L o g ( 1 x z ) z < ε , 0 < x < δ .
Thus, if P n < δ , then one can infer that
t 0 t ξ ^ υ ( τ ) ( τ ) τ t 0 t ξ ^ υ P n ( τ ) ( τ ) τ = i = 1 n t i 1 t i ξ ^ υ ( τ ) ( τ ) ξ ^ υ P n ( τ ) ( τ ) τ υ P n ( t i ) < δ t i 1 t i ξ ^ υ ( τ ) ( τ ) ξ ^ υ P n ( τ ) ( τ ) τ υ P n ( t i ) < δ t i 1 t i ξ ^ υ ( τ ) ( τ ) z τ + t i 1 t i z ξ ^ υ P n ( τ ) ( τ ) τ < 2 ε t t 0 .
This shows that
lim P n 0 e ^ z , P n ( t , t 0 ) = e ^ z ( t , t 0 )
holds, which means that Conclusion (i) holds.
It follows from Lemma 2 that h ^ α , P n ( t , t 0 ) 0 applies for any α > 1 . Then, through using Definition 7 and Lemma 5, we have
lim P n 0 h ^ α , P n ( t , t 0 ) = lim P n 0 L , t 0 , P n 1 { 1 z α + 1 } ( t ) = lim P n 0 1 2 π i c i c + i 1 z α + 1 e ^ z , P n ( t , t 0 ) d z = 1 2 π i c i c + i 1 z α + 1 e ^ z ( t , t 0 ) d z = L , t 0 1 { 1 z α + 1 } ( t ) = h ^ α ( t , t 0 ) 0 ,
which shows that Conclusion (ii) holds.
The proof is thus completed. □
Definition 14.
Let α ( 0 , 1 ] , f ( t ) be a real function defined on Ω. Then, f ( t ) is said to be α-analytic at t 0 if there exists an interval N + ( t 0 ) = [ t 0 , t 0 + δ ] T such that, for all t N + ( t 0 ) , f ( t ) can be expressed as n = 1 c n h ^ n α 1 ( t , t 0 ) ( c n R ) .
Since D h ^ α + 1 ( t , t 0 ) = h ^ α ( t , t 0 ) and h ^ α ( t , t 0 ) 0 , for any α > 1 , we can obtain h ^ α ( t , t 0 ) , which increases for any α > 0 . This implies that there exists a ρ > 0 , such that the series n = 1 c n h ^ n α 1 ( t , t 0 ) is absolutely convergent for any t [ t 0 , t 0 + ρ ) T and is divergent for any t > t 0 + ρ . We called ρ the convergence radius of the series.
Definition 15.
A point t 0 is said to be an α-ordinary point of the equation L n α y ( t ) = 0 if the functions a k ( t ) ( k = 0 , 1 , , n 1 ) are α-analytic at t 0 .
The following properties are clear.
Property 2.
Let α ( 0 , 1 ] . If f ( t ) is α-analytic at t 0 , with the convergence radius ρ ( ρ > 0 ) , then
D , t 0 α f ( t ) = D , t 0 α n = 1 c n h ^ n α 1 ( t , t 0 ) = n = 1 c n D , t 0 α h ^ n α 1 ( t , t 0 ) .
Property 3.
(see also [24,25]) Let α ( 0 , 1 ] , t t 0 . Then,
D , t 0 α h ^ β ( t , t 0 ) = h ^ β α ( t , t 0 ) , β > α 1 .
and
D , t 0 α C h ^ β ( t , t 0 ) = h ^ β α ( t , t 0 ) , β > α 1 .

5.1. Solutions around an Ordinary Point of a Fractional Differential Equation of Order α

In this section, we analyze the existence of solutions for the equation
L α y ( t ) : = y ( α ( t ) + p ( t ) y ( t ) = 0
around an α -ordinary point t 0 Ω , with p ( t ) defined in Ω . We separately consider the cases where y ( α ( t ) represents the Riemann–Liouville and Caputo fractional derivatives of an order α of the function y ( t ) . Since t 0 is an α -ordinary point, p ( t ) can be expressed as follows:
p ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) ,
where this series is convergent for t [ t 0 , t 0 + ρ ) T , with ρ > 0 .
Theorem 12.
Let α ( 0 , 1 ] and t 0 Ω be an α-ordinary point for the equation
L α y ( t ) : = D , t 0 α y ( t ) + p ( t ) y ( t ) = 0 .
Then, there exists a unique function
y ( t ) = n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) ,
which is the solution to Equation (37) for t ( t 0 , t 0 + ρ ) T and satisfies the initial condition a 0 = lim t t 0 y h ^ α ( · , t 0 ) ( t ) .
Proof. 
We sought a solution of (37) as follows:
y ( t ) = n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) .
If y ( t ) is given by (38), then when using (34), (36) and D , t 0 α h ^ α 1 ( t , t 0 ) = 0 (see Formula (85) in [24]), we have
D , t 0 α y ( t ) = D , t 0 α n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) = D , t 0 α a 0 h ^ α 1 ( t , t 0 ) + n = 1 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) = n = 1 a n h ^ n α 1 ( t , t 0 ) = n = 0 a n + 1 h ^ ( n + 1 ) α 1 ( t , t 0 )
and
p ( t ) y ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) = n = 1 k = 0 n 1 p n k h ^ ( n k ) α 1 ( t , t 0 ) a k h ^ ( k + 1 ) α 1 ( t , t 0 ) = n = 1 k = 0 n 1 p n k a k h ^ ( n + 1 ) α 1 ( t , t 0 ) .
Through substituting (39) and (40) into (37), we can obtain
D , t 0 α y ( t ) + p ( t ) y ( t ) = n = 0 a n + 1 h ^ ( n + 1 ) α 1 ( t , t 0 ) + n = 1 k = 0 n 1 p n k a k h ^ ( n + 1 ) α 1 ( t , t 0 ) = 0 .
Thus, we have the following recurrence formula, which allows us to express a n ( n > 0 ) in terms of a 0 as follows:
a 1 = 0 , a n + 1 = k = 0 n 1 p n k a k , n = 1 , 2 , · · · .
When the following equation does not apply
y h ^ α ( · , t 0 ) = a 0 + n = 2 a n h ^ n α ( t , t 0 ) ,
then we have a 0 = lim t t 0 y h ^ α ( · , t 0 ) ( t ) .
If we select a n with the above formula, then (38) is a solution of (37). The proof is thus completed. □
Theorem 13.
Let α ( 0 , 1 ] and t 0 Ω be an α-ordinary point for equation
L α y ( t ) : = D , t 0 α C y ( t ) + p ( t ) y ( t ) = 0 ,
where p ( t ) is defined in Ω . Then, there exists a unique function y ( t ) = n = 0 a n h ^ n α ( t , t 0 ) , which is the solution to Equation (41) for t ( t 0 , t 0 + ρ ) T and satisfies the initial condition y ( t 0 ) = a 0 .
Proof. 
We sought the series representation for the solution in the form
y ( t ) = n = 0 a n h ^ n α ( t , t 0 ) .
Using (34) and (36), we obtained
D , t 0 α C y ( t ) = D , t 0 α C n = 0 a n h ^ n α ( t , t 0 ) = D , t 0 α C a 0 + n = 1 a n h ^ n α ( t , t 0 ) = n = 1 a n h ^ ( n 1 ) α ( t , t 0 ) = n = 0 a n + 1 h ^ n α ( t , t 0 )
and
p ( t ) y ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) n = 0 a n h ^ n α ( t , t 0 ) = n = 1 k = 0 n 1 p n k h ^ ( n k ) α 1 ( t , t 0 ) a k h ^ k α ( t , t 0 ) = n = 1 k = 0 n 1 p n k a k h ^ n α ( t , t 0 ) .
Through substituting (43) and (44) into (41), we obtained
D , t 0 α C y ( t ) + p ( t ) y ( t ) = n = 0 a n + 1 h ^ n α ( t , t 0 ) + n = 1 k = 0 n 1 p n k a k h ^ n α ( t , t 0 ) = 0 .
Thus, we obtained the following recurrence formula, which allows us to express a n ( n > 0 ) , in terms of a 0 , as follows:
a 1 = 0 , a n + 1 = k = 0 n 1 p n k a k . n = 1 , 2 , · · · .
If we select a n via the above formula, then (42) is a solution of (41). The proof is thus completed. □

5.2. Solutions around an Ordinary Point of a Fractional Differential Equation of Order 2 α

In this section, we shall consider the solutions around an α -ordinary point t 0 Ω to equation
L 2 α y ( t ) : = y ( 2 α ( t ) + p ( t ) y ( α ( t ) + q ( t ) y ( t ) = 0 ,
where p ( t ) and q ( t ) are defined on Ω , and y ( 2 α and y ( α represent the Riemann–Liouville or Caputo sequential derivatives of order 2 α and α of the function y ( t ) .
Since t 0 is an α -ordinary point of (45), then
p ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) ( t [ t 0 , t 0 + ρ 1 ] T , ρ 1 > 0 )
and
q ( t ) = n = 1 q n h ^ n α 1 ( t , t 0 ) ( t [ t 0 , t 0 + ρ 2 ] T , ρ 2 > 0 ) .
Theorem 14.
Let α ( 0 , 1 ] and t 0 Ω be an α-ordinary point of equation
L 2 α y ( t ) : = D , t 0 2 α y ( t ) + p ( t ) D , t 0 α y ( t ) + q ( t ) y ( t ) = 0 .
Then, there exists a unique function that is formed as
y ( t ) = n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 )
for t ( t 0 , t 0 + ρ ) T with ρ = min { ρ 1 , ρ 2 } , which is a solution to Equation (48) and satisfies the following initial conditions:
lim t t 0 y h ^ α ( · , t 0 ) ( t ) = a 0
and
lim t t 0 D , t 0 α y h ^ α ( · , t 0 ) ( t ) = a 1 .
Proof. 
We sought the solution to Equation (48) in the form
y ( t ) = n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) .
Through calculating D , t 0 α y ( t ) and D , t 0 2 α y ( t ) and taking (34) into account, we obtained
D , t 0 α y ( t ) = D , t 0 α a 0 h ^ α 1 ( t , t 0 ) + n = 1 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) = n = 1 a n h ^ n α 1 ( t , t 0 )
and
D , t 0 2 α y ( t ) = D , t 0 α n = 1 a n h ^ n α 1 ( t , t 0 ) = D , t 0 α a 1 h ^ α 1 + n = 2 a n h ^ n α 1 ( t , t 0 ) = n = 2 a n h ^ ( n 1 ) α 1 ( t , t 0 ) = n = 0 a n + 2 h ^ ( n + 1 ) α 1 ( t , t 0 ) .
Through using (34), (46) and (47), we obtained
p ( t ) D , t 0 α y ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) n = 1 a n h ^ n α 1 ( t , t 0 ) = n = 2 k = 1 n 1 p n k h ^ ( n k ) α 1 ( t , t 0 ) a k h ^ k α 1 ( t , t 0 ) = n = 2 k = 1 n 1 p n k a k h ^ n α 1 ( t , t 0 ) = n = 1 k = 1 n p n + 1 k a k h ^ n α 1 ( t , t 0 )
and
q ( t ) y ( t ) = n = 1 q n h ^ n α 1 ( t , t 0 ) n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) = n = 1 k = 0 n 1 q n k h ^ ( n k ) α 1 ( t , t 0 ) a k h ^ ( k + 1 ) α 1 ( t , t 0 ) = n = 1 k = 0 n 1 q n k a k h ^ ( n + 1 ) α 1 ( t , t 0 ) .
When substituting these values into Equation (48), we have
D , t 0 2 α y ( t ) + p ( t ) D , t 0 α y ( t ) + q ( t ) y ( t ) = n = 0 a n + 2 h ^ ( n + 1 ) α 1 ( t , t 0 ) + n = 1 k = 1 n p n + 1 k a k h ^ n α 1 ( t , t 0 ) + n = 1 k = 0 n 1 q n k a k h ^ ( n + 1 ) α 1 ( t , t 0 ) = 0 .
Thus, we arrived at the following recurrence formula for the coefficients a n :
a n + 2 = k = 0 n p n k a k + 1 + q n k a k , a 2 = 0 .
Through using these a n , we can obtain
y ( t ) = n = 0 a n h ^ ( n + 1 ) α 1 ( t , t 0 ) ,
which is a solution of (48). The proof is thus completed. □
Theorem 15.
Let α ( 0 , 1 ] and t 0 Ω be an α-ordinary point of equation
L 2 α y ( t ) : = D , t 0 2 α C y ( t ) + p ( t ) D , t 0 α C y ( t ) + q ( t ) y ( t ) = 0 .
Then, for any given number a 0 , a 1 R , there exists a unique function that is formed as
y ( t ) = n = 0 a n h ^ n α ( t , t 0 )
for t ( t 0 , t 0 + ρ ) T with ρ = min { ρ 1 , ρ 2 } , which is a solution to Equation (49) and satisfies the following initial conditions:
a 2 = 0 , a n + 2 = k = 0 n 1 p n k a k + 1 + q n k a k , n = 1 , 2 , · · · .
Proof. 
We sought a solution to Equation (49) in the form y ( t ) = n = 0 a n h ^ n α ( t , t 0 ) .
When calculating D , t 0 α C y ( t ) and D , t 0 2 α C y ( t ) and taking (34) into account, we can obtain
D , t 0 α C y ( t ) = D , t 0 α C a 0 + n = 1 a n h ^ n α ( t , t 0 ) = n = 1 a n h ^ ( n 1 ) α ( t , t 0 )
and
D , t 0 2 α C y ( t ) = D , t 0 α C a 1 + n = 2 a n h ^ ( n 1 ) α ( t , t 0 ) = n = 2 a n h ^ ( n 2 ) α ( t , t 0 ) = n = 0 a n + 2 h ^ n α ( t , t 0 ) .
When using (34), (46) and (47), we obtain
p ( t ) D , t 0 α C y ( t ) = n = 1 p n h ^ n α 1 ( t , t 0 ) n = 1 a n h ^ ( n 1 ) α ( t , t 0 ) = n = 2 k = 1 n 1 p n k h ^ ( n k ) α 1 ( t , t 0 ) a k h ^ ( k 1 ) α ( t , t 0 ) = n = 2 k = 1 n 1 p n k a k h ^ ( n 1 ) α ( t , t 0 ) = n = 2 k = 0 n 2 p n k 1 a k + 1 h ^ ( n 1 ) α ( t , t 0 ) = n = 1 k = 0 n 1 p n k a k + 1 h ^ n α ( t , t 0 )
and
q ( t ) y ( t ) = n = 1 q n h ^ n α 1 ( t , t 0 ) n = 0 a n h ^ n α ( t , t 0 ) = n = 1 k = 0 n 1 q n k h ^ ( n k ) α ( t , t 0 ) a k h ^ k α ( t , t 0 ) = n = 1 k = 0 n 1 q n k a k h ^ n α ( t , t 0 ) .
This shows that (50) is a solution of (49) if and only if
D , t 0 2 α C y ( t ) + p ( t ) D , t 0 α C y ( t ) + q ( t ) y ( t ) = n = 0 a n + 2 h ^ n α ( t , t 0 ) + n = 1 k = 0 n 1 p n k a k + 1 h ^ n α ( t , t 0 ) + n = 1 k = 0 n 1 q n k a k h ^ n α ( t , t 0 ) = a 2 + n = 1 k = 0 n 1 p n k a k + 1 + q n k a k h ^ n α ( t , t 0 ) = 0 .
Thus, we arrived at the following recurrence formula for the coefficients a n :
a 2 = 0 , a n + 2 = k = 0 n 1 p n k a k + 1 + q n k a k , n = 1 , 2 , · · · .
When using the above a n , we can obtain (50), which is a solution of (49). The proof is thus completed. □

6. Applications and Conclusions

Fractional differential equations have a wide range of practical applications. For example, in reviewing Equation (20), which we discussed earlier, when t T = R , Equation (20) is the fractional-order equation of continuous variables. When t T = Z , then Equation (20) degenerates into the usual difference equation. When T selects different closed sets in R , our equation corresponds to different fractional-order equations or difference equations with different step lengths. Therefore, it can be seen that the results we obtained have a wide range of applications. In addition, fractional-order equations are also used in many specific fields in real life, including population models, infectious disease models, heat conduction models, wave equations, etc.
In recent years, an increasing number of researchers have focused on fractional-order infectious disease models [39,40,41,42]. The global epidemic in recent years has made the research on fractional-order COVID-19 [43,44,45] an exceptionally hot topic. The traveling wave solution characterizes the problem of the speed of transmission of infectious diseases during their spread. In recent years, the discussion of the spread and speed of transmission of infectious diseases has become a hot topic, especially since the outbreak of COVID-19, which has once again sounded an alarm for humanity. Inspired by [46,47], we present the following wave equation system, which corresponds to a fractional-order SIR infectious disease model:
c D α S ( t ) = d 1 D 2 α S ( t ) + Λ t β S ( t ) I ( t ) N ( t ) μ S ( t ) , c D α I ( t ) = d 2 D 2 α I ( t ) + β S ( t ) I ( t ) N ( t ) + ρ R ( t ) ( μ + ϕ ) I ( t ) , c D α R ( t ) = d 3 D 2 α R ( t ) + ϕ I ( t ) ρ R ( t ) μ R ( t ) , N ( t ) = S ( t ) + I ( t ) + R ( t ) ,
where t T , 0 < α 1 ; S represents susceptible individuals; I represents infected individuals; R represents recovered individuals; Λ represents total recruitment scale; c represents disease propagation speed; d 1 , d 2 , d 3 represent diffusion coefficient; β represents the infection rate; ρ represents the relapse rate; μ represents the natural mortality rate; ϕ represents the cure rate; and the total population is represented by N. If we choose d 1 = d 2 = d 3 = d , via adding the first three equations of (51), we can obtain
c D α S ( t ) + I ( t ) + R ( t ) = d D 2 α S ( t ) + I ( t ) + R ( t ) + Λ t μ S ( t ) + I ( t ) + R ( t ) .
By further organizing the fourth equation in (51), it can be concluded that
d D 2 α N ( t ) c D α N ( t ) μ N ( t ) = Λ t .
It is evident that the form of (52) is completely consistent with (20). Next, we used the relevant results from Section 4 to solve (52). With (21), we obtained the characteristic polynomial of (52)
P 2 λ = λ 2 c d λ μ d .
When P 2 λ = 0 , we can obtain the eigenvalues of (52) as follows:
λ 1 , 2 = c d ± 1 d c 2 + 4 d μ 2 .
According to Theorem 8, we can obtain the solution to the wave Equation (52)
N ( t ) = ( n = 0 c d + 1 d c 2 + 4 d μ 2 n I , t 0 . ( n + 1 ) α Λ ( t ) ) · ( n = 0 c d 1 d c 2 + 4 d μ 2 n I , t 0 . ( n + 1 ) α Λ ( t ) ) .
From the practical significance of infectious diseases, N ( t ) characterizes the speed of disease transmission in the overall population. We find from (53) that the value of N ( t ) only depends on the four parameters c , d , μ and Λ . Furthermore, (53) established clear relationships between the four variables, providing a mathematical basis for further quantitative research on the speed and control of disease transmission. Moreover, T is a general time scale, and by selecting different types of time scales, it can be used for the study of various intermittent, seasonal and discontinuous long-term epidemics. Due to space limitations, we will provide further research in a separate article.
The theory of fractional-order differential equations is a new topic providing many directions for further research. In this paper, we presented a general theory for sequential differential equations of a fractional order on time scales. We have given the explicit solution, in the case of constant coefficients, for both the homogeneous and the non-homogeneous problem using the ∇-Mittag-Leffler function, Laplace transform method, operational method and operational decomposition method. In addition, we proved the fact that the nabla generalized power function is positive, which is notable as this is also an open problem proposed in [37], which is also important for the studying of fractional-order differential equations on time scales. By introducing the concept of α -analyticity on time scales and in using the Laplace transform method, we studied the existence of a new class of fractional, sequential differential equations with convolution-type variable coefficients. We hope these results will provide a general and unified theory to study the continuous and discrete sequential differential equations of a fractional order on time scales.

Author Contributions

Methodology, J.Z.; Formal analysis, C.-C.Z. and J.Z.; Writing—original draft, J.Z.; Writing—review & editing, C.-C.Z. and J.Z.; Funding acquisition, C.-C.Z. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Jiangsu Province, China (grant no. BK20190578) and by the Jiangsu Province Colleges and Universities Undergraduate Scientific Research Innovative Program.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

We wish to thank the student Wu Ling for her hard work on the first draft of this article.

Conflicts of Interest

The authors declare that there are no conflicts of interests regarding the publication of this article.

References

  1. Fan, W.; Cheng, H. A novel numerical inverse technique for multi-parameter time fractional radially symmetric anomalous diffusion problem with initial singularity. Comput. Math. Appl. 2024, 158, 95–101. [Google Scholar] [CrossRef]
  2. Li, Z.; Zhang, Z. Stabilization Control for a Class of Fractional-Order HIV-1 Infection Model with Time Delays. Axioms 2023, 12, 695. [Google Scholar] [CrossRef]
  3. Malik, M.; Sajid, M.; Kumar, V. Controllability of singular dynamic systems on time scales. Asian J. Control 2022, 24, 2771–2777. [Google Scholar] [CrossRef]
  4. Zeng, C.; Chen, Y.Q. Optimal random search, fractional dynamics and fractional calculus. Fract. Calc. Appl. Anal. 2014, 17, 321–332. [Google Scholar] [CrossRef]
  5. Agarwal, R.P. Certain fractional q-integrals and q-derivatives. Proc. Camb. Philos. Soc. 1969, 66, 365–370. [Google Scholar] [CrossRef]
  6. Diaz, R.; Osler, T.J. Differences of fractional order. Math. Comput. 1974, 28, 185–202. [Google Scholar] [CrossRef]
  7. Guo, G.C.; Baleanu, D. Chaos synchronization of the discrete fractional logistic map. Signal Process. 2014, 102, 96–99. [Google Scholar]
  8. Guo, G.C.; Baleanu, D. Discrete chaos in fractional delayed logistic maps. Nonlinear Dyn. 2015, 80, 1697–1703. [Google Scholar]
  9. Guo, G.C.; Baleanu, D. Discrete fractional logistic map and its chaos. Nonlinear Dyn. 2014, 75, 283–287. [Google Scholar]
  10. Guo, G.C.; Baleanu, D.; Zeng, S.D. Discrete chaos in fractional sine and standard maps. Phys. Lett. A 2014, 378, 484–487. [Google Scholar]
  11. Atici, F.M.; Eloe, P.W. A Transform Method in Discrete Fractional Calculus. Inter. J. Diff. Equ. 2007, 2, 165–176. [Google Scholar]
  12. Atici, F.M.; Eloe, P.W. Initial value problems in discrete fractional calculus. Proc. Am. Math. Soc. 2009, 137, 981–989. [Google Scholar] [CrossRef]
  13. Atici, F.M.; Şengül, S. Modeling with fractional difference equations. J. Math. Anal. Appl. 2010, 369, 1–9. [Google Scholar] [CrossRef]
  14. Holm, M. Sum and Difference Compositions in Discrete Fractional Calculus. CUBO Math. J. 2011, 13, 153–184. [Google Scholar] [CrossRef]
  15. Anastassiou, G.A. Nabla discrete fractional calculus and nabla inequalities. Math. Comput. Model. 2010, 51, 562–571. [Google Scholar] [CrossRef]
  16. Atici, F.M.; Eloe, P.W. Discrecte fractional caculus with the nabla operator. Electron. J. Qual. Theory Differ. Equ. Spec. Ed. I 2009, 2009, 1–12. [Google Scholar]
  17. Hilger, S. Analysis on measure chains a unified approach to continuous and discrete calculas. Results Math. 1990, 18, 18–56. [Google Scholar] [CrossRef]
  18. Williams, P.A. Unifying Fractional Calculus with Time Scales. Doctoral Thesis, The University of Melbourne, Melbourne, Australia, 2012. [Google Scholar]
  19. Anastassiou, G.A. Foundations of nabla fractional calculus on time scales and inequalities. Comput. Math. Appl. 2010, 59, 3750–3762. [Google Scholar] [CrossRef]
  20. Bastos, N. Fractional Calculus on Time Scales. Doctoral Thesis, The University of Aveiro, Aveiro, Portugal, 2012. [Google Scholar]
  21. Anastassiou, G.A. Principles of delta fractional calculus on time scales and inequalities. Math. Comput. Model. 2010, 52, 556–566. [Google Scholar] [CrossRef]
  22. Bastos, N.; Ferreira, R.; Torres, D. Discrete-time fractional variational problems. Sign. Proc. 2011, 91, 513–524. [Google Scholar] [CrossRef]
  23. Holm, M.T. The Laplace transform in discrete fractional calculus. Comput. Math. Appl. 2011, 62, 1591–1601. [Google Scholar] [CrossRef]
  24. Wu, L.; Zhu, J. Fractional Cauchy problem with Riemann-Liouville derivative on time scales. Abstr. Appl. Anal. 2013, 2013, 795701. [Google Scholar] [CrossRef]
  25. Zhu, J.; Wu, L. Fractional Cauchy Problem with Caputo Nabla Derivative on Time Scales. Abstr. Appl. Anal. 2015, 2015, 486054. [Google Scholar] [CrossRef]
  26. Zhu, J.; Zhu, Y. Fractional Cauchy Problem with Riemann-Liouville Fractional Delta Derivative on Time Scales. Abstr. Appl. Anal. 2013, 2013, 401596. [Google Scholar] [CrossRef]
  27. Duque, C.; Leiva, H.; Tridane, A. Relative asymptotic equivalence of dynamic equations on time scales. Adv. Contin. Discret. Models 2022, 2022, 4. [Google Scholar] [CrossRef]
  28. Fang, J.J.; Mou, D.S.; Zhang, H.C.; Wang, Y.Y. Discrete fractional soliton dynamics of the fractional Ablowitz-Ladik model. Optik 2021, 228, 166186. [Google Scholar] [CrossRef]
  29. Kisela, T. On dynamical systems with Nabla half derivative on time scales. Mediterr. J. Math. 2020, 17, 187. [Google Scholar] [CrossRef]
  30. Liu, M.; Dong, H.; Fang, Y.; Dong, H. Modelling and analysis of dynamic systems on time-space scales and application in burgers equation. J. Appl. Anal. Comput. 2022, 12, 2555–2577. [Google Scholar] [CrossRef]
  31. Messina, E.; Raffoul, Y.; Vecchio, A. Qualitative analysis of dynamic equations on time scales using Lyapunov functions. Diff. Equ. Appl. 2022, 14, 215–226. [Google Scholar] [CrossRef]
  32. Wu, G.C.; Baleanu, D.; Zeng, S.D.; Deng, Z.G. Discrete fractional diffusion equation. Nonlinear Dyn. 2015, 80, 281–286. [Google Scholar] [CrossRef]
  33. Gogoi, B.; Saha, U.K.; Hazarika, B. Existence of solution of a nonlinear fractional dynamic equation with initial and boundary conditions on time scales. J. Anal. 2024, 32, 85–102. [Google Scholar] [CrossRef]
  34. Gogoi, B.; Hazarika, B.; Saha, U.K.; Tikare, S. Periodic boundary value problems for fractional dynamic equations on time scales. Results Math. 2023, 78, 228. [Google Scholar] [CrossRef]
  35. Gogoi, B.; Saha, U.K.; Hazarika, B.; Torres, D.F.; Ahmad, H. Nabla Fractional Derivative and Fractional Integral on Time Scales. Axioms 2021, 10, 317. [Google Scholar] [CrossRef]
  36. Morsy, A.; Nisar, K.S.; Ravichandran, C.; Anusha, C. Sequential fractional order Neutral functional Integro differential equations on time scales with Caputo fractional operator over Banach spaces. AIMS Math. 2023, 8, 5934–5949. [Google Scholar] [CrossRef]
  37. Kisela, T. Power functions and essentials of frational calculus on isolated time scales. Adv. Diff. Equ. 2013, 2013, 259. [Google Scholar] [CrossRef]
  38. Bohner, M.; Peterson, A. Advances in Dynamic Equations on Time Scales; Birkhäuser: Boston, MA, USA, 2003. [Google Scholar]
  39. Akdim, K.; Ez-Zetouni, A.; Zahid, M. The influence of awareness campaigns on the spread of an infectious disease: Aqualitative analysis of a fractional epidemic model. Model. Earth Syst. Environ. 2022, 8, 1311–1319. [Google Scholar] [CrossRef] [PubMed]
  40. Bounkaicha, C.; Allali, K. Dynamics of a time fractional order spatio-temporal SIR with vaccination and temporary immunity. Partial Diff. Equ. Appl. Math. 2023, 7, 100524. [Google Scholar] [CrossRef]
  41. Paul, S.; Mahata, A.; Mukherjee, S.; Roy, B. Dynamics of SIQR epidemic model with fractional order derivative. Partial Diff. Equ. Appl. Math. 2022, 5, 100216. [Google Scholar] [CrossRef] [PubMed]
  42. Sadki, M.; Harroudi, S.; Allali, K. Fractional-order SIR epidemic model with treatment cure rate. Partial Diff. Equ. Appl. Math. 2023, 8, 100593. [Google Scholar] [CrossRef]
  43. Barros, L.C.; Lopes, M.M.; Pedro, F.S.; Esmi, E.; Santos, J.P.; Sánchez, D.E. The memory effect on fractional calculus: An application in the spread of COVID-19. Comput. Appl. Math. 2021, 40, 72. [Google Scholar] [CrossRef]
  44. Hamou, A.A.; Azroul, E.; Alaoui, A.L. Fractional Model and Numerical Algorithms for Predicting COVID-19 with Isolation and Quarantine Strategies. Int. J. Appl. Comput. Math. 2021, 7, 142. [Google Scholar] [CrossRef] [PubMed]
  45. Khalaf, S.L.; Kadhim, M.S.; Khudair, A.R. Studying of COVID-19 fractional model: Stability analysis. Partial Diff. Equ. Appl. Math. 2023, 7, 100470. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, C.C.; Li, W.T.; Yang, F.Y. Traveling waves in a nonlocal dispersal SIRH model with relapse. Comput. Math. Appl. 2017, 73, 1707–1723. [Google Scholar] [CrossRef]
  47. Zhu, C.C.; Li, W.T.; Yang, F.Y. Traveling waves of a reaction-diffusion SIRQ epidemic model with relapse. J. Appl. Anal. Comput. 2017, 7, 147–171. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, C.-C.; Zhu, J. Fractional-Order Sequential Linear Differential Equations with Nabla Derivatives on Time Scales. Axioms 2024, 13, 447. https://doi.org/10.3390/axioms13070447

AMA Style

Zhu C-C, Zhu J. Fractional-Order Sequential Linear Differential Equations with Nabla Derivatives on Time Scales. Axioms. 2024; 13(7):447. https://doi.org/10.3390/axioms13070447

Chicago/Turabian Style

Zhu, Cheng-Cheng, and Jiang Zhu. 2024. "Fractional-Order Sequential Linear Differential Equations with Nabla Derivatives on Time Scales" Axioms 13, no. 7: 447. https://doi.org/10.3390/axioms13070447

APA Style

Zhu, C. -C., & Zhu, J. (2024). Fractional-Order Sequential Linear Differential Equations with Nabla Derivatives on Time Scales. Axioms, 13(7), 447. https://doi.org/10.3390/axioms13070447

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