Next Article in Journal
Limit Property of an L2-Normalized Solution for an L2-Subcritical Kirchhoff-Type Equation with a Variable Exponent
Next Article in Special Issue
Energy-Conserving Explicit Relaxed Runge–Kutta Methods for the Fractional Nonlinear Schrödinger Equation Based on Scalar Auxiliary Variable Approach
Previous Article in Journal
On the Čech-Completeness of the Space of τ-Smooth Idempotent Probability Measures
Previous Article in Special Issue
A Note on a Fractional Extension of the Lotka–Volterra Model Using the Rabotnov Exponential Kernel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Investigation of Fractional Estimations in Liouville–Caputo Operators of Mixed Order with Applications

by
Pshtiwan Othman Mohammed
1,2,*,
Alina Alb Lupas
3,*,
Ravi P. Agarwal
4,
Majeed A. Yousif
5,
Eman Al-Sarairah
6,7 and
Mohamed Abdelwahed
8
1
Department of Mathematics, College of Education, University of Sulaimani, Sulaymaniyah 46001, Iraq
2
Research and Development Center, University of Sulaimani, Sulaymaniyah 46001, Iraq
3
Department of Mathematics and Computer Science, University of Oradea, 410087 Oradea, Romania
4
Department of Mathematics and Systems Engineering, Florida Institute of Technology, Melbourne, FL 32901, USA
5
Department of Mathematics, College of Education, University of Zakho, Duhok 42001, Iraq
6
Department of Mathematics, Khalifa University of Science and Technology, Abu Dhabi P.O. Box 127788, United Arab Emirates
7
Department of Mathematics, Al-Hussein Bin Talal University, P.O. Box 20, Ma’an 71111, Jordan
8
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Axioms 2024, 13(8), 570; https://doi.org/10.3390/axioms13080570
Submission received: 10 July 2024 / Revised: 9 August 2024 / Accepted: 19 August 2024 / Published: 21 August 2024
(This article belongs to the Special Issue Fractional Calculus - Theory and Applications II)

Abstract

:
In this study, to approximate nabla sequential differential equations of fractional order, a class of discrete Liouville–Caputo fractional operators is discussed. First, some special functions are re-called that will be useful to make a connection with the proposed discrete nabla operators. These operators exhibit inherent symmetrical properties which play a crucial role in ensuring the consistency and stability of the method. Next, a formula is adopted for the solution of the discrete system via binomial coefficients and analyzing the Riemann–Liouville fractional sum operator. The symmetry in the binomial coefficients contributes to the precise approximation of the solutions. Based on this analysis, the solution of its corresponding continuous case is obtained when the step size p 0 tends to 0. The transition from discrete to continuous domains highlights the symmetrical nature of the fractional operators. Finally, an example is shown to testify the correctness of the presented theoretical results. We discuss the comparison of the solutions of the operators along with the numerical example, emphasizing the role of symmetry in the accuracy and reliability of the numerical method.

1. Introduction

Since the 1990s, fractional calculus (in both continuous and discrete versions) has become an outstanding topic in applied mathematics research. The research work on fractional calculus has continued to deepen with the advance in technology (see earlier studies [1,2,3,4,5,6]). In the last few decades, many authors have studied fractional differential equations and fractional systems formed by Riemann–Liouville (R-L) or Liouville–Caputo (L-C) operators. In addition, these studies focus on the existence and uniqueness of their solutions including initial and boundary conditions (see [7,8,9,10,11]).
Discrete fractional operators arising from general fractional calculus arise in a wide range of applications and they have been of continued interest in applied mathematics for over a century. In very recent years, discrete fractional operators have been widely used to discretize fractional order models for both linear and nonlinear systems of equations. These discretizations are often expressed via fractional differences and sums (see, e.g., [5,12]). Numerous applications of discrete fractional equations may also be found in modeling mathematical issues, such as mathematical analysis [13,14], mathematical physics [15,16,17], uncertainty theory [18,19], stability analysis [20,21], and monotonicity and positivity analyses [22,23,24,25,26].
Based on the current research articles, it is known that effective and important operators for analyzing the stability, existence, and uniqueness of fractional order systems are R-L and L-C fractional operators (see, for example, [27,28,29,30]). Of relevant interest are the steady states of such systems, where non-integer orders may appear, playing an important role in the behavior of the solutions of the R-L or L-C problems (see, for example, [31,32,33,34,35,36] and the references therein).
Recently, a fractional difference technique was used by Mozyrska et al. [37] to solve a fractional differential problem of delta L-C type with specific initial value conditions:
Δ ( α 1 ) p 0 L C p 0 α w ( n p 0 ) = y ( n p 0 + ( β 1 ) p 0 ) , Δ ( β 1 ) p 0 L C p 0 β y ( n p 0 ) = f n p 0 , w ( n p 0 + ( α 1 ) p 0 ) ,
with
Δ ( α 1 ) p 0 L C p 0 α w ( 0 ) = w 0 , w ( ( α 1 ) p 0 ) = w 1 ,
where α , β ( 0 , 1 ] , p 0 > 0 , w 0 and w 1 are constant vectors in R .
However, the basic concepts of solving nabla L-C fractional problems still lack adequate research. Thus, in this article, we present a numerical solution of a nabla L-C fractional system by using the nabla fractional difference technique.
The innovative contents of this article are summarized below. Section 2 briefly reviews the RL and L-C fractional operators and presents some useful existing concepts. Section 3 contains the main contribution of this article, where an L-C approximation algorithm based on fractional derivatives is introduced, and the foundation of its solution and uniqueness is examined by using the Lipschitz condition (LiC). Numerical problems are implemented to illustrate the efficiency and accuracy of the proposed scheme in Section 4.

