Open AccessArticle

The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator

Iqbal M. Batiha

^1,2

Shameseddin Alshorm

^1,*

Abdallah Al-Husban

³,

Rania Saadeh

⁴

Gharib Gharib

⁴ and

Shaher Momani

^2,5

Department of Mathematics, Al Zaytoonah University of Jordan, Amman 11733, Jordan

Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman P.O. Box 346, United Arab Emirates

Department of Mathematics, Faculty of Science and Technology, Irbid National University, Irbid P.O. Box 2600, Jordan

⁴

Department of Mathematics, Faculty of Science, Zarqa University, Zarqa 13110, Jordan

⁵

Department of Mathematics, The University of Jordan, Amman 11942, Jordan

Author to whom correspondence should be addressed.

Symmetry 2023, 15(4), 938; https://doi.org/10.3390/sym15040938

Submission received: 1 January 2023 / Revised: 20 February 2023 / Accepted: 21 February 2023 / Published: 19 April 2023

(This article belongs to the Special Issue Symmetry in Numerical Analysis and Computational Fluid Dynamics)

Download Versions Notes

Abstract

In this paper, we aim to present a novel n-point composite fractional formula for approximating a Riemann–Liouville fractional integral operator. With the use of the definite fractional integral’s definition coupled with the generalized Taylor’s formula, a novel three-point central fractional formula is established for approximating a Riemann–Liouville fractional integrator. Such a new formula, which emerges clearly from the symmetrical aspects of the proposed numerical approach, is then further extended to formulate an n-point composite fractional formula for approximating the same operator. Several numerical examples are introduced to validate our findings.

Keywords:

Richardson extrapolation; Riemann–Liouville fractional derivative and integral; Lagrange interpolating polynomial; Caputo derivative

1. Introduction

Fractional calculus has many uses in the mathematical modeling of chemical phenomena, physics, technical, and economics. It has contributed in a significant way in developing many topics and implementations in applied mathematics, although there are different definitions of the fractional-order operators, see [1,2,3,4,5,6]. In recent times, many evolutions in the theory of fractional calculus have been investigated to be employed in many fields of science and engineering. In particular, the fundamental of fractional calculus has been approved as an illustrious mathematical facility used to describe many actual applications [7,8]. As a result, various fractional-order differentiators and integrators have been established and approved by many researchers. It is important to highlight that there are two main operators for fractional-order differentiators; the first one is Caputo’s derivative operator with a power law function of convolution of a given function related to a local derivative, whereas the other one is the Riemann–Liouville derivative operator with a power law kernel of convolution [9]. In light of the various views of many mathematicians, the Caputo fractional-order differentiator has confirmed that it is more satisfactory for several real applications than that of the Riemann–Liouville derivative operator [10]. This is due to its suitability for using the assumed initial conditions when the fractional derivatives are taken [11]. Regardless of the best operator among the two former operators, the Riemann–Liouville integrator represents an inverse operator for both. This is because the Caputo differentiator is just a modification for the Riemann–Liouville differentiator [9,11].

The fractional-order integrator supposes that different constructs are not compatible and not constantly equivalent with each other. Actually, the fractional-order integrator is commonly employed for expressing an indefinite integral. In former research, there were only two endeavors attempted to establish a generalization of the fundamental theorem in fractional calculus [12,13]. In particular, different forms of such theorem were outlined in [12], whereas the authors in [13] proposed a more directed scheme into implementation than that of the first reference, but also the fractional definite integral definition had not been yet handled at both. Nor did the formulations, which were introduced in [14,15], outline the definition of the fractional definite integral. It looks as though nobody has formulated such a definition [16,17] except M. Ortigueira and J. Machado in their work in [18], whereby their fractional definite integral definition was set and the fundamental theorem of fractional calculus was consequently analyzed.

In light of the former discussion and based on what is usually carried out in the classical numerical analysis, we aim in this work to derive a novel formula enabling one to approximate a definite fractional integral called the n-point composite fractional formula for approximating a Riemann–Liouville fractional integrator. Thanks to O. Manuel and J. Machado’s definition of the definite fractional integral coupled with the generalized Taylor’s formula, a new three-point central fractional formula is first derived for approximating the Riemann–Liouville fractional integral operator. Then, this formula is extended to n composite points for formulating our main result in this work.

The remainder of this manuscript is constructed in the following manner. Section 2 aims to recollect some essentials and definitions regarding fractional calculus. Section 3 intends to illustrate the primary results of this work, whereby it will exemplify how we will derive the aimed n-point composite fractional formula for approximating the Riemann–Liouville fractional integrator. Section 4 provides two numerical examples that validate our theoretical outcomes, followed by the last section that outlines the concluding remarks of this work.

2. Preliminaries

