# Determination of Coefficients of High-Order Schemes for Riemann-Liouville Derivative

Authors: Rifang Wu; Hengfei Ding; Changpin Li
Journal: The Scientific World Journal
 (2014)

---

## Abstract

Although there have existed some numerical algorithms for the fractional differential equations, developing high-order methods (i.e., with convergence order greater than or equal to 2) is just the beginning. Lubich has ever proposed the high-order schemes when he studied the fractional linear multistep methods, where he constructed thepth order schemes (
p
=
2
,
3
,
4
,
5
,
6
) for the αth order Riemann-Liouville integral and αth order Riemann-Liouville derivative. In this paper, we study such a problem and develop recursion formulas to compute these coefficients in the higher-order schemes. The coefficients of higher-order schemes (
p
=
7,8
,
9,10
) are also obtained. We first find that these coefficients are oscillatory, which is similar to Runge’s phenomenon. So, they are not suitable for numerical calculations. Finally, several numerical examples are implemented to testify the efficiency of the numerical schemes for p
=
3
,
…
,
6.


---

## Body

## 1. Introduction

Loosely speaking, fractional calculus is often regarded as the generalization of classical calculus. From the viewpoint of rigorous mathematics, this is not the case; see [1]. Now, fractional calculus has been successfully applied in the fields of chemistry, physics, finance, signal processing, bioengineering, and control. For details, see [2–11] and references cited therein.Although there are more than six kinds of fractional derivatives, the commonly used derivatives are Riemann-Liouville, Grünwald-Letnikov, and Caputo ones. In this paper, we focus on Riemann-Liouville derivative. Under suitable conditions, the Riemann-Liouville derivative can be discretized by the discrete form of the Grünwald-Letnikov one.Some numerical approximate formulas for fractional calculus have been proposed [4, 12–24]. It is worth mentioning that the fractional linear multistep (also high-order) methods for the Riemann-Liouville integrals and derivatives were firstly proposed in [16]. The high-order numerical methods for Caputo derivatives were firstly constructed in [15]. The high-order algorithms for Riesz derivative was firstly considered by Ding et al. [25].In the following, we give some basic definitions, notations, and properties of the fractional calculus [1, 3, 7, 26].Definition 1.
Letf be defined on the interval [
a
,
b
] and α
>
0. Then, the left Riemann-Liouville integral of order α is defined as
(1)
D
a
,
x
-
α
f
(
x
)
=
1
Γ
(
α
)
∫
a
x
‍
f
(
t
)
(
x
-
t
)
1
-
α
d
t
,
where Γ
(
·
) is Euler’s gamma function.Definition 2.
Letf be defined on the interval (
a
,
b
] and α
>
0 and let n be the smallest integer greater than α
(
n
-
1
≤
α
<
n
). Then, the left Riemann-Liouville derivative of order α is defined by
(2)
D
RL
a
,
x
α
f
(
x
)
=
D
n
D
a
,
x
-
(
n
-
α
)
f
(
x
)
=
1
Γ
(
n
-
α
)
(
d
d
x
)
n
∫
a
x
‍
f
(
t
)
(
x
-
t
)
α
-
n
+
1
d
t
.Definition 3.
Letf be defined on the interval [
a
,
b
] and α
>
0 and let n be the smallest integer greater than α
(
n
-
1
<
α
≤
n
). Then, the left Caputo derivative of order α is defined as follows:
(3)
D
C
a
,
x
α
f
(
x
)
=
D
a
,
x
-
(
n
-
α
)
D
n
f
(
x
)
=
1
Γ
(
n
-
α
)
∫
a
x
‍
f
(
n
)
(
t
)
(
x
-
t
)
α
-
n
+
1
d
t
.Definition 4.
Letf be defined on the interval [
a
,
b
] and α
>
0. Then, the left Grünwald-Letnikov derivative of order α is defined as
(4)
D
GL
a
,
x
α
f
(
x
)
=
lim
⁡
h
→
0
m
h
=
x
-
a
1
h
α
∑
ℓ
=
0
m
‍
(
-
1
)
ℓ
(
α
ℓ
)
f
(
x
-
ℓ
h
)
.Lemma 5.
Suppose thatf
(
x
) is differentiable in the senses of both Caputo and Riemann-Liouville. Then,
(5)
D
C
a
,
x
α
f
(
x
)
=
D
R
L
a
,
x
α
f
(
x
)
-
∑
k
=
0
n
-
1
‍
f
(
k
)
(
a
)
Γ
(
k
+
1
-
α
)
(
x
-
a
)
k
-
α
,
where n
-
1
<
α
<
n
∈
ℤ
+.Lemma 6.
Letf
∈
C
n
[
a
,
b
]; then the finite Grünwald-Letnikov derivative
(6)
D
~
G
L
a
,
x
α
f
(
x
)
=
1
h
α
∑
ℓ
=
0
[
(
x
-
a
)
/
h
]
‍
(
-
1
)
ℓ
(
α
ℓ
)
f
(
x
-
ℓ
h
)
yields a first-order approximation for the Riemann-Liouville derivative D
R
L
a
,
x
α
f
(
x
) if f
(
a
+
)
=
0; that is,
(7)
D
R
L
a
,
x
α
f
(
x
)
=
D
~
G
L
a
,
x
α
f
(
x
)
+
𝒪
(
h
)
.
But when f
(
a
+
)
≠
0, then one has
(8)
D
R
L
a
,
x
α
f
(
x
)
=
D
~
G
L
a
,
x
α
f
(
x
)
+
𝒪
(
h
)
+
𝒪
(
f
(
a
+
)
)
.The present paper is organized as follows. We divide Section2 into three subsections. In Section 2.1, we introduce the recursion formulas of coefficients of orders 1 and 2, which are usually used, and some properties are given. In Section 2.2, coefficients of orders from 3 to 6 are presented for reference. And coefficients of orders from 7 to 10 are displayed in Section 2.3 which are oscillatory. Such a phenomenon is similar to Runge’s phenomenon. So, the high-order (more than 6th-order) schemes seem not to be suitable for numerical calculations, which is the same as the case of ordinary differential equation. In Section 3, some numerical experiments are carried out to support the computational schemes. Section 4 concludes this paper.


## 2. The Determination of the Coefficients

Firstly, we divide the given interval[
a
,
b
] into
(9)
Λ
:
a
=
x
0
<
x
1
<
⋯
<
x
M
=
b
and x
m
=
a
+
m
h, in which h
=
(
b
-
a
)
/
M, m
=
0,1
,
…
,
M.In most situations, we naturally use the following formula to approximate the Riemann-Liouville derivative:(10)
D
RL
a
,
x
α
f
(
x
m
)
=
1
h
α
∑
ℓ
=
0
m
‍
ϖ
1
,
ℓ
(
α
)
f
(
x
m
-
ℓ
)
+
𝒪
(
h
)
+
𝒪
(
f
(
a
+
)
)
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
m
=
0,1
,
…
,
M
,
where
(11)
ϖ
1
,
ℓ
(
α
)
=
(
-
1
)
ℓ
(
α
ℓ
)
=
(
-
1
)
ℓ
Γ
(
α
+
1
)
Γ
(
ℓ
+
1
)
Γ
(
α
-
ℓ
+
1
)
,
hhhhhhhhhhhhhhhhhhhhhh
ℓ
=
0,1
,
…
,
m
,
are the binomial coefficients. And they have the following recurrence relationships:
(12)
ϖ
1,0
(
α
)
=
1
,
ϖ
1
,
ℓ
(
α
)
=
(
1
-
1
+
α
ℓ
)
ϖ
1
,
ℓ
-
1
(
α
)
,
ℓ
=
1,2
,
…
,
m
.Obviously,ϖ
1
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
,
m
) are just the first m
+
1 coefficients of Taylor series of the expansion of the following function:
(13)
W
1
(
α
)
(
z
)
=
(
1
-
z
)
α
=
∑
ℓ
=
0
∞
‍
(
-
1
)
ℓ
(
α
ℓ
)
z
ℓ
=
∑
ℓ
=
0
∞
‍
ϖ
1
,
ℓ
(
α
)
z
ℓ
.We can see that formula (10) has only the first-order accuracy if f
(
a
+
)
=
0. Therefore, to seek high accurate numerical methods for the fractional derivatives is of great importance. In [16], Lubich firstly proposed numerical schemes of orders 2, 3, 4, 5, and 6 Riemann-Liouvile derivative, named fractional linear multistep formulas. Here, it must be mentioned that the fractional linear multistep method is different from the usual linear multistep method. The former is of varied steps. That is to say, the value of the mth step x
m relies on the preceding step values x
0
,
x
1
,
…
,
x
m
-
1, which means that the number of multisteps is increasing, while the latter is of fixed number of multisteps. Under the homogeneous initial value conditions, that is, f
(
k
)
(
a
+
)
=
0
(
k
=
0,1
,
…
,
p
-
1
), the fractional linear multistep scheme for the Riemann-Liouville derivative has the following form:
(14)
D
RL
a
,
x
α
f
(
x
m
)
=
1
h
α
∑
ℓ
=
0
m
‍
ϖ
p
,
ℓ
(
α
)
f
(
x
m
-
ℓ
)
+
𝒪
(
h
p
)
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
p
=
2
,
…
,
10
.Ifp
≥
2, the corresponding numerical methods are often called the high-order methods. The coefficients ϖ
p
,
ℓ
(
α
) in the above equations are those of the Taylor series expansions of the corresponding generating functions W
p
(
α
)
(
z
); that is,
(15)
W
p
(
α
)
(
z
)
=
∑
ℓ
=
0
∞
‍
ϖ
p
,
ℓ
(
α
)
z
ℓ
,
p
=
2
,
…
,
10
.
The corresponding generating functions for p
=
2
,
…
,
6 are given as follows [16]:
(16)
W
2
(
α
)
(
z
)
=
(
3
2
-
2
z
+
1
2
z
2
)
α
,
W
3
(
α
)
(
z
)
=
(
11
6
-
3
z
+
3
2
z
2
-
1
3
z
3
)
α
,
W
4
(
α
)
(
z
)
=
(
25
12
-
4
z
+
3
z
2
-
4
3
z
3
+
1
4
z
4
)
α
,
W
5
(
α
)
(
z
)
=
(
137
60
-
5
z
+
5
z
2
-
10
3
z
3
+
5
4
z
4
-
1
5
z
5
)
α
,
W
6
(
α
)
(
z
)
=
(
147
60
-
6
z
+
15
2
z
2
-
20
3
z
3
+
15
4
z
4
h
h
h
-
6
5
z
5
+
1
6
z
6
)
α
.
Using the similar method, we list the generating functions for p
=
7
,
…
,
10,
(17)
W
7
(
α
)
(
z
)
=
(
363
140
-
7
z
+
21
2
z
2
-
35
3
z
3
+
35
4
z
4
h
h
h
-
21
5
z
5
+
7
6
z
6
-
1
7
z
7
)
α
,
W
8
(
α
)
(
z
)
=
(
761
280
-
8
z
+
14
z
2
-
56
3
z
3
+
35
2
z
4
h
h
k
-
56
5
z
5
+
14
3
z
6
-
8
7
z
7
+
1
8
z
8
)
α
,
W
9
(
α
)
(
z
)
=
(
7129
2520
-
9
z
+
18
z
2
-
28
z
3
+
63
2
z
4
-
126
5
z
5
h
h
h
+
14
z
6
-
36
7
z
7
+
9
8
z
8
-
1
9
z
9
)
α
,
W
10
(
α
)
(
z
)
=
(
14762
5040
-
10
z
+
45
2
z
2
-
40
z
3
+
105
2
z
4
-
252
5
z
5
h
h
h
+
35
z
6
-
120
7
z
7
+
45
8
z
8
-
10
9
z
9
+
1
10
z
10
)
α
.From the above introduction, the key question is how to compute the coefficientsϖ
p
,
ℓ
(
α
), p
=
2
,
…
,
10. As far as we know, there are three methods to compute these coefficients, in which one way is to use the fast Fourier transform method [19].Another way of computing the coefficientsϖ
p
,
ℓ
(
α
), p
=
2
,
…
,
10, is by the automatic differentiation techniques [13]
(18)
ϖ
p
,
ℓ
(
α
)
=
1
ℓ
u
0
∑
j
=
0
ℓ
-
1
‍
[
α
(
ℓ
-
j
)
-
j
]
ϖ
p
,
j
(
α
)
u
ℓ
-
j
,
p
=
2
,
…
,
10
.
Here, the values u
k
(
k
=
0,1
,
…
) in the above formula denote the Taylor expansion coefficients of the generating functions W
p
(
α
)
(
z
)
(
p
=
2
,
…
,
10
) of the classical linear multistep methods.The third method is to use the expansion of series, where the coefficients are not given in the form of recurrence relationships but the exact expressions; see the Appendix. Such analytical expressions are useful for theoretical analysis, such as stability and convergence. Here, we introduce the second method.


