Next Article in Journal
An Integrated Decision-Making Approach for Green Supplier Selection in an Agri-Food Supply Chain: Threshold of Robustness Worthiness
Previous Article in Journal
On the Stability of Convection in a Non-Newtonian Vertical Fluid Layer in the Presence of Gold Nanoparticles: Drug Agent for Thermotherapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Riemann—Liouville and Caputo Fractional Forward Difference Monotonicity Analysis

by
Pshtiwan Othman Mohammed
1,*,
Thabet Abdeljawad
2,3,4,* and
Faraidun Kadir Hamasalh
1
1
Department of Mathematics, College of Education, University of Sulaimani, Sulaymaniyah 46001, Kurdistan Region, Iraq
2
Department of Mathematics and General Sciences, Prince Sultan University, P.O. Box 66833, Riyadh 11586, Saudi Arabia
3
Department of Medical Research, China Medical University, Taichung 40402, Taiwan
4
Department of Computer Science and Information Engineering, Asia University, Taichung 41354, Taiwan
*
Authors to whom correspondence should be addressed.
Submission received: 13 May 2021 / Revised: 3 June 2021 / Accepted: 3 June 2021 / Published: 7 June 2021

Abstract

:
Monotonicity analysis of delta fractional sums and differences of order υ ( 0 , 1 ] on the time scale h Z are presented in this study. For this analysis, two models of discrete fractional calculus, Riemann–Liouville and Caputo, are considered. There is a relationship between the delta Riemann–Liouville fractional h-difference and delta Caputo fractional h-differences, which we find in this study. Therefore, after we solve one, we can apply the same method to the other one due to their correlation. We show that y ( z ) is υ -increasing on M a + υ h , h , where the delta Riemann–Liouville fractional h-difference of order υ of a function y ( z ) starting at a + υ h is greater or equal to zero, and then, we can show that y ( z ) is υ -increasing on M a + υ h , h , where the delta Caputo fractional h-difference of order υ of a function y ( z ) starting at a + υ h is greater or equal to 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) for each z M a + h , h . Conversely, if y ( a + υ h ) is greater or equal to zero and y ( z ) is increasing on M a + υ h , h , we show that the delta Riemann–Liouville fractional h-difference of order υ of a function y ( z ) starting at a + υ h is greater or equal to zero, and consequently, we can show that the delta Caputo fractional h-difference of order υ of a function y ( z ) starting at a + υ h is greater or equal to 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) on M a , h . Furthermore, we consider some related results for strictly increasing, decreasing, and strictly decreasing cases. Finally, the fractional forward difference initial value problems and their solutions are investigated to test the mean value theorem on the time scale h Z utilizing the monotonicity results.

1. Introduction

Fractional differentiation and integration have opened many new doors for researchers in recent decades due to their wide and novel applicability in many fields of science including mathematical analysis, technology, and engineering (see [1,2,3,4,5,6,7]). Many techniques are used to deal with these new differential and integral operators; for instance, some researchers used analytical techniques including Laplace transform, spline interpolation, Green function, Crank–Nicolson approximation method, method of separation of variable, and many others to derive exact solutions to linear differential or integral equations (see [8,9,10,11,12,13,14]). Using the fixed-point technique, some researchers provided the conditions under which differential and integral equations have unique solutions. Some others provided numerical schemes that could be used to solve numerically differential and integral equations with fractional order. Very recently, fractional differentiation and integration found application in image processing, where the fractional kernel is used to remove noise in a given image.
On the other hand, fractional operators were employed in fuzzy theory. In fact, so far, researchers have developed a new class of differential and integral equations called fuzzy fractional differential and integral equations. This topic is highly regarded as its applications are found in many fields too. We point out that fractional differential and integral operators can be represented differently in continuous form, discrete form and discretized form. Discrete fractional calculus has been the focus of many researchers in recent years. For recent research on this topic, we advise the readers to refer to [15,16,17,18,19,20,21,22,23,24,25,26,27].
The aim of this study is to investigate the υ -monotonicity analysis on h-discrete delta fractional models in the sense of Riemann–Liouville ( RL ) and Caputo fractional operators on the time scale h Z . The remainder of our article is structured as follows. In Section 2, we provide some notations and make some preparations. In Section 3, the monotonicity results and some corollaries are presented. In Section 4, some results related to the RL fractional forward difference equation ( RL FFDE ) (Section 4.1) and the Caputo fractional forward difference equation (Caputo FFDE ) (Section 4.2) are prepared. Additionally, we discuss the mean value theorem ( MVT ) later as an application of our monotonicity results. The paper is concluded in Section 5.

2. Preliminaries