In this part, some essential definitions and necessary preliminaries in connection with fractional calculus are recalled. This actually paves the way to our main results later on.

Definition 1

([19,20]). The fractional Riemann–Liouville integrator of a function

h (t)

of order

μ > 0

can be expressed as:

\begin{matrix} J^{μ} h (t) = \frac{1}{Γ (μ)} \int_{0}^{t} h (s) {(t - s)}^{μ - 1} d s, \end{matrix}

(1)

where

t > 0

and

μ > 0

In what follows, we recall certain properties of the fractional Riemann–Liouville integral operator for completeness:

\begin{matrix} 1) & J^{0} h (t) = h (t) . \end{matrix}

(2)

\begin{matrix} 2) & J^{μ} t^{γ} = \frac{Γ (γ + 1)}{Γ (μ + γ + 1)} t^{μ + γ}, γ \geq - 1 . \end{matrix}

(3)

\begin{matrix} 3) & J^{μ} J^{β} h (t) = J^{β} J^{μ} h (t) μ, β \geq 0 . \end{matrix}

(4)

\begin{matrix} 4) & J^{μ} J^{β} h (t) = J^{μ + β} h (t) μ, β \geq 0 . \end{matrix}

(5)

Definition 2

([19,20]). The Caputo fractional differentiator of a function

h (t)

of order

μ > 0

can be expressed as:

\begin{matrix} D_{t}^{μ} h (t) = \frac{1}{Γ (m - μ)} \int_{0}^{t} {(t - s)}^{m - μ - 1} h^{(m)} (s) d s, m - 1 < μ \leq m, \end{matrix}

(6)

where

m \in N

and

t > 0

In the following content, we list some properties of the Caputo differentiator [19,20]:

\begin{matrix} 1) & D_{t}^{μ} c = 0, where c is constant . \end{matrix}

(7)

\begin{matrix} 2) & D_{t}^{μ} t^{ρ} = \frac{Γ (ρ + 1)}{Γ (ρ - μ + 1)} t^{ρ - μ}, where ρ > μ - 1 . \end{matrix}

(8)

\begin{matrix} 3) & D_{t}^{μ} (μ h (t) + ω g (t)) = μ D_{t}^{μ} (h (t)) + ω D_{t}^{μ} (g (t)), \end{matrix}

(9)

where

μ

and

ω

are constant. In the same regard, we report below some other properties related to the composition between the previous two operators [19,20]:

D_{t}^{α} J^{α} h (t) = h (t),

(10)

and

J^{α} D_{t}^{α} h (t) = h (t) - \sum_{i = 1}^{n} h^{i} (0^{+}) \frac{t^{i}}{i!},

(11)

where

t > 0

and

n - 1 < α \leq n

such that

n \in N

Definition 3

([19,20]). The Riemann–Liouville fractional differentiator of a function

h (t)

of order

μ > 0

might be outlined in terms of the Riemann–Liouville fractional integrator as:

D^{μ} h (t) = D^{m} [J^{ρ} h (t)],

(12)

where

ρ = m - μ

0 < ρ < 1

and m is the smallest integral greater than μ.

Hereinafter, we recall two highly significant results that help us in deriving the main results of this work; the first one referred to M. Ortigueira and J. Machado, who established a proper formula to find the exact values of given definite fractional integrals [18], while the other one referred to Z. Odibat and S. Momani, who provide a generalization to the well-known Taylor theorem in its fractional case [20].

Definition 4

([18]). The definite fractional integral of the function f of order α is given by:

J_{a}^{α} f (x) = \int_{a}^{b} f^{(- α + 1)} (x) . d x = \int_{a}^{b} D_{a}^{- α + 1} f (x) d x,

(13)

where

- \infty < a < b < \infty

and

α - 1 < n \leq α

such that

n \in N

Theorem 1

([20]). (Generalized Taylor’s Theorem) Suppose that

D_{t}^{k α} f (x) \in C^{n + 1} (a, b]

for

k = 0, 1, \dots, n + 1,

where

0 < α \leq 1

. Then, the function f can be expanded about

x = x_{0}

as follows:

\begin{matrix} f (x) = \sum_{i = 0}^{n} \frac{{(x - x_{0})}^{i α}}{Γ (i α + 1)} D_{t}^{i α} f (x_{0}) + \frac{{(x - x_{0})}^{(n + 1) α}}{Γ ((n + 1) α + 1)} D_{t}^{(n + 1) α} f (ξ), \end{matrix}

(14)

where

0 < ξ < x

and

x \in (a, b]

3. Main Results

In this section, we aim to present the main results of this work. However, before we begin, we should first recall a very significant result related to the Caputo fractional differentiator. In fact, this result has been recently derived in [21] for the purpose of approximating such a differentiator,

D_{t}^{α}

, where

0 < α \leq 1