2. Preliminaries

In this section, we present a comprehensive introduction to R-L fractional operators of order α > 0 and some related properties. Throughout the article, we suppose that x 1 ( n p 0 ) = n p 0 p 0 and x 2 ( n p 0 ) = n p 0 + p 0 .
Definition 1
(see [1,38]). The R-L integral operator is defined by
I a α f ( t ) : = 1 Γ ( α ) a t ( t r ) α 1 f ( r ) d r for a < t ,
and its discrete version is given by
a p 0 α f ( t ) : = p 0 Γ ( α ) r = a p 0 + 1 t p 0 t + p 0 r p 0 p 0 α 1 ¯ f ( r p 0 ) for t N a + p 0 , p 0 ,
where N a , p 0 and t p 0 α ¯ are given by
N a , p 0 : = { a , a + p 0 , a + 2 p 0 , } ,
t p 0 α ¯ = p 0 α Γ t p 0 + α Γ t p 0 = Γ ( α + 1 ) t p 0 + α 1 t p 0 1 p 0 α ,
For t p 0 , α R such that neither t p 0 + α nor t p 0 is a pole of Γ. Furthermore, we recall the binomial formula
β α : = Γ ( β + 1 ) Γ ( α + 1 ) Γ ( β α + 1 ) .
Remark 1.
By considering Definition 1, it can be noted that
a p 0 α f ( t ) = p 0 α Γ ( α ) s = 0 n Γ ( n s + α ) Γ ( n s + 1 ) f ( a + s p 0 + p 0 ) = p 0 α s = 0 n n s ( α ) 1 n s f ( a + s p 0 + p 0 ) = p 0 α s 0 = 0 n ( 1 ) s 0 α s 0 f ( a + p 0 s 0 p 0 ) ,
For t = a + ( n + α ) p 0 , n N 0 , and it is seen that (see [39])
s 0 α 1 s 0 = ( 1 ) s 0 α s 0 ,
The following definitions present the R-L fractional operators associated with (1) and (2), respectively.
Definition 2
(see [1,23,40]). For α [ 0 , 1 ) , the R-L fractional derivative is defined by
D a R L α f ( t ) = 1 Γ ( 1 α ) d d t a t ( t r ) α f ( r ) d r for t > a ,
and the nabla R-L fractional difference is defined by
a RL p 0 α f ( t ) = p 0 Γ ( α ) r = a p 0 + 1 t p 0 ( t + p 0 r p 0 ) p 0 α 1 ¯ f ( r p 0 ) ,
For t N a + m p 0 , p 0 , where m 1 < α < m , m N 1 . It is worth noting that
t p 0 α ¯ = Γ t p 0 + α Γ t p 0 p 0 α ,
such that t p 0 α ¯ 0 when Γ t p 0 is undefined.
Specifically, the nabla difference operator can be expressed as follows:
p 0 f ( t ) = 1 p 0 f ( t ) f ( t p 0 ) for t N a + p 0 , p 0 .
Below, we recall the corresponding L-C derivative and difference associated with (1) and (2), respectively.
Definition 3
(see [1,23,40] [Lemma 2.3]). Assume that α [ 0 , 1 ) . Then, the L-C fractional derivative can be written as
D a L C α f ( t ) = a t ( t r ) α Γ ( 1 α ) d d r f ( r ) d r = lim p 0 0 1 p 0 α j = 0 t a p 0 + 1 ( 1 ) j α j f ( t j p 0 ) ( t a ) α Γ ( 1 α ) f ( a ) .
And the nabla Liouville–Caputo p 0 fractional difference can be written as
a L C p 0 α f ( t ) = p 0 Γ ( α ) r = a p 0 + 1 t p 0 ( t + p 0 r p 0 ) p 0 α 1 ¯ f ( r p 0 ) ( t a ) p 0 α ¯ Γ ( 1 α ) f ( a ) ,
for t in N a + p 0 , p 0 .
Lemma 1
(see [23] [Property 2.1]). The following power difference formula can be deduced:
p 0 t p 0 α ¯ = α t p 0 α 1 ¯ ,
for non-negative values of α , p 0 , and t in N 0 , p 0 .
Lemma 2
(see [5] [Theorem 3.93]). For α > 0 and β R , one can have
a 0 p 0 α ( t a 0 ) p 0 β ¯ = Γ ( β + 1 ) Γ ( β + 1 + α ) ( t a 0 ) p 0 β + α ¯ ,
and
a 0 R L p 0 α ( t a 0 ) p 0 β ¯ = Γ ( β + 1 ) Γ ( β + 1 α ) ( t a 0 ) p 0 β α ¯ ,
where t N a 0 , p 0 , and β + α , β α 0 .
Lemma 3
(see [5] [Theorem 3.107]). Let f be defined on N a 0 + p 0 , p 0 with α , h , β > 0 . Then,
a 0 p 0 α a 0 p 0 β f ( t ) = a 0 p 0 ( α + β ) f ( t ) = a 0 p 0 β a 0 p 0 α f ( t ) ,
where t N a 0 , p 0 .
Lemma 4.
Suppose that f is defined on N a 0 , p 0 , and 0 < α 1 , p 0 > 0 . Then, we have
a 0 p 0 α a 0 L C p 0 α f ( t ) = f ( t ) f ( a 0 ) .
Proof. 
This proof can be deduced from the following identity (see [38] [Proposition 6]):
a 0 p 0 α a 0 R L p 0 α f ( t ) = f ( t ) ,
using Definition (11) and identity (12). □
Next, we consider binomial special functions, which are defined in [41].
Definition 4
(see [41]). For α , β > 0 , we define ψ κ , s and ψ ˜ κ , s as follows:
ψ κ , s ( n p 0 ) : = n κ + α κ + β s n κ p 0 α κ + β s , when n N κ , 0 , when n N κ ,
and
ψ ˜ κ , s ( n p 0 ) : = n + μ 1 n p 0 β , for n N 0 , 0 , for n N 0 ,
for n , κ , s N 0 , α , β > 0 and μ = α κ + β s .
Remark 2
(see [41]). The above special function has some major properties:
(i)
ψ 0 , 0 ( n p 0 ) = 1 ;
(ii)
ψ 1 , 0 ( n p 0 ) = 0 p 0 α 1 ( n p 0 ) = n + α 1 n 1 p 0 α
and ψ 0 , 1 ( n p 0 ) = 0 p 0 β 1 ( n p 0 ) = n + β 1 n 1 p 0 β ;
(iii)
ψ 2 , 0 ( n p 0 ) = 0 p 0 α 1 ( n p 0 p 0 ) = n + α 2 n 2 p 0 α ;
(iv)
ψ κ , s ( n p 0 ) = ( n κ + 1 ) p 0 α κ + β s ¯ Γ ( α κ + β s + 1 ) .
(v)
For 0 < α , β 1 , and p 0 > 0 , we have
0 p 0 α ψ κ , s ( x 1 ( n p 0 ) ) = ψ κ + 1 , s ( n p 0 ) ,
and
0 p 0 β ψ κ , s ( n p 0 ) = ψ κ , s + 1 ( n p 0 ) .
(vi)
It can be noted, for n = 0 , that
ψ κ , s ( 0 ) = 0 ,
for κ > 0 .