Related concepts regarding the discrete fractional operators used in the current article are shown in this section.
Definition 1 (see [20,21,22]) Let f be defined on the time scale h Z ; then, the forward h-difference operator is given by
Δ h f ( z ) = f ( σ h ( z ) ) f ( z ) h z M a , h : = { a , a + h , a + 2 h , } ,
and the backward h-difference operator is given by
h f ( z ) = f ( z ) f ( ρ h ( z ) ) h ( z M a , h ) ,
where σ h ( z ) : = z + h and ρ h ( z ) : = z h .
Definition 2 (see [20,21,22]). Let υ > 0 and f : M a , h R with a starting point a. Then, the delta left RL fractional h-sum of order υ is given by
Δ a h υ f ( z ) = 1 Γ ( υ ) r = a h z h υ z σ h ( r h ) h ( υ 1 ) f ( r h ) h ( z M a + υ h , h ) ,
and for a function f : b , h M : = { , b 2 h , b h , b } R with an end point b, the delta right RL fractional h-sum of order υ is given by
h Δ b υ f ( z ) = 1 Γ ( υ ) r = z h + υ b h r h σ h ( z ) h ( υ 1 ) f ( r h ) h ( z b υ h , h M ) ,
where the h-falling factorial function z h ( υ ) is defined by
z h ( υ ) = h υ Γ z h + 1 Γ z h + 1 υ ( z , υ R ) ,
and we use the convention that division at a pole yields 0.
Definition 3 (see [22]). For 0 υ < 1 , the delta left RL fractional h-difference of order υ is defined by
Δ a R L h υ f ( z ) = Δ h Δ a h ( 1 υ ) f ( z ) = 1 Γ ( 1 υ ) Δ h r = a h z h + υ 1 z σ h ( r h ) h ( υ ) f ( r h ) h ( z M a + ( 1 υ ) h , h ) ,
and the delta right RL fractional h-difference of order υ is defined by
Δ h R L b υ f ( z ) = Δ h h Δ b ( 1 υ ) f ( z ) = 1 Γ ( 1 υ ) Δ h r = z h + 1 υ b h ρ h ( z ) r h h ( υ ) f ( r h ) h ( z b ( 1 υ ) h , h M ) .
Definition 4 (see [22]). For 0 υ < 1 , the delta left Caputo fractional h-difference of order υ is defined by
Δ a C h υ f ( z ) = Δ a h ( 1 υ ) Δ h f ( z ) = 1 Γ ( 1 υ ) r = a h z h + υ 1 z σ h ( r h ) h ( υ ) ( Δ h f ) ( r h ) h ( z M a + ( 1 υ ) h , h ) ,
and the delta right Caputo fractional h-difference of order υ is defined by
Δ h C b υ f ( z ) = h Δ b ( 1 υ ) ( Δ h f ) ( z ) = 1 Γ ( 1 υ ) r = z h + 1 υ b h ρ h ( z ) r h h ( υ ) ( Δ h f ) ( r h ) h ( z b ( 1 υ ) h , h M ) .
In the following lemma, we show that z h ( υ ) increases on M 0 , h .
Lemma 1.
Let υ > 0 and h > 0 , then Δ h z h ( υ ) = υ z h ( υ 1 ) . Moreover, z h ( υ ) increases on M 0 , h .
Proof. 
From Definition 1 and Equation (3), we have
Δ h z h ( υ ) = ( z + h ) h ( υ ) z h ( υ ) h = 1 h h υ Γ z h + 2 Γ z h + 2 υ h υ Γ z h + 1 Γ z h + 1 υ = h υ 1 Γ z h + 1 Γ z h + 1 υ z h + 1 z h + 1 υ 1 = υ h υ 1 Γ z h + 1 Γ z h + 1 ( υ 1 ) = υ z h ( υ 1 ) .
Since υ , h > 0 , it follows that
Δ h z h ( υ ) = ( z + h ) h ( υ ) z h ( υ ) h = υ z h ( υ 1 ) 0 ,
which implies that ( z + h ) h ( υ ) z h ( υ ) , and this completes the proof. □
The following theorem can be seen as an equivalence definition to the delta RL fractional h-differences.
Theorem 1.
Let 0 < υ < 1 . Then, the delta left and delta right RL fractional h-differences of order υ defined on M a + ( 1 υ ) h , h and b ( 1 υ ) h , h M , respectively, are defined by
Δ a R L h υ f ( z ) = 1 Γ ( υ ) r = a h z h + υ ( z σ h ( r h ) ) h ( υ 1 ) f ( r h ) h ,
and
Δ h R L b υ f ( z ) = 1 Γ ( υ ) r = z h υ b h ( ρ h ( z ) r h ) h ( υ 1 ) f ( r h ) h .
Proof. 
From Definitions 1 and 3, we have
Δ a R L h υ f ( z ) = 1 Γ ( 1 υ ) Δ h r = a h z h + υ 1 z σ h ( r h ) h ( υ ) f ( r h ) h = 1 h Γ ( 1 υ ) r = a h z h + υ ( z + h σ h ( r h ) ) h ( υ ) f ( r h ) h r = a h z h + υ 1 z σ h ( r h ) h ( υ ) f ( r h ) h = 1 h Γ ( 1 υ ) r = a h z h + υ 1 ( z + h σ h ( r h ) ) h ( υ ) f ( r h ) h r = a h z h + υ 1 z σ h ( r h ) h ( υ ) f ( r h ) h + 1 h Γ ( 1 υ ) ( z + h σ h ( z + υ h ) ) h ( υ ) f ( z + υ h ) h = 1 Γ ( 1 υ ) r = a h z h + υ 1 Δ h z σ h ( r h ) h ( υ ) f ( r h ) h + ( υ h ) h ( υ ) h Γ ( 1 υ ) f ( z + υ h ) h = 1 ( υ ) Γ ( υ ) r = a h z h + υ 1 Δ h z σ h ( r h ) h ( υ ) f ( r h ) h + h 1 υ f ( z + υ h ) h .
Then by using Lemma 1, we get
Δ a R L h υ f ( z ) = 1 Γ ( υ ) r = a h z h + υ 1 z σ h ( r h ) h ( υ 1 ) f ( r h ) h + h υ 1 f ( z + υ h ) h = 1 Γ ( υ ) r = a h z h + υ z σ h ( r h ) h ( υ 1 ) f ( r h ) h ,
which is the required Equation (4). In the same manner as Equation (4), we can prove Equation (5), and thus, the proof is completed. □
A relationship between the delta RL fractional and delta Caputo fractional h-differences are presented in the following proposition.
Proposition 1.
Let 0 υ < 1 , then
Δ a C h υ f ( z ) = Δ a R L h υ f ( z ) 1 Γ ( 1 υ ) ( z a ) h ( υ ) f ( a )
for z M a + ( 1 υ ) h , h , and
Δ h C b υ f ( z ) = Δ h R L b υ f ( z ) 1 Γ ( 1 υ ) ( b z ) h ( υ ) f ( b )
for z b ( 1 υ ) h , h M .
Proof. 
From Definition 4 and the fact that z σ h ( r h ) h ( υ ) = 0 at r = z h + υ , we have
Δ a C h υ f ( z ) = Δ a h ( 1 υ ) Δ h f ( z ) = Δ a h ( 1 υ ) f ( z + h ) f ( z ) h = 1 h Γ ( 1 υ ) r = a + h h z h + υ ( z + h σ h ( r h ) ) h ( υ ) f ( r h ) h 1 h Γ ( 1 υ ) r = a h z h + υ 1 z σ h ( r h ) h ( υ ) f ( r h ) h = 1 h Γ ( 1 υ ) r = a h + 1 z h + υ ( z + h σ h ( r h ) ) h ( υ ) f ( r h ) h 1 h Γ ( 1 υ ) r = a h z h + υ z σ h ( r h ) h ( υ ) f ( r h ) h
= 1 h Γ ( 1 υ ) r = a h z h + υ ( z + h σ h ( r h ) ) h ( υ ) f ( r h ) h 1 Γ ( 1 υ ) ( z a ) h ( υ ) f ( a ) 1 h Γ ( 1 υ ) r = a h z h + υ z σ h ( r h ) h ( υ ) f ( r h ) h = 1 Γ ( 1 υ ) r = a h z h + υ Δ h z σ h ( r h ) h ( υ ) f ( r h ) h 1 Γ ( 1 υ ) ( z a ) h ( υ ) f ( a ) .
By using Lemma 1 and Theorem 1, we get
Δ a C h υ f ( z ) = 1 Γ ( υ ) r = a h z h + υ z σ h ( r h ) h ( υ 1 ) f ( r h ) h 1 Γ ( 1 υ ) ( z a ) h ( υ ) f ( a ) = Δ a R L h υ f ( z ) 1 Γ ( 1 υ ) ( z a ) h ( υ ) f ( a ) ,
which is the required Equation (6). Using the same technique used for Equation (6), we can prove Equation (7), and thus, the proof is completed. □
In the following lemma, we prove and modify a power rule that appeared in ([22], Lemma 4). We state and prove the modified result in a simpler way as follows:
Lemma 2.
Let υ > 0 , μ > 1 and h > 0 , then
Δ a + μ h h υ ( z a ) h ( μ ) = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) ( z a ) h ( υ + μ )
for z M a + ( μ + υ ) h , h .
Proof. 
Following ([4], Lemma 1), we have
Δ a h + μ υ z h a h ( μ ) = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) z h a h ( υ + μ )
for z M a + ( μ + υ ) h , h . Calculating both sides of Equation (9), we get
LHS : = Δ a h + μ υ z h a h ( μ ) = 1 Γ ( υ ) r = a h + μ z h υ z h σ h ( r h ) ( υ 1 ) r a h ( μ ) = 1 Γ ( υ ) r = a h + μ z h υ z r h h h ( υ 1 ) r h a h ( μ ) = 1 Γ ( υ ) r = a h + μ z h υ Γ z r h h h + 1 Γ z r h h h + 1 ( υ 1 ) · Γ r h a h + 1 Γ r h a h + 1 μ = 1 Γ ( υ ) r = a h + μ z h υ ( z r h h ) h ( υ 1 ) h υ 1 · ( r h a ) h ( μ ) h μ = 1 Γ ( υ ) r = a h + μ z h υ 1 h υ + μ z σ h ( r h ) h ( υ 1 ) ( r h a ) h ( μ ) h ,
and
RHS : = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) z h a h ( υ + μ ) = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) Γ z a h + 1 Γ z a h + 1 ( υ + μ ) = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) ( z a ) h ( υ + μ ) h υ + μ .
Substituting the LHS and RHS results into Equation (9), we get
1 Γ ( υ ) r = a h + μ z h υ 1 h υ + μ z σ h ( r h ) h ( υ 1 ) ( r h a ) h ( μ ) h = Γ ( μ + 1 ) Γ ( μ + 1 + υ ) ( z a ) h ( υ + μ ) h υ + μ .
Multiplying by a positive constant h υ + μ on both sides of the equality, we get the desired Equation (8). □