. It has been derived in light of using a similar manner used in [22], and it has been called the modified three-point fractional formula for approximating the Caputo fractional differentiator.

Theorem 2

([21]). Suppose that

f \in C^{3} [a, b]

and

x_{0}, x_{1}, x_{2}

are distinct points in the interval

[a, b]

such that

a = x_{0} < x_{1} = x_{0} + h < x_{2} = x_{0} + 2 h = b

with

h > 0

. Then, the modified three-point fractional formula for approximating the Caputo fractional differentiator can be given by:

\begin{matrix} D_{t}^{α} f (x) = & \frac{x^{2 - α}}{h^{2} Γ (3 - α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) \\ - \frac{x^{1 - α}}{2 h^{2} Γ (2 - α)} (f (x_{0}) (x_{1} + x_{2}) - 2 f (x_{1}) (x_{0} + x_{2}) + f (x_{2}) (x_{0} + x_{1})) \\ + \frac{f^{(3)} (ξ)}{6} (\frac{6}{Γ (4 - α)} x^{3 - α} - \frac{2 (x_{0} + x_{1} + x_{2})}{Γ (3 - α)} x^{2 - α} + \frac{(x_{0} x_{1} + x_{0} x_{2} + x_{1} x_{2})}{Γ (2 - α)} x^{1 - α}), \end{matrix}

(15)

for a generally unknown

ξ \in (a, b)

, where

x \in [a, b]

In consequence of the previous result, we can outline the next result that could help one to approximate the Caputo fractional differentiator

D_{t}^{2 α}

, where

0 < α \leq 1

Corollary 1.

Suppose that

f \in C^{3} [a, b]

and

x_{0}, x_{1}, x_{2}

are distinct points in the interval

[a, b]

such that

a = x_{0} < x_{1} = x_{0} + h < x_{2} = x_{0} + 2 h = b

with

h > 0

. Let

0 < α \leq 1

, then the modified three-point fractional formula for approximating Caputo fractional differentiator

D_{t}^{2 α}

can be given by:

D_{t}^{2 α} f (x) = \frac{x^{2 - 2 α}}{h^{2} Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) + \frac{f^{(3)} (ξ)}{6} (\frac{6 x^{3 - 2 α}}{Γ (4 - 2 α)} - \frac{2 (x_{0} + x_{1} + x_{2})}{Γ (3 - 2 α) x^{2 - 2 α}}) .

(16)

for a generally unknown

ξ \in (a, b)

, where

x \in [a, b]

Proof.

In order to prove this result, we simply apply all properties of the Caputo differentiator reported in (7) on (15). In particular, one can operate

D_{t}^{α}

twice again into the modified three-point fractional formula (15) to obtain the desired result. □

Remark 1.

It is worth mentioning that the well-known classical central three-point formulae reported for approximating the first and the second derivatives in [22] can be obtained easily by just substituting

α = 1

into (15) and (16), respectively.

In the following content, we aim to derive a new formula for approximating the Riemann–Liouville fractional integral operator

J^{α}

, which would be, from now on, called the three-point central fractional formula for approximating Riemann–Liouville fractional integrator. This will be carried out by applying on the generalized Taylor’s Theorem 1 coupled with considering definite fractional integral’s Definition 4.

Theorem 3.

Let

D_{t}^{n α} f \in C^{4} [a, b]

for

n = 0, 1, 2, 3, 4

and

0 < α \leq 1

. Suppose

a = x_{0} < x_{1} = x_{0} + h < x_{2} = x_{0} + 2 h = b

with

h > 0

. Then, the three-point central fractional formula about

x = x_{1}

for approximating a Riemann–Liouville fractional integrator is given by:

\begin{matrix} J^{α} f (x) & = 2 h f (x_{1}) + \frac{2 h^{3 α} x_{1}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) \\ + \frac{2 h^{3 α}}{6 Γ (3 α + 1)} f^{(3)} (ξ) (\frac{6 x_{1}^{(3 - 2 α)}}{Γ (4 - 2 α)} - \frac{2 (x_{0} + x_{1} + x_{2})}{Γ (3 - 2 α)} x_{1}^{(2 - 2 α)}) \\ + \frac{2 h^{5 α}}{Γ (5 α + 1)} D_{t}^{4 α} f (ξ), \end{matrix}

(17)

for a generally unknown

ξ \in (a, b)

, where

x \in [a, b]

Proof.

In order to prove this result, we first assume that all the above assumptions hold. Then, by expanding the function f about

x = x_{1}

using the generalized Taylor Theorem 1, we obtain:

\begin{matrix} f (x) = & f (x_{1}) + D_{t}^{α} f (x_{1}) \frac{{(x - x_{1})}^{α}}{Γ (α + 1)} + D_{t}^{2 α} f (x_{1}) \frac{{(x - x_{1})}^{2 α}}{Γ (2 α + 1)} \\ + D_{t}^{3 α} f (x_{1}) \frac{{(x - x_{1})}^{3 α}}{Γ (3 α + 1)} + D_{t}^{4 α} f (ξ) \frac{{(x - x_{1})}^{4 α}}{Γ (4 α + 1)}, \end{matrix}

(18)

for some

ξ \in (a, b)

. Now, by applying the Riemann–Liouville fractional integral operator to both sides of the above equality, we obtain:

\begin{matrix} J^{α} f (x) = & J^{α} f (x_{1}) + \frac{D_{t}^{α} f (x_{1})}{Γ (α + 1)} J^{α} [{(x - x_{1})}^{α}] + \frac{D_{t}^{2 α} f (x_{1})}{Γ (2 α + 1)} J^{α} [{(x - x_{1})}^{2 α}] \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (3 α + 1)} J^{α} [{(x - x_{1})}^{3 α}] + \frac{D_{t}^{4 α} f (ξ)}{Γ (4 α + 1)} J^{α} [{(x - x_{1})}^{4 α}] . \end{matrix}