### 2.1. Properties ofϖ
1
,
ℓ
(
α
) and ϖ
2
,
ℓ
(
α
)

As far as we know, the coefficients of the 1st-order and 2nd-order schemes are often used. Here, we present their properties for reference. We know thatϖ
1
,
ℓ
(
α
) are the binomial coefficients in formula (12). They have the following properties.Property 1.
The coefficientsϖ
1
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 0
<
α
<
1 have the following properties:(1)
ϖ
1,0
(
α
)
=
1, ϖ
1
,
ℓ
(
α
)
=
(
-
1
)
ℓ
Γ
(
α
+
1
)
/
(
Γ
(
ℓ
+
1
)
Γ
(
α
-
ℓ
+
1
)
)
<
0, ℓ
≥
1;
(2)
ϖ
1
,
ℓ
-
1
(
α
)
≤
ϖ
1
,
ℓ
(
α
), ℓ
≥
2;
(3)
∑
ℓ
=
0
∞
‍
ϖ
1
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
1
,
ℓ
(
α
)
=
0.Parallelly, the properties of the coefficients of the 1st-order schemes for1
<
α
<
2 are given as follows.Property  1′. The coefficients ϖ
1
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 1
<
α
<
2 have the following properties:(1)
ϖ
1,0
(
α
)
=
1, ϖ
1,1
=
-
α
<
0, ϖ
1
,
ℓ
(
α
)
=
(
-
1
)
ℓ
Γ
(
α
+
1
)
/
(
Γ
(
ℓ
+
1
)
Γ
(
α
-
ℓ
+
1
)
)
>
0, ℓ
≥
2;
(2)
ϖ
1
,
ℓ
-
1
(
α
)
≥
ϖ
1
,
ℓ
(
α
), ℓ
≥
3;
(3)
∑
ℓ
=
0
∞
‍
ϖ
1
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
1
,
ℓ
(
α
)
=
0.The above Properties 1 and 1′ are known to us. Next, we present the coefficientsϖ
2
,
ℓ
(
α
) which are somewhat complex. Let i
=
ℓ
-
j; formula (18) can be written as
(19)
ϖ
p
,
ℓ
(
α
)
=
1
ℓ
u
0
∑
i
=
1
ℓ
‍
(
α
i
-
ℓ
+
i
)
ϖ
p
,
ℓ
-
i
(
α
)
u
i
,
p
=
2
,
…
,
10
.
In case p
=
2, the corresponding generating function is W
2
(
α
)
(
z
). So, u
0
=
3
/
2, u
1
=
-
2, u
2
=
1
/
2, u
i
=
0, i
=
3,4
,
…, and ϖ
2
,
ℓ
(
α
) are coefficients of the Taylor series expansion of  W
2
(
α
)
(
z
); it is easy to get
(20)
ϖ
2,0
(
α
)
=
(
3
2
)
α
.
From (19), we have
(21)
ϖ
2,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
2,0
(
α
)
u
1
=
-
4
3
α
ϖ
2,0
(
α
)
,
ϖ
2
,
ℓ
(
α
)
=
1
ℓ
u
0
[
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
u
1
+
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
u
2
]
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
,
hhhhhhhhhhhhhhhhhhhhhhhhhhhh
ℓ
=
3,4
,
…
.By the tedious calculations, one has the following result.Property 2.
The coefficientsϖ
2
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 0
<
α
<
1 have the following properties:(1)
ϖ
2,0
(
α
)
=
(
3
/
2
)
α
>
0, ϖ
2,1
(
α
)
=
-
(
4
/
3
)
α
ϖ
2,0
(
α
)
<
0;
(22)
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
≤
0
for
0
<
α
≤
5
8
,
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
>
0
for
5
8
<
α
<
1
;
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
≤
0
for
0
<
α
≤
7
8
,
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
>
0
for
7
8
<
α
<
1
;
ϖ
2
,
ℓ
(
α
)
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
<
0
,
ℓ
=
4,5
,
…
(2)
ϖ
2
,
ℓ
-
1
(
α
)
≤
ϖ
2
,
ℓ
(
α
), ℓ
≥
5;
(3)
∑
ℓ
=
0
∞
‍
ϖ
2
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
2
,
ℓ
(
α
)
=
0.Similarly, the coefficients of the 2nd scheme for1
<
α
<
2 have the following properties.Property  2′. The coefficients ϖ
2
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 1
<
α
<
2 have the following properties;(1)
ϖ
2,0
(
α
)
=
(
3
/
2
)
α
>
0, ϖ
2,1
(
α
)
=
-
(
4
/
3
)
α
ϖ
2,0
(
α
)
<
0;
(23)
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
>
0
;
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
<
0
;
ϖ
2
,
ℓ
(
α
)
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
>
0
,
ℓ
=
4,5
,
…
(2)
ϖ
2
,
ℓ
-
1
(
α
)
≥
ϖ
2
,
ℓ
(
α
), ℓ
≥
5;
(3)
∑
ℓ
=
0
∞
‍
ϖ
2
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
2
,
ℓ
(
α
)
=
0.The proofs of Properties2 and 2′ can refer to [4].


### 2.2. Determination ofϖ
p
,
ℓ
(
α
)
,
p
=
3,4
,
5,6

In this subsection, we present the recurrence relationships of the coefficientsϖ
p
,
ℓ
(
α
)
(
p
=
3,4
,
5,6
) for reference.(1) Whenp
=
3, then u
0
=
11
/
6, u
1
=
-
3, u
2
=
3
/
2, u
3
=
-
1
/
3, u
i
=
0, i
=
4,5
,
…,
(24)
ϖ
3,0
(
α
)
=
(
11
6
)
α
,
ϖ
3,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
3,0
(
α
)
u
1
=
-
18
11
α
ϖ
3,0
(
α
)
,
ϖ
3,2
(
α
)
=
1
2
u
0
[
(
α
-
2
+
1
)
ϖ
3,1
(
α
)
u
1
+
(
2
α
-
2
+
2
)
ϖ
3,0
(
α
)
u
2
]
=
3
11
[
-
3
(
α
-
1
)
ϖ
3,1
(
α
)
+
3
α
ϖ
3,0
(
α
)
]
,
ϖ
3
,
ℓ
(
α
)
=
1
ℓ
u
0
[
(
α
-
ℓ
+
1
)
ϖ
3
,
ℓ
-
1
(
α
)
u
1
+
(
2
α
-
ℓ
+
2
)
ϖ
3
,
ℓ
-
2
(
α
)
u
2
+
(
3
α
-
ℓ
+
3
)
ϖ
3
,
ℓ
-
3
(
α
)
u
3
]
=
6
11
ℓ
[
-
3
(
α
-
ℓ
+
1
)
ϖ
3
,
ℓ
-
1
(
α
)
+
3
2
(
2
α
-
ℓ
+
2
)
ϖ
3
,
ℓ
-
2
(
α
)
h
h
h
h
h
h
-
1
3
(
3
α
-
ℓ
+
3
)
ϖ
3
,
ℓ
-
3
(
α
)
]
,
ℓ
=
3,4
,
…
.(2) Whenp
=
4, then u
0
=
25
/
12, u
1
=
-
4, u
2
=
3, u
3
=
-
4
/
3, u
4
=
1
/
4, u
i
=
0, i
=
5,6
,
…,
(25)
ϖ
4,0
(
α
)
=
(
25
12
)
α
,
ϖ
4,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
4,0
(
α
)
u
1
=
-
48
25
α
ϖ
4,0
(
α
)
,
ϖ
4,2
(
α
)
=
6
25
[
-
4
(
α
-
1
)
ϖ
4,1
(
α
)
+
6
α
ϖ
4,0
(
α
)
]
,
ϖ
4,3
(
α
)
=
4
25
[
-
4
(
α
-
2
)
ϖ
4,2
(
α
)
+
3
(
2
α
-
1
)
ϖ
4,1
(
α
)
-
4
α
ϖ
4,0
(
α
)
]
,
ϖ
4
,
ℓ
(
α
)
=
12
25
ℓ
[
1
4
-
4
(
α
-
ℓ
+
1
)
ϖ
4
,
ℓ
-
1
(
α
)
+
3
(
2
α
-
ℓ
+
2
)
ϖ
4
,
ℓ
-
2
(
α
)
hhhhhh
-
4
3
(
3
α
-
ℓ
+
3
)
ϖ
4
,
ℓ
-
3
(
α
)
hhhhhh
+
1
4
(
4
α
-
ℓ
+
4
)
ϖ
4
,
ℓ
-
4
(
α
)
]
,
ℓ
=
4,5
,
…
.(3) Whenp
=
5, then u
0
=
137
/
60, u
1
=
-
5, u
2
=
5, u
3
=
-
10
/
3, u
4
=
5
/
4, u
5
=
-
1
/
5, u
i
=
0, i
=
6,7
,
…,
(26)
ϖ
5,0
(
α
)
=
(
137
60
)
α
,
ϖ
5,1
(
α
)
=
-
300
137
α
ϖ
5,0
(
α
)
,
ϖ
5,2
(
α
)
=
30
137
[
-
5
(
α
-
1
)
ϖ
5,1
(
α
)
+
10
α
ϖ
5,0
(
α
)
]
,
ϖ
5,3
(
α
)
=
20
137
[
-
5
(
α
-
2
)
ϖ
5,2
(
α
)
+
5
(
2
α
-
1
)
ϖ
5,1
(
α
)
-
10
α
ϖ
5,0
(
α
)
]
,
ϖ
5,4
(
α
)
=
15
137
[
10
3
-
5
(
α
-
3
)
ϖ
5,3
(
α
)
+
5
(
2
α
-
2
)
ϖ
5,2
(
α
)
h
h
h
h
-
10
3
(
3
α
-
1
)
ϖ
5,1
(
α
)
+
5
α
ϖ
5,0
(
α
)
]
,
ϖ
5
,
ℓ
(
α
)
=
60
137
ℓ
[
1
5
-
5
(
α
-
ℓ
+
1
)
ϖ
5
,
ℓ
-
1
(
α
)
+
5
(
2
α
-
ℓ
+
2
)
ϖ
5
,
ℓ
-
2
(
α
)
hhhhh
-
10
3
(
3
α
-
ℓ
+
3
)
ϖ
5
,
ℓ
-
3
(
α
)
hhhhh
+
5
4
(
4
α
-
ℓ
+
4
)
ϖ
5
,
ℓ
-
4
(
α
)
hhhhh
-
1
5
(
5
α
-
ℓ
+
5
)
ϖ
5
,
ℓ
-
5
(
α
)
]
,
ℓ
=
5,6
,
…
.(4) Whenp
=
6, then u
0
=
147
/
60, u
1
=
-
6, u
2
=
15
/
2, u
3
=
-
20
/
3, u
4
=
15
/
4, u
5
=
-
6
/
5, u
6
=
1
/
6, u
i
=
0, i
=
7,8
,
…,
(27)
ϖ
6,0
(
α
)
=
(
147
60
)
α
,
ϖ
6,1
(
α
)
=
-
360
147
α
ϖ
6,0
(
α
)
,
ϖ
6,2
(
α
)
=
30
147
[
-
6
(
α
-
1
)
ϖ
6,1
(
α
)
+
15
α
ϖ
6,0
(
α
)
]
,
ϖ
6,3
(
α
)
=
20
147
[
-
6
(
α
-
2
)
ϖ
6,2
(
α
)
+
15
2
(
2
α
-
1
)
ϖ
6,1
(
α
)
-
20
α
ϖ
6,0
(
α
)
]
,
ϖ
6,4
(
α
)
=
15
147
[
-
6
(
α
-
3
)
ϖ
6,3
(
α
)
+
15
2
(
2
α
-
2
)
ϖ
6,2
(
α
)
h
h
h
h
-
20
3
(
3
α
-
1
)
ϖ
6,1
(
α
)
+
15
α
ϖ
6,0
(
α
)
15
2
]
,
ϖ
6,5
(
α
)
=
12
147
[
-
6
(
α
-
4
)
ϖ
6,4
(
α
)
+
15
2
(
2
α
-
3
)
ϖ
6,3
(
α
)
h
h
h
h
-
20
3
(
3
α
-
2
)
ϖ
6,2
(
α
)
+
15
4
(
4
α
-
1
)
ϖ
6,1
(
α
)
h
h
h
h
-
6
α
ϖ
6,0
(
α
)
15
2
]
,
ϖ
6
,
ℓ
(
α
)
=
60
147
ℓ
[
-
6
(
α
-
ℓ
+
1
)
ϖ
6
,
ℓ
-
1
(
α
)
+
15
2
(
2
α
-
ℓ
+
2
)
ϖ
6
,
ℓ
-
2
(
α
)
hhhhh
-
20
3
(
3
α
-
ℓ
+
3
)
ϖ
6
,
ℓ
-
3
(
α
)
hhhhh
+
15
4
(
4
α
-
ℓ
+
4
)
ϖ
6
,
ℓ
-
4
hhhhh
-
6
5
(
5
α
-
ℓ
+
5
)
ϖ
6
,
ℓ
-
5
(
α
)
hhhhh
+
1
6
(
6
α
-
ℓ
+
6
)
ϖ
6
,
ℓ
-
6
]
,
ℓ
=
6,7
,
…
.In Figures1, 2, 3, and 4, we display the coefficients ϖ
p
,
ℓ
(
α
)
(
p
=
3,4
,
5,6
) for different α; it can be seen that ϖ
p
,
ℓ
(
α
)
→
0 when ℓ
→
∞, which coincides with the convergence [16].Figure 1
The values of the coefficientsϖ
3
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 2
The values of the coefficientsϖ
4
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 3
The values of the coefficientsϖ
5
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 4
The values of the coefficientsϖ
6
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Remark 7.
In [15], the high-order schemes for Caputo derivative were firstly derived. Here, one can get another way to construct the high-order numerical algorithms for Caputo derivatives. If the homogeneous initial value conditions are satisfied, one has the following numerical schemes due to Lemma 5:
(28)
D
C
a
,
x
α
f
(
x
m
)
=
1
h
α
∑
ℓ
=
0
m
‍
ϖ
p
,
ℓ
(
α
)
f
(
x
m
-
ℓ
)
+
𝒪
(
h
p
)
,
h
h
h
h
h
h
h
h
h
h
h
h
p
=
2,3
,
4,5
,
6
.