3. The Monotonicity Results

This section illustrates the monotonicity of a discrete function. Monotonicity analysis of discrete functions defined on M 1 0 was originally introduced in [27], and there is extensive literature on monotonicity analysis techniques and its extensions on M a , h ; for example, see [22,23,26].
Definition 5 (see [22,23,26]). Let 0 < h , υ 1 , and y : M a , h R be a function satisfying y ( a ) 0 . Then, y is called an υ-increasing function on M a , h if
y ( z + h ) υ y ( z ) z M a , h .
Observe that, if y ( z ) is increasing on M a , h , then y ( z + h ) y ( z ) for all z M a , h , and thus, y ( z ) is υ-increasing on M a , h .
Definition 6 (see [22,23,26]). Let 0 < h , υ 1 , and y : M a , h R be a function satisfying y ( a ) 0 . Then, y is called an υ-decreasing function on M a , h if
y ( z + h ) υ y ( z ) z M a , h .
Observe that, if y ( z ) is decreasing on M a , h , then y ( z + h ) y ( z ) for all z M a , h , and thus, y ( z ) is υ-decreasing on M a , h .
Remark 1.
Note that, if υ = 1 in Definition 5, then the increasing and υ-increasing concepts coincide and that, if υ = 1 in Definition 6, then the decreasing and υ-decreasing concepts coincide.
To provide motivation for the above monotonicity definitions, we prove a few fundamental results of the discrete RL and Caputo fractional operators.
Theorem 2.
Let y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . Suppose that Δ a + υ h R L h υ y ( z ) 0 for 0 < h 1 , 0 < υ < 1 and z M a + h , h . Then, y ( z ) is υ-increasing on M a + υ h , h .
Proof. 
From the assumption and proof of Theorem 1, we have
Δ a + υ h R L h υ y ( z ) = h υ y ( z + υ h ) + 1 Γ ( υ ) r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) y ( r h ) h 0 .
For z = a , we see that
Δ a + υ h R L h υ y ( a ) = h υ y ( a + υ h ) 0 ,
and, thus, y ( a + υ h ) 0 .
For z = a + h , we see that
Δ a + υ h R L h υ y ( a + h ) = h υ y ( a + υ h + h ) + 1 Γ ( υ ) ( υ h ) h ( υ 1 ) y ( a + υ h ) h = h υ y ( a + υ h + h ) υ h υ y ( a + υ h ) 0 ,
and, thus, y ( a + υ h + h ) υ y ( a + υ h ) since h υ > 0 and y ( a + υ h ) 0 .
Now, inductively, we show that
y ( z + h ) υ y ( z ) , z M a + υ h , h .
Suppose that y ( k + h ) υ y ( k ) , k < z provided that k , z M a + υ h , h . Then, we have to show that y ( z + h ) υ y ( z ) .
Replace z by z + h in Equation (10) to get
h υ y ( z + υ h + h ) + 1 Γ ( υ ) r = a h + υ z h + υ ( z + h σ h ( r h ) ) h ( υ 1 ) y ( r h ) h = h υ y ( z + υ h + h ) υ Γ ( 1 υ ) r = a h + υ z h + υ ( z + h σ h ( r h ) ) h ( υ 1 ) y ( r h ) h 0 .
It follows that
h υ y ( z + υ h + h ) υ Γ ( 1 υ ) [ ( z + h σ h ( a + υ h ) ) h ( υ 1 ) y ( a + υ h ) h + ( z + h σ h ( a + υ h + h ) ) h ( υ 1 ) y ( a + υ h + h ) h + + ( z + h σ h ( z + υ h ) ) h ( υ 1 ) y ( z + υ h ) h ] 0 ,
or equivalently,
h υ y ( z + υ h + h ) υ Γ ( 1 υ ) [ ( z + h σ h ( a + υ h ) ) h ( υ 1 ) y ( a + υ h ) h + ( z + h σ h ( a + υ h + h ) ) h ( υ 1 ) y ( a + υ h + h ) h + + ( z + h σ h ( z + υ h ) ) h ( υ 1 ) y ( z + υ h ) h ] υ Γ ( 1 υ ) ( z + h σ h ( z + υ h ) ) h ( υ 1 ) y ( z + υ h ) h = υ Γ ( 1 υ ) h υ 1 Γ ( 1 υ ) y ( z + υ h ) h = υ h υ Γ ( 1 υ ) y ( z + υ h ) .
This implies that y ( z + υ h + h ) y ( z + υ h ) , and this completes the proof. □
Corollary 1.
Let y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . Suppose that
a + υ h C Δ h υ y ( z ) 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) ,
for 0 < h 1 , 0 < υ < 1 , and z M a + h , h . Then, y ( z ) is υ-increasing on M a + υ h , h .
Proof. 
The proof follows from Proposition 1 and Theorem 2. □
Theorem 3.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . If y ( z ) is increasing on M a + υ h , h , then
Δ a + υ h R L h υ y ( z ) 0 z M a , h .
Proof. 
For each z M a + h , h , we have to show that
Δ a + υ h R L h υ y ( z ) = 1 Γ ( υ ) r = a h + υ z h + υ z σ h ( r h ) h ( υ 1 ) y ( r h ) h 0 .
From Equation (10) with z = a , we have
Δ a + υ h R L h υ y ( a ) = h υ y ( a + υ h ) 0 ,
since h υ > 0 and y ( a + υ h ) 0 . Assume that Δ a + υ h R L h υ y ( i ) 0 , i < z ; then, we have to show that Δ a + υ h R L h υ y ( z ) 0 .
From Equation (10), we see that
Δ a + υ h R L h υ y ( z ) = h υ y ( z + υ h ) + 1 Γ ( υ ) r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) y ( r h ) h = h υ y ( z + υ h ) υ Γ ( 1 υ ) r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) y ( r h ) h = h υ y ( z + υ h ) υ h υ y ( z + υ h h ) υ Γ ( 1 υ ) r = a h + υ z h + υ 2 z σ h ( r h ) h ( υ 1 ) y ( r h ) h .
Write r = a h + υ + , = 0 , 1 , . . . , z h a h + υ 2 . We see that
υ z σ h ( r h ) h ( υ 1 ) Γ ( 1 υ ) · h = υ h υ Γ ( 1 υ ) Γ z a h υ Γ z a h + 1 = υ h υ Γ ( 1 υ ) Γ m υ Γ m + 1 , denoting m = z a h , = h υ Γ ( υ ) Γ m υ Γ m + 1 .
Since z a h + υ 2 , we get m > 2 υ > υ for υ ( 0 , 1 ) and, hence, m > m υ > 0 . Additionally, h υ > 0 , and from the graph of the gamma function, we see that Γ ( υ ) < 0 for υ ( 0 , 1 ) . Therefore,
υ z σ h ( r h ) h ( υ 1 ) Γ ( 1 υ ) · h > 0 .
From the assumption that y ( z ) increases, then y ( z + υ h ) y ( z + υ h h ) y ( a + υ h ) for each z M a , h . It follows that
y ( z + υ h h ) y ( r h ) , r = a h + υ , a h + υ , . . . , z h + υ 2 .
Considering Equations (11)–(13), it follows that
Δ a + υ h R L h υ y ( z ) = h υ y ( z + υ h ) υ h υ y ( z + υ h h ) υ Γ ( 1 υ ) r = a h + υ z h + υ 2 z σ h ( r h ) h ( υ 1 ) y ( z + υ h h ) h = h υ y ( z + υ h ) υ Γ ( 1 υ ) r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) y ( z + υ h h ) h = h υ y ( z + υ h ) ( υ h ) h ( υ ) Γ ( 1 υ ) y ( z + υ h h ) + ( υ h ) h ( υ ) Γ ( 1 υ ) y ( z + υ h h ) υ h Γ ( 1 υ ) y ( z + υ h h ) r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) = h υ y ( z + υ h ) y ( z + υ h h ) + 1 Γ ( 1 υ ) y ( z + υ h h ) ( υ h ) h ( υ ) υ h r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) .
Since y ( z ) increases, it follows that
Δ a + υ h R L h υ y ( z ) 1 Γ ( 1 υ ) y ( z + υ h h ) ( υ h ) h ( υ ) υ h r = a h + υ z h + υ 1 z σ h ( r h ) h ( υ 1 ) = 1 Γ ( 1 υ ) y ( z + υ h h ) [ ( υ h ) h ( υ ) υ h ( z a υ h h ) h ( υ 1 ) υ h ( z a υ h 2 h ) h ( υ 1 ) υ h ( υ h ) h ( υ 1 ) ] = h υ y ( z + υ h h ) Γ ( 1 υ ) [ Γ ( 1 υ ) Γ ( 1 ) υ Γ z a h υ Γ z a h + 1 υ Γ z a h υ 1 Γ z a h υ Γ ( 2 υ ) Γ ( 3 ) υ Γ ( 1 υ ) Γ ( 2 ) ] = h υ y ( z + υ h h ) Γ ( 1 υ ) [ Γ ( 1 υ ) Γ ( 1 ) υ Γ ( 1 υ ) Γ ( 2 ) υ Γ ( 2 υ ) Γ ( 3 ) υ Γ ( 3 υ ) Γ ( 4 ) υ Γ z a h υ Γ z a h + 1 ]
= h υ y ( z + υ h h ) Γ ( 1 υ ) [ Γ ( 1 υ ) Γ ( 2 ) ( 1 υ ) υ Γ ( 2 υ ) Γ ( 3 ) υ Γ ( 3 υ ) Γ ( 4 ) υ Γ z a h υ Γ z a h + 1 ]    = h υ y ( z + υ h h ) Γ ( 1 υ ) Γ ( 2 υ ) Γ ( 3 ) ( 2 υ ) υ Γ ( 3 υ ) Γ ( 4 ) υ Γ z a h υ Γ z a h + 1              = h υ y ( z + υ h h ) Γ ( 1 υ ) Γ ( 3 υ ) Γ ( 4 ) ( 3 υ ) υ Γ z a h υ Γ z a h + 1 .
Continue this process to get
Δ a + υ h R L h υ y ( z ) h υ y ( z + υ h h ) Γ ( 1 υ ) Γ z a h υ Γ z a h υ Γ z a h υ Γ z a h + 1 = h υ y ( z + υ h h ) Γ ( 1 υ ) Γ z a h υ Γ z a h 1 υ 1 z a h = h υ y ( z + υ h h ) Γ ( 1 υ ) z a h υ Γ z a h υ z a h Γ z a h = h υ y ( z + υ h h ) Γ ( 1 υ ) Γ z a h υ + 1 Γ z a h + 1 = h υ y ( z + υ h h ) Γ ( 1 υ ) Γ z a υ h h + 1 Γ z a υ h h + 1 ( υ ) = ( z ( a + υ h ) ) h ( υ ) Γ ( 1 υ ) y ( z + υ h h ) 0 ,
which completes the proof. □
Corollary 2.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . If y ( z ) increases on M a + υ h , h , then
a + υ h C Δ h υ y ( z ) 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) z M a , h .
Proof. 
The proof follows from Proposition 1 and Theorem 3. □
Theorem 4.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) > 0 . If y ( z ) strictly increases on M a + υ h , h , then
Δ a + υ h R L h υ y ( z ) > 0 z M a , h .
Proof. 
The proof is similar to Theorem 3, so it is omitted. □
Corollary 3.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) > 0 . If y ( z ) strictly increases on M a + υ h , h , then
a + υ h C Δ h υ y ( z ) > 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) z M a , h .
Theorem 5.
Let y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . Suppose that Δ a + υ h R L h υ y ( z ) 0 for 0 < h 1 , 0 < υ < 1 and z M a + h , h . Then, y ( z ) is υ-decreasing on M a + υ h , h .
Proof. 
Let g : M a + υ h , h R be a function such that g ( z ) = y ( z ) . Therefore,
Δ a + υ h R L h υ g ( z ) = Δ a + υ h R L h υ ( y ) ( z ) = Δ a + υ h R L h υ y ( z ) 0 .
Apply Theorem 2 to g ( z ) , and thus, the proof is completed. □
Corollary 4.
Let y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . Suppose that a + υ h C Δ h υ y ( z ) 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) for 0 < h 1 , 0 < υ < 1 and z M a + h , h . Then, y ( z ) is υ-decreasing on M a + υ h , h .
Theorem 6.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . If y ( z ) decreases on M a + υ h , h , then
Δ a + υ h R L h υ y ( z ) 0 z M a , h .
Proof. 
The proof is obtained by applying Theorem 4 to g ( z ) = y ( z ) . □
Corollary 5.
Let 0 < h 1 , 0 < υ < 1 , and y : M a + υ h , h R be a function satisfying y ( a + υ h ) 0 . If y ( z ) decreases on M a + υ h , h , then
a + υ h C Δ h υ y ( z ) 1 Γ ( 1 υ ) ( z ( a + υ h ) ) h ( υ ) y ( a + υ h ) z M a , h .