(19)

With the help of using (13) coupled with considering

ρ = - α + 1

, we obtain:

\begin{matrix} J^{α} f (x) = & f (x_{1}) \int_{a}^{b} d x + \frac{D_{t}^{α} f (x_{1})}{Γ (α + 1)} \int_{a}^{b} D^{ρ} {(x - x_{1})}^{α} d x + \frac{D_{t}^{2 α} f (x_{1})}{Γ (2 α + 1)} \int_{a}^{b} D^{ρ} {(x - x_{1})}^{2 α} d x \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (3 α + 1)} \int_{a}^{b} D^{ρ} {(x - x_{1})}^{3 α} d x + \frac{D_{t}^{4 α} f (ξ)}{Γ (4 α + 1)} \int_{a}^{b} D^{ρ} {(x - x_{1})}^{4 α} d x . \end{matrix}

(20)

This consequently yields:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{D_{t}^{α} f (x_{1})}{Γ (α + 1)} \frac{Γ (α + 1)}{Γ (α + 1 - ρ)} \int_{a}^{b} {(x - x_{1})}^{α - ρ} d x \\ + \frac{D_{t}^{2 α} f (x_{1})}{Γ (2 α + 1)} \frac{Γ (2 α + 1)}{Γ (2 α + 1 - ρ)} \int_{a}^{b} {(x - x_{1})}^{2 α - ρ} d x \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (3 α + 1)} \frac{Γ (3 α + 1)}{Γ (3 α + 1 - ρ)} \int_{a}^{b} {(x - x_{1})}^{3 α - ρ} d x \\ + \frac{D_{t}^{4 α} f (ξ)}{Γ (4 α + 1)} \frac{Γ (4 α + 1)}{Γ (4 α + 1 - ρ)} \int_{a}^{b} {(x - x_{1})}^{4 α - ρ} d x . \end{matrix}

(21)

By simplifying the above equality, we obtain:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{D_{t}^{α} f (x_{1})}{Γ (2 α)} \int_{a}^{b} {(x - x_{1})}^{2 α - 1} d x + \frac{D_{t}^{2 α} f (x_{1})}{Γ (3 α)} \int_{a}^{b} {(x - x_{1})}^{3 α - 1} d x \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (4 α)} \int_{a}^{b} {(x - x_{1})}^{4 α - 1} d x + \frac{D_{t}^{4 α} f (ξ)}{Γ (5 α)} \int_{a}^{b} {(x - x_{1})}^{5 α - 1} d x, \end{matrix}

(22)

which immediately means:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{D_{t}^{α} f (x_{1})}{Γ (2 α)} \frac{{(x - x_{1})}^{2 α}}{(2 α)} |_{x_{0}}^{x_{2}} + \frac{D_{t}^{2 α} f (x_{1})}{Γ (3 α)} \frac{{(x - x_{1})}^{3 α}}{(3 α)} |_{x_{0}}^{x_{2}} \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (4 α)} \frac{{(x - x_{1})}^{4 α}}{(4 α)} |_{x_{0}}^{x_{2}} + \frac{D_{t}^{4 α} f (ξ)}{Γ (5 α)} \frac{{(x - x_{1})}^{5 α}}{(5 α)} |_{x_{0}}^{x_{2}} . \end{matrix}

(23)

This implies:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{D_{t}^{α} f (x_{1})}{Γ (2 α + 1)} ({(h)}^{2 α} - {(- h)}^{2 α}) + \frac{D_{t}^{2 α} f (x_{1})}{Γ (3 α + 1)} ({(h)}^{3 α} - {(- h)}^{3 α}) \\ + \frac{D_{t}^{3 α} f (x_{1})}{Γ (4 α + 1)} ({(h)}^{4 α} - {(- h)}^{4 α}) + \frac{D_{t}^{4 α} f (ξ)}{Γ (5 α + 1)} ({(h)}^{5 α} - {(- h)}^{5 α}) . \end{matrix}