3. Solution of Difference Systems

In this section, our target is to approximate the following system of difference equations:
0 L C p 0 α w ( x 1 ( n p 0 ) ) = y ( x 2 ( n p 0 ) ) , 0 L C p 0 β y ( n p 0 ) = f x 2 ( n p 0 ) , w ( x 1 ( n p 0 ) ) ,
subject to the conditions
0 L C p 0 α w ( 2 p 0 ) = w 0 , w ( p 0 ) = w 1 ,
where α , β ( 0 , 1 ] , p 0 > 0 , w 0 and w 1 are constant vectors in R .
Theorem 1.
System (18) with (19) has the solution
w ( x 1 ( n p 0 ) ) = w 1 + ψ 2 , 0 ( n p 0 ) w 0 + 0 p 0 α g ( x 1 ( n p 0 ) ) ,
where g ( n p 0 ) = 0 p 0 β F ( x 2 ( n p 0 ) ) and F ( x 2 ( n p 0 ) ) = f x 2 ( n p 0 ) , w ( x 1 ( n p 0 ) ) , for n N 0 .
Proof. 
By considering Lemma 4 with t = x 2 ( n p 0 ) , x 1 ( n p 0 ) , one can have
0 p 0 β 0 L C p 0 β y ( x 2 ( n p 0 ) ) = y ( x 2 ( n p 0 ) ) y ( 0 ) = 0 L C p 0 α w ( x 1 ( n p 0 ) ) w 0 ,
and
0 p 0 α 0 L C p 0 α w ( x 1 ( n p 0 ) ) = w ( x 1 ( n p 0 ) ) w 1 .
It follows from (21) that
0 L C p 0 α w ( x 1 ( n p 0 ) ) = w 0 + 0 p 0 β F ( x 2 ( n p 0 ) ) ,
and it follows from this, Remark 2 (ii), and (22) that
w ( x 1 ( n p 0 ) ) = w 1 + 0 p 0 α w 0 + 0 p 0 β F ( x 2 ( n p 0 ) ) ( x 1 ( n p 0 ) ) = w 1 + 0 p 0 α 1 ( x 1 ( n p 0 ) ) w 0 + 0 p 0 α g ( x 1 ( n p 0 ) ) = w 1 + n 2 + α n 2 p 0 α w 0 + 0 p 0 α g ( n p 0 ) = w 1 + ψ 2 , 0 ( n p 0 ) w 0 + 0 p 0 α g ( x 1 ( n p 0 ) ) ,
where g ( n p 0 ) = 0 p 0 β F ( x 2 ( n p 0 ) ) and F ( x 2 ( n p 0 ) ) = f x 2 ( n p 0 ) , w ( x 1 ( n p 0 ) ) , for n N 0 . Hence, the proof is completed. □
Remark 3.
Following [42] [Theorem 4.1], we can deduce that system (18) has a unique solution (i.e., (20) is unique) such that y and f satisfy the Lipschitz condition, i.e.,
y ( n 1 ) y ( n 2 ) L 1 n 1 n 2 ,
and
f ( n , w ) f ( n 2 , w ˜ ) L 2 w w ˜ ,
for some positive constant L 1 and L 2 .
In the following theorem, we examine the solution of a particular case of system (18).
Theorem 2.
Systems (18) and (19) have the following solution:
w ( x 1 ( n p 0 ) ) = κ = 0 A κ ψ κ , κ w 1 + ψ κ + 2 , κ w 0 ( n p 0 ) ,
where f x 1 ( n p 0 ) , w ( x 1 ( n p 0 ) ) = A w ( x 1 ( n p 0 ) ) , and n N 0 .
Proof. 
First, we try to investigate (23). In the case of n = 0 , it follows from (20) that
w ( p 0 ) = w 1 + ψ 2 , 0 ( 0 ) w 0 + 0 p 0 α g ( p 0 ) = w 1 ,
which is equivalent to
w ( p 0 ) = ψ 0 , 0 w 1 + ψ 2 , 0 w 0 ( 0 ) + κ = 1 A κ ψ κ , κ w 1 + ψ κ + 1 , κ w 0 ( 0 ) = 0 ,
according to Remark 2 (vi) and the convention that
0 p 0 α g ( p 0 ) = p 0 Γ ( α ) r = 1 1 r p 0 p 0 α 1 ¯ g ( r p 0 ) = 0 .
Next, for n > 0 , we define
w m + 1 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + 0 p 0 α g m ( n p 0 ) ,
where g m ( n p 0 ) = 0 p 0 β F m ( x 2 ( n p 0 ) ) and F m ( x 2 ( n p 0 ) ) = A w m ( x 1 ( n p 0 ) ) , for m N 0 .
Then, we see that
w 1 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + 0 p 0 α g 0 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + A w 1 0 p 0 α 0 p 0 β 1 ( n p 0 ) ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + A w 1 0 p 0 α ψ 0 , 1 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + A ψ 1 , 1 ( n p 0 ) w 1 ,
where we use Remark 2 (iii), (v), and
g 0 ( n p 0 ) = 0 p 0 β F 0 ( x 2 ( n p 0 ) ) = A w 0 ( x 1 ( n p 0 ) ) 0 p 0 β 1 ( x 2 ( n p 0 ) ) = A w 1 ψ 0 , 1 ( n p 0 ) .
Using the same technique together with using Remark 2 (iii) and (v), we can deduce
w 2 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + 0 p 0 α g 1 ( x 1 ( n p 0 ) ) = ψ 0 , 0 ( n p 0 ) w 1 + ψ 2 , 0 ( n p 0 ) w 0 + A ψ 1 , 1 ( n p 0 ) w 1 + ψ 3 , 1 ( n p 0 ) w 0 + A 2 ψ 2 , 2 ( n p 0 ) w 1 ,
where
g 1 ( n p 0 ) = 0 p 0 β F 1 ( x 2 ( n p 0 ) ) = A w 1 0 p 0 β ψ 0 , 0 ( n p 0 ) + w 0 0 p 0 β ψ 2 , 0 ( n p 0 ) + A w 1 0 p 0 β ψ 1 , 1 ( n p 0 ) = A ψ 0 , 1 ( n p 0 ) w 1 + ψ 2 , 1 ( n p 0 ) w 0 + A ψ 1 , 2 ( n p 0 ) w 1 .
By the same process, we can obtain the desired result in (23). □
The following lemmas are useful to obtain our main result.
Lemma 5.
If 0 < α 1 , then the following summation formula can be obtained:
lim p 0 0 t p 0 a α ¯ = ( t a ) α ,
for p 0 > 0 , t > a , where a R , t p 0 = a + p 0 + n p 0 , and n = t a p 0 + 1 .
Proof. 
Consider
t p 0 a α ¯ = Γ ( n + 1 α ) Γ ( n + 1 ) p 0 α = Γ ( n + 1 α ) Γ ( n + 1 ) t p 0 a α ( n + 1 ) α .
As p 0 0 (or n ), it follows that
lim p 0 0 t p 0 a α ¯ = lim p 0 0 Γ ( n + 1 α ) Γ ( n + 1 ) t p 0 a α ( n + 1 ) α = lim p 0 0 t p 0 a α · lim n Γ ( n + 1 α ) Γ ( n + 1 ) ( n + 1 ) α ( t a ) α ,
where it is seen that
lim p 0 0 t p 0 a = lim p 0 0 p 0 + ( t a ) + p 0 = t a ,
and the asymptotic formula is (see [39])
lim n Γ ( n + 1 α ) Γ ( n + 1 ) ( n + 1 ) α = 1 + O 1 n + 1 1 .
This proves the result. □
Lemma 6.
If α ( 0 , 1 ] , then the following summation formula holds:
a R L p 0 α w ( t p 0 ) = p 0 α j = 0 n ( 1 ) j α j w t p 0 j p 0 ,
for t > a and p 0 > 0 with a R , t p 0 = a + p 0 + n p 0 , and n N 0 .
Proof. 
Considering Definition (8), one can have
a R L p 0 α w ( t p 0 ) = 1 Γ ( α ) m = 0 t a p 0 1 t p 0 a m p 0 p 0 α 1 ¯ p 0 w ( a + p 0 + m p 0 ) = p 0 m = 0 n n p 0 + p 0 m p 0 p 0 α 1 ¯ Γ ( α ) w ( a + p 0 + m p 0 ) = p 0 α m = 0 n Γ ( n m α ) Γ ( n m + 1 ) Γ ( α ) w ( a + p 0 + m p 0 ) = p 0 α m = 0 n n m α 1 n m w ( a + p 0 + m p 0 ) = p 0 α j = 0 n j α 1 j w ( a + x 2 ( n p 0 ) j p 0 ) = p 0 α j = 0 n ( 1 ) j α j w ( a + x 2 ( n p 0 ) j p 0 ) = p 0 α j = 0 n ( 1 ) j α j w ( t p 0 j p 0 ) ,
where we used (6). Thus, the proof is completed. □
Theorem 3.
If 0 < α 1 , w is a continuous function, and w is integrable on [ a , Λ ] with Λ > 0 , then, for p 0 > 0 and t > a with a R , one can have
D a L C α w ( t ) = lim p 0 0 a L C p 0 α w ( t p 0 ) ,
where t p 0 = a + p 0 + n p 0 and n N 0 .
Proof. 
By using (25) in (11), we can deduce
a L C p 0 α w ( t p 0 ) = p 0 α j = 0 n ( 1 ) j α j w t p 0 j p 0 t p 0 a p 0 α ¯ Γ ( 1 α ) w ( a ) ,
for t p 0 = a + p 0 + n p 0 . As p 0 0 on both sides of (26), we obtain the result by using Lemma 5. □
Proposition 1.
If α ( 0 , 1 ] and p 0 > 0 , then the solution of
D 0 L C α w ( t ) = f t , w ( t , w ( 0 ) = w 0 ,
can be estimated with the solution of
a L C α w ˜ ( t ) = f t , w ˜ ( t ) , w ˜ ( a ) = w 0 ,
via the following limit:
lim p 0 0 w ˜ t p 0 = w ( t ) ,
where t p 0 = a + p 0 + n p 0 and n = ( t p 0 ) / p 0 + 1 .
Proposition 2.
If α ( 0 , 1 ] and p 0 > 0 , then the solution of
D 0 L C β D 0 L C α w ( t ) = f t , w ( t ) , w ( 0 ) = w 0 ,
can be approximated with the solutions of (18) and (19) via the limit
lim p 0 0 w ˜ t p 0 = w ( t ) ,
where n and t p 0 are given in Proposition 1.

4. Numerical Tests

In this section, we will present the result of the numerical experiment in [40] to demonstrate the effectiveness of the proposed technique.
For this reason, we consider
D 0 L C α w ( t ) = w ( t ) with w ( 0 ) = 1 ,
where its solution is given in [40] by
w ( t ) = E t α : = κ = 0 t α κ Γ ( α κ + 1 ) .
To compare the efficiency of the methods, we report, in Figure 1 and Figure 2, the actual computational results ( w ) compared with the numerical results ( w ˜ ) in the cases of α = 1 2 and t = 0 , 0.1 , , 1 for different values of p 0 and Λ .
Furthermore, we know that
D 0 L C α D 0 L C α w ( t ) = w ( t ) , D 0 L C α w ( 0 ) = 1 , w ( 0 ) = 1 ,
where its exact solution is w ( t ) = E t α .
As a result, by considering Figure 1 and Figure 2 and Table 1 and Table 2, we conclude that the exact solution can be obtained from w ˜ such that p 0 is sufficiently small. This implies the applicability of Theorem 3 where p 0 0 .

5. Concluding Remarks and Future Works

To summarize, we investigated the complete solution of a sequential fractional differential problem of Liouville–Caputo type including certain initial value conditions,
D 0 L C β D 0 L C α w ( t ) = f t , w ( t ) ,
with the condition w ( 0 ) = w 1 , by analyzing the corresponding system of the fractional difference,
0 L C p 0 α w ˜ ( x 1 ( n p 0 ) ) = y ( x 2 ( n p 0 ) ) , 0 L C p 0 α y ( n h + h ) = f x 1 ( n p 0 ) , w ˜ ( x 1 ( n p 0 ) ) ,
subject to the conditions
0 L C p 0 α w ˜ ( p 0 ) = w ˜ 0 , w ˜ ( 2 p 0 ) = w ˜ a ,
via the following limitation:
lim p 0 0 w ˜ t p 0 = w ( t ) ,
where t p 0 and n are defined in Proposition 1. We used the expressions ψ κ , s and ψ ˜ κ , s in the solution of the discrete fractional difference system as stated in Definition 4.
Furthermore, the numerical example shows that the new reconstruction method can obtain a very high order of accuracy; when the value of p 0 is small enough (or Λ is large enough), the approximation results are in good agreement with the exact one.
Concerning future works, we plan to extend the numerical schemes presented here to apply them to other types of discrete fractional systems including exponential and Mittag–Leffler systems in kernels to establish similar results (visit [22,38] to see these operators).

Author Contributions

Conceptualization, R.P.A.; Formal analysis, M.A.Y.; Funding acquisition, A.A.L.; Investigation, E.A.-S.; Methodology, P.O.M.; Project administration, R.P.A.; Software, M.A.Y. and E.A.-S.; Supervision, P.O.M.; Visualization, M.A.; Writing—original draft, M.A.; Writing—review and editing, A.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

The publication of this research was supported by the University of Oradea, Romania.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors acknowledge Researchers Supporting Project number (RSP2024R136), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  2. Cinar, M.; Secer, A.; Ozisik, M.; Bayram, M. Derivation of optical solitons of dimensionless Fokas-Lenells equation with perturbation term using Sardar sub-equation method. Opt. Quant. Electron 2022, 54, 402. [Google Scholar]
  3. Ehsan, H.; Abbas, M.; Nazir, T.; Mohammed, P.O.; Chorfi, N.; Baleanu, D. Efficient analytical algorithms to study Fokas dynamical models involving M-truncated derivative. Qual. Theory Dyn. Syst. 2024, 23, 49. [Google Scholar]
  4. Dos Santos, J.P.C.; Arjunan, M.M.; Cuevas, C. Existence results for fractional neutral integro-differential equations with state-dependent delay. Comput. Math. Appl. 2011, 62, 1275–1283. [Google Scholar]
  5. Goodrich, C.S.; Peterson, A.C. Discrete Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  6. Korpinar, Z.; Inc, M.; Baleanu, D.; Bayram, M. Theory and application for the time fractional Gardner equation with Mittag-Leffler kernel. J. Taibah Univ. Sci. 2019, 13, 813–819. [Google Scholar]
  7. Hashemi, M.S.; Inc, M.; Bayram, M. Symmetry properties and exact solutions of the time fractional Kolmogorov-Petrovskii-Piskunov equation. Rev. Mex. Física 2019, 65, 529–535. [Google Scholar]
  8. Chadha, A.; Pandey, D.N. Existence of mild solutions for a fractional equation with state-dependent delay via resolvent operators. Nonlinear Stud. 2015, 22, 71–85. [Google Scholar]
  9. Phuong, N.D.; Sakar, F.M.; Etemad, S.; Rezapour, S. A novel fractional structure of a multi-order quantum multi-integro-differential problem. Adv. Differ. Equ. 2020, 2020, 633. [Google Scholar]
  10. Rezapour, S.; Imran, A.; Hussain, A.; Martínez, F.; Etemad, S.; Kaabar, M.K.A. Condensing Functions and Approximate Endpoint Criterion for the Existence Analysis of Quantum Integro-Difference FBVPs. Symmetry 2021, 13, 469. [Google Scholar] [CrossRef]
  11. Etemad, S.; Matar, M.M.; Ragusa, M.A.; Rezapour, S. Tripled Fixed Points and Existence Study to a Tripled Impulsive Fractional Differential System via Measures of Noncompactness. Mathematics 2022, 10, 25. [Google Scholar]
  12. Atici, F.; Uyanik, M. Analysis of discrete fractional operators. Appl. Anal. Discret. Math. 2015, 9, 139–149. [Google Scholar]
  13. Atici, F.M.; Atici, M.; Belcher, M.; Marshall, D. A new approach for modeling with discrete fractional equations. Fund. Inform. 2017, 151, 313–324. [Google Scholar]
  14. Atici, F.; Sengul, S. Modeling with discrete fractional equations. J. Math. Anal. Appl. 2010, 369, 1–9. [Google Scholar]
  15. Goodrich, C.S. On discrete sequential fractional boundary value problems. J. Math. Anal. Appl. 2012, 385, 111–124. [Google Scholar]
  16. Noureen, R.; Naeem, M.N.; Baleanu, D.; Mohammed, P.O.; Almusawa, M.Y. Application of trigonometric B-spline functions for solving Caputo time fractional gas dynamics equation. AIMS Math. 2023, 8, 25343–25370. [Google Scholar]
  17. Silem, A.; Wu, H.; Zhang, D.-J. Discrete rogue waves and blow-up from solitons of a nonisospectral semi-discrete nonlinear Schrödinger equation. Appl. Math. Lett. 2021, 116, 107049. [Google Scholar]
  18. Liu, Y.; Liu, H.; Zhu, Y. An approach for numerical solutions of Caputo-Hadamard uncertain fractional differential equations. Fractal Fract. 2022, 6, 693. [Google Scholar] [CrossRef]
  19. Lu, Q.; Zhu, Y. Comparison theorems and distributions of solutions to uncertain fractional difference equations. J. Comput. Appl. Math. 2020, 376, 112884. [Google Scholar]
  20. Chen, C.R.; Bohner, M.; Jia, B.G. Ulam-Hyers stability of Caputo fractional difference equations. Math. Meth. Appl. Sci. 2019, 42, 7461–7470. [Google Scholar]
  21. Lizama, C. The Poisson distribution, abstract fractional difference equations, and stability. Proc. Am. Math. Soc. 2017, 145, 3809–3827. [Google Scholar]
  22. Mohammed, P.O.; Abdeljawad, T.; Hamasalh, F.K. On discrete delta Caputo-Fabrizio fractional operators and monotonicity analysis. Fractal Fract. 2021, 5, 116. [Google Scholar] [CrossRef]
  23. Liu, X.; Du, F.; Anderson, D.; Jia, B. Monotonicity results for nabla fractional h-difference operators. Math. Meth. Appl. Sci. 2021, 44, 1207–1218. [Google Scholar]
  24. Abdeljawad, T.; Baleanu, D. Monotonicity analysis of a nabla discrete fractional operator with discrete Mittag-Leffler kernel. Chaos Solit. Fractals 2017, 116, 1–5. [Google Scholar]
  25. Goodrich, C.S.; Lizama, C. Positivity, monotonicity, and convexity for convolution operators. Discret. Contin. Dyn. Syst. 2020, 40, 4961–4983. [Google Scholar]
  26. Wu, G.-C.; Baleanu, D.; Zeng, S.-D. Discrete chaos in fractional sine and standard maps. Phys. Lett. A 2014, 378, 484–487. [Google Scholar]
  27. He, J.W.; Zhang, L.; Zhou, Y.; Ahmad, B. Existence of solutions for fractional difference equations via topological degree methods. Adv. Differ. Equ. 2018, 2018, 153. [Google Scholar]
  28. Huang, L.-L.; Park, J.H.; Wu, G.-C.; Mo, Z.-W. Variable-order fractional discrete-time recurrent neural networks. J. Comput. Appl. Math. 2020, 370, 112633. [Google Scholar]
  29. Mohammed, P.O.; Baleanu, D.; Al-Sarairah, E.; Abdeljawad, T.; Chorfi, N. Theoretical and numerical computations of convexity analysis for fractional differences using lower boundedness. Fractals 2023, 31, 2340183. [Google Scholar]
  30. Mohammed, P.O.; Almusawa, M.Y. On analysing discrete sequential operators of fractional order and their monotonicity results. AIMS Math. 2023, 8, 12872–12888. [Google Scholar]
  31. Segi Rahmat, M.R.; Md Noorani, M.S. Caputo type fractional difference operator and its application on discrete time scales. Adv. Differ. Equ. 2015, 2015, 160. [Google Scholar]
  32. Almusawa, M.Y.; Mohammed, P.O. Approximation of sequential fractional systems of Liouville-Caputo type by discrete delta difference operators. Chaos Solit. Fractals 2023, 176, 114098. [Google Scholar]
  33. Fečkan, M.; Pospíšil, M.; Danca, M.-F.; Wang, J. Caputo delta weakly fractional difference equations. Fract. Calc. Appl. Anal. 2022, 25, 2222–2240. [Google Scholar]
  34. Srivastava, H.M.; Mohammed, P.O.; Ryoo, C.S.; Hamed, Y.S. Existence and uniqueness of a class of uncertain Liouville-Caputo fractional difference equations. J. King Saud Univ. Sci. 2021, 33, 101497. [Google Scholar]
  35. Mohammed, P.O.; Goodrich, C.S.; Srivastava, H.M.; Al-Sarairah, E.; Hamed, Y.S. A study of monotonicity analysis for the delta and nabla discrete fractional operators of the Liouville-Caputo family. Axioms 2023, 12, 114. [Google Scholar] [CrossRef]
  36. Mozyrska, D.; Girejko, E. Overview of Fractional h-difference Operators. In Advances in Harmonic Analysis and Operator Theory. Operator Theory: Advances and Applications; Almeida, A., Castro, L., Speck, F.O., Eds.; Birkhäuser: Basel, Switzerland, 2013; Volume 229. [Google Scholar]
  37. Mozyrska, D.; Girejko, E.; Wyrwas, M. Fractional nonlinear systems with sequential operators. Cent. Eur. J. Phys. 2013, 11, 1295–1303. [Google Scholar]
  38. Abdeljawad, T. Different type kernel h-fractional differences and their fractional h–sums. Chaos Solit. Fractals 2018, 116, 146–156. [Google Scholar]
  39. Podlubny, I. Fractional Differential Equations; AP: New York, NY, USA, 1999. [Google Scholar]
  40. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equation; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  41. Girejko, E.; Mozyrska, D. Positivity of fractional discrete systems with sequential h-differences. In Proceedings of the 2012 IEEE 4th International Conference on Nonlinear Science and Complexity (NSC), Budapest, Hungary, 6–11 August 2012; pp. 137–140. [Google Scholar]
  42. Wang, M.; Jia, B.; Chen, C.; Zhu, X.; Du, F. Discrete fractional Bihari inequality and uniqueness theorem of solutions of nabla fractional difference equations with non-Lipschitz nonlinearities. Appl. Math. Comput. 2020, 367, 125118. [Google Scholar]
Figure 1. Numerical results for different values of p 0 .
Figure 1. Numerical results for different values of p 0 .
Axioms 13 00570 g001
Figure 2. Numerical results for different values of p 0 .
Figure 2. Numerical results for different values of p 0 .
Axioms 13 00570 g002
Table 1. Numerical results in case of p 0 = 0.1 .
Table 1. Numerical results in case of p 0 = 0.1 .
E t 1 2 Numerical Result | Error |
1.0000000000000001.0038044584697540.003804458469754
1.4867633976736791.4924416140809310.005678216407252
1.7990172441881771.7997757870838070.000758542895631
2.1076992038372702.1082387050239360.000539501186666
2.4300431414976622.4353511170277520.005307975530090
2.7742859576700092.7820776299710300.007791672301020
3.1462130322103353.1555531390526270.009340106842292
3.5508026836468893.5521017457316270.001299062084737
3.9928358341927083.9985240708014300.005688236608722
4.4771848107957354.4818787172063170.004693906410582
5.0089800807622835.0090991014572960.000119020695013
Table 2. Numerical results in case of p 0 = 0.01 .
Table 2. Numerical results in case of p 0 = 0.01 .
E t 1 2 Numerical Result | Error |
1.0000000000000001.0007112157804340.000711215780434
1.4867633976736791.4869851444076960.000221746734017
1.7990172441881771.7991346618390320.000117417650856
2.1076992038372702.1079958797104880.000296675873218
2.4300431414976622.4303619197995880.000318778301926
2.7742859576700092.7747101244297230.000424166759714
3.1462130322103353.1467208904949960.000507858284661
3.5508026836468893.5508881994439800.000085515797090
3.9928358341927083.9930983164274060.000262482234698
4.4771848107957354.4779858254185050.000801014622770
5.0089800807622835.0090093010398450.000029220277562
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

Mohammed, P.O.; Lupas, A.A.; Agarwal, R.P.; Yousif, M.A.; Al-Sarairah, E.; Abdelwahed, M. Theoretical Investigation of Fractional Estimations in Liouville–Caputo Operators of Mixed Order with Applications. Axioms 2024, 13, 570. https://doi.org/10.3390/axioms13080570

AMA Style

Mohammed PO, Lupas AA, Agarwal RP, Yousif MA, Al-Sarairah E, Abdelwahed M. Theoretical Investigation of Fractional Estimations in Liouville–Caputo Operators of Mixed Order with Applications. Axioms. 2024; 13(8):570. https://doi.org/10.3390/axioms13080570

Chicago/Turabian Style

Mohammed, Pshtiwan Othman, Alina Alb Lupas, Ravi P. Agarwal, Majeed A. Yousif, Eman Al-Sarairah, and Mohamed Abdelwahed. 2024. "Theoretical Investigation of Fractional Estimations in Liouville–Caputo Operators of Mixed Order with Applications" Axioms 13, no. 8: 570. https://doi.org/10.3390/axioms13080570

APA Style

Mohammed, P. O., Lupas, A. A., Agarwal, R. P., Yousif, M. A., Al-Sarairah, E., & Abdelwahed, M. (2024). Theoretical Investigation of Fractional Estimations in Liouville–Caputo Operators of Mixed Order with Applications. Axioms, 13(8), 570. https://doi.org/10.3390/axioms13080570

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