4. Fractional Forward Difference Initial Value Problem and Mean Value Theorem

In this section, we move on from monotonicity analysis to the MVT . Having established the monotonicity analysis for the discrete RL and Caputo fractional operators, we now obtain the MVT using those discrete monotonicity results.

4.1. Establishing the Riemann–Liouville case

The following is the main results to start off the MVT here.
Lemma 3
([2,3]). Let y : M a , h R and 0 < h 1 , 0 < υ < 1 , then
Δ a h υ Δ h y ( z ) = Δ h Δ a h υ y ( z ) ( z a ) h ( υ 1 ) Γ ( υ ) y ( a ) .
Lemma 4
([2,3]). Let y : M a , h R and 0 < h 1 , 0 < υ < 1 , then
a + β h Δ h υ Δ a h β y ( z ) = Δ a h ( υ + β ) y ( z ) = Δ a + υ h h β Δ a h υ y ( z ) , z M a + ( υ + β ) h , h .
Theorem 7.
Let y : M a , h R and 0 < h 1 , 0 < υ < 1 , then
Δ a + h h υ Δ a + υ h R L h υ y ( z ) = y ( z ) h 1 υ Γ ( υ ) ( z ( a + h ) ) h ( υ 1 ) y ( a + υ h ) .
Proof. 
Let g ( z ) : = Δ a + h h υ Δ a + υ h R L h υ y ( z ) , then we have
g ( z ) = Δ a + h h υ Δ a + υ h R L h υ y ( z ) = Δ a + h h υ Δ h Δ a + υ h h ( 1 υ ) y ( z ) = Δ a + h h υ Δ h f ( z ) , where f ( z ) = Δ a + υ h R L h ( 1 υ ) y ( z ) .
Apply Lemma 3 but replace a by a + h to get
g ( z ) = Δ h Δ a + h h υ f ( z ) ( z ( a + h ) ) h ( υ 1 ) Γ ( υ ) f ( a + h ) .
Calculating f ( a + h ) using Definition 2, we get
f ( a + h ) = Δ a + υ h h ( 1 υ ) y ( a + h ) = 1 Γ ( 1 υ ) r = a h + υ a + h h ( 1 υ ) ( a + h σ h ( r h ) ) h ( υ 1 ) y ( r h ) h = 1 Γ ( 1 υ ) a + h σ h ( a + υ h ) h ( υ ) y ( a + υ h ) h = 1 Γ ( 1 υ ) ( υ h ) h ( υ ) y ( a + υ h ) h = h Γ ( 1 υ ) h υ Γ ( 1 υ ) Γ ( 1 ) y ( a + υ h ) = h 1 υ y ( a + υ h ) .
Substituting this and the value of f ( z ) in Equation (14), we get
g ( z ) = Δ h Δ a + h h υ Δ a + υ h h ( 1 υ ) y ( z ) ( z ( a + h ) ) h ( υ 1 ) Γ ( υ ) h 1 υ y ( a + υ h ) .
Applying Lemma 4 using β = 1 υ , we get
g ( z ) = Δ h Δ a + h h ( υ + 1 υ ) y ( z ) h 1 υ Γ ( υ ) ( z ( a + h ) ) h ( υ 1 ) y ( a + υ h ) = y ( z ) h 1 υ Γ ( υ ) ( z ( a + h ) ) h ( υ 1 ) y ( a + υ h ) .
This completes the proof. □
Consider RL FFDE :
Δ a + υ h R L h υ y ( z ) = f ( z , y ( z ) ) z M a + h , h
with the initial condition
Δ a + υ h R L h ( 1 υ ) y ( z ) | z = a + h = h 1 υ y ( a + υ h ) = c 1 ,
where a R , h , υ ( 0 , 1 ) , and c 1 is a constant.
Theorem 8.
y is a solution of RL FFDE in Equation (15) with the initial condition in Equation (16) if and only if it has the representation
y ( z ) = ( z ( a + h ) ) h ( υ 1 ) Γ ( υ ) c 1 + Δ a + h h υ f ( z , y ( z ) ) .
Proof. 
This proof follows from Theorem 7. □
According to Theorem 7, we can write
Δ a + h h υ Δ a + υ h R L h υ y ( z ) = y ( z ) N h ( υ , z , a ) y ( a + υ h ) ,
where N h ( υ , z , a ) = h 1 υ Γ ( υ ) ( z ( a + h ) ) h ( υ 1 ) .
It is worthwhile to mention that analyzing the monotonicity property of such fractional difference operators is useful for better understanding the qualitative properties of solutions of different discrete fractional dynamic equations. The monotonicity properties are part of the basics of the discrete fractional calculus, where for example, we used them in this article to prove a discrete fractional version of MVT . Then, it is of interest to obtain the following MVT for RL FFDE in Equation (17).
Theorem 9.
Let f and g be functions defined on I : = M a + h , h M b h , h = { a + h , a + 2 h , a + 3 h , , b 2 h , b h , b } , where b = a + k h for some k M . Assume that g strictly increases, g ( a + υ h ) > 0 , and 0 < υ < 1 , 0 < h 1 . Then, there exist s 1 , s 2 I such that
Δ a + υ h R L h υ f ( s 1 ) Δ a + υ h R L h υ g ( s 1 ) f ( b ) N h ( υ , b , a ) f ( a + υ h ) g ( b ) N h ( υ , b , a ) g ( a + υ h ) a + υ h R L Δ h υ f ( s 2 ) a + υ h R L Δ h υ g ( s 2 ) .
Proof. 
First, we prove that the denominators in the inequality in Equation (18) are not zero. Since g strictly increases with g ( a + υ h ) > 0 , then by Theorem 4, we have
Δ a + υ h R L h υ g ( z ) > 0 z I .
Apply the fractional sum operator Δ a + h h υ on both sides in Equation (19) to get
Δ a + h h υ Δ a + υ h R L h υ g ( z ) > Δ a + h h υ ( 0 ) .
This together with Equation (17) imply that
g ( z ) N h ( υ , z , a ) g ( a + υ h ) > 0 z I ,
and particularly g ( b ) N h ( υ , b , a ) g ( a + υ h ) > 0 .
To end the proof, we use the contradiction technique: we assume that Equation (18) is not true. Then, either
f ( b ) N h ( υ , b , a ) f ( a + υ h ) g ( b ) N h ( υ , b , a ) g ( a + υ h ) < Δ a + υ h R L h υ f ( z ) Δ a + υ h R L h υ g ( z ) z I ,
or
f ( b ) N h ( υ , b , a ) f ( a + υ h ) g ( b ) N h ( υ , b , a ) g ( a + υ h ) > Δ a + υ h R L h υ f ( z ) Δ a + υ h R L h υ g ( z ) z I .
By means of Equation (19), we can multiply both sides of Equation (20) by a positive constant, and thus, we get
f ( b ) N h ( υ , b , a ) f ( a + υ h ) g ( b ) N h ( υ , b , a ) g ( a + υ h ) Δ a + υ h R L h υ g ( z ) < Δ a + υ h R L h υ f ( z ) z I .
Applying the fractional sum operator Δ a + h h υ on both sides of this inequality, we get
f ( b ) N h ( υ , b , a ) f ( a + υ h ) g ( b ) N h ( υ , b , a ) g ( a + υ h ) Δ a + h h υ Δ a + υ h R L h υ g ( z ) < Δ a + h h υ Δ a + υ h R L h υ f ( z ) .
This at z = b and Equation (17) imply
f ( b ) N h ( υ , b , a ) f ( a + υ h ) < f ( b ) N h ( υ , b , a ) f ( a + υ h ) ,
and this implies that f ( b ) < f ( b ) , which is a contradiction. In the similar manner, the inequality in Equation (21) leads to a contradiction. The proof follows. □