### 2.3. Determination ofϖ
p
,
ℓ
(
p
=
7,8
,
9,10
)

In this subsection, we present the recursion formulas ofϖ
p
,
ℓ
(
p
=
7,8
,
9,10
) for reference.(1) Whenp
=
7, then u
0
=
363
/
140, u
1
=
-
7, u
2
=
21
/
2, u
3
=
-
35
/
3, u
4
=
35
/
4, u
5
=
-
21
/
5, u
6
=
7
/
6, u
7
=
-
1
/
7, and u
i
=
0, i
=
8,9
,
…. The coefficients are given as follows:
(29)
ϖ
7,0
(
α
)
=
(
363
140
)
α
,
ϖ
7,1
(
α
)
=
-
980
363
α
ϖ
7,0
(
α
)
,
ϖ
7,2
(
α
)
=
70
363
[
-
7
(
α
-
1
)
ϖ
7,1
(
α
)
+
21
α
ϖ
7,0
(
α
)
]
,
ϖ
7,3
(
α
)
=
140
1089
[
-
7
(
α
-
2
)
ϖ
7,2
(
α
)
+
21
2
(
2
α
-
1
)
ϖ
7,1
(
α
)
-
35
α
ϖ
7,0
(
α
)
]
,
ϖ
7,4
(
α
)
=
35
363
[
35
3
-
7
(
α
-
3
)
ϖ
7,3
(
α
)
+
21
(
α
-
1
)
ϖ
7,2
(
α
)
h
h
h
h
-
35
3
(
3
α
-
1
)
ϖ
7,1
(
α
)
+
35
α
ϖ
7,0
(
α
)
]
,
ϖ
7,5
(
α
)
=
28
363
[
-
7
(
α
-
4
)
ϖ
7,4
(
α
)
+
21
2
(
2
α
-
3
)
ϖ
7,3
(
α
)
h
h
h
h
-
35
3
(
3
α
-
2
)
ϖ
7,2
(
α
)
+
35
4
(
4
α
-
1
)
ϖ
7,1
(
α
)
h
h
h
h
-
21
α
ϖ
7,0
(
α
)
21
2
]
,
ϖ
7,6
(
α
)
=
70
1089
[
21
5
-
7
(
α
-
5
)
ϖ
7,5
(
α
)
+
21
(
α
-
2
)
ϖ
7,4
(
α
)
h
h
h
h
h
-
35
(
α
-
1
)
ϖ
7,3
(
α
)
+
35
2
(
2
α
-
1
)
ϖ
7,2
(
α
)
h
h
h
h
h
-
21
5
(
5
α
-
1
)
ϖ
7,1
(
α
)
+
7
α
ϖ
7,0
(
α
)
]
,
ϖ
7
,
ℓ
=
140
363
ℓ
[
-
7
(
α
-
ℓ
+
1
)
ϖ
7
,
ℓ
-
1
(
α
)
+
21
2
(
2
α
-
ℓ
+
2
)
ϖ
7
,
ℓ
-
2
(
α
)
h
h
h
h
h
-
35
3
(
3
α
-
ℓ
+
3
)
ϖ
7
,
ℓ
-
3
(
α
)
h
h
h
h
h
+
35
4
(
4
α
-
ℓ
+
4
)
ϖ
7
,
ℓ
-
4
(
α
)
h
h
h
h
h
-
21
5
(
5
α
-
ℓ
+
5
)
ϖ
7
,
ℓ
-
5
(
α
)
h
h
h
h
h
+
7
6
(
6
α
-
ℓ
+
6
)
ϖ
7
,
ℓ
-
6
(
α
)
h
h
h
h
h
-
1
7
(
7
α
-
ℓ
+
7
)
ϖ
7
,
ℓ
-
7
(
α
)
]
,
ℓ
=
7,8
,
…
.(2) Whenp
=
8, then u
0
=
761
/
280, u
1
=
-
8, u
2
=
14, u
3
=
-
56
/
3, u
4
=
35
/
2, u
5
=
-
56
/
5, u
6
=
14
/
3, u
7
=
-
8
/
7, u
8
=
1
/
8, and u
i
=
0, i
=
9,10
,
…. The coefficients are given as follows:
(30)
ϖ
8,0
(
α
)
=
(
761
280
)
α
,
ϖ
8,1
(
α
)
=
-
2240
761
α
ϖ
8,0
(
α
)
,
ϖ
8,2
(
α
)
=
140
761
[
-
8
(
α
-
1
)
ϖ
8,1
(
α
)
+
28
α
ϖ
8,0
(
α
)
]
,
ϖ
8,3
(
α
)
=
280
2283
[
-
8
(
α
-
2
)
ϖ
8,2
(
α
)
h
h
h
h
h
+
14
(
2
α
-
1
)
ϖ
8,1
(
α
)
-
56
α
ϖ
8,0
(
α
)
]
,
ϖ
8,4
(
α
)
=
70
761
[
56
3
-
8
(
α
-
3
)
ϖ
8,3
(
α
)
+
28
(
α
-
1
)
ϖ
8,2
(
α
)
h
h
h
h
-
56
3
(
3
α
-
1
)
ϖ
8,1
(
α
)
+
70
α
ϖ
8,0
(
α
)
]
,
ϖ
8,5
(
α
)
=
56
761
[
56
3
-
8
(
α
-
4
)
ϖ
8,4
(
α
)
+
14
(
2
α
-
3
)
ϖ
8,3
(
α
)
h
h
h
h
-
56
3
(
3
α
-
2
)
ϖ
8,2
(
α
)
+
35
2
(
4
α
-
1
)
ϖ
8,1
(
α
)
h
h
h
h
-
56
α
ϖ
8,0
(
α
)
56
3
]
,
ϖ
8,6
(
α
)
=
140
2283
[
56
5
-
8
(
α
-
5
)
ϖ
8,5
(
α
)
+
28
(
α
-
2
)
ϖ
8,4
(
α
)
h
h
h
h
h
-
56
(
α
-
1
)
ϖ
8,3
(
α
)
+
35
(
2
α
-
1
)
ϖ
8,2
(
α
)
h
h
h
h
h
-
56
5
(
5
α
-
1
)
ϖ
8,1
(
α
)
+
28
α
ϖ
8,0
(
α
)
]
,
ϖ
8,7
(
α
)
=
40
761
[
56
3
-
8
(
α
-
6
)
ϖ
8,6
(
α
)
+
14
(
2
α
-
5
)
ϖ
8,5
(
α
)
h
h
h
h
-
56
3
(
3
α
-
4
)
ϖ
8,4
(
α
)
+
35
2
(
4
α
-
3
)
ϖ
8,3
(
α
)
h
h
h
h
-
56
5
(
5
α
-
2
)
ϖ
8,2
(
α
)
+
14
3
(
6
α
-
1
)
ϖ
8,1
(
α
)
h
h
h
h
-
8
α
ϖ
8,0
(
α
)
56
3
]
,
ϖ
8
,
ℓ
(
α
)
=
280
761
ℓ
[
1
8
-
8
(
α
-
ℓ
+
1
)
ϖ
8
,
ℓ
-
1
(
α
)
h
h
h
h
h
+
14
(
2
α
-
ℓ
+
2
)
ϖ
8
,
ℓ
-
2
(
α
)
h
h
h
h
h
-
56
3
(
3
α
-
ℓ
+
3
)
ϖ
8
,
ℓ
-
3
(
α
)
h
h
h
h
h
+
35
2
(
4
α
-
ℓ
+
4
)
ϖ
8
,
ℓ
-
4
(
α
)
h
h
h
h
h
-
56
5
(
5
α
-
ℓ
+
5
)
ϖ
8
,
ℓ
-
5
(
α
)
h
h
h
h
h
+
14
3
(
6
α
-
ℓ
+
6
)
ϖ
8
,
ℓ
-
6
(
α
)
h
h
h
h
h
-
8
7
(
7
α
-
ℓ
+
7
)
ϖ
8
,
ℓ
-
7
(
α
)
h
h
h
h
h
+
1
8
(
8
α
-
ℓ
+
8
)
ϖ
8
,
ℓ
-
8
(
α
)
]
,
ℓ
=
8,9
,
…
.(3) Whenp
=
9, then u
0
=
7129
/
2520, u
1
=
-
9, u
2
=
18, u
3
=
-
28, u
4
=
63
/
2, u
5
=
-
126
/
5, u
6
=
14, u
7
=
-
36
/
7, u
8
=
9
/
8, and u
9
=
-
1
/
9, u
i
=
0, i
=
10,11
,
…. The coefficients are displayed as follows:
(31)
ϖ
9,0
(
α
)
=
(
7129
2520
)
α
,
ϖ
9,1
(
α
)
=
-
22680
7129
α
ϖ
9,0
(
α
)
,
ϖ
9,2
(
α
)
=
1260
7129
[
-
9
(
α
-
1
)
ϖ
9,1
(
α
)
+
36
α
ϖ
9,0
(
α
)
]
,
ϖ
9,3
(
α
)
=
840
7129
[
-
9
(
α
-
2
)
ϖ
9,2
(
α
)
+
18
(
2
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
84
α
ϖ
9,0
(
α
)
]
,
ϖ
9,4
(
α
)
=
630
7129
[
-
9
(
α
-
3
)
ϖ
9,3
(
α
)
+
36
(
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
28
(
3
α
-
1
)
ϖ
9,1
(
α
)
+
126
α
ϖ
9,0
(
α
)
]
,
ϖ
9,5
(
α
)
=
504
7129
[
63
2
-
9
(
α
-
4
)
ϖ
9,4
(
α
)
+
18
(
2
α
-
3
)
ϖ
9,3
(
α
)
h
h
h
h
h
-
28
(
3
α
-
2
)
ϖ
9,2
(
α
)
+
63
2
(
4
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
126
α
ϖ
9,0
(
α
)
63
2
]
,
ϖ
9,6
(
α
)
=
420
7129
[
126
5
-
9
(
α
-
5
)
ϖ
9,5
(
α
)
+
36
(
α
-
2
)
ϖ
9,4
(
α
)
h
h
h
h
h
-
84
(
α
-
1
)
ϖ
9,3
(
α
)
+
63
(
2
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
1
)
ϖ
9,1
(
α
)
+
84
α
ϖ
9,0
(
α
)
]
,
ϖ
9,7
(
α
)
=
360
7129
[
63
2
-
9
(
α
-
6
)
ϖ
9,6
(
α
)
+
18
(
2
α
-
5
)
ϖ
9,5
(
α
)
h
h
h
h
h
-
28
(
3
α
-
4
)
ϖ
9,4
(
α
)
+
63
2
(
4
α
-
3
)
ϖ
9,3
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
2
)
ϖ
9,2
(
α
)
+
14
(
6
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
36
α
ϖ
9,0
(
α
)
63
2
]
,
ϖ
9,8
(
α
)
=
315
7129
[
36
7
-
9
(
α
-
7
)
ϖ
9,7
(
α
)
+
36
(
α
-
3
)
ϖ
9,6
(
α
)
h
h
h
h
h
-
28
(
3
α
-
5
)
ϖ
9,5
(
α
)
+
126
(
α
-
1
)
ϖ
9,4
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
3
)
ϖ
9,3
(
α
)
+
28
(
3
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
36
7
(
7
α
-
1
)
ϖ
9,1
(
α
)
+
72
α
ϖ
9,0
(
α
)
]
,
ϖ
9
,
ℓ
(
α
)
=
2520
7129
ℓ
[
-
9
(
α
-
ℓ
+
1
)
ϖ
9
,
ℓ
-
1
(
α
)
+
18
(
2
α
-
ℓ
+
2
)
ϖ
9
,
ℓ
-
2
(
α
)
1
9
h
h
h
h
h
h
-
28
(
3
α
-
ℓ
+
3
)
ϖ
9
,
ℓ
-
3
(
α
)
h
h
h
h
h
h
+
63
2
(
4
α
-
ℓ
+
4
)
ϖ
9
,
ℓ
-
4
(
α
)
h
h
h
h
h
h
-
126
5
(
5
α
-
ℓ
+
5
)
ϖ
9
,
ℓ
-
5
(
α
)
h
h
h
h
h
h
+
14
(
6
α
-
ℓ
+
6
)
ϖ
9
,
ℓ
-
6
(
α
)
h
h
h
h
h
h
-
36
7
(
7
α
-
ℓ
+
7
)
ϖ
9
,
ℓ
-
7
(
α
)
h
h
h
h
h
h
+
9
8
(
8
α
-
ℓ
+
8
)
ϖ
9
,
ℓ
-
8
(
α
)
h
h
h
h
h
h
-
1
9
(
9
α
-
ℓ
+
9
)
ϖ
9
,
ℓ
-
9
(
α
)
]
,
ℓ
=
9,10
,
…
.(4) Whenp
=
10, then u
0
=
7381
/
2520, u
1
=
-
10, u
2
=
45
/
2, u
3
=
-
40, u
4
=
105
/
2, u
5
=
-
252
/
5, u
6
=
35, u
7
=
-
120
/
7, u
8
=
45
/
8, u
9
=
-
10
/
9, u
10
=
1
/
10, and u
i
=
0, i
=
11,12
,
…. The coefficients are given as follows:
(32)
ϖ
10,0
(
α
)
=
(
7381
2520
)
α
,
ϖ
10,1
(
α
)
=
-
25200
7381
α
ϖ
10,0
(
α
)
,
ϖ
10,2
(
α
)
=
1260
7381
[
-
10
(
α
-
1
)
ϖ
10,1
(
α
)
+
45
α
ϖ
10,0
(
α
)
]
,
ϖ
10,3
(
α
)
=
840
7381
[
-
10
(
α
-
2
)
ϖ
10,2
(
α
)
+
45
2
(
2
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
120
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,4
(
α
)
=
630
7381
[
-
10
(
α
-
3
)
ϖ
10,3
(
α
)
+
45
(
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
40
(
3
α
-
1
)
ϖ
10,1
(
α
)
+
210
α
ϖ
10,0
(
α
)
]
,
ϖ
10,5
(
α
)
=
504
7381
[
-
10
(
α
-
4
)
ϖ
10,4
(
α
)
+
45
2
(
2
α
-
3
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
40
(
3
α
-
2
)
ϖ
10,2
(
α
)
+
105
2
(
4
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
252
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,6
(
α
)
=
420
7381
[
252
5
-
10
(
α
-
5
)
ϖ
10,5
(
α
)
+
45
(
α
-
2
)
ϖ
10,4
(
α
)
h
h
h
h
h
-
120
(
α
-
1
)
ϖ
10,3
(
α
)
+
105
(
2
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
1
)
ϖ
10,1
(
α
)
+
210
α
ϖ
10,0
(
α
)
]
,
ϖ
10,7
(
α
)
=
360
7381
[
-
10
(
α
-
6
)
ϖ
10,6
(
α
)
+
45
2
(
2
α
-
5
)
ϖ
10,5
(
α
)
h
h
h
h
h
-
40
(
3
α
-
4
)
ϖ
10,4
(
α
)
+
105
2
(
4
α
-
3
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
2
)
ϖ
10,2
(
α
)
+
35
(
6
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
120
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,8
(
α
)
=
315
7381
[
120
7
-
10
(
α
-
7
)
ϖ
10,7
(
α
)
+
45
(
α
-
3
)
ϖ
10,6
(
α
)
h
h
h
h
h
-
40
(
3
α
-
5
)
ϖ
10,5
(
α
)
+
210
(
α
-
1
)
ϖ
10,4
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
3
)
ϖ
10,3
(
α
)
+
70
(
3
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
120
7
(
7
α
-
1
)
ϖ
10,1
(
α
)
+
45
α
ϖ
10,0
(
α
)
]
,
ϖ
10,9
(
α
)
=
280
7381
[
-
10
(
α
-
8
)
ϖ
10,8
(
α
)
+
45
2
(
2
α
-
7
)
ϖ
10,7
(
α
)
h
h
h
h
h
-
120
(
α
-
2
)
ϖ
10,6
(
α
)
+
105
2
(
4
α
-
5
)
ϖ
10,5
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
4
)
ϖ
10,4
(
α
)
+
105
(
2
α
-
1
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
120
7
(
7
α
-
2
)
ϖ
10,2
(
α
)
+
45
8
(
8
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
10
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10
,
ℓ
(
α
)
=
2520
7381
ℓ
[
1
10
-
10
(
α
-
ℓ
+
1
)
ϖ
10
,
ℓ
-
1
(
α
)
h
h
h
h
h
h
+
45
2
(
2
α
-
ℓ
+
2
)
ϖ
10
,
ℓ
-
2
(
α
)
h
h
h
h
h
h
-
40
(
3
α
-
ℓ
+
3
)
ϖ
10
,
ℓ
-
3
(
α
)
h
h
h
h
h
h
+
105
2
(
4
α
-
ℓ
+
4
)
ϖ
10
,
ℓ
-
4
(
α
)
h
h
h
h
h
h
-
252
5
(
5
α
-
ℓ
+
5
)
ϖ
10
,
ℓ
-
5
(
α
)
h
h
h
h
h
h
+
35
(
6
α
-
ℓ
+
6
)
ϖ
10
,
ℓ
-
6
(
α
)
h
h
h
h
h
h
-
120
7
(
7
α
-
ℓ
+
7
)
ϖ
10
,
ℓ
-
7
(
α
)
h
h
h
h
h
h
+
45
8
(
8
α
-
ℓ
+
8
)
ϖ
10
,
ℓ
-
8
(
α
)
h
h
h
h
h
h
-
10
9
(
9
α
-
ℓ
+
9
)
ϖ
10
,
ℓ
-
9
(
α
)
h
h
h
h
h
h
+
1
10
(
10
α
-
ℓ
+
10
)
ϖ
10,10
-
ℓ
(
α
)
]
,
f
f
f
f
f
f
f
f
f
f
f
f
f
f
ℓ
=
10,11
,
…
.Next, we plot the figures of coefficients ofϖ
p
,
ℓ
(
α
)
(
p
=
7,8
,
9,10
) to show the evolutions with ℓ. From Figures 5, 6, 7, and 8, we can see that these coefficients are violently oscillatory such that the approximation behaves like Runge’s phenomenon, which is similar to the case of ordinary differential equation. So, to seek high-order (≥7th-order) schemes by this form seems not to be appropriate.Figure 5
The values of the coefficientsϖ
7
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 6
The values of the coefficientsϖ
8
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 7
The values of the coefficientsϖ
9
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 8
The values of the coefficientsϖ
10
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.


## 2.1. Properties ofϖ
1
,
ℓ
(
α
) and ϖ
2
,
ℓ
(
α
)

As far as we know, the coefficients of the 1st-order and 2nd-order schemes are often used. Here, we present their properties for reference. We know thatϖ
1
,
ℓ
(
α
) are the binomial coefficients in formula (12). They have the following properties.Property 1.
The coefficientsϖ
1
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 0
<
α
<
1 have the following properties:(1)
ϖ
1,0
(
α
)
=
1, ϖ
1
,
ℓ
(
α
)
=
(
-
1
)
ℓ
Γ
(
α
+
1
)
/
(
Γ
(
ℓ
+
1
)
Γ
(
α
-
ℓ
+
1
)
)
<
0, ℓ
≥
1;
(2)
ϖ
1
,
ℓ
-
1
(
α
)
≤
ϖ
1
,
ℓ
(
α
), ℓ
≥
2;
(3)
∑
ℓ
=
0
∞
‍
ϖ
1
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
1
,
ℓ
(
α
)
=
0.Parallelly, the properties of the coefficients of the 1st-order schemes for1
<
α
<
2 are given as follows.Property  1′. The coefficients ϖ
1
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 1
<
α
<
2 have the following properties:(1)
ϖ
1,0
(
α
)
=
1, ϖ
1,1
=
-
α
<
0, ϖ
1
,
ℓ
(
α
)
=
(
-
1
)
ℓ
Γ
(
α
+
1
)
/
(
Γ
(
ℓ
+
1
)
Γ
(
α
-
ℓ
+
1
)
)
>
0, ℓ
≥
2;
(2)
ϖ
1
,
ℓ
-
1
(
α
)
≥
ϖ
1
,
ℓ
(
α
), ℓ
≥
3;
(3)
∑
ℓ
=
0
∞
‍
ϖ
1
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
1
,
ℓ
(
α
)
=
0.The above Properties 1 and 1′ are known to us. Next, we present the coefficientsϖ
2
,
ℓ
(
α
) which are somewhat complex. Let i
=
ℓ
-
j; formula (18) can be written as
(19)
ϖ
p
,
ℓ
(
α
)
=
1
ℓ
u
0
∑
i
=
1
ℓ
‍
(
α
i
-
ℓ
+
i
)
ϖ
p
,
ℓ
-
i
(
α
)
u
i
,
p
=
2
,
…
,
10
.
In case p
=
2, the corresponding generating function is W
2
(
α
)
(
z
). So, u
0
=
3
/
2, u
1
=
-
2, u
2
=
1
/
2, u
i
=
0, i
=
3,4
,
…, and ϖ
2
,
ℓ
(
α
) are coefficients of the Taylor series expansion of  W
2
(
α
)
(
z
); it is easy to get
(20)
ϖ
2,0
(
α
)
=
(
3
2
)
α
.
From (19), we have
(21)
ϖ
2,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
2,0
(
α
)
u
1
=
-
4
3
α
ϖ
2,0
(
α
)
,
ϖ
2
,
ℓ
(
α
)
=
1
ℓ
u
0
[
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
u
1
+
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
u
2
]
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
,
hhhhhhhhhhhhhhhhhhhhhhhhhhhh
ℓ
=
3,4
,
…
.By the tedious calculations, one has the following result.Property 2.
The coefficientsϖ
2
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 0
<
α
<
1 have the following properties:(1)
ϖ
2,0
(
α
)
=
(
3
/
2
)
α
>
0, ϖ
2,1
(
α
)
=
-
(
4
/
3
)
α
ϖ
2,0
(
α
)
<
0;
(22)
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
≤
0
for
0
<
α
≤
5
8
,
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
>
0
for
5
8
<
α
<
1
;
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
≤
0
for
0
<
α
≤
7
8
,
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
>
0
for
7
8
<
α
<
1
;
ϖ
2
,
ℓ
(
α
)
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
<
0
,
ℓ
=
4,5
,
…
(2)
ϖ
2
,
ℓ
-
1
(
α
)
≤
ϖ
2
,
ℓ
(
α
), ℓ
≥
5;
(3)
∑
ℓ
=
0
∞
‍
ϖ
2
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
2
,
ℓ
(
α
)
=
0.Similarly, the coefficients of the 2nd scheme for1
<
α
<
2 have the following properties.Property  2′. The coefficients ϖ
2
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for 1
<
α
<
2 have the following properties;(1)
ϖ
2,0
(
α
)
=
(
3
/
2
)
α
>
0, ϖ
2,1
(
α
)
=
-
(
4
/
3
)
α
ϖ
2,0
(
α
)
<
0;
(23)
ϖ
2,2
(
α
)
=
α
(
8
α
-
5
)
9
(
3
2
)
α
>
0
;
ϖ
2,3
(
α
)
=
4
α
(
1
-
α
)
(
8
α
-
7
)
81
(
3
2
)
α
<
0
;
ϖ
2
,
ℓ
(
α
)
=
2
3
ℓ
[
-
2
(
α
-
ℓ
+
1
)
ϖ
2
,
ℓ
-
1
(
α
)
+
1
2
(
2
α
-
ℓ
+
2
)
ϖ
2
,
ℓ
-
2
(
α
)
]
>
0
,
ℓ
=
4,5
,
…
(2)
ϖ
2
,
ℓ
-
1
(
α
)
≥
ϖ
2
,
ℓ
(
α
), ℓ
≥
5;
(3)
∑
ℓ
=
0
∞
‍
ϖ
2
,
ℓ
(
α
)
=
0, lim
⁡
ℓ
→
∞
⁡
ϖ
2
,
ℓ
(
α
)
=
0.The proofs of Properties2 and 2′ can refer to [4].


## 2.2. Determination ofϖ
p
,
ℓ
(
α
)
,
p
=
3,4
,
5,6

In this subsection, we present the recurrence relationships of the coefficientsϖ
p
,
ℓ
(
α
)
(
p
=
3,4
,
5,6
) for reference.(1) Whenp
=
3, then u
0
=
11
/
6, u
1
=
-
3, u
2
=
3
/
2, u
3
=
-
1
/
3, u
i
=
0, i
=
4,5
,
…,
(24)
ϖ
3,0
(
α
)
=
(
11
6
)
α
,
ϖ
3,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
3,0
(
α
)
u
1
=
-
18
11
α
ϖ
3,0
(
α
)
,
ϖ
3,2
(
α
)
=
1
2
u
0
[
(
α
-
2
+
1
)
ϖ
3,1
(
α
)
u
1
+
(
2
α
-
2
+
2
)
ϖ
3,0
(
α
)
u
2
]
=
3
11
[
-
3
(
α
-
1
)
ϖ
3,1
(
α
)
+
3
α
ϖ
3,0
(
α
)
]
,
ϖ
3
,
ℓ
(
α
)
=
1
ℓ
u
0
[
(
α
-
ℓ
+
1
)
ϖ
3
,
ℓ
-
1
(
α
)
u
1
+
(
2
α
-
ℓ
+
2
)
ϖ
3
,
ℓ
-
2
(
α
)
u
2
+
(
3
α
-
ℓ
+
3
)
ϖ
3
,
ℓ
-
3
(
α
)
u
3
]
=
6
11
ℓ
[
-
3
(
α
-
ℓ
+
1
)
ϖ
3
,
ℓ
-
1
(
α
)
+
3
2
(
2
α
-
ℓ
+
2
)
ϖ
3
,
ℓ
-
2
(
α
)
h
h
h
h
h
h
-
1
3
(
3
α
-
ℓ
+
3
)
ϖ
3
,
ℓ
-
3
(
α
)
]
,
ℓ
=
3,4
,
…
.(2) Whenp
=
4, then u
0
=
25
/
12, u
1
=
-
4, u
2
=
3, u
3
=
-
4
/
3, u
4
=
1
/
4, u
i
=
0, i
=
5,6
,
…,
(25)
ϖ
4,0
(
α
)
=
(
25
12
)
α
,
ϖ
4,1
(
α
)
=
1
u
0
(
α
-
1
+
1
)
ϖ
4,0
(
α
)
u
1
=
-
48
25
α
ϖ
4,0
(
α
)
,
ϖ
4,2
(
α
)
=
6
25
[
-
4
(
α
-
1
)
ϖ
4,1
(
α
)
+
6
α
ϖ
4,0
(
α
)
]
,
ϖ
4,3
(
α
)
=
4
25
[
-
4
(
α
-
2
)
ϖ
4,2
(
α
)
+
3
(
2
α
-
1
)
ϖ
4,1
(
α
)
-
4
α
ϖ
4,0
(
α
)
]
,
ϖ
4
,
ℓ
(
α
)
=
12
25
ℓ
[
1
4
-
4
(
α
-
ℓ
+
1
)
ϖ
4
,
ℓ
-
1
(
α
)
+
3
(
2
α
-
ℓ
+
2
)
ϖ
4
,
ℓ
-
2
(
α
)
hhhhhh
-
4
3
(
3
α
-
ℓ
+
3
)
ϖ
4
,
ℓ
-
3
(
α
)
hhhhhh
+
1
4
(
4
α
-
ℓ
+
4
)
ϖ
4
,
ℓ
-
4
(
α
)
]
,
ℓ
=
4,5
,
…
.(3) Whenp
=
5, then u
0
=
137
/
60, u
1
=
-
5, u
2
=
5, u
3
=
-
10
/
3, u
4
=
5
/
4, u
5
=
-
1
/
5, u
i
=
0, i
=
6,7
,
…,
(26)
ϖ
5,0
(
α
)
=
(
137
60
)
α
,
ϖ
5,1
(
α
)
=
-
300
137
α
ϖ
5,0
(
α
)
,
ϖ
5,2
(
α
)
=
30
137
[
-
5
(
α
-
1
)
ϖ
5,1
(
α
)
+
10
α
ϖ
5,0
(
α
)
]
,
ϖ
5,3
(
α
)
=
20
137
[
-
5
(
α
-
2
)
ϖ
5,2
(
α
)
+
5
(
2
α
-
1
)
ϖ
5,1
(
α
)
-
10
α
ϖ
5,0
(
α
)
]
,
ϖ
5,4
(
α
)
=
15
137
[
10
3
-
5
(
α
-
3
)
ϖ
5,3
(
α
)
+
5
(
2
α
-
2
)
ϖ
5,2
(
α
)
h
h
h
h
-
10
3
(
3
α
-
1
)
ϖ
5,1
(
α
)
+
5
α
ϖ
5,0
(
α
)
]
,
ϖ
5
,
ℓ
(
α
)
=
60
137
ℓ
[
1
5
-
5
(
α
-
ℓ
+
1
)
ϖ
5
,
ℓ
-
1
(
α
)
+
5
(
2
α
-
ℓ
+
2
)
ϖ
5
,
ℓ
-
2
(
α
)
hhhhh
-
10
3
(
3
α
-
ℓ
+
3
)
ϖ
5
,
ℓ
-
3
(
α
)
hhhhh
+
5
4
(
4
α
-
ℓ
+
4
)
ϖ
5
,
ℓ
-
4
(
α
)
hhhhh
-
1
5
(
5
α
-
ℓ
+
5
)
ϖ
5
,
ℓ
-
5
(
α
)
]
,
ℓ
=
5,6
,
…
.(4) Whenp
=
6, then u
0
=
147
/
60, u
1
=
-
6, u
2
=
15
/
2, u
3
=
-
20
/
3, u
4
=
15
/
4, u
5
=
-
6
/
5, u
6
=
1
/
6, u
i
=
0, i
=
7,8
,
…,
(27)
ϖ
6,0
(
α
)
=
(
147
60
)
α
,
ϖ
6,1
(
α
)
=
-
360
147
α
ϖ
6,0
(
α
)
,
ϖ
6,2
(
α
)
=
30
147
[
-
6
(
α
-
1
)
ϖ
6,1
(
α
)
+
15
α
ϖ
6,0
(
α
)
]
,
ϖ
6,3
(
α
)
=
20
147
[
-
6
(
α
-
2
)
ϖ
6,2
(
α
)
+
15
2
(
2
α
-
1
)
ϖ
6,1
(
α
)
-
20
α
ϖ
6,0
(
α
)
]
,
ϖ
6,4
(
α
)
=
15
147
[
-
6
(
α
-
3
)
ϖ
6,3
(
α
)
+
15
2
(
2
α
-
2
)
ϖ
6,2
(
α
)
h
h
h
h
-
20
3
(
3
α
-
1
)
ϖ
6,1
(
α
)
+
15
α
ϖ
6,0
(
α
)
15
2
]
,
ϖ
6,5
(
α
)
=
12
147
[
-
6
(
α
-
4
)
ϖ
6,4
(
α
)
+
15
2
(
2
α
-
3
)
ϖ
6,3
(
α
)
h
h
h
h
-
20
3
(
3
α
-
2
)
ϖ
6,2
(
α
)
+
15
4
(
4
α
-
1
)
ϖ
6,1
(
α
)
h
h
h
h
-
6
α
ϖ
6,0
(
α
)
15
2
]
,
ϖ
6
,
ℓ
(
α
)
=
60
147
ℓ
[
-
6
(
α
-
ℓ
+
1
)
ϖ
6
,
ℓ
-
1
(
α
)
+
15
2
(
2
α
-
ℓ
+
2
)
ϖ
6
,
ℓ
-
2
(
α
)
hhhhh
-
20
3
(
3
α
-
ℓ
+
3
)
ϖ
6
,
ℓ
-
3
(
α
)
hhhhh
+
15
4
(
4
α
-
ℓ
+
4
)
ϖ
6
,
ℓ
-
4
hhhhh
-
6
5
(
5
α
-
ℓ
+
5
)
ϖ
6
,
ℓ
-
5
(
α
)
hhhhh
+
1
6
(
6
α
-
ℓ
+
6
)
ϖ
6
,
ℓ
-
6
]
,
ℓ
=
6,7
,
…
.In Figures1, 2, 3, and 4, we display the coefficients ϖ
p
,
ℓ
(
α
)
(
p
=
3,4
,
5,6
) for different α; it can be seen that ϖ
p
,
ℓ
(
α
)
→
0 when ℓ
→
∞, which coincides with the convergence [16].Figure 1
The values of the coefficientsϖ
3
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 2
The values of the coefficientsϖ
4
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 3
The values of the coefficientsϖ
5
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 4
The values of the coefficientsϖ
6
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Remark 7.
In [15], the high-order schemes for Caputo derivative were firstly derived. Here, one can get another way to construct the high-order numerical algorithms for Caputo derivatives. If the homogeneous initial value conditions are satisfied, one has the following numerical schemes due to Lemma 5:
(28)
D
C
a
,
x
α
f
(
x
m
)
=
1
h
α
∑
ℓ
=
0
m
‍
ϖ
p
,
ℓ
(
α
)
f
(
x
m
-
ℓ
)
+
𝒪
(
h
p
)
,
h
h
h
h
h
h
h
h
h
h
h
h
p
=
2,3
,
4,5
,
6
.


## 2.3. Determination ofϖ
p
,
ℓ
(
p
=
7,8
,
9,10
)

In this subsection, we present the recursion formulas ofϖ
p
,
ℓ
(
p
=
7,8
,
9,10
) for reference.(1) Whenp
=
7, then u
0
=
363
/
140, u
1
=
-
7, u
2
=
21
/
2, u
3
=
-
35
/
3, u
4
=
35
/
4, u
5
=
-
21
/
5, u
6
=
7
/
6, u
7
=
-
1
/
7, and u
i
=
0, i
=
8,9
,
…. The coefficients are given as follows:
(29)
ϖ
7,0
(
α
)
=
(
363
140
)
α
,
ϖ
7,1
(
α
)
=
-
980
363
α
ϖ
7,0
(
α
)
,
ϖ
7,2
(
α
)
=
70
363
[
-
7
(
α
-
1
)
ϖ
7,1
(
α
)
+
21
α
ϖ
7,0
(
α
)
]
,
ϖ
7,3
(
α
)
=
140
1089
[
-
7
(
α
-
2
)
ϖ
7,2
(
α
)
+
21
2
(
2
α
-
1
)
ϖ
7,1
(
α
)
-
35
α
ϖ
7,0
(
α
)
]
,
ϖ
7,4
(
α
)
=
35
363
[
35
3
-
7
(
α
-
3
)
ϖ
7,3
(
α
)
+
21
(
α
-
1
)
ϖ
7,2
(
α
)
h
h
h
h
-
35
3
(
3
α
-
1
)
ϖ
7,1
(
α
)
+
35
α
ϖ
7,0
(
α
)
]
,
ϖ
7,5
(
α
)
=
28
363
[
-
7
(
α
-
4
)
ϖ
7,4
(
α
)
+
21
2
(
2
α
-
3
)
ϖ
7,3
(
α
)
h
h
h
h
-
35
3
(
3
α
-
2
)
ϖ
7,2
(
α
)
+
35
4
(
4
α
-
1
)
ϖ
7,1
(
α
)
h
h
h
h
-
21
α
ϖ
7,0
(
α
)
21
2
]
,
ϖ
7,6
(
α
)
=
70
1089
[
21
5
-
7
(
α
-
5
)
ϖ
7,5
(
α
)
+
21
(
α
-
2
)
ϖ
7,4
(
α
)
h
h
h
h
h
-
35
(
α
-
1
)
ϖ
7,3
(
α
)
+
35
2
(
2
α
-
1
)
ϖ
7,2
(
α
)
h
h
h
h
h
-
21
5
(
5
α
-
1
)
ϖ
7,1
(
α
)
+
7
α
ϖ
7,0
(
α
)
]
,
ϖ
7
,
ℓ
=
140
363
ℓ
[
-
7
(
α
-
ℓ
+
1
)
ϖ
7
,
ℓ
-
1
(
α
)
+
21
2
(
2
α
-
ℓ
+
2
)
ϖ
7
,
ℓ
-
2
(
α
)
h
h
h
h
h
-
35
3
(
3
α
-
ℓ
+
3
)
ϖ
7
,
ℓ
-
3
(
α
)
h
h
h
h
h
+
35
4
(
4
α
-
ℓ
+
4
)
ϖ
7
,
ℓ
-
4
(
α
)
h
h
h
h
h
-
21
5
(
5
α
-
ℓ
+
5
)
ϖ
7
,
ℓ
-
5
(
α
)
h
h
h
h
h
+
7
6
(
6
α
-
ℓ
+
6
)
ϖ
7
,
ℓ
-
6
(
α
)
h
h
h
h
h
-
1
7
(
7
α
-
ℓ
+
7
)
ϖ
7
,
ℓ
-
7
(
α
)
]
,
ℓ
=
7,8
,
…
.(2) Whenp
=
8, then u
0
=
761
/
280, u
1
=
-
8, u
2
=
14, u
3
=
-
56
/
3, u
4
=
35
/
2, u
5
=
-
56
/
5, u
6
=
14
/
3, u
7
=
-
8
/
7, u
8
=
1
/
8, and u
i
=
0, i
=
9,10
,
…. The coefficients are given as follows:
(30)
ϖ
8,0
(
α
)
=
(
761
280
)
α
,
ϖ
8,1
(
α
)
=
-
2240
761
α
ϖ
8,0
(
α
)
,
ϖ
8,2
(
α
)
=
140
761
[
-
8
(
α
-
1
)
ϖ
8,1
(
α
)
+
28
α
ϖ
8,0
(
α
)
]
,
ϖ
8,3
(
α
)
=
280
2283
[
-
8
(
α
-
2
)
ϖ
8,2
(
α
)
h
h
h
h
h
+
14
(
2
α
-
1
)
ϖ
8,1
(
α
)
-
56
α
ϖ
8,0
(
α
)
]
,
ϖ
8,4
(
α
)
=
70
761
[
56
3
-
8
(
α
-
3
)
ϖ
8,3
(
α
)
+
28
(
α
-
1
)
ϖ
8,2
(
α
)
h
h
h
h
-
56
3
(
3
α
-
1
)
ϖ
8,1
(
α
)
+
70
α
ϖ
8,0
(
α
)
]
,
ϖ
8,5
(
α
)
=
56
761
[
56
3
-
8
(
α
-
4
)
ϖ
8,4
(
α
)
+
14
(
2
α
-
3
)
ϖ
8,3
(
α
)
h
h
h
h
-
56
3
(
3
α
-
2
)
ϖ
8,2
(
α
)
+
35
2
(
4
α
-
1
)
ϖ
8,1
(
α
)
h
h
h
h
-
56
α
ϖ
8,0
(
α
)
56
3
]
,
ϖ
8,6
(
α
)
=
140
2283
[
56
5
-
8
(
α
-
5
)
ϖ
8,5
(
α
)
+
28
(
α
-
2
)
ϖ
8,4
(
α
)
h
h
h
h
h
-
56
(
α
-
1
)
ϖ
8,3
(
α
)
+
35
(
2
α
-
1
)
ϖ
8,2
(
α
)
h
h
h
h
h
-
56
5
(
5
α
-
1
)
ϖ
8,1
(
α
)
+
28
α
ϖ
8,0
(
α
)
]
,
ϖ
8,7
(
α
)
=
40
761
[
56
3
-
8
(
α
-
6
)
ϖ
8,6
(
α
)
+
14
(
2
α
-
5
)
ϖ
8,5
(
α
)
h
h
h
h
-
56
3
(
3
α
-
4
)
ϖ
8,4
(
α
)
+
35
2
(
4
α
-
3
)
ϖ
8,3
(
α
)
h
h
h
h
-
56
5
(
5
α
-
2
)
ϖ
8,2
(
α
)
+
14
3
(
6
α
-
1
)
ϖ
8,1
(
α
)
h
h
h
h
-
8
α
ϖ
8,0
(
α
)
56
3
]
,
ϖ
8
,
ℓ
(
α
)
=
280
761
ℓ
[
1
8
-
8
(
α
-
ℓ
+
1
)
ϖ
8
,
ℓ
-
1
(
α
)
h
h
h
h
h
+
14
(
2
α
-
ℓ
+
2
)
ϖ
8
,
ℓ
-
2
(
α
)
h
h
h
h
h
-
56
3
(
3
α
-
ℓ
+
3
)
ϖ
8
,
ℓ
-
3
(
α
)
h
h
h
h
h
+
35
2
(
4
α
-
ℓ
+
4
)
ϖ
8
,
ℓ
-
4
(
α
)
h
h
h
h
h
-
56
5
(
5
α
-
ℓ
+
5
)
ϖ
8
,
ℓ
-
5
(
α
)
h
h
h
h
h
+
14
3
(
6
α
-
ℓ
+
6
)
ϖ
8
,
ℓ
-
6
(
α
)
h
h
h
h
h
-
8
7
(
7
α
-
ℓ
+
7
)
ϖ
8
,
ℓ
-
7
(
α
)
h
h
h
h
h
+
1
8
(
8
α
-
ℓ
+
8
)
ϖ
8
,
ℓ
-
8
(
α
)
]
,
ℓ
=
8,9
,
…
.(3) Whenp
=
9, then u
0
=
7129
/
2520, u
1
=
-
9, u
2
=
18, u
3
=
-
28, u
4
=
63
/
2, u
5
=
-
126
/
5, u
6
=
14, u
7
=
-
36
/
7, u
8
=
9
/
8, and u
9
=
-
1
/
9, u
i
=
0, i
=
10,11
,
…. The coefficients are displayed as follows:
(31)
ϖ
9,0
(
α
)
=
(
7129
2520
)
α
,
ϖ
9,1
(
α
)
=
-
22680
7129
α
ϖ
9,0
(
α
)
,
ϖ
9,2
(
α
)
=
1260
7129
[
-
9
(
α
-
1
)
ϖ
9,1
(
α
)
+
36
α
ϖ
9,0
(
α
)
]
,
ϖ
9,3
(
α
)
=
840
7129
[
-
9
(
α
-
2
)
ϖ
9,2
(
α
)
+
18
(
2
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
84
α
ϖ
9,0
(
α
)
]
,
ϖ
9,4
(
α
)
=
630
7129
[
-
9
(
α
-
3
)
ϖ
9,3
(
α
)
+
36
(
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
28
(
3
α
-
1
)
ϖ
9,1
(
α
)
+
126
α
ϖ
9,0
(
α
)
]
,
ϖ
9,5
(
α
)
=
504
7129
[
63
2
-
9
(
α
-
4
)
ϖ
9,4
(
α
)
+
18
(
2
α
-
3
)
ϖ
9,3
(
α
)
h
h
h
h
h
-
28
(
3
α
-
2
)
ϖ
9,2
(
α
)
+
63
2
(
4
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
126
α
ϖ
9,0
(
α
)
63
2
]
,
ϖ
9,6
(
α
)
=
420
7129
[
126
5
-
9
(
α
-
5
)
ϖ
9,5
(
α
)
+
36
(
α
-
2
)
ϖ
9,4
(
α
)
h
h
h
h
h
-
84
(
α
-
1
)
ϖ
9,3
(
α
)
+
63
(
2
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
1
)
ϖ
9,1
(
α
)
+
84
α
ϖ
9,0
(
α
)
]
,
ϖ
9,7
(
α
)
=
360
7129
[
63
2
-
9
(
α
-
6
)
ϖ
9,6
(
α
)
+
18
(
2
α
-
5
)
ϖ
9,5
(
α
)
h
h
h
h
h
-
28
(
3
α
-
4
)
ϖ
9,4
(
α
)
+
63
2
(
4
α
-
3
)
ϖ
9,3
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
2
)
ϖ
9,2
(
α
)
+
14
(
6
α
-
1
)
ϖ
9,1
(
α
)
h
h
h
h
h
-
36
α
ϖ
9,0
(
α
)
63
2
]
,
ϖ
9,8
(
α
)
=
315
7129
[
36
7
-
9
(
α
-
7
)
ϖ
9,7
(
α
)
+
36
(
α
-
3
)
ϖ
9,6
(
α
)
h
h
h
h
h
-
28
(
3
α
-
5
)
ϖ
9,5
(
α
)
+
126
(
α
-
1
)
ϖ
9,4
(
α
)
h
h
h
h
h
-
126
5
(
5
α
-
3
)
ϖ
9,3
(
α
)
+
28
(
3
α
-
1
)
ϖ
9,2
(
α
)
h
h
h
h
h
-
36
7
(
7
α
-
1
)
ϖ
9,1
(
α
)
+
72
α
ϖ
9,0
(
α
)
]
,
ϖ
9
,
ℓ
(
α
)
=
2520
7129
ℓ
[
-
9
(
α
-
ℓ
+
1
)
ϖ
9
,
ℓ
-
1
(
α
)
+
18
(
2
α
-
ℓ
+
2
)
ϖ
9
,
ℓ
-
2
(
α
)
1
9
h
h
h
h
h
h
-
28
(
3
α
-
ℓ
+
3
)
ϖ
9
,
ℓ
-
3
(
α
)
h
h
h
h
h
h
+
63
2
(
4
α
-
ℓ
+
4
)
ϖ
9
,
ℓ
-
4
(
α
)
h
h
h
h
h
h
-
126
5
(
5
α
-
ℓ
+
5
)
ϖ
9
,
ℓ
-
5
(
α
)
h
h
h
h
h
h
+
14
(
6
α
-
ℓ
+
6
)
ϖ
9
,
ℓ
-
6
(
α
)
h
h
h
h
h
h
-
36
7
(
7
α
-
ℓ
+
7
)
ϖ
9
,
ℓ
-
7
(
α
)
h
h
h
h
h
h
+
9
8
(
8
α
-
ℓ
+
8
)
ϖ
9
,
ℓ
-
8
(
α
)
h
h
h
h
h
h
-
1
9
(
9
α
-
ℓ
+
9
)
ϖ
9
,
ℓ
-
9
(
α
)
]
,
ℓ
=
9,10
,
…
.(4) Whenp
=
10, then u
0
=
7381
/
2520, u
1
=
-
10, u
2
=
45
/
2, u
3
=
-
40, u
4
=
105
/
2, u
5
=
-
252
/
5, u
6
=
35, u
7
=
-
120
/
7, u
8
=
45
/
8, u
9
=
-
10
/
9, u
10
=
1
/
10, and u
i
=
0, i
=
11,12
,
…. The coefficients are given as follows:
(32)
ϖ
10,0
(
α
)
=
(
7381
2520
)
α
,
ϖ
10,1
(
α
)
=
-
25200
7381
α
ϖ
10,0
(
α
)
,
ϖ
10,2
(
α
)
=
1260
7381
[
-
10
(
α
-
1
)
ϖ
10,1
(
α
)
+
45
α
ϖ
10,0
(
α
)
]
,
ϖ
10,3
(
α
)
=
840
7381
[
-
10
(
α
-
2
)
ϖ
10,2
(
α
)
+
45
2
(
2
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
120
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,4
(
α
)
=
630
7381
[
-
10
(
α
-
3
)
ϖ
10,3
(
α
)
+
45
(
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
40
(
3
α
-
1
)
ϖ
10,1
(
α
)
+
210
α
ϖ
10,0
(
α
)
]
,
ϖ
10,5
(
α
)
=
504
7381
[
-
10
(
α
-
4
)
ϖ
10,4
(
α
)
+
45
2
(
2
α
-
3
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
40
(
3
α
-
2
)
ϖ
10,2
(
α
)
+
105
2
(
4
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
252
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,6
(
α
)
=
420
7381
[
252
5
-
10
(
α
-
5
)
ϖ
10,5
(
α
)
+
45
(
α
-
2
)
ϖ
10,4
(
α
)
h
h
h
h
h
-
120
(
α
-
1
)
ϖ
10,3
(
α
)
+
105
(
2
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
1
)
ϖ
10,1
(
α
)
+
210
α
ϖ
10,0
(
α
)
]
,
ϖ
10,7
(
α
)
=
360
7381
[
-
10
(
α
-
6
)
ϖ
10,6
(
α
)
+
45
2
(
2
α
-
5
)
ϖ
10,5
(
α
)
h
h
h
h
h
-
40
(
3
α
-
4
)
ϖ
10,4
(
α
)
+
105
2
(
4
α
-
3
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
2
)
ϖ
10,2
(
α
)
+
35
(
6
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
120
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10,8
(
α
)
=
315
7381
[
120
7
-
10
(
α
-
7
)
ϖ
10,7
(
α
)
+
45
(
α
-
3
)
ϖ
10,6
(
α
)
h
h
h
h
h
-
40
(
3
α
-
5
)
ϖ
10,5
(
α
)
+
210
(
α
-
1
)
ϖ
10,4
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
3
)
ϖ
10,3
(
α
)
+
70
(
3
α
-
1
)
ϖ
10,2
(
α
)
h
h
h
h
h
-
120
7
(
7
α
-
1
)
ϖ
10,1
(
α
)
+
45
α
ϖ
10,0
(
α
)
]
,
ϖ
10,9
(
α
)
=
280
7381
[
-
10
(
α
-
8
)
ϖ
10,8
(
α
)
+
45
2
(
2
α
-
7
)
ϖ
10,7
(
α
)
h
h
h
h
h
-
120
(
α
-
2
)
ϖ
10,6
(
α
)
+
105
2
(
4
α
-
5
)
ϖ
10,5
(
α
)
h
h
h
h
h
-
252
5
(
5
α
-
4
)
ϖ
10,4
(
α
)
+
105
(
2
α
-
1
)
ϖ
10,3
(
α
)
h
h
h
h
h
-
120
7
(
7
α
-
2
)
ϖ
10,2
(
α
)
+
45
8
(
8
α
-
1
)
ϖ
10,1
(
α
)
h
h
h
h
h
-
10
α
ϖ
10,0
(
α
)
45
2
]
,
ϖ
10
,
ℓ
(
α
)
=
2520
7381
ℓ
[
1
10
-
10
(
α
-
ℓ
+
1
)
ϖ
10
,
ℓ
-
1
(
α
)
h
h
h
h
h
h
+
45
2
(
2
α
-
ℓ
+
2
)
ϖ
10
,
ℓ
-
2
(
α
)
h
h
h
h
h
h
-
40
(
3
α
-
ℓ
+
3
)
ϖ
10
,
ℓ
-
3
(
α
)
h
h
h
h
h
h
+
105
2
(
4
α
-
ℓ
+
4
)
ϖ
10
,
ℓ
-
4
(
α
)
h
h
h
h
h
h
-
252
5
(
5
α
-
ℓ
+
5
)
ϖ
10
,
ℓ
-
5
(
α
)
h
h
h
h
h
h
+
35
(
6
α
-
ℓ
+
6
)
ϖ
10
,
ℓ
-
6
(
α
)
h
h
h
h
h
h
-
120
7
(
7
α
-
ℓ
+
7
)
ϖ
10
,
ℓ
-
7
(
α
)
h
h
h
h
h
h
+
45
8
(
8
α
-
ℓ
+
8
)
ϖ
10
,
ℓ
-
8
(
α
)
h
h
h
h
h
h
-
10
9
(
9
α
-
ℓ
+
9
)
ϖ
10
,
ℓ
-
9
(
α
)
h
h
h
h
h
h
+
1
10
(
10
α
-
ℓ
+
10
)
ϖ
10,10
-
ℓ
(
α
)
]
,
f
f
f
f
f
f
f
f
f
f
f
f
f
f
ℓ
=
10,11
,
…
.Next, we plot the figures of coefficients ofϖ
p
,
ℓ
(
α
)
(
p
=
7,8
,
9,10
) to show the evolutions with ℓ. From Figures 5, 6, 7, and 8, we can see that these coefficients are violently oscillatory such that the approximation behaves like Runge’s phenomenon, which is similar to the case of ordinary differential equation. So, to seek high-order (≥7th-order) schemes by this form seems not to be appropriate.Figure 5
The values of the coefficientsϖ
7
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 6
The values of the coefficientsϖ
8
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 7
The values of the coefficientsϖ
9
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.Figure 8
The values of the coefficientsϖ
10
,
ℓ
(
α
)
(
ℓ
=
0,1
,
…
) for α
=
0.3.


## 3. Numerical Examples

In order to verify the reasonability of the coefficients forp
=
3,4
,
5,6, we give the following two numerical examples. These numerical results show that the coefficients are efficient.Example 8.
Consider the function(
x
)
=
x
q
(
q
=
3,4
,
5,6
), x
∈
[
0,1
]. The numerical absolute error and convergence order at x
=
1 by higher-order difference scheme (14) with different p
(
p
=
3,4
,
5,6
) are shown in Tables 1, 2, 3, and 4.Table 1
The numerical results of Example8 with q
=
3 by formula (14) with p
=
3.
α
h
The absolute error
The convergence order
0.2
1
/
10
2.645451
e
-
004
—
1
/
20
3.262235
e
-
005
3.0196
1
/
40
4.051498
e
-
006
3.0093
1
/
80
5.048441
e
-
007
3.0045
0.4
1
/
10
4.245164
e
-
004
—
1
/
20
5.164722
e
-
005
3.0391
1
/
40
6.373660
e
-
006
3.0185
1
/
80
7.917542
e
-
007
3.0090
0.6
1
/
10
4.388473
e
-
004
—
1
/
20
5.266513
e
-
005
3.0588
1
/
40
6.457975
e
-
006
3.0277
1
/
80
7.997522
e
-
007
3.0135
0.8
1
/
10
2.902946
e
-
004
—
1
/
20
3.435614
e
-
005
3.0789
1
/
40
4.185974
e
-
006
3.0369
1
/
80
5.167862
e
-
007
3.0179Table 2
The numerical results of Example8 with q
=
4 by formula (14) with p
=
4.
α
h
The absolute error
The convergence order
0.2
1
/
10
8.565957
e
-
005
—
1
/
20
5.244691
e
-
006
4.0297
1
/
40
3.248631
e
-
007
4.0130
1
/
80
2.021606
e
-
008
4.0063
0.4
1
/
10
1.389875
e
-
004
—
1
/
20
8.345704
e
-
006
4.0578
1
/
40
5.123066
e
-
007
4.0578
1
/
80
3.174360
e
-
008
4.0125
0.6
1
/
10
1.451794
e
-
004
—
1
/
20
8.553695
e
-
006
4.0851
1
/
40
5.203490
e
-
007
4.0390
1
/
80
3.210299
e
-
008
4.0187
0.8
1
/
10
9.693830
e
-
005
—
1
/
20
5.608540
e
-
006
4.1114
1
/
40
3.381060
e
-
007
4.0521
1
/
80
2.076976
e
-
008
4.0249Table 3
The numerical results of Example8 with q
=
5 by formula (14) with p
=
5.
α
h
The absolute error
The convergence order
0.2
1
/
5
1.243084
e
-
003
—
1
/
10
3.601600
e
-
005
5.1091
1
/
20
1.098356
e
-
006
5.0352
1
/
40
3.392476
e
-
008
5.0169
0.4
1
/
5
2.128798
e
-
003
—
1
/
10
5.880181
e
-
005
5.1780
1
/
20
1.756897
e
-
006
5.0648
1
/
40
5.363365
e
-
008
5.0337
0.6
1
/
5
2.304179
e
-
003
—
1
/
10
6.169164
e
-
005
5.2230
1
/
20
1.810078
e
-
006
5.0910
1
/
40
5.461289
e
-
008
5.0507
0.8
1
/
5
1.548320
e
-
003
—
1
/
10
4.138297
e
-
005
5.2255
1
/
20
1.193106
e
-
006
5.1162
1
/
40
3.557536
e
-
008
5.0677Table 4
The numerical results of Example8 with q
=
6 by formula (14) with p
=
6.
α
h
The absolute error
The convergence order
0.2
1
/
5
1.366727
e
-
003
—
1
/
10
1.834110
e
-
005
6.2195
1
/
20
2.831669
e
-
007
6.0173
1
/
40
4.352325
e
-
009
6.0237
0.4
1
/
5
2.458586
e
-
003
—
1
/
10
2.979151
e
-
005
6.3668
1
/
20
4.561725
e
-
007
6.0292
1
/
40
6.997415
e
-
009
6.0266
0.6
1
/
5
2.796453
e
-
003
—
1
/
10
3.106344
e
-
005
6.4922
1
/
20
4.745505
e
-
007
6.0325
1
/
40
7.487801
e
-
009
5.9859
0.8
1
/
5
81.968603
e
-
003
—
1
/
10
2.092182
e
-
005
6.5560
1
/
20
3.168932
e
-
007
6.0449
1
/
40
4.769865
e
-
009
6.0539Example 9.
Let us consider a fractional ordinary differential equation(33)
D
RL
0
,
x
α
y
(
x
)
=
Γ
(
8
-
α
)
Γ
(
8
-
2
α
)
x
7
-
2
α
-
2880
Γ
(
7
-
α
)
x
6
-
α
+
Γ
(
6
+
α
)
20
x
5
,
with initial values
(34)
D
RL
0
,
x
α
-
1
y
(
x
)
|
x
=
0
=
0
,
α
∈
(
0,1
)
.The exact solution of the above equation is given by(35)
y
(
x
)
=
x
7
-
α
-
4
x
6
+
6
x
5
+
α
.At this moment, we use the numerical formula (14) with different order p to solve this equation. The absolute error and numerical convergence order are listed in Tables 5, 6, 7, and 8.Table 5
The numerical results of Example9 by formula (14) with p
=
3.
α
h
The absolute error
The convergence order
0.2
1
/
20
8.448207
e
-
004
—
1
/
40
1.101436
e
-
004
2.9393
1
/
80
1.406092
e
-
005
2.9696
1
/
160
1.776240
e
-
006
2.9848
0.4
1
/
20
2.742719
e
-
003
—
1
/
40
3.594722
e
-
004
2.9317
1
/
80
4.599387
e
-
005
2.9664
1
/
160
5.816133
e
-
006
2.9816
0.6
1
/
20
6.698506
e
-
003
—
1
/
40
8.826119
e
-
004
2.9240
1
/
80
1.132122
e
-
004
2.9628
1
/
160
1.433358
e
-
005
2.9816
0.8
1
/
20
1.458663
e
-
002
—
1
/
40
1.931852
e
-
003
2.9166
1
/
80
2.484219
e
-
004
2.9591
1
/
160
3.149142
e
-
005
2.9798Table 6
The numerical results of Example9 by formula (14) with p
=
4.
α
h
The absolute error
The convergence order
0.2
1
/
10
7.061963
e
-
004
—
1
/
20
4.697130
e
-
005
3.9102
1
/
40
3.067961
e
-
006
3.9364
1
/
80
1.966873
e
-
007
3.9633
0.4
1
/
10
2.626111
e
-
003
—
1
/
20
1.735473
e
-
004
3.9195
1
/
40
1.117528
e
-
005
3.9569
1
/
80
7.093827
e
-
007
3.9776
0.6
1
/
10
7.072031
e
-
003
—
1
/
20
4.720003
e
-
004
3.9053
1
/
40
3.043142
e
-
005
3.9552
1
/
80
1.931056
e
-
006
3.9781
0.8
1
/
10
1.660047
e
-
002
—
1
/
20
1.125103
e
-
003
3.8831
1
/
40
7.302226
e
-
005
3.9456
1
/
80
4.647896
e
-
006
3.9737Table 7
The numerical results of Example9 by formula (14) with p
=
5.
α
h
The absolute error
The convergence order
0.2
1
/
5
3.255980
e
-
004
—
1
/
10
2.611815
e
-
005
3.6400
1
/
20
1.709675
e
-
006
3.9333
1
/
40
7.187111
e
-
008
4.5722
0.4
1
/
5
4.502359
e
-
003
—
1
/
10
1.361298
e
-
004
5.0476
1
/
20
4.930325
e
-
006
4.7872
1
/
40
1.689483
e
-
007
4.8670
0.6
1
/
5
1.726870
e
-
002
—
1
/
10
4.858555
e
-
004
5.1515
1
/
20
1.493401
e
-
005
5.0239
1
/
40
4.654620
e
-
007
5.0038
0.8
1
/
5
4.809699
e
-
002
—
1
/
10
1.421046
e
-
003
5.0809
1
/
20
4.352885
e
-
005
5.0288
1
/
40
1.348082
e
-
006
5.0130Table 8
The numerical results of Example9 by formula (14) with p
=
6.
α
h
The absolute error
The convergence order
0.2
1
/
5
1.624059
e
-
003
—
1
/
10
1.004444
e
-
005
7.3371
1
/
20
4.333263
e
-
007
4.5348
1
/
40
7.431670
e
-
009
5.8656
0.4
1
/
5
3.378083
e
-
003
—
1
/
10
1.682553
e
-
006
10.9713
1
/
20
3.570178
e
-
007
2.2366
1
/
40
7.098342
e
-
009
5.6524
0.6
1
/
5
4.899244
e
-
003
—
1
/
10
1.910223
e
-
005
8.0027
1
/
20
3.435689
e
-
008
9.1189
1
/
40
2.811996
e
-
009
3.6109
0.8
1
/
5
4.824629
e
-
003
—
1
/
10
3.070796
e
-
005
7.2957
1
/
20
3.129379
e
-
007
6.6166
1
/
40
3.973945
e
-
009
6.2992From the numerical results presented in Tables1–8, we can see that the coefficients of the fractional linear multistep method for p
=
3,4
,
5,6 are efficient.The coefficients ofϖ
p
,
ℓ
(
α
)
(
p
=
7,8
,
9,10
) are violently oscillatory which may not be suitable for numerical calculations. Now, we take an example to show this.Example 10.
Consider the functionf
(
x
)
=
x
3, x
∈
[
0,1
]. Consider
(36)
D
RL
0
,
x
α
x
3
=
Γ
(
4
)
Γ
(
4
-
α
)
x
3
-
α
.
We use the scheme (14) with p
=
7,8
,
9,10 to numerically compute D
RL
0
,
x
α
x
3. See Figures 9, 10, 11, and 12. From these figures, we can see that the results are not numerically stable, which can be regarded as fractional Runge’s phenomenon. So, it is not necessary to derive more than 6th-order schemes for Riemann-Liouville derivative by generating function method.Figure 9
The exact value ofD
RL
0,1
α
x
3 and numerical solution of p
=
7 for α
=
0.5 and step length of h
=
1
/
120.Figure 10
The exact value ofD
RL
0,1
α
x
3 and numerical solution of p
=
8 for α
=
0.5 and step length of h
=
1
/
120.Figure 11
The exact value ofD
RL
0,1
α
x
3 and numerical solution of p
=
9 for α
=
0.5 and step length of h
=
1
/
120.Figure 12
The exact value ofD
RL
0,1
α
x
3 and numerical solution of p
=
10 for α
=
0.5 and step length of h
=
1
/
120.


## 4. Conclusion

In this paper, we propose recursion formulas to compute the coefficients of the fractional linear multistep schemes. The numerical experiments have been carried out to support the derived numerical schemes. Here, we should note that thepth order (
p
≥
7
) schemes are not stable.



---

*Source: 402373-2014-04-15.xml*