(24)

One might observe that when

α = 1

, the above equality will be defined if

{(- h)}^{2 α} = {(h)}^{2 α}

{(- h)}^{3 α} = {- (h)}^{3 α}

{(- h)}^{4 α} = {(h)}^{4 α}

, and

{(- h)}^{5 α} = {- (h)}^{5 α}

. Based on this observation, we can assert the following equation:

\begin{matrix} J^{α} f (x) = 2 h f (x_{1}) + \frac{2 h^{3 α}}{Γ (3 α + 1)} D_{t}^{2 α} f (x_{1}) + \frac{2 h^{5 α}}{Γ (5 α + 1)} D_{t}^{4 α} f (ξ) . \end{matrix}

(25)

Now, by substituting the modified three-point fractional formula for approximating Caputo fractional differentiator

D_{t}^{2 α}

reported in (16) into (25), we obtain:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{2 h^{3 α}}{Γ (3 α + 1)} [\frac{x_{1}^{2 - 2 α}}{Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) \\ + \frac{f^{(3)} (ξ)}{6} (\frac{6 x_{1}^{(3 - 2 α)}}{Γ (4 - 2 α)} - \frac{2 (x_{0} + x_{1} + x_{2})}{Γ (3 - 2 α)} x_{1}^{(2 - 2 α)})] + \frac{2 h^{5 α}}{Γ (5 α + 1)} D_{t}^{4 α} f (ξ), \end{matrix}

(26)

for some

ξ \in (a, b)

. Simplifying the above equation yields:

\begin{matrix} J^{α} f (x) = & 2 h f (x_{1}) + \frac{2 h^{3 α} x_{1}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) \\ + \frac{2 h^{3 α}}{6 Γ (3 α + 1)} f^{(3)} (ξ) (\frac{6 x_{1}^{(3 - 2 α)}}{Γ (4 - 2 α)} - \frac{2 (x_{0} + x_{1} + x_{2})}{Γ (3 - 2 α)} x_{1}^{(2 - 2 α)}) + \frac{2 h^{5 α}}{Γ (5 α + 1)} D^{4 α} f (ξ), \end{matrix}

(27)

which represents the aimed desired result. □

In light of the previous discussions, one could notice that the error term reported in (27) is too high when one wants to estimate it. From this point of view and to obtain more accuracy in regard to approximating the Riemann–Liouville fractional integrator, we introduce in what followed, one of the main results for this work, the n-point composite fractional formula for approximating the Riemann–Liouville fractional integrator. This would be achieved by tracking the same procedure used to obtain the n-point composite Simpson/trapezoidal formulae [22].

Corollary 2.

Let

D_{t}^{k α} f \in C^{4} [a, b]

for

k = 0, 1, 2, 3, 4

and

0 < α \leq 1

. Suppose

a = x_{0} < x_{1} = x_{0} + h < x_{2} = x_{0} + 2 h < \dots < x_{n} = x_{0} + n h = b

with

h > 0

and

n \geq 2

. Then, the n-point composite fractional formula for approximating the Riemann–Liouville fractional integrator is given by:

\begin{matrix} J^{α} f (x) & \approx 2 h \sum_{i = 0}^{\frac{n - 2}{2}} f (x_{2 i + 1}) + \frac{2 h^{3 α}}{Γ (3 α + 1) Γ (3 - 2 α)} \prod_{i = 0}^{\frac{n - 2}{2}} x_{2 i + 1} (f (x_{2 i}) - 2 f (x_{2 i + 1}) + f (x_{2 i + 2})), \end{matrix}

(28)

for some

ξ \in (a, b)

, where

x \in [a, b]

Proof.

To prove this result, one may rewrite

J^{α} f (x)

as:

J^{α} f (x) = J_{x_{0}}^{α} f (x) + J_{x_{2}}^{α} f (x) + \dots + J_{x_{n - 2}}^{α} f (x) .

(29)

With the use of (1), the above equality can be re-expressed as:

\begin{matrix} J^{α} f (x) & = \frac{1}{Γ (α)} \int_{a}^{x} {(x - t)}^{α - 1} f (t) d t \\ = \frac{1}{Γ (α)} \int_{x_{0}}^{x} {(x - t)}^{α - 1} f (t) + \frac{1}{Γ (α)} \int_{x_{2}}^{x} {(x - t)}^{α - 1} f (t) + \dots + \frac{1}{Γ (α)} \int_{x_{n - 2}}^{x} {(x - t)}^{α - 1} f (t) . \end{matrix}

(30)