4.2. Establishing the Caputo Case

Let us discuss the following results.
Lemma 5.
For any a R and h , υ > 0 , we have
Δ a + h h υ z ( a + υ h ) h ( υ ) = Γ ( 1 υ ) 1 h 1 υ Γ ( υ ) z ( a + h ) h ( υ 1 ) ( z M a + υ h , h ) .
Proof. 
Let f be a function defined on M a , h ; then, from the Definition 2, we have
Δ a + h h υ f ( z ) = 1 Γ ( υ ) r = a h + 1 z h υ z σ h ( r h ) h ( υ 1 ) f ( r h ) h = Δ a h υ f ( z ) h Γ ( υ ) z ( a + h ) h ( υ 1 ) f ( a ) .
Considering f ( z ) = z ( a + υ h ) h ( υ ) in Equation (22), we get
Δ a + h h υ z ( a + υ h ) h ( υ ) = Δ a h υ z ( a + υ h ) h ( υ ) h Γ ( υ ) z ( a + h ) h ( υ 1 ) υ h ) h ( υ )
Then, by applying Lemma 2 and the identity (3), it follows that
Δ a + h h υ z ( a + υ h ) h ( υ ) = Γ ( 1 υ ) Γ ( 1 υ + υ ) z ( a + υ h ) h ( 0 ) h Γ ( υ ) z ( a + h ) h ( υ 1 ) h υ Γ 1 υ h h Γ 1 υ h h + υ = Γ ( 1 υ ) h 1 υ Γ ( υ ) z ( a + h ) h ( υ 1 ) Γ ( 1 υ ) .
Hence, the proof is complete. □
Theorem 10.
Let y : M a , h R and 0 < h 1 , 0 < υ < 1 ; then,
Δ a + h h υ Δ a + υ h C h υ y ( z ) = y ( z ) y ( a + υ h ) .
Proof. 
We use Proposition 1 in order to proceed
Δ a + υ h C h υ y ( z ) = Δ a + υ h R L h υ y ( z ) 1 Γ ( 1 υ ) z ( a + υ h ) h ( υ ) y ( a + υ h ) .
Taking Δ a + h h υ on both sides of Equation (24) and then using Theorem 7 and Lemma 5, we get
Δ a + h h υ Δ a + υ h C h υ y ( z ) = Δ a + h h υ Δ a + υ h R L h υ y ( z ) y ( a + υ h ) Γ ( 1 υ ) Δ a + h h υ z ( a + υ h ) h ( υ ) = y ( z ) h 1 υ Γ ( υ ) ( z ( a + h ) ) h ( υ 1 ) y ( a + υ h ) y ( a + υ h ) Γ ( 1 υ ) · Γ ( 1 υ ) 1 h 1 υ Γ ( υ ) z ( a + h ) h ( υ 1 ) = y ( z ) y ( a + υ h ) ,
and after some simplifications, the result is obtained. □
Consider the following Caputo FFDE :
Δ a + υ h C h υ y ( z ) = f ( z , y ( z ) ) z M a + h , h
with the initial condition
Δ a + υ h R L h ( 1 υ ) y ( z ) | z = a + h = h 1 υ y ( a + υ h ) = h 1 υ c 2 ,
where a R , h , υ ( 0 , 1 ) , and c 2 is a constant.
Theorem 11.
y is a solution of Caputo FFDE in Equation (25) with the initial condition in Equation (26) if and only if it has the representation
y ( z ) = c 2 + Δ a + h h υ f ( z , y ( z ) ) .
Proof. 
The proof follows from Theorem 10 directly. □
Remark 2.
In the case of Caputo, MVT for Caputo FFDE in Equation (23),
Δ a + υ h C h υ f ( s 1 ) Δ a + υ h C h υ g ( s 1 ) f ( b ) f ( a + υ h ) g ( b ) g ( a + υ h ) a + υ h C Δ h υ f ( s 2 ) a + υ h C Δ h υ g ( s 2 ) ,
does not hold, where the assumptions of Theorem 9 are supposed to be given. The reason for this is that we do not know whether Δ a + υ h C h υ g ( z ) > 0 when, by assumption and Corollary 3, we have
Δ a + υ h C h υ g ( z ) > 1 Γ ( 1 υ ) z ( a + υ h ) h ( υ ) g ( a + υ h ) z I < 0 ,
since the three quantities 1 Γ ( 1 υ ) , z ( a + υ h ) h ( υ ) , and g ( a + υ h ) are all positive for υ ( 0 , 1 ) and z I .

