survey of the history of the fractional calculus. 1. Introduction and Preliminaries. The concept of the differentiation operator D = d/dx is familiar to all who.
Comm. Korean Math. Soc. 12 (1997), No. 2, pp. 457–466
A FORMULA RELATED TO FRACTIONAL CALCULUS JUNESANG CHOI ABSTRACT. We give a formula related to fractional calculus which may be useful in solving some fractional linear differential equations. We also give a brief survey of the history of the fractional calculus.
1. Introduction and Preliminaries The concept of the differentiation operator D = d/d x is familiar to all who have studied the elementary calculus. And for suitable functions f , the nth derivative of f , namely D n f (x) = d n f (x)/d x n is well-defined provided that n is a positive integer. In 1695, L’Hˆopital inquired of Leibniz what meaning could be ascribed to D n f if n were a fraction. Since that time the fractional calculus has drawn the attention of many famous mathematicians, such as Euler, Laplace, Fourier, Abel, Liouville, Riemann, and Laurent. But it was not until 1884 that the theory of generalized operators achieved a level in its development suitable as a point of departure for the modern mathematician. By then the theory had been extended to include operators D ν , where ν could be rational or irrational, positive or negative, real or complex. Thus the name fractional calculus became somewhat of a misnomer. A better description might be differentiation and integration to an arbitrary order. However, by tradition this theory is referred to as the fractional calculus. For more and detailed historical survey, see [1], [4], [6]. As noted in [4], in the period 1975 to the present, about 400 papers have been published relating to the fractional calculus. The fractional calculus have been used in many fields of science and engineering. We give some references used in mathematics [2], [3], [5]. More than one version of the fractional integral exists. We introduce only one version here. Received September 25, 1996. Revised October 16, 1996. 1991 Mathematics Subject Classification: Primary 26A33; Secondary 34A30. Key words and phrases: fractional calculus.
458
Junesang Choi
Let Re(ν) > 0 and let f be piecewise continuous on (0, ∞) and integrable on any finite subinterval of [0, ∞). Then for t > 0 we call
−ν 0 Dt
(1.1)
Z
1 f (t) = 0(ν)
t
(t − ξ )ν−1 f (ξ )dξ
0
the Riemann-Liouville fractional integral of f of order ν and 0 is the wellknown Gamma function. We also denote 0 Dt−ν as D −ν . In particular setting f (t) = eat (a constant) in (1.1) yields
(1.2)
D
1 e = 0(ν)
−ν at
Z
t
(t − ξ )ν−1 eaξ dξ, ν > 0.
0
If we make the change of variable x = t − ξ, (1.2) becomes (1.3)
D
eat e = 0(ν)
−ν at
Z
t
x ν−1 e−ax d x, ν > 0,
0
which is not an elementary function. For Re(ν) > 0 the incomplete gamma function γ ∗ (ν, t) may be defined as (1.4)
1 γ (ν, t) = 0(ν)t ν ∗
Z
t
ξ ν−1 e−ξ dξ.
0
Thus we may write (1.2) as (1.5)
D −ν eat = t ν eat γ ∗ (ν, at).
Since the right-hand side of (1.5) is the fractional integral of an exponential, it is not surprising that this function frequently arises in the study of the fractional calculus. We shall call it E t (ν, a), (1.6)
E t (ν, a) = t ν eat γ ∗ (ν, at).
Now the concept of the fractional derivative is introduced. Let f satisfy conditions of the fractional integral (1.1) and let µ > 0. Let m be the smallest
A Formula related to fractional calculus
459
integer that exceeds µ. Then the fractional derivative of f of order µ is defined as (1.7)
D µ f (t) = D m [D −ν f (t)], µ > 0, t > 0
(if it exists) where ν = m − µ > 0. Suppose that f (t) = eat . Then the fractional derivative of eat of order µ is (1.8)
D µ eat = D m [D −ν eat ],
where µ, ν and m have the same meaning as above. With the help of (1.6), we readily obtain (1.9)
D m E t (ν, a) = E t (ν − m, a) = E t (−µ, a)
since µ = m − ν. Thus we conclude that (1.10)
D µ eat = E t (−µ, a), t > 0.
Now we introduce the concept of fractional differential equations. As a first attempt define a fractional differential equation, let rm , rm−1 , · · · , r0 be a strictly deceasing sequence of nonnegative numbers. Then if b1 , b2 , · · · , bm are constants, rm (1.11) D + b1 D rm−1 + · · · + bm D r0 y(t) = 0 is a candidate. But even this equation is a little too complex. We shall impose the additional requirement that the r j be rational numbers. Thus if q is the least common multiple of the denominators of the nonzero r j , we may write (1.11) as nν (1.12) D + a1 D (n−1)ν + · · · + an D 0 y(t) = 0, t ≥ 0, where (1.13)
ν=
1 . q
If q = 1 , then ν = 1, and (1.12) is simply an ordinary differential equation. We shall call (1.12) a fractional linear differential equation with constant coefficients of order (n, q), or more briefly, a fractional differential equation of order (n, q). In this note we are aiming at getting a formula which can be occurred in solving the differential equation (1.12).
460
Junesang Choi
2. Explicit representation of the solution of (1.12) Let P(x) = x n + a1 x n−1 + · · · + an
(2.1)
be the corresponding indicial polynomial to (1.12). Let (2.2) y1 (t) = L−1 P −1 (s ν ) , where L−1 denotes the inverse Laplace transform. Then if N is the smallest integer with the property that N ≥ nν, y1 (t), y2 (t), . . . , y N (t), where yj +1 (t) = D j y1 (t),
j = 0, 1, . . . , N − 1
are N linearly independent solutions of (1.12) which is well-known. We introduce the Laplace transform method to solve the equation (1.12). Let α1 , . . . , αn with αi 6= α j for i 6= j be the zeros of P(x) and let (2.3)
A−1 m = D P(αm ),
m = 1, 2, . . . , n.
But (2.4)
P −1 (s ν ) =
sν
A1 A2 An + ν +···+ ν , − α1 s − α2 s − αn
and we can readily obtain (2.5)
L
−1
1 ν s −a
=
q X
a j −1 E t ( j ν − 1, a q ).
j =1
If we make the change of dummy index of summation k = q − j in the sum of (2.5), (2.6)
L
−1
1 ν s −a
=
q−1 X k=0
a q−k−1 E t (−kν, a q ).
A Formula related to fractional calculus
461
Using (2.6) we see the inverse Laplace transform of (2.4) yields y1 (t) =
(2.7)
n X
Am
m=1
q−1 X
αmq−k−1 E t (−kν, αmq )
k=0
is a solution of (1.12). Hence y2 , y3 . . . , y N of the solutions of (1.12) may be constructed from (2.7). We begin our method of solving (1.12) by recalling that D pν E t (−kν, a) = E t (−(k + p)ν, a)
(2.8)
provided that kν < 1, which it is for k = 0, 1, . . . , q − 1. If we define e(t) as e(t) =
(2.9)
q−1 X
α q−k−1 E t (−kν, α q )
k=0
(where for the moment α is an arbitrary constant), then (2.8) implies that D ν e(t) = αe(t).
(2.10)
For p a positive integer greater than 1, D e(t) = α e(t) + pν
(2.11)
p
p−1 p−1−k −1−kν X α t k=1
0(−kν)
.
Formula (2.11) is also valid for p = 0 or 1, since in these cases the sum in (2.11) is vacuous. Thus if we write n X (2.12) P(D ν )e(t) = an− p D pν e(t), a0 = 1, p=0
(2.11) implies that (2.13)
ν
P(D )e(t) = P(α)e(t) +
n X p=2
an− p
p−1 p−1−k −1−kν X α t k=1
0(−kν)
.
462
Junesang Choi
Now if α1 , . . . , αn are the distinct zeros of P(x), and if ej (t) =
(2.14)
q−1 X
q−k−1
αj
q
E t (−kν, α j ),
k=0
then certainly, from (2.13), " n # n X X ν P(D ) cm em (t) = cm P(αm )em (t) m=1
(2.15)
+
m=1
n X
cm
m=1
n X
an− p
p=2
p−1 p−1−k −1−kν X αm t k=1
0(−kν)
for any arbitrary constants c1 , c2 , . . . , cn . But since α1 , . . . , αn are the roots of P(x) = 0, the first term on the right-hand side of (2.15) vanishes and # " n p−1 −1−kν X X t P(D ν ) cm em (t) = an− p cm αmp−1−k . 0(−kν) m=1 p=2 k=1 m=1 "
(2.16)
n X
#
n X
Thus if we can choose the cm such that the right-hand side of (2.16) vanishes, then (2.17)
y1 (t) =
n X
cm em (t)
m=1
will be a solution of (1.12). Ignoring the trivial solution c1 = c2 = · · · = cn = 0, we see that if we let cm = Am , where the Am are given by (2.3), then we know that the sum in brackets in (2.16) is zero. Thus (2.18)
y1 (t) =
n X
Am em (t)
m=1
is a solution of (1.12). But this is precisely (2.7). Note that in the latter method, the role of (2.11) is important.
A Formula related to fractional calculus
463
If the indicial polynomial P(x) in (2.1) has a double zero at α, then e(t)∗e(t) is used to solve (1.12), where e(t)∗e(t) is the convolution of e(t) with itself. We also have an explicit representation of e(t) ∗ e(t) in terms of the E t functions, namely, e(t) ∗ e(t) =L−1 (s ν − α)−2 =
(2.19)
q−1 q−1 X X
α 2q− j −k−2 {t E t (−( j + k)ν, α q )
j =0 k=0
+ ( j + k)ν E t (1 − ( j + k)ν, α q )}. Finally we observe that an explicit representation of D pν e(t) ∗ e(t) is as important as that of D pν e(t) in (2.11) in solving (1.12) if the corresponding indicial polynomial P(x) has a double zero at α. Fortunately we can give an explicit representation of D pν e(t) ∗ e(t): For p a positive integer, we have D pν e(t) ∗ e(t) =α p e(t) ∗ e(t) + pα p−1 e(t) (2.20)
+
p−1 j −1 X X
α p−k−2
j =0 k=0
t −kν−1 . 0(−kν)
3. Proof of (2.20) First we give some elementary identities needed to our purpose without proof : (3.1)
(3.2)
(3.3)
(3.4)
D µ E t (λ, a) = E t (λ − µ, a). D µ t E t (ω, c) = t E t (ω − µ, c) + µE t (ω − µ + 1, c), ω > −2. E t (−1, a) = a E t (0, a).
E t (ν, a) = a E t (ν + 1, a) +
tν . 0(ν + 1)
464
Junesang Choi
For convenience, let F( pν) = D pν e(t) ∗ e(t). With the aid of (3.1) and (3.2) we obtain F( pν) =
q−1 q−1 X X
α 2q− j −k−2 {D pν t E t (−( j + k)ν, α q )
j =0 k=0
+ ( j + k)ν D pν E t (1 − ( j + k)ν, α q )} =
q−1 q−1 X X
α 2q− j −k−2 {t E t (−( j + k + p)ν, α q )
j =0 k=0
+ pν E t (1 − ( j + k + p)ν, α q ) + ( j + k)ν E t (1 − ( j + k + p)ν, α q )}. Setting j 0 = j + p in the resulting equation and dropping the prime on j 0 , we get X X
q+ p−1 q−1
F( pν) =
j=p
(3.5)
=
I ( p, q, j, k)
k=0
q−1 q−1 X X j =0 k=0
−
X X
q+ p−1 q−1
I ( p, q, j, k) +
p−1 q−1 X X
j =q
I ( p, q, j, k)
k=0
I ( p, q, j, k),
j =0 k=0
where I ( p, q, j, k) =α 2q− j −k+ p−2 t E t (−( j + k)ν, α q ) +( j + k)ν E t (1 − ( j + k)ν, α q ) . Setting j 0 = j − q in the second sum of (3.5) and dropping the prime on j 0
A Formula related to fractional calculus
465
yields F( pν) =α p e(t) ∗ e(t) +
(3.6)
p−1 q−1 X X
α q− j −k+ p−2 {t E t (−( j + k)ν − 1, α q )
j =0 k=0
+ (( j + k)ν + 1)E t (−( j + k)ν, α q ) − α q t E t (−( j + k)ν, α q ) − α q ( j + k)ν E t (1 − ( j + k)ν, α q )}. Using (3.4) we find that t −( j +k)ν , 0(−( j + k)ν) (3.7) ( j + k)ν E t (−( j + k)ν, α q ) = ( j + k)να q E t (1 − ( j + k)ν, α q ) t E t (−( j + k)ν − 1, α q ) = α q t E t (−( j + k)ν, α q ) + t −( j +k)ν . − 0(−( j + k)ν) Putting (3.7) into (3.6) leads to (3.8)
F( pν) = α e(t) ∗ e(t) + q
q−1 p−1 X X
α q− j −k+ p−2 E t (−( j + k)ν, α q ).
k=0 j =0
Setting j + k = k 0 in the summation part of (3.8) and dropping the prime on k 0 , we obtain (3.9) p−1 q−1 p−1 q+ j −1 X X X X q− j −k+ p−2 q α E t (−( j + k)ν, α ) = α q−k+ p−2 E t (−kν, α q ). j =0 k=0
j =0
k= j
Now we also separate the inner sum in the right side of (3.9) into three parts: q+ j −1
X
(3.10)
α
q−k+ p−2
E t (−kν, α ) = q
k= j
q−1 X k=0
q+ j −1
+
X k=q
α q−k+ p−2 E t (−kν, α q )
α
q−k+ p−2
E t (−kν, α ) − q
j −1 X k=0
α q−k+ p−2 E t (−kν, α q ).
466
Junesang Choi
Setting k − q = k 0 in the second sum of the right side of (3.10) and dropping the prime on k 0 with the aid of (3.4) and (2.9) yields q+ j −1
(3.11)
X
α q−k+ p−2 E t (−kν, α q ) = α p−1 e(t) +
k= j
j −1 X
α −k+ p−2
k=0
t −kν−1 . 0(−kν)
Finally putting (3.11) into (3.9) and considering (3.8) leads to the desired result (2.20). We conclude this note by giving some special cases of (2.20). Setting p = 1, 2, 3 in (2.20) with 1/ 0(0) = 0, we obtain D ν e(t) ∗ e(t) = αe(t) ∗ e(t) + e(t), (3.12)
D 2ν e(t) ∗ e(t) = α 2 e(t) ∗ e(t) + 2αe(t), D 3ν e(t) ∗ e(t) = α 3 e(t) ∗ e(t) + 3α 2 e(t) +
t −ν−1 . 0(−ν)
References [1] B. Ross, Editor, Fractional Calaclus and its Applications, LNM 457, Springer-Verlag, 1975. [2] B. Ross, Serenpidity in Mathematics or How one is led to discover that..., Amer. Math. Monthly 90 (1983), 562–566. [3] J. A. Dur´an, S. L. Kalla and H. M. Srivastava, Fractional calculus and the sums of certain families of infinite series, J. M. A. A. 190 (1995), 738–754. [4] K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley snd Sons, Inc., 1993. [5] K. Nishimoto, Unification of the integrals and derivatives (A serenpidity in fractional calculus), Journal of Fractional Calculus 6 (1994), 1–14. [6] S. G. Samko, A. A. Kilbas and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach Science Publishers, 1993.
Department of Mathematics Dongguk University Kyongju 780-714, Korea