Now, due to the approximate version of the three-point central fractional formula for approximating the Riemann–Liouville fractional integrator has the form:

\begin{matrix} J^{α} f (x) \approx 2 h f (x_{1}) + \frac{2 h^{3 α} x_{1}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})), \end{matrix}

(31)

then equality (30) with the use of (31) can yield:

\begin{matrix} J^{α} f (x) & \approx [2 h f (x_{1}) + \frac{2 h^{3 α} x_{1}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{0}) - 2 f (x_{1}) + f (x_{2}))] \\ + [2 h f (x_{3}) + \frac{2 h^{3 α} x_{3}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{2}) - 2 f (x_{3}) + f (x_{4}))] \\ + [2 h f (x_{5}) + \frac{2 h^{3 α} x_{5}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{4}) - 2 f (x_{5}) + f (x_{6}))] \\ + \\ ⋮ \\ + [2 h f (x_{n - 1}) + \frac{2 h^{3 α} x_{n - 1}^{2 - 2 α}}{Γ (3 α + 1) Γ (3 - 2 α)} (f (x_{n - 2}) - 2 f (x_{n - 1}) + f (x_{n}))] . \end{matrix}

(32)

Immediately, simplifying the above equation gives:

\begin{matrix} J^{α} f (x) = & 2 h (f (x_{1}) + f (x_{3}) + f (x_{5}) + \dots + f (x_{n - 1})) + \frac{2 h^{3 α}}{Γ (3 α + 1) Γ (3 - 2 α)} \\ \times [ & x_{1}^{2 - 2 α} (f (x_{0}) - 2 f (x_{1}) + f (x_{2})) \\ + x_{3}^{2 - 2 α} (f (x_{2}) - 2 f (x_{3}) + f (x_{4})) \\ + x_{5}^{2 - 2 α} (f (x_{4}) - 2 f (x_{5}) + f (x_{6})) \\ + \\ ⋮ \\ + x_{n - 1}^{2 - 2 α} (f (x_{n - 2}) - 2 f (x_{n - 1}) + f (x_{n}))] . \end{matrix}

(33)

In other words, we have the following composite formula:

\begin{matrix} J^{α} f (x) & \approx 2 h \sum_{i = 0}^{\frac{n - 2}{2}} f (x_{2 i + 1}) + \frac{2 h^{3 α}}{Γ (3 α + 1) Γ (3 - 2 α)} \prod_{i = 0}^{\frac{n - 2}{2}} x_{2 i + 1} (f (x_{2 i}) - 2 f (x_{2 i + 1}) + f (x_{2 i + 2})) . \end{matrix}

(34)

□

4. Numerical Examples

In this part, we intend to verify the effectiveness of our final proposed formula; the n-point composite fractional formula for approximating the Riemann–Liouville integrator for two specific functions. Tables are utilized to display and compare the gained findings.

Example 1.

Let us consider the main function

f (x) = 2 x^{3} + 8 x

. Suppose the fractional-order has the values

α = 1, 0.75, 0.5

and the interval is

[0, 2]

with

n = 10

. Then, applying Riemann–Liouville integrator generates the following cases:

Case 1: When $α = 1$ , we obtain:

$\begin{matrix} J^{1} f (x) & = 0.5 x^{4} + 4 x^{2} . \end{matrix}$

(35)
Case 2: When $α = 0.75$ , we obtain:

$\begin{matrix} J^{0.75} f (x) & = 0.7234 x^{3.75} + 6.0182 x^{1.75} . \end{matrix}$

(36)
Case 3: When $α = 0.5$ , we obtain:

$\begin{matrix} J^{0.5} f (x) & = 1.0316 x^{3.5} + 7.7549 x^{1.5} . \end{matrix}$

(37)

The results from the above expressions (35)–(37) are compared with the approximate values generated by using the n-point composite fractional formula. Accordingly, a numerical comparison is generated for different values of α and reported in Table 1.

It is worth mentioning that the average elapsed time (or the CPU time) that is required to give the approximate value for each of the above Riemann–Liouville integrators using the n-point composite fractional formula is 0.096549 s, which is regarded too short in comparison with finding their exact values analytically.

Example 2.

Consider now the main function is

f (x) = s i n (x) = \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{(2 k + 1)!} x^{2 k + 1}

, and suppose the fractional-order has the values

α = 1, 0.65, 0.25

. The interval here is also chosen as

[0, 2]

with

n = 10

. Then, by applying the Riemann–Liouville integrator, we obtain the following cases:

Case 1: When $α = 1$ , we have:

$\begin{matrix} J^{1} f (x) & = \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{(2 k + 1)!} \frac{Γ (2 k + 2)}{Γ (2 k + 3)} x^{2 k + 2} . \end{matrix}$

(38)
Case 2: When $α = 0.65$ , we have:

$\begin{matrix} J^{0.65} f (x) & = \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{(2 k + 1)!} \frac{Γ (2 k + 2)}{Γ (2 k + 2.65)} x^{2 k + 1.65} . \end{matrix}$

(39)
Case 3: When $α = 0.45$ , we have:

$\begin{matrix} J^{0.45} f (x) & = \sum_{k = 0}^{\infty} \frac{{(- 1)}^{k}}{(2 k + 1)!} \frac{Γ (2 k + 2)}{Γ (2 k + 2.45)} x^{2 k + 1.45} . \end{matrix}$

(40)

Consequently, the results from applying the above three expressions (38)–(40) are compared with the approximate values generated by using the n-point composite fractional formula. As a result, another numerical comparison is generated for different values of α and reported in Table 2.

Herein, the average elapsed time required to give the approximate value for each of the above Riemann–Liouville integrators using the n-point composite fractional formula is 0.099873 s. This time is also regarded as too short in comparison with finding their exact values analytically.

In light of the previous examples, one might notice that the proposed n-point composite fractional formula can provide a good approximation for the Riemann–Liouville fractional integral operator as compared with the provided exact values. This inference can lead us in the close future to apply such a formula in several applications, especially those applications that are related to differential equations.

5. Conclusions

In this paper, a new numerical formula called the n-point composite fractional formula for approximating Riemann–Liouville fractional integrators has been successfully established. In particular, this formula has been derived based on establishing another novel formula called the three-point central fractional formula for approximating Riemann–Liouville fractional integrators. Two main notions have been used for attaining these results; the generalized Taylor theorem and the definite fractional integral. Several numerical applications, including solving fractional differential equations, partial fractional differential equations and others, have been left to future considerations.

Author Contributions

Conceptualization, S.A.; Methodology, I.M.B., S.A., R.S. and S.M.; Software, S.A. and G.G.; Validation, S.A.; Formal analysis, S.A. and A.A.-H.; Investigation, S.A.; Data curation, S.A.; Supervision, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Aho, A.V.; Hopcroft, J.E.; Ullman, J.D. The Design and Analysis of Computer Algorithms; Addison-Wesley: Reading, MA, USA, 1974. [Google Scholar]
Batiha, I.M.; Obeidat, A.; Alshorm, S.; Alotaibi, A.; Alsubaie, H.; Momani, S.; Albdareen, M.; Zouidi, F.; Eldin, S.M.; Jahanshahi, H. A Numerical Confirmation of a Fractional-Order COVID-19 Model’s Efficiency. Symmetry 2022, 14, 2583. [Google Scholar] [CrossRef]
Batiha, I.M.; El-Khazali, R.; Ababneh, O.Y.; Ouannas, A.; Batyha, R.M.; Momani, S. Optimal design of PI^ρD^μ-controller for artificial ventilation systems for COVID-19 patients. Aims Math. 2023, 8, 657–675. [Google Scholar] [CrossRef]
Williams, W.K.; Vijayakumar, V.; Nisar, K.S.; Shukla, A. Atangana–Baleanu Semilinear Fractional Differential Inclusions with Infinite Delay: Existence and Approximate Controllability. ASME J. Comput. Nonlinear Dynam. 2023, 18, 021005. [Google Scholar] [CrossRef]
Bose, C.S.V.; Udhayakumar, R. Analysis on the Controllability of Hilfer Fractional Neutral Differential Equations with Almost Sectorial Operators and Infinite Delay via Measure of Noncompactness. Qual. Theory Dyn. Syst. 2023, 22, 22. [Google Scholar] [CrossRef]
Dineshkumar, C.; Udhayakumar, R.; Vijayakumar, V.; Shukla, A.; Nisar, K.S. Discussion on the Approximate Controllability of Nonlocal Fractional Derivative by Mittag-Leffler Kernel to Stochastic Differential Systems. Qual. Theory Dyn. Syst. 2023, 22, 27. [Google Scholar] [CrossRef]
Dababneh, A.; Djenina, N.; Ouannas, A.; Grassi, G.; Batiha, I.M.; Jebril, I.H. A New Incommensurate Fractional-Order Discrete COVID-19 Model with Vaccinated Individuals Compartment. Fractal Fract. 2022, 6, 456. [Google Scholar] [CrossRef]
Ouannas, A.; Batiha, I.M.; Bekiros, S.; Liu, J.; Jahanshahi, H.; Aly, A.A.; Alghtani, A.H. Synchronization of the Glycolysis Reaction-Diffusion Model via Linear Control Law. Entropy 2021, 23, 1516. [Google Scholar] [CrossRef] [PubMed]
Albadarneh, R.B.; Batiha, I.M.; Adwai, A.; Tahat, N.; Alomari, A.K. Numerical approach of Riemann-Liouville fractional derivative operator. Int. J. Electr. Comput. Eng. 2021, 11, 5367–5378. [Google Scholar] [CrossRef]
Bezziou, M.; Dahmani, Z.; Jebril, I.; Kaid, M. Caputo-hadamard approach applications: Solvability for an integro-differential problem of lane and emden type. J. Math. Comput. Sci. 2021, 11, 1629–1649. [Google Scholar]
Albadarneh, R.B.; Batiha, I.; Alomari, A.K.; Tahat, N. Numerical approach for approximating the Caputo fractional-order derivative operator. Aims Math. 2021, 6, 12743–12756. [Google Scholar] [CrossRef]
Grigoletto, E.C.; de Oliveira, E.C. Fractional Versions of the Fundamental Theorem of Calculus. Appl. Math. 2013, 4, 23–33. [Google Scholar] [CrossRef]
Tarasov, V.E. Fractional vector calculus and fractional Maxwell’s equations. Ann. Phys. 2008, 323, 2756–2778. [Google Scholar] [CrossRef]
Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Gordon and Breach Science Publishers: London, UK, 1993. [Google Scholar]
Machado, J.T.; Mainardi, F.; Kiryakova, V. Fractional Calculus: Quo Vadimus, 495–526. Fract. Calc. Appl. Anal. 2015, 18. [Google Scholar] [CrossRef]
Machado, J.T.; Mainardi, F.; Kiryakova, V.; Atanacković, T. Fractional Calculus: D’où venos-nous? Que sommes-Nous Où allons-Nous? Fract. Calc. Appl. Anal. 2016, 19, 1074–1104. [Google Scholar] [CrossRef]
Ortigueira, M.; Machado, J. Fractional Definite Integral. Fractal Fract. 2017, 1, 2. [Google Scholar] [CrossRef]
Batiha, I.M.; El-Khazali, R.; AlSaedi, A.; Momani, S. The General Solution of Singular Fractional-Order Linear Time-Invariant Continuous Systems with Regular Pencils. Entropy 2018, 20, 400. [Google Scholar] [CrossRef] [PubMed]
Odibat, M.Z.; Momani, S. An algorithm for the numerical solution of differential equations of fractional order. J. Appl. Math. Inform. 2008, 26, 15–27. [Google Scholar]
Batiha; Alshorm, I.M.S.; Ouannas, A.; Momani, S.; Ababneh, O.Y.; Albdareen, M. Modified Three-Point Fractional Formulas with Richardson Extrapolation. Mathematics 2022, 10, 3489. [Google Scholar] [CrossRef]
Burden, R.L.; Faires, J.D. Numerical Analysis, 9th ed.; Brooks/Cole: Pacific Grove, CA, USA, 2011. [Google Scholar]