5. Conclusions

In this paper, a discrete υ -monotonicity analysis for discrete functions defined on M a , h in the framework of the discrete RL fractional sums, and RL and Caputo fractional differences on the time scale h Z were successfully studied. The relation between delta RL and delta Caputo fractional h-differences was established. The RL FFDE and Caputo FFDE were considered, and their solutions were discussed. Utilizing the monotonicity results and RL FFDE solution, the MVT was presented. However, MVT for Caputo FFDE was not valid using its corresponding monotonicity results and Caputo FFDE .
A further extension of this study includes improving our findings to study other classes of discrete fractional sums and differences, including those of exponential kernel Δ a C F R h υ y ( z ) or Mittag–Liffler kernel Δ a A B R h υ y ( z ) , defined in [22]. Fortunately, some research is in progress in this area.

Author Contributions

Conceptualization, P.O.M. and T.A.; data curation, F.K.H.; formal analysis, P.O.M. and T.A.; funding acquisition, T.A.; investigation, P.O.M.; methodology, F.K.H.; project administration, P.O.M. and T.A.; software, P.O.M. and T.A.; supervision, T.A. and F.K.H.; validation, F.K.H.; visualization, F.K.H.; writing—original draft, P.O.M.; writing—review and editing, T.A. and F.K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The second author would like to thank Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group number RG-DES-2017-01-17.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  2. Daftardar-Gejji, V. Fractional Calculus and Fractional Differential Equations; Springer: Singapore, 2019. [Google Scholar]
  3. 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]
  4. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin, Germany, 2010. [Google Scholar]
  5. Wu, J.; Yuan, J.; Gao, W. Analysis of fractional factor system for data transmission in SDN. Appl. Math. Nonlinear Sci. 2019, 4, 191–196. [Google Scholar] [CrossRef] [Green Version]
  6. Kurt, A.; Şenol, M.; Tasbozan, O.; Chand, M. Two reliable methods for the solution of fractional coupled Burgers’ equation arising as a model of Polydispersive sedimentation. Appl. Math. Nonlinear Sci. 2019, 4, 523–534. [Google Scholar] [CrossRef] [Green Version]
  7. Touchent, K.A.; Hammouch, Z.; Mekkaoui, T. A modified invariant subspace method for solving partial differential equations with non-singular kernel fractional derivatives. Appl. Math. Nonlinear Sci. 2020, 5, 35–48. [Google Scholar] [CrossRef]
  8. Hamasalh, F.K.; Mohammed, P.O. Generalized quartic fractional spline approximation function with applications. IJOPCM 2015, 8, 67–80. [Google Scholar] [CrossRef]
  9. Srivastava, H.M.; Mohammed, P.O.; Guirao, J.L.G.; Hamed, Y.S. Some higher-degree Lacunary fractional splines in the approximation of fractional differential equations. Symmetry 2021, 13, 422. [Google Scholar] [CrossRef]
  10. Onal, M.; Esen, A. A Crank-Nicolson approximation for the time fractional Burgers equation. Appl. Math. Nonlinear Sci. 2020, 5, 177–184. [Google Scholar] [CrossRef]
  11. Martinez, M.; Mohammed, P.O.; Valdes, J.E.N. Non-conformable fractional Laplace transform. Kragujev. J. Math. 2022, 46, 341–354. [Google Scholar]
  12. Wang, Y. The Green’s function of a class of two-term fractional differential equation boundary value problem and its applications. Adv. Differ. Equ. 2020, 2020, 80. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, S.; Hong, S. Variable separation method for a nonlinear time fractional partial differential equation with forcing term. J. Comput. Appl. Math. 2018, 339, 297–305. [Google Scholar] [CrossRef]
  14. Tarasov, V. Handbook of Fractional Calculus with Applications; Applications in Physics, Part A; De Gruyter: Boston, MA, USA; Berlin, Germany, 2019; Volume 4. [Google Scholar]
  15. Goodrich, C.; Peterson, A. Discrete Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  16. Goodrich, C. Existence of a positive solution to a system of discrete fractional boundary value problems. Appl. Math. Comput. 2011, 217, 4740–4753. [Google Scholar] [CrossRef]
  17. Atici, F.; Eloe, P. A transform method in discrete fractional calculus. Internat. J. Differ. Equ. 2007, 2, 165–176. [Google Scholar]
  18. Atici, F.; Eloe, P. Initial value problems in discrete fractional calculus. Proc. Am. Math. Soc. 2009, 137, 981–989. [Google Scholar] [CrossRef]
  19. Mohammed, P.O. A generalized uncertain fractional forward difference equations of Riemann-Liouville type. J. Math. Res. 2019, 11, 43–50. [Google Scholar] [CrossRef]
  20. Abdeljawad, T. On delta and nabla caputo fractional differences and dual identities. Discret. Dyn. Nat. Soc. 2013, 2013, 12. [Google Scholar] [CrossRef]
  21. Abdeljawad, T.; Jarad, F.; Atangana, A.; Mohammed, P.O. On a new type of fractional difference operators on h-step isolated time scales. J. Fract. Calc. Nonlinear Sys. 2021, 1, 46–74. [Google Scholar]
  22. Abdeljawad, T. Different type kernel h–fractional differences and their fractional h–sums. Chaos Solit. Fract. 2018, 116, 146–156. [Google Scholar] [CrossRef]
  23. Mohammed, P.O.; Hamasalh, F.K.; Abdeljawad, T. Difference monotonicity analysis on discrete fractional operators with discrete generalized Mittag-Leffler kernels. Adv. Differ. Equ. 2021, 2021, 213. [Google Scholar] [CrossRef]
  24. Mohammed, P.O.; Abdeljawad, T. Discrete generalized fractional operators defined using h-discrete Mittag-Leffler kernels and applications to AB fractional difference systems. Math. Meth. Appl. Sci. 2020, 1–26. [Google Scholar] [CrossRef]
  25. Abdeljawad, T.; Atici, F. On the Definitions of Nabla Fractional Operators. Abstr. Appl. Anal. 2012, 2012, 406757. [Google Scholar] [CrossRef]
  26. Suwan, I.; Owies, S.; Abdeljawad, T. Monotonicity results for h-discrete fractional operators and application. Adv. Differ. Equ. 2018, 2018, 207. [Google Scholar] [CrossRef]
  27. Atici, F.; Uyanik, M. Analysis of discrete fractional operators. Appl. Anal. Discr. Math. 2015, 9, 139–149. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mohammed, P.O.; Abdeljawad, T.; Hamasalh, F.K. On Riemann—Liouville and Caputo Fractional Forward Difference Monotonicity Analysis. Mathematics 2021, 9, 1303. https://0-doi-org.brum.beds.ac.uk/10.3390/math9111303

AMA Style

Mohammed PO, Abdeljawad T, Hamasalh FK. On Riemann—Liouville and Caputo Fractional Forward Difference Monotonicity Analysis. Mathematics. 2021; 9(11):1303. https://0-doi-org.brum.beds.ac.uk/10.3390/math9111303

Chicago/Turabian Style

Mohammed, Pshtiwan Othman, Thabet Abdeljawad, and Faraidun Kadir Hamasalh. 2021. "On Riemann—Liouville and Caputo Fractional Forward Difference Monotonicity Analysis" Mathematics 9, no. 11: 1303. https://0-doi-org.brum.beds.ac.uk/10.3390/math9111303

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