Table 1. A numerical comparison between the exact and approximate values.

$n ∖ α$	0.5	0.75	1
2	22.27105816166	26.16734499344	23.68537600000
4	22.69349110102	26.26331428545	23.69932800000
6	24.03916177943	26.50011822893	23.72599466666
8	27.18307909882	26.96769964223	23.77049600000
10	28.49632000000	26.79576811488	23.84000000000
Exact:	28.69344858658	26.46463660779	24

Table 2. A numerical comparison between the exact and approximate values.

$n ∖ α$	0.45	0.65	1
2	0.988434759968	1.109574041883	1.416657631095
4	1.046664086456	1.129294584780	1.418163344357
6	1.098873318034	1.271316766203	1.420407266983
8	1.156966874475	1.233598909270	1.423035132930
10	1.293240352453	1.306986260607	1.425632059945
Exact:	1.273563291144	1.362300353262	1.416146836547

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.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Batiha, I.M.; Alshorm, S.; Al-Husban, A.; Saadeh, R.; Gharib, G.; Momani, S. The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator. Symmetry 2023, 15, 938. https://doi.org/10.3390/sym15040938

AMA Style

Batiha IM, Alshorm S, Al-Husban A, Saadeh R, Gharib G, Momani S. The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator. Symmetry. 2023; 15(4):938. https://doi.org/10.3390/sym15040938

Chicago/Turabian Style

Batiha, Iqbal M., Shameseddin Alshorm, Abdallah Al-Husban, Rania Saadeh, Gharib Gharib, and Shaher Momani. 2023. "The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator" Symmetry 15, no. 4: 938. https://doi.org/10.3390/sym15040938

APA Style

Batiha, I. M., Alshorm, S., Al-Husban, A., Saadeh, R., Gharib, G., & Momani, S. (2023). The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator. Symmetry, 15(4), 938. https://doi.org/10.3390/sym15040938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

The n-Point Composite Fractional Formula for Approximating Riemann–Liouville Integrator

Abstract

1. Introduction

2. Preliminaries

3. Main Results

4. Numerical Examples

5. Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI