9611

Full download ebooks at ebookmeta.
com
Fractional Differential Equations - An Approach

via Fractional Derivatives 1st Edition Bangti Jin
https://ebookmeta.com/product/fractional-differential-
equations-an-approach-via-fractional-derivatives-1st-
edition-bangti-jin/
OR CLICK BUTTON
DOWLOAD NOW
Download more ebook from https://ebookmeta.com

More products digital (pdf, epub, mobi) instant
download maybe you interests ...
Fractional Stochastic Differential Equations -

Applications to Covid-19 Modeling 1st Edition Abdon
Atangana
https://ebookmeta.com/product/fractional-stochastic-differential-
equations-applications-to-covid-19-modeling-1st-edition-abdon-
atangana/
Numerical Methods for Fractal-Fractional Differential

Equations and Engineering 1st Edition Muhammad Altaf
Khan
https://ebookmeta.com/product/numerical-methods-for-fractal-
fractional-differential-equations-and-engineering-1st-edition-
muhammad-altaf-khan/
Methods of Mathematical Modelling Fractional

Differential Equations Mathematics and its Applications
1st Edition Harendra Singh (Editor)
https://ebookmeta.com/product/methods-of-mathematical-modelling-
fractional-differential-equations-mathematics-and-its-
applications-1st-edition-harendra-singh-editor/
Natures Patterns and the Fractional Calculus Fractional

Calculus in Applied Sciences and Engineering 1st
Edition West
https://ebookmeta.com/product/natures-patterns-and-the-
fractional-calculus-fractional-calculus-in-applied-sciences-and-
engineering-1st-edition-west/
Fractional Order Systems Ivo Petráš Editor
https://ebookmeta.com/product/fractional-order-systems-ivo-
petras-editor/
Differential Equations: A Linear Algebra Approach 1st

Edition Anindya Dey
https://ebookmeta.com/product/differential-equations-a-linear-
algebra-approach-1st-edition-anindya-dey/
Stochastic Models for Fractional Calculus Mark M.

Meerschaert
https://ebookmeta.com/product/stochastic-models-for-fractional-
calculus-mark-m-meerschaert/
Fractional Inequalities In Banach Algebras 1st Edition

George A. Anastassiou
https://ebookmeta.com/product/fractional-inequalities-in-banach-
algebras-1st-edition-george-a-anastassiou/
An Elementary Course on Partial Differential Equations

Aftab Alam
https://ebookmeta.com/product/an-elementary-course-on-partial-
differential-equations-aftab-alam/
Applied Mathematical Sciences
Bangti Jin
Fractional
Differential
Equations
An Approach via Fractional Derivatives
Volume 206
Series Editors
Anthony Bloch, Department of Mathematics, University of Michigan, Ann Arbor,
MI, USA
C. L. Epstein, Department of Mathematics, University of Pennsylvania,
Philadelphia, PA, USA
Alain Goriely, Department of Mathematics, University of Oxford, Oxford, UK
Leslie Greengard, New York University, New York, NY, USA
Advisory Editors
J. Bell, Center for Computational Sciences and Engineering, Lawrence Berkeley
National Laboratory, Berkeley, CA, USA
P. Constantin, Department of Mathematics, Princeton University, Princeton, NJ,
USA
R. Durrett, Department of Mathematics, Duke University, Durham, CA, USA
R. Kohn, Courant Institute of Mathematical Sciences, New York University,
New York, NY, USA
R. Pego, Department of Mathematical Sciences, Carnegie Mellon University,
Pittsburgh, PA, USA
L. Ryzhik, Department of Mathematics, Stanford University, Stanford, CA, USA
A. Singer, Department of Mathematics, Princeton University, Princeton, NJ, USA
A. Stevens, Department of Applied Mathematics, University of Münster, Münster,
Germany
S. Wright, Computer Sciences Department, University of Wisconsin, Madison, WI,
USA
Founding Editors
F. John, New York University, New York, NY, USA
J. P. LaSalle, Brown University, Providence, RI, USA
L. Sirovich, Brown University, Providence, RI, USA
The mathematization of all sciences, the fading of traditional scientific boundaries,
the impact of computer technology, the growing importance of computer modeling
and the necessity of scientific planning all create the need both in education and
research for books that are introductory to and abreast of these developments. The
purpose of this series is to provide such books, suitable for the user of mathematics,
the mathematician interested in applications, and the student scientist. In particular,
this series will provide an outlet for topics of immediate interest because of the
novelty of its treatment of an application or of mathematics being applied or lying
close to applications. These books should be accessible to readers versed in
mathematics or science and engineering, and will feature a lively tutorial style, a
focus on topics of current interest, and present clear exposition of broad appeal.
A compliment to the Applied Mathematical Sciences series is the Texts in Applied
Mathematics series, which publishes textbooks suitable for advanced undergraduate
and beginning graduate courses.
More information about this series at http://www.springer.com/series/34

Bangti Jin
Fractional Differential
Equations
An Approach via Fractional Derivatives
123
Bangti Jin
Department of Computer Science
University College London
London, UK
ISSN 0066-5452 ISSN 2196-968X (electronic)

ISBN 978-3-030-76042-7 ISBN 978-3-030-76043-4 (eBook)
https://doi.org/10.1007/978-3-030-76043-4
Mathematics Subject Classification: 26A33, 33E12, 35R11, 34A08, 35R30
© The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature
Switzerland AG 2021
This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether
the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of
illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and
transmission or information storage and retrieval, electronic adaptation, computer software, or by similar
or dissimilar methodology now known or hereafter developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this
publication does not imply, even in the absence of a specific statement, that such names are exempt from
the relevant protective laws and regulations and therefore free for general use.
The publisher, the authors and the editors are safe to assume that the advice and information in this
book are believed to be true and accurate at the date of publication. Neither the publisher nor the
authors or the editors give a warranty, expressed or implied, with respect to the material contained
herein or for any errors or omissions that may have been made. The publisher remains neutral with regard
to jurisdictional claims in published maps and institutional affiliations.
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface
Fractional differential equations (FDES), i.e., differential equations involving

fractional-order derivatives, have received much recent attention in engineering,
physics, biology and mathematics, due to their extraordinary modeling capability
for describing certain anomalous transport phenomena observed in real world.
However, the relevant mathematical theory of FDES is still far from complete when
compared with the more established integer-order counterparts.
The main modeling tools, i.e., fractional integrals and fractional derivatives,
have a history nearly as old as calculus itself. In 1695, in a letter to Leibniz,
L’Hospital asked about the possible meaning of a half derivative. Leibniz in a
letter dated September 30, 1695—now deemed as the birthday of fractional cal-
culus—replied “It will lead to a paradox, from which one day useful consequences
will be drawn.” The first definitions of general fractional derivatives, e.g.,
Riemann-Liouville fractional integral and derivatives, that supported rigorous
analysis were proposed in the nineteenth century. By the middle of the twentieth
century, fractional calculus was already an almost fully developed field. The
practical applications of these operators, as conjectured by Leibniz, were realized
only much later. Recent experimental observations of “fractional” diffusion have
led to a flurry of studies in physics and engineering. By now FDES have been
established as the lingua franca in the study of anomalous diffusion processes.
This motivates a mathematical study of differential equations containing one or
more fractional derivatives. This is the primary purpose of this textbook: we give a
self-contained introduction to basic mathematical theory of FDES, e.g.,
well-posedness (existence, uniqueness and regularity of the solutions) and other
analytic properties, e.g., comparison principle, maximal principle and analyticity,
and contrast them with the integer counterparts. The short answer is that there is
much in common, but with some significant and important differences and addi-
tional technical challenges. One distinct feature is that the solution usually con-
tains weak singularity, even if the problem data are smooth, due to the nonlocality
of the operators and weak singularity of the associated kernels. These new features
have raised many outstanding challenges for relevant fields, e.g., mathematical
analysis, numerical analysis, inverse problems and optimal control, which have
v
vi Preface
witnessed exciting developments in recent years. We will also brief ly touch upon
these topics to give a f lavor, illustrating distinct influences of the nonlocality of
fractional operators.
The term FDE is a big umbrella for any differential equations involving one or
more fractional derivatives. This textbook only considers the models that involve
only one fractional derivative in either spatial or temporal variable, mostly of the
Djrbashian-Caputo type, and only brief ly touches upon that of the Riemann-
Liouville type. This class of FDES has direct physical meaning and allows
initial/boundary values as for the classical integer-order differential equations, and
thus has been predominant in the mathematical modeling of many practical
applications.
The book consists of seven chapters and one appendix, which are organized as
follows. Chapter 1 briefly describes continuous time random walk, which shows
that the probability density function of certain stochastic process satisfies a FDE.
Chapters 2 and 3 describe, respectively, fractional calculus (various fractional
integral and derivatives), and two prominent special functions, i.e., Mittag-Leffler
function and Wright function. These represent the main mathematical tools for
FDES, and thus Chapters 2 and 3 build the foundation for the study of FDES in
Chapters 4–7. Chapters 4–5 are devoted to initial and boundary value problems for
fractional ODEs, respectively. Chapters 6–7 describe mathematical theory for
time-fractional diffusion (subdiffusion) in Hilbert spaces and Hölder spaces,
respectively. In the appendix, we recall basic facts of function spaces, integral
transforms and fixed point theorems that are extensively used in the book. Note that
Chapters 4–7 are organized in a manner such that they only loosely depend on each
other and the reader can directly proceed to the chapter of interest and study each
chapter independently, after going through Chapters 2–3 (and the respective parts
of the appendix).
The selected materials focus mainly on solution theory, especially the issues of
existence, uniqueness and regularity, which are also important in other related
areas, e.g., numerical analysis. Needless to say, these topics can represent only a
few small tips of current exciting developments within the FDE community, and
many important topics and results have to be left out due to page limitation. Also
we do not aim for a comprehensive treatment and describing the results in the most
general form, and instead we only state the results that are commonly used in order
to avoid unnecessary technicalities. (Actually each chapter can be expanded into a
book itself!) Whenever known, references for the stated results are provided, and
the pointers to directly relevant further extensions are also given. Nonetheless, the
reference list is nowhere close to complete and clearly biased by the author’s
knowledge.
There are many excellent books available, especially at the research level,
focusing on various aspects of fractional calculus or FDEs, notably [SKM93] on
Riemann-Liouville fractional integral/derivatives, [Pod99, KST06] on fractional
derivatives and ODEs, [Die10] on ODEs with Djrbashian-Caputo fractional derivative.
Our intention is to take a broader scope than the more focused research monographs
and to provide a relatively self-contained gentle introduction to current
Preface vii
mathematical theory, especially distinct features when compared with the

integer-order counterparts. Besides, it contains discussions on numerical algorithms
and inverse problems, which have so far been only scarcely discussed in textbooks
on fdes, despite their growing interest in the community. Thus, it may serve as an
introductory reading for young researchers entering the field. It can also be used as
the textbook for a single semester course. Many exercises (of different degree of
difficulty) are provided at the end of each chapter, and quite a few of them give
further extensions to the results stated in the book. The typical audience would be
senior undergraduate and graduate students with a solid background in real analysis
and basic theory of ordinary and partial differential equations, especially Chapters 6
–7.
The book project was initiated in 2015, and but stopped in Summer 2016. The
writing was resumed in later 2019, and the scope and breadth have grown sig-
nificantly, in order to properly reflect several important developments in the last
few years. The majority of the writing was done during the COVID 19 pandemic, and
partly during an extended research stay in 2020 at Department of Mathematics,
The Chinese University of Hong Kong, when the book was being finalized. The
hospitality is greatly appreciated. The writing has benefited enormously from the
collaborations with several researchers on the topic of FDES during the past decade,
especially Raytcho Lazarov, Buyang Li, Yikan Liu, Joseph Pasciak, William
Rundell, Yubin Yan and Zhi Zhou. Further several researchers have kindly pro-
vided constructive comments on the book at different stages. In particular, Yubin
Yan and Lele Yuan kindly proofread several chapters of the book and gave many
useful comments, which have helped eliminate many typos.
Last but not least, I would like to take this opportunity to thank my family for
their constant support over the writing of the book, which has taken much of the
time that I should have spent with them.
London, UK Bangti Jin

March 2021
Contents
Part I Preliminaries
1 Continuous Time Random Walk . . . . . . . ....... . . . . . . . . . . . . . . 3
1.1 Random Walk on a Lattice . . . . . . . . ....... . . . . . . . . . . . . . . 3
1.2 Continuous Time Random Walk . . . . ....... . . . . . . . . . . . . . . 6
1.3 Simulating Continuous Time Random Walk . . . . . . . . . . . . . . . . 16
2 Fractional Calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.1 Gamma Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2 Riemann-Liouville Fractional Integral . . . . . . . . . . . . . . . . . . . . . 21
2.3 Fractional Derivatives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.3.1 Riemann-Liouville fractional derivative . . . . . . . . . . . . . . . 28
2.3.2 Djrbashian-Caputo fractional derivative . . . . . . . . . . . . . . . 41
2.3.3 Grünwald-Letnikov fractional derivative . . . . . . . . . . . . . . 51
3 Mittag-Leffler and Wright Functions . . . . . . . . . . . . . . . . . . . . . . . . 59

3.1 Mittag-Leffler Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.1.1 Basic analytic properties . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.1.2 Mittag-Leffler function Efi;1 ðxÞ . . . . . . . . . . . . . . . . . . . . 68
3.2 Wright Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
3.2.1 Basic analytic properties . . . . . . . . . . . . . . . . . . . . . . . . . . 75
3.2.2 Wright function Wq;l ðxÞ . . . . . . . . . . . . . . . . . . . . . . . . 79
3.3 Numerical Algorithms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
3.3.1 Mittag-Leffler function Efi;fl ðzÞ . . . . . . . . . . . . . . . . . . . . . 84
3.3.2 Wright function Wq;l ðxÞ . . . . . . . . . . . . . . . . . . . . . . . . . . 89
ix
x Contents
Part II Fractional Ordinary Differential Equations
4 Cauchy Problem for Fractional ODEs . . . . . . . . . . . . . . . . . . . . . . . 97

4.1 Gronwall’s Inequalities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
4.2 ODEs with a Riemann-Liouville Fractional Derivative . . . . . . . . . 104
4.3 ODEs with a Djrbashian-Caputo Fractional Derivative . . . . . . . . . 113
5 Boundary Value Problem for Fractional ODEs . . . . . . . . . . . . . . . . 137

5.1 Green’s Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
5.1.1 Riemann-Liouville case . . . . . . . . . . . . . . . . . . . . . . . . . . 137
5.1.2 Djrbashian-Caputo case . . . . . . . . . . . . . . . . . . . . . . . . . . 144
5.2 Variational Formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
5.2.1 One-sided fractional derivatives . . . . . . . . . . . . . . . . . . . . 150
5.2.2 Two-sided mixed fractional derivatives . . . . . . . . . . . . . . . 158
5.3 Fractional Sturm-Liouville Problem . . . . . . . . . . . . . . . . . . . . . . . 165
5.3.1 Riemann-Liouville case . . . . . . . . . . . . . . . . . . . . . . . . . . 166
5.3.2 Djrbashian-Caputo case . . . . . . . . . . . . . . . . . . . . . . . . . . 168
Part III Time-Fractional Diffusion
6 Subdiffusion: Hilbert Space Theory . . . . . . . . . . . . . . . . . . . . . . . . . 175

6.1 Existence and Uniqueness in an Abstract Hilbert Space . . . . . . . . 175
6.2 Linear Problems with Time-Independent Coefficients . . . . . . . . . . 185
6.2.1 Solution representation . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
6.2.2 Existence, uniqueness and regularity . . . . . . . . . . . . . . . . . 195
6.3 Linear Problems with Time-Dependent Coefficients . . . . . . . . . . . 208
6.4 Nonlinear Subdiffusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
6.4.1 Lipschitz nonlinearity . . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
6.4.2 Allen-Cahn equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
6.4.3 Compressible Navier-Stokes problem . . . . . . . . . . . . . . . . 229
6.5 Maximum Principles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
6.6 Inverse Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 247
6.6.1 Backward subdiffusion . . . . . . . . . . . . . . . . . . . . . . . . . . . 248
6.6.2 Inverse source problems . . . . . . . . . . . . . . . . . . . . . . . . . . 251
6.6.3 Determining fractional order . . . . . . . . . . . . . . . . . . . . . . . 256
6.6.4 Inverse potential problem . . . . . . . . . . . . . . . . . . . . . . . . . 258
6.7 Numerical Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 267
6.7.1 Convolution quadrature . . . . . . . . . . . . . . . . . . . . . . . . . . 267
6.7.2 Piecewise polynomial interpolation . . . . . . . . . . . . . . . . . . 269
Contents xi
7 Subdiffusion: Hölder Space Theory . . . . . . . . . . . . . . . . . . . . . . . . . 277

7.1 Fundamental Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 277
7.1.1 Fundamental solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . 277
7.1.2 Fractional h-functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 284
7.2 Hölder Regularity in One Dimension . . . . . . . . . . . . . . . . . . . . . . 289
7.2.1 Subdiffusion in R . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 290
7.2.2 Subdiffusion in R þ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 296
7.2.3 Subdiffusion on bounded intervals . . . . . . . . . . . . . . . . . . 302
7.3 Hölder Regularity in Multi-Dimension . . . . . . . . . . . . . . . . . . . . . 306
7.3.1 Subdiffusion in Rd . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 306
7.3.2 Subdiffusion in Rdþ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 311
7.3.3 Subdiffusion on bounded domains . . . . . . . . . . . . . . . . . . . 325
A Mathematical Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 331

A.1 AC Spaces and Hölder Spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . 331
A.1.1 AC spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 331
A.1.2 Hölder spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 332
A.2 Sobolev Spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 333
A.2.1 Lebesgue spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 333
A.2.2 Sobolev spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 334
A.2.3 Fractional Sobolev spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . 336
A.2.4 H_ ðXÞ spaces . . . . . . . . . . . .
s
. . . . . . . . . . . . . . . . . . . . . . . . . 337
A.2.5 Bochner spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 339
A.3 Integral Transforms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 340
A.3.1 Laplace transform . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 340
A.3.2 Fourier transform. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 341
A.4 Fixed Point Theorems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 342
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 345
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 367
Acronyms
CðzÞ Gamma function

C;u contour in complex plane
hfi , hfi fractional h functions
AC n ðDÞ space of absolutely continuous functions
Bðfi; flÞ Beta function
Cþ the set fz 2 C : Rz [ 0g
C the set fz 2 C : RðzÞ\0g
Cðfi; flÞ binomial coefficient
C k;c ðXÞ Hölder spaces
C fi left-sided Djrbashian-Caputo fractional derivative of order fi
a Dx
C fi right-sided Djrbashian-Caputo fractional derivative of order fi
x Db
C fi regularized left-sided Djrbashian-Caputo fractional derivative of
a Dx
order fi
R fi left-sided Riemann-Liouville fractional derivative of order fi
a Dx
R fi right-sided Riemann-Liouville fractional derivative of order fi
x Db
GL fi left-sided Grunwald-Letnikov fractional derivative
a Dx
E fi;fl ðzÞ two-parameter Mittag-Leffler function
E expectation
2F 1 ða; b; c; xÞ hypergeometric function
F Fourier transform
Gfi , Gfi fundamental solutions for subdiffusion
H_ s ðXÞ (fractional order) Sobolev space
H s ðXÞ Sobolev spaces
fi
aI x left-sided Riemann-Liouville fractional integral of order fi
fi
xIb right-sided Riemann-Liouville fractional integral of order fi
L Laplace transform
Ml Mainardi function
N the set of natural numbers
N0 the set of nonnegative integers, i.e., N ¼ N [ f0g
xiii
xiv Acronyms
Rþ the set ð0; 1Þ

R the set ð1; 0Þ
W q;l ðzÞ Wright function
W k;p ðXÞ Sobolev spaces
@ fit left-sided Djrbashian-Caputo fractional derivative of order fi
(in time)
R fi left-sided Riemann-Liouville fractional derivative of order fi
@t
(in time)
ae almost everywhere
erf error function
erfc complementary error function
BDF backward differentiation formula
BVP boundary value problem
CQ convolution quadrature
CTRW continuous time random walk
FDE fractional differential equation
ODE ordinary differential equation
PDE partial differential equation
PDF probability density function
Part I
Preliminaries
Chapter 1
Continuous Time Random Walk
The classical diffusion equation
∂t u(x, t) − κΔu(x, t) = 0 (1.1)
is widely used to describe transport phenomena observed in nature, where u denotes

the concentration of a substance, κ is the diffusion coefficient, and (1.1) describes how
the concentration evolves over time. The model (1.1) can be derived from a purely
macroscopic argument, based on conservation of mass, i.e., ∂t u + ∇ · J = 0 (J is the
flux), and Fick’s first law of diffusion, i.e., J = −κ∇u. Alternatively, following the
work of Albert Einstein in 1905 [Ein05], one can also derive it from the underlying
stochastic process, under a Brownian motion assumption on the particle movement.
Then the probability density function (pdf) p(x, t) of the particle satisfies (1.1). In
this chapter, with continuous time random walk, we show that the pdf p(x, t) of the
particle satisfies fractional differential equations (fdes) of the form
C α
0 Dt p(x, t) = κα ∂xx
2
p(x, t),
κμ μ μ
∂t p(x, t) = 2 (1 − β)−∞RDx p(x, t) + (1 + β)RxD∞ p(x, t) ,
α μ
where 0 < α < 1, 1 < μ < 2, −1 ≤ β ≤ 1, and the notation C 0 Dt and −∞D x ,
R
etc. denote Djrbashian-Caputo / Riemann-Liouville fractional derivatives, etc. in

time / space; see Definitions 2.2 and 2.3 in Chapter 2 for the precise definitions.
In physics, these models are used to describe so-called fractional kinetics [MK00,
MK04, ZDK15], and will represent the main objects of this book.
1.1 Random Walk on a Lattice
Perhaps the simplest stochastic approach to derive equation (1.1) is to consider the
random walk framework on a lattice, which is also known as a Brownian walk. At
each time step (with a time step size Δt), the walker randomly jumps to one of its
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2021 3
B. Jin, Fractional Differential Equations, Applied Mathematical Sciences 206,
https://doi.org/10.1007/978-3-030-76043-4_1
4 1 Continuous Time Random Walk
four nearest neighboring sites on a square lattice (with a space step Δx); see Fig. 1.1
for a schematic illustration.
Fig. 1.1 Random walk on a square lattice, starting from the circle.
In the one-dimensional case, such a process can be modeled by
p j (t + Δt) = 12 p j−1 (t) + 12 p j+1 (t),
where the index j denotes the position on the lattice (grid), and p j (t) the probability
that the walker is at grid j at time t. It relates the probability of being at position j at
time t + Δt to that of the two adjacent sites j ± 1 at time t. The factor 12 expresses the
directional isotropy of the jumps: the jumps to the left and right are equally likely.
A rearrangement gives
p j (t + Δt) − p j (t) 1 p j−1 (t) − 2p j (t) + p j+1 (t)
= .
(Δx)2 2 (Δx)2
The right-hand side is the central finite difference approximation of the second-order
2 p(x, t) at grid j. In the limit Δt, Δx → 0+ , if the pdf p(x, t) is smooth
derivative ∂xx
in x and t, discarding the higher order terms in Δt and Δx in the following Taylor
expansions
p(x, t + Δt) = p(x, t) + Δt∂t p(x, t) + O((Δt)2 ),

(Δx)2 2
p(x ± Δx, t) = p(x, t) ± Δx∂x p(x, t) + 2 ∂xx p(x, t) + O((Δx)3 )
leads to
∂t p(x, t) = κ∂xx
2
p(x, t). (1.2)
The limit is taken such that the quotient
κ= lim (2Δt)−1 (Δx)2 (1.3)

Δx→0+,Δt→0+
is a positive constant, and κ is called the diffusion coefficient, which connects the
spatial and time scales.
1.1 Random Walk on a Lattice 5
Equation (1.2) can also be regarded as a consequence of the central limit theorem.
Suppose that the jump length Δx has a pdf given by λ(x) so that for a < b,
∫b
P(a < Δx < b) = a λ(x)dx. Then Fourier transform gives
∫ ∞ ∫ ∞
=
λ(ξ) e−iξ x λ(x)dx = 1 − iξ x − 12 ξ 2 x 2 + . . . λ(x)dx
−∞ −∞
= 1 − iξ μ1 − 12 ξ 2 μ2 + . . . ,
∫∞
where μ j is the jth moment μ j = −∞ x j λ(x)dx, provided that these moments do
exist, which holds if λ(x) decays sufficiently fast to zero as x → ±∞. Further, assume
that the pdf λ(x) is normalized and even, i.e., μ1 = 0, μ2 = 1 and μ3 = 0. Then
= 1 − 1 ξ 2 + O(ξ 4 ).
λ(ξ) 2
Denote by Δxi the jump length taken at the ith step. In the random walk model,
the steps Δx1 , Δx2 , . . . are independent. The sum of the independent and identically
distributed (i.i.d.) random variables Δxn (according to the pdf with a rescaled λ(x))
xn = Δx1 + Δx2 + . . . + Δxn
gives the position of the walker after n steps. It is also a random variable. Now we
recall a standard result on the sum of two independent random variables, where ∗
denotes convolution of f and g.
Theorem 1.1 If X and Y are independent random variables with a pdf given by f
and g, respectively, then the sum Z = X + Y has a pdf f ∗ g.
n,
By Theorem 1.1, the random variable Xn has a Fourier transform pn (ξ) = (λ(ξ))
and the normalized sum n− 2 xn has the Fourier transform
1
1 n
pn n− 2 ξ = 1 − 2n ξ + O(n−2 ) .
1 2
ξ2
Taking the limit as n → ∞ gives p(ξ) = e− 2 and inverting the Fourier transform
x2
gives the standard Gaussian distribution p(x) = (2π)− 2 e− 2 , cf. Example A.4 in the
1
appendix. This is precisely the central limit theorem, asserting that the long-term
average behavior of i.i.d. random variables is Gaussian. One requirement for the
whole procedure to work is the finiteness of the second moment μ2 of λ(x).
Now we interpret xn as the particle position after n time steps, by correlating the
time step size Δt with the variance of Δx according to the ansatz (1.3). This can be
easily achieved by rescaling the variance of λ(x) to 2κt. Then by the scaling rule for
the Fourier transform, pn (ξ) is given by
pn (n− 2 ξ) = (1 − n−1 κtξ 2 + O(n−2 ))n

1
and taking limit as n → ∞ gives p(ξ) = e−ξ κt . By inverse Fourier transform, the
2
pdf being at certain position x and time t is governed by the standard diffusion model
(1.2), and is given by the Gaussian

x2
p(x, t) = (4πκt)− 2 e− 4κ t .
1
It is the fundamental solution, i.e., the solution p(x, t) of (1.2) with the initial condi-
tion p(x, 0) = δ(x), the Dirac function concentrated at x = 0. Note that at any fixed
time∫t > 0, p(x, t) is a Gaussian distribution in x, with mean zero and variance 2κt,
∞
i.e., −∞ x 2 p(x, t)dx = 2κt, which scales linearly with the time t. This represents one
distinct feature of normal diffusion processes.
1.2 Continuous Time Random Walk
There are several stochastic models for deriving differential equations involving a
fractional-order derivative. We describe the continuous time random walk framework
(ctrw) due to Montroll and Weiss [MW65]. ctrw generalizes the random walk
model in Section 1.1, in which the length of each jump and the time elapsed between
two successive jumps follow a given pdf. We assume that these two random variables
are independent, though in theory one can allow correlation in order to achieve more
flexible modeling. In one spatial dimension, the picture is as follows: a walker moves
along the x-axis, starting at a position x0 at time t0 = 0. At time t1 , the walker jumps
to x1 , then at time t2 jumps to x2 , and so on. We assume that the temporal and spatial
increments Δtn = tn −tn−1 and Δxn = xn − xn−1 are i.i.d. random variables, following
pdfs ψ(t) and λ(x), respectively, known as the waiting time distribution and jump
length distribution, respectively. Namely, the probability of Δtn lying in any interval
[a, b] ⊂ R+ and Δxn lying in any interval [c, d] ⊂ R are given by
∫ b ∫ d
P(a < Δtn < b) = ψ(t)dt and P(c < Δxn < d) = λ(x)dx,
a c
respectively. Now the goal is to determine the probability that the walker lies in a
given spatial interval at time t. For given pdfs ψ and λ, the position x of the walker
can be regarded as a step function of t.
Example 1.1 Suppose that the distribution ψ(t) is exponential with a parameter τ > 0,
t
i.e., ψ(t) = τ −1 e− τ for t ∈ R+ , and the jump length distribution λ(x) is Gaussian
x2
−
with mean zero and variance σ 2 , i.e., λ(x) = (2πσ 2 )− 2 e
1
2σ 2 for x ∈ R. Then the
waiting time Δtn and jump length Δxn satisfy
E[Δtn ] = τ, E[Δtn2 ] = τ 2, E[Δxn ] = 0, E[Δxn2 ] = σ 2,
where E denotes taking expectation with respect to the underlying distribution. The
position x(t) of the walker is a step function (in time t). Two sample trajectories of
ctrw with τ = 1 and σ = 1 are given in Fig. 1.2.
1.2 Continuous Time Random Walk 7
20 20
0 0
x
x
-20 -20
0 20 40 60 80 100 0 20 40 60 80 100
t t
Fig. 1.2 Two realizations of 1D ctrw, with exponential waiting time distribution (with τ = 1) and
Gaussian jump length distribution (with σ = 1), starting from x0 = 0.
Now we derive the pdf of the total waiting time tn after n steps and the total jump
length xn − x0 after n steps. The main tools in the derivation are Laplace and Fourier
transforms in Appendices A.3.1 and A.3.2, respectively. We denote by ψn (t) the pdf
of tn = Δt1 + Δt2 + . . . + Δtn , and ψ1 = ψ. By Theorem 1.1, we have
∫ t
ψn (t) = (ψn−1 ∗ ψ)(t) = ψn−1 (s)ψ(t − s)ds.
0
Then the characteristic

∫∞ function ψn (z) of ψn (t) (i.e., its Laplace transform) ψn (z) =
L[ψn ](z) = 0 e−zt ψn (t)dt, by the convolution rule for Laplace transform in (A.8),
is given by ψn (z) = (ψ(z))n . Let Ψ(t) denote the survival probability, i.e., the
probability of the walker not jumping within a time t (or equivalently, the probability
of remaining stationary for at least a duration t). Then
∫ ∞ ∫ t
Ψ(t) = ψ(s)ds = 1 − ψ(s)ds, 0 < t < ∞.
t 0
The characteristic function for the survival probability Ψ is
Ψ(z) = z−1 − z −1 ψ(z).
Then the probability of taking exactly n steps up to time t is given by

∫ t
χn (t) = ψn (s)Ψ(t − s)ds = (ψn ∗ Ψ)(t).
0
The convolution rule for Laplace transform yields
χn (z) = ψn (z)Ψ(z) = ψ(z)n z −1 (1 − ψ(z)).
Next we derive the pdf of the jump length xn − x0 after n steps. We denote by
λn (x) the pdf of the random variable xn − x0 = Δx1 + Δx2 + . . . + Δxn . Appealing
to Theorem 1.1 again yields
∫ ∞
λn (x) = (λn−1 ∗ λ)(x) = λn−1 (y)λ(x − y)dy,
−∞
n (ξ) is given by
and by the convolution rule for Fourier transform, λ
λ n,
n (ξ) = (λ(ξ)) n ≥ 0.
We denote by p(x, t) the pdf of the walker at the position x at time t. Since χn (t)
is the probability of taking n steps up to time t,
∞
p(x, t) = λn (x) χn (t).
n=0
We denote the Fourier-Laplace transform of p(x, t) by

∫ ∞ ∫ ∞
p(ξ, z) = LF [p](ξ, z) = e−zt e−iξ x p(x, t)dxdt.
0 −∞
Then p(ξ, z) is given by

∞ ∞
p(ξ, z) = n (ξ) χn (z) = 1 − ψ(z)
λ ψ(z)]n .
[λ(ξ) (1.4)
n=0
z n=0
Since λ and ψ are pdfs,

∫ ∞ ∫ ∞
=
λ(0) λ(x)dx = 1 and ψ(0) = ψ(t)dt = 1. (1.5)
−∞ 0
ψ(z)| < 1 so the series in (1.4) is absolutely convergent:

If ξ 0 or z > 0, then | λ(ξ)
∞
ψ(z)]n = (1 − λ(ξ)
[λ(ξ) ψ(z))−1,
n=0
and hence we obtain the following fundamental relation for the pdf p(x, t) in the
Laplace-Fourier domain
1 − ψ(z) 1
p(ξ, z) = . (1.6)
z ψ(z)
1 − λ(ξ)
t
Example 1.2 For the waiting time distribution ψ(t) = τ −1 e− τ and the jump length
2
− x2
distribution λ(x) = (2πσ 2 )− 2 e in Example 1.1, we have ψ(z) = (1 + τz)−1 and
1
2σ
σ2 ξ 2
= e− 2 . Thus, the Fourier-Laplace transform p(ξ, z) of the pdf p(x, t) of the
λ(ξ)
walker at position x (relative to x0 ) at time t is given by
σ2 ξ 2
p(ξ, z) = τ(1 + τz − e− 2 )−1 .
Different types of ctrw processes can be classified according to the characteristic

waiting time and jump length variance ς 2 defined by
∫ ∞ ∫ ∞
=: E[Δtn ] = tψ(t)dt and ς 2 =: E[(Δxn )2 ] = x 2 λ(x)dx,
0 −∞
being finite or diverging. Below we shall discuss the following three scenarios: (i)
Finite and ς 2 , (ii) Diverging and finite ς 2 and (iii) Diverging ς 2 and finite .
(i) Finite characteristic waiting time and jump length variance. When both
characteristic waiting time and jump length variance ς are finite, ctrw does not
lead to anything new. Specifically, assume that the pdfs ψ(t) and λ(x) are normalized:
∫ ∞ ∫ ∞ ∫ ∞
tψ(t)dt = 1, xλ(x)dx = 0, x 2 λ(x)dx = 1. (1.7)
0 −∞ −∞
These conditions are satisfied after suitable rescaling if the waiting time pdf ψ(t) has
a finite mean, and the jump length pdf λ(x) has a finite variance. Recall the relation

(1.5): ψ(0) = 1 = λ(0). Since
∫ ∞
dk ψ
= e−zt (−t)k ψ(t)dt,
dz k 0
using the normalization condition (1.7), we have

∫ ∞
ψ (0) = − tψ(t)dt = −1.
0
Similarly, we deduce
∫ ∞ ∫ ∞

λ (0) = −i xλ(x)dx = 0, (0) = −
λ x 2 λ(x)dx = −1.
−∞ −∞
Next, for τ > 0 and σ > 0, let the random variables Δtn and Δxn follow the
rescaled pdfs
ψτ (t) = τ −1 ψ(τ −1 t) and λσ (x) = σ −1 λ(σ −1 x). (1.8)
Now consider the limit of p(ξ, z; σ, τ) as τ, σ → 0+ . Then simple computation shows

=: E[Δtn ] = τ, E[Δxn ] = 0, and ς 2 =: E[Δxn2 ] = σ 2 . By the scaling rules for
Fourier and Laplace transforms,
ψτ (z) = ψ(τz) and σ (ξ) = λ(σξ).

λ
Consequently,
1 − ψ(τz) 1
p(ξ, z; σ, τ) = .
z
1 − ψ(τz)λ(σξ)
The Taylor expansion of ψ(z) around z = 0 yields
ψ(z) = ψ(0) + ψ (0)z + . . . = 1 − z + O(z 2 ) as z → 0.

(0) = 0
Further we assume that the pdf λ(x) is even, i.e., λ(−x) = λ(x). Then λ

and the Taylor expansion of λ(ξ) around ξ = 0 is given by
= λ(0)
λ(ξ) +λ (0)ξ 2 + . . . = 1 − 1 ξ 2 + O(ξ 4 ).
(0)ξ + 1 λ
2 2
Using the last two relations, simple algebraic manipulation gives

τ 1 + O(τz)
p(ξ, z; σ, τ) = . (1.9)
τz + 2σ ξ
1 2 2 1 + O(τz + σ 2 ξ 2 )
Now in the formula (1.9), we let σ, τ → 0+ while keeping the relation

limσ→0+,τ→0+ (2τ)−1 σ 2 = κ, for some κ > 0, and obtain
τ 1
p(ξ, z) = lim = . (1.10)
τz + 12 σ 2 ξ 2 z + κξ 2
The scaling relation here for κ is identical to that in the random walk framework, cf.
(1.3). Upon inverting the Laplace transform, we find
∫
1 ezt
dz = e−κ ξ t ,
2
p(ξ, t) =
2πi C z + κξ 2
and then inverting the Fourier transform gives the familiar Gaussian pdf p(x, t) =
x2
(4πκt)− 2 e− 4κ t . Last we verify that the pdf p(x, t) satisfies (1.2):
1
LF [∂t p − κ∂xx
2
p](ξ, z) = z p(ξ, z) − p(ξ, 0) + κξ 2 p(ξ, z)
= (z + κξ 2 ) p(ξ, z) − p(ξ, 0) = 1 − p(ξ, 0) = 0,
where the last step follows by p̃(ξ, 0) = δ̃(ξ) = 1. Hence, the pdf p(x, t) for x(t) − x0
at the position x (relative to x0 ) of the walker at time t indeed satisfies (1.2). The
ctrw framework recovers the classical diffusion model, as long as the waiting time
pdf ψ(t) has a finite mean and the jump length pdf λ(x) has finite first and second
moments. Further, the displacement variance is given by
∫ ∞
μ2 (t) = E p(t) [(x(t) − x0 )2 ] = x 2 p(x, t)dx = 2κt,
−∞
which grows linearly with the time t. One striking thing is that any pair of pdfs with
finite characteristic waiting time and jump length variance ς 2 lead to the same
result, to the lowest order.
(ii) Divergent mean waiting time. Now we consider the situation where the charac-
teristic waiting time diverges, but the jump length variance ς 2 is finite. This occurs
for example when the particle might be trapped in a potential well, and it takes a
long time for the particle to leave the well. To model such phenomena, we employ
a heavy-tailed (i.e., the tail is not exponentially bounded) waiting time pdf with the
asymptotic behavior
ψ(t) ∼ at −1−α, as t → ∞ (1.11)
for some α ∈ (0, 1), where a > 0 is a constant. One example of such pdf is
ψ(t) = α(1 + t)−1−α . The (asymptotic) power law decay in (1.11) is heavy tailed and
allows occasional very large waiting time between consecutive jumps. Again, the
specific form of ψ(t) is irrelevant, and only the power law decay at large time matters.
The parameter α determines the asymptotic decay of the pdf ψ(t): the closer is α to
zero, the slower the decay is and more likely the long waiting time is.
Now the question is whether such occasional long waiting time will change
completely the dynamics of the stochastic process∫in the large time. For the power
∞
law decay, the mean waiting time is divergent, i.e., 0 tψ(t)dt = +∞. So one cannot
assume the normalizing condition on the pdf ψ in (1.7), and the above analysis
breaks
∫ ∞down. Nonetheless,∫ in this part, the assumption on λ(x) remains unchanged,
∞
i.e., −∞ xλ(x)dx = 0 and −∞ x 2 λ(x)dx = 1.
15
10
5
x
-5
0 1000 2000 3000
t
(a) (b)
Fig. 1.3 ctrw with a power law waiting time pdf. Panel (a) gives one sample trajectory with
waiting time pdf ψ(t) = α(1 + t)−1−α , α = 34 , and standard Gaussian jump length pdf; Panel (b) is
a zoom-in of the first cluster of Panel (a).
In Fig. 1.3, we show one sample trajectory of ctrw with a power law waiting time
pdf ψ(t) = α(1 + t)−1−α , α = 34 , and the standard Gaussian jump length pdf, and an
enlarged version of the initial time interval. One observes clearly the occasional but
large waiting time appearing at different time scales, cf. Fig. 1.3(b). This behavior
is dramatically different from the case of finite mean waiting time in Fig. 1.2.
The following result bounds the Laplace transform ψ(z) for z → 0. The notation
Γ(z) denotes the Gamma function; see (2.1) in Section 2.1 for its definition.
Theorem 1.2 Let α ∈ (0, 1), and a > 0. If ψ(t) = at −1−α + O(t −2−α ) as t → ∞, then
ψ(z) = 1 − Bα z α + O(z) as z → 0
with Bα = aα−1 Γ(1 − α).

Proof By the assumption on ψ, there exists t0 > 0 such that
|t 1+α ψ(t) − a| ≤ ct −1, ∀t0 ≤ t < ∞ (1.12)
and we consider 0 < z < t0−1 so that zt0 < 1. Since ψ(0) = 1, we have
∫ ∞ ∫ t0
1 − ψ(z) = (1 − e−zt )ψ(t)dt = (1 − e−zt )(ψ(t) − at −1−α )dt
0 0
∫ ∞ ∫ ∞ 3
−zt −1−α
+ (1 − e )(ψ(t) − at )dt + (1 − e−zt )at −1−α dt =: Ii .
t0 0 i=1
Since 0 ≤ 1 − e−s ≤ s for all s ∈ [0, 1], the term I1 can be bounded by
∫ t0 ∫ t0
|I1 | ≤ ztψ(t)dt + ztat −1−α dt
0 0
∫ ∞ ∫ t0
≤ zt0 ψ(t)dt + az t −α dt ≤ cα,t0 z.
0 0
For the term I2 , using (1.12) and changing variables s = zt, we deduce
∫ ∞
|I2 | ≤ (1 − e−zt )ct −2−α dt
t0
∫ 1 ∫ ∞
≤ cz 1+α s−1−α ds + s−2−α ds
t0 z 1
−1 −α
≤ cz 1+α
(α ((t0 z) − 1) + (1 + α)−1 ) ≤ cα,t0 z.
Last, for the term I3 , by changing variables s = zt, we have

∫ ∞ ∫ ∞
−zt −1−α α 1 − e−s
I3 = (1 − e )at dt = az ds.
0 0 s1+α
Integration by parts and the identity (2.2) below give
∫ ∞ ∫ ∞
1 − e−s −1
ds = α e−s s−α ds = α−1 Γ(1 − α).
0 s1+α 0
Combining the preceding estimates yields the desired inequality.
Now we introduce the following rescaled pdfs for the incremental waiting time
Δtn and jump length Δxn : ψτ (t) = τ −1 ψ(τ −1 t) and λσ (x) = σ −1 λ(σ −1 x). By (1.6),
p(ξ, z; σ, τ) is given by
1 − ψ(τz) 1
p(ξ, z; σ, τ) = .
z
1 − ψ(τz)λ(σξ)
Under the power law waiting time pdf, from Theorem 1.2 we deduce ψ(z) = 1 −
= 1 − 1 ξ 2 + O(ξ 4 ) as ξ → 0. Hence,
Bα z α + O(z) as z → 0 and λ(ξ) 2
Bα τ α z α−1 1 + O(τ 1−α z 1−α )

p(ξ, z; σ, τ) = . (1.13)
Bα τ α z α + 12 σ 2 ξ 2 1 + O(τ 1−α z 1−α + τ α z α + σξ)
Once again, in (1.13), we find the Fourier-Laplace transform p(ξ, z) by letting

σ, τ → 0+ , while keeping limσ→0+,τ→0+ (2Bα τ α )−1 σ 2 = κα for some κα > 0. Thus
we obtain
Bα τ α z α−1 z α−1
p(ξ, z) = lim p(ξ, z; σ, τ) = lim = . (1.14)
Bα τ α z α + 12 σ 2 ξ 2 z α + κα ξ 2
Formally, we recover the formula (1.10) by setting α = 1.

Last, we invert the Fourier-Laplace transform p(ξ, z) using the Mittag-Leffler
function Eα,1 (z) (cf. (3.1) in Chapter 3) and Wright function Wρ,μ (z) (cf. (3.24)
in Chapter 3). Using the Laplace transform formula of the Mittag-Leffler function
Eα,1 (z) in Lemma 3.2 below, we deduce
p(ξ, t) = Eα,1 (−κα t α ξ 2 ),
and then by the Fourier transform of the Wright function Wρ,μ (z) in Proposition 3.4,
we obtain an explicit expression of the pdf p(x, t)
p(x, t) = (4κα t α )− 2 W− α2 ,1− α2 (−(κα t α )− 2 |x|).

1 1
With α = 1, this formula recovers the familiar Gaussian density. Then Lemma 2.9
implies that the pdf p(x, t) satisfies
C α
LF 0 Dt p(x, t) − κα ∂xx
2
p(x, t) = z α p(ξ, z) − z α−1 p(ξ, 0) + κα ξ 2 p(ξ, z)
= (z α + κα ξ 2 ) p(ξ, z) − z α−1 p(ξ, 0) ≡ 0,
since p(ξ, 0) =
δ(ξ) = 1. Thus the pdf p(x, t) satisfies the following time-fractional
diffusion equation with a Djrbashian-Caputo fractional derivative C α
0 Dt p in time t
(see Definition 2.3 in Chapter 2 for the precise definition)
C α
0 Dt p(x, t) − κα ∂xx
2
p(x, t) = 0, 0 < t < ∞, −∞ < x < ∞.
It can be rewritten using a Riemann-Liouville fractional derivative 0RDt1−α p as
∂t p(x, t) = κα 0RDt1−α ∂xx

2
p(x, t).
In this book, we focus on the approach based on the Djrbashian-Caputo

∫∞ one.
Now we compute the mean square displacement μ2 (t) = −∞ x 2 p(x, t) dx for
t > 0. To derive an explicit formula, we resort to the Laplace transform
∫ ∞
d2
μ2 (z) = x 2 p(x, z)dx = − p(ξ, z)|ξ=0
−∞ dξ 2
d2
=− (z + κα z 1−α ξ 2 )−1 |ξ=0 = 2κα z −1−α,
dξ 2
which upon inverse Laplace transform yields μ2 (t) = 2Γ(1 + α)−1 κα t α . Thus the
mean square displacement grows only sublinearly with the time t, which, at large
time t, is slower than that in normal diffusion, and as α → 1− , it recovers the formula
for normal diffusion. Such a diffusion process is called subdiffusion or anomalously
slow diffusion in the literature. Subdiffusion has been observed in a large number
of physical applications, e.g., column experiments [HH98], thermal diffusion in
fractal domains [Nig86], non-Fickian transport in geological formations [BCDS06],
and protein transport within membranes [Kou08]. The mathematical theory for the
subdiffusion model will be discussed in detail in Chapters 6 and 7.
(iii) Diverging jump length variance. Last we turn to the case of diverging jump
length variance ς 2 , but finite characteristic waiting time . To be specific, we assume
an exponential waiting time ψ(t), and that the jump length follows a (possibly asym-
metric) Lévy distribution. See the monograph [Nol20] for an extensive treatment of
univariate stable distributions and [GM98] for the use of stable distributions within
fractional calculus.
The most convenient way to define Lévy random variables is via their character-
istic function (Fourier transform)
−|ξ | μ (1 − iβ sign(ξ) tan μπ

2 ), μ 1,
μ, β) =
ln λ(ξ;
−|ξ |(1 + iβ sign(ξ) π2 ln |ξ |), μ = 1,
where μ ∈ (0, 2] and β ∈ [−1, 1]. The parameter μ determines the rate at which the
tails of the pdf taper off. When μ = 2, it recovers a Gaussian distribution, irrespective
of the value of β, whereas with μ = 1 and β = 0, the stable density is identical to
the Cauchy distribution, i.e., λ(x) = (π(1 + x 2 ))−1 . The parameter β determines the
degree of asymmetry of the distribution. When β is negative (respectively positive),
the distribution is skewed to the left (respectively right). When β = 0, the expression
μ) = e− |ξ | μ .
simplifies to λ(ξ;
Generally, the pdf of a Lévy stable distribution is not available in closed form.
However, it is possible to compute the density numerically [Nol97] or to simulate
random variables from such stable distributions, cf. Section 1.3. Nonetheless, fol-
lowing the proof of Theorem 1.2, when 1 < μ < 2, one can obtain an inverse power
law asymptotic [Nol20, Theorem 1.2, p. 13]
λ(x) ∼ aμ,β |x| −1−μ as |x| → ∞ (1.15)
for some aμ,β > 0, from which it follows directly that the jump length variance
diverges, i.e., ∫ ∞
x 2 λ(x) dx = ∞. (1.16)
−∞
In general, the pth moment of a stable random variable is finite if and only if p < μ.
50
40
30
x
20
10
0
0 50 100 150 200 250
t
(a) (b)
Fig. 1.4 ctrw with a Lévy jump length pdf. Panel (a) is one sample trajectory with an exponential
waiting time pdf ψ(t) = e−t , and a Lévy jump length pdf with μ = 32 and β = 0. Panel (b) is a
zoom-in of Panel (a).
For ctrw, one especially relevant case is λ(x; μ, β) with μ ∈ (1, 2] for the jump
length distribution. In Fig. 1.4, we show one sample trajectory of ctrw with a
Lévy jump length distribution. Due to the asymptotic property (1.15) of the jump
length pdf, very long jumps may occur with a much higher probability than for an
exponentially decaying pdf like a Gaussian. The scaling nature of the Lévy jump
length pdf leads to clusters along the trajectory, i.e., local motion is occasionally
interrupted by long sojourns, at all length scales.
Below we derive the diffusion limit using the scaling technique, by appealing
to the scaled pdfs ψ∫τ (t) and λσ (x), for some τ, σ > 0. The pdf ψ is assumed to
∞
be normalized, i.e., 0 tψ(t)dt = 1. Using Taylor expansion, for small ξ, we have
(μ 1)
μπ
e−λ(ξ;μ,β) = 1 − |ξ | μ 1 − iβ sign(ξ) tan + O(|ξ | 2μ ).
2
Hence, for μ 1, we have
ψτ (z) = 1 − τz + O(τ 2 z 2 ) as τ → 0,

σ (ξ) = 1 − σ μ |ξ | μ 1 − iβ sign(ξ) tan μπ + O(σ 2μ |ξ | 2μ ) as σ → 0.
λ
2
These two identities and Fourier-Laplace convolution formulas yield
1
p(ξ, z) = lim p(ξ, z; σ, τ) = , (1.17)
z + κμ |ξ | μ 1 − iβ sign(ξ) tan μπ
2
where the limit is taken under the condition κμ = limσ→0+,τ→0+ τ −1 σ μ for some
κμ > 0, which represents the diffusion coefficient.
Next we invert the Fourier transform of p(ξ, z) and obtain
μ (1−iβ sign(ξ) tan μ π

p(ξ, t) = e−κμ |ξ | 2 )t
,
which is the characteristic function of a Lévy stable distribution, for which an explicit
expression is generally unavailable. This can be regarded as a generalized central
limit theorem for stable distributions.
One can verify by Fourier and Laplace inversion that the pdf p(x, t) satisfies
κ μ μ
∂t p(x, t) = 2μ (1 − β)−∞RDx p(x, t) + (1 + β)Rx D∞ p(x, t) ,
in view of Theorem 2.8 below. In the symmetric case β = 0, it simplifies to

μ
∂t p(x, t) = −κμ (−Δ)R2 p(x, t),
μ
where (−Δ)R2 denotes the Riesz fractional operator
μ
μ μ
−(−Δ)R2 f = 12 ( −∞RDx f + R
x D∞ f ),
μ μ
where the notation −∞RDx and R x D∞ denote the left-sided and right-sided Riemann-
Liouville fractional derivatives, cf. Definition 2.2 in Chapter 2 for the precise def-
initions. This operator can be viewed as the one-dimensional version of fractional
Laplacian (see the survey [Gar19] for a nice introduction to fractional Laplacian).
It recovers the Gaussian density as μ → 2. Using (1.15), we deduce the following
power law asymptotic as |x| → ∞
p(x, t) ∼ κμ t|x| −1−μ, μ < 2.

∫∞
Due to this property, the mean squared displacement diverges −∞ |x| 2 p(x, t) dx =
∞. That is, the particle following such a random walk model spreads faster than
normal diffusion, and it is commonly known as superdiffusion in the literature. It
has been proposed for several applications, including solute transport [BWM00]
and searching patterns of animals (e.g., albatross, bacteria, plankton, jackals, spider
monkey) [RFMM+04]. The divergent mean squares displacement has caused some
controversy in practice, and various modifications have been proposed. Boundary
value problems (bvps) for superdiffusion are involved, since the long jumps make the
very definition of a proper boundary condition actually quite intricate. In Chapter 5,
we present mathematical theory for 1D stationary superdiffusion models.
1.3 Simulating Continuous Time Random Walk
To simulate ctrw, one needs to generate random numbers for a given pdf, e.g.,
power laws and Lévy stable distributions. The task of generating random numbers
for an arbitrary pdf is often highly nontrivial, especially in high-dimensional spaces.
There are a number of possible methods, and we describe only the transformation
1.3 Simulating Continuous Time Random Walk 17
method, which is simple and easy to implement. For more advanced methods, we
refer to the monograph [Liu08].
The transformation method requires that we have access to random numbers
uniformly distributed in the unit interval 0 ≤ r ≤ 1, which can be generated by any
standard pseudo-random number generators. Now suppose that p(x) is a continuous
pdf from which we wish to draw samples. The pdfs p(x) and p(r) (the uniform
distribution on [0, 1]) are related by
dr dr
p(x) = p(r) = ,
dx dx
where the second equality follows because p(r) = 1 over the interval [0, 1]. Integrat-
ing both sides with respect to x, we obtain the complementary cumulative distribution
function F(x) in terms of r:
∫ ∞ ∫ 1
F(x) = p(x ) dx = dr = 1 − r,
x r
or equivalent x = F −1 (1 − r), where F −1

denotes the inverse function of the function
F(x). For many pdfs, this function can be evaluated in closed form.
Example 1.3 In this example we illustrate the method on the power law density
ψ(x) = α(1 + x)−1−α . One can compute directly
∫ ∞ ∫ ∞
α 1
F(x) = ψ(x ) dx = dx = ,
x x (1 + x ) 1+α (1 + x)α
and the inverse function F −1 (x) is given by F −1 (x) = x − α − 1. This last formula
1
provides an easy way to generate random variables from the power law density ψ(t),
needed for simulating the subdiffusion model.
Generating samples from the Lévy stable distribution λ(x; μ, β) is more delicate.
The main challenge lies in the fact that there are no analytic expressions for the
inverse F −1 , and thus the transformation method is not easy to apply. Two well-
known exceptions are the Gaussian (μ = 2) and Cauchy (μ = 1, β = 0). One popular
algorithm is from [CMS76, WW95]. For μ ∈ (0, 2] and β ∈ [−1, 1], it generates a
random variable X ∼ λ(x; μ, β) as in Algorithm 1.
Exercises
Exercise 1.1 Show formula (1.13).

Exercise 1.2 Prove the asymptotic (1.15). (Hint: [Nol20, Section 3.6])
Exercise 1.3 Show the jump length variance formula (1.16).

Algorithm 1 Generating stable random variables.

1: Generate a random variable V uniformly distributed on (− π2 , π
2) and an independent exponen-
tial random variable W with mean 1.
2: For μ 1, compute
1−μ
sin(μ(V + Bμ, β )) cos(V − μ(V + Bμ, β )) μ
X = Cμ, β 1
,
(cos V ) μ W
where
arctan(β tan μπ μπ 2μ1
2 )
Bμ, β = and Cμ, β = 1 + β 2 tan2 .
μ 2
3: For μ = 1, compute

2 π W cos V
X= + βV tan V − β ln π .
π 2 2 + βV
Exercise 1.4 Prove that the pth moment of a stable random variable is finite if and
only if p < μ.
Exercise 1.5 Derive the representation (1.17).
Exercise 1.6 Develop an algorithm for generating random variables from the Cauchy
distribution λ(x) = (π(1 + x 2 ))−1 and implement it.
Exercise 1.7 Implement the algorithm for generating random variables from Lévy
stable distribution, with μ = 1.5, and β = 0.
Exercise 1.8 Derive the diffusion limit for the ctrw with a power law waiting time
and Lévy jump length distribution, and the corresponding governing equation.
Chapter 2
Fractional Calculus
In this chapter, we describe basic properties of fractional-order integrals and deriva-

tives, including Riemann-Liouville fractional integral and derivative, Djrbashian-
Caputo and Grünwald-Letnikov fractional derivatives. They will serve as the main
modeling tools in fdes.
2.1 Gamma Function
This function will appear in nearly every formula in this book! The usual way to
define the Gamma function Γ(z) is
∫ ∞
Γ(z) = t z−1 e−t dt, (z) > 0. (2.1)
0
It is analytic in the right half complex plane C+ := {z ∈ C : (z) > 0} and

integration by parts shows the following recursion formula
Γ(z + 1) = zΓ(z). (2.2)
Since Γ(1) = 1, (2.2) implies that for any positive integer n ∈ N, Γ(n + 1) = n!.
The following reflection formula [AS65, 6.1.17, p. 256]
π
Γ(z)Γ(1 − z) = , 0 < (z) < 1, (2.3)
sin(πz)
and Legendre’s duplication formula [AS65, 6.1.18, p. 256]
√
Γ(z)Γ(z + 12 ) = π21−2z Γ(2z) (2.4)
are very useful in practice. For example, one only needs to approximate √ Γ(z) for
(z) ≥ 12 and uses (2.3) for (z) < 12 . It also gives directly Γ( 12 ) = π.
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2021 19

B. Jin, Fractional Differential Equations, Applied Mathematical Sciences 206,
https://doi.org/10.1007/978-3-030-76043-4_2
20 2 Fractional Calculus
It is possible to continue Γ(z) analytically into the left half complex plane C− :=
{z ∈ C : (z) < 0}. This is often done by representing the reciprocal Gamma
1
function Γ(z) as an infinite product (cf. Exercise 2.2)
1 n−z
= lim z(z + 1) . . . (z + n), ∀z ∈ C. (2.5)
Γ(z) n→∞ n!
1
Hence, Γ(z) is an entire function of z with zeros at z = 0, −1, −2, . . ., and accordingly,
Γ(z) has poles at z = 0, −1, . . .. The integral representation for Γ(z)
1
due to Hermann
Hankel [Han64] is very useful. It is often derived using Laplace transform (cf.
Appendix A.3). Substituting t = su in (2.1) yields
∫ ∞
Γ(z)
= uz−1 e−su du,
sz 0
which can be regarded as the Laplace transform of uz−1 for a fixed z ∈ C. Hence
uz−1 can be interpreted as an inverse Laplace transform
∫
z−1 −1 Γ(z) 1 Γ(z)
u =L = esu z ds,
sz 2πi C s
where C is any deformed Bromwich contour that winds around the negative real axis
R− = (−∞, 0) in the anticlockwise sense. Now substituting ζ = su yields
∫
1 1
= ζ −z e ζ dζ . (2.6)
Γ(z) 2πi C
The integrand ζ −z has a branch cut on R− but is analytic elsewhere so that (2.6) is
independent of the contour C. The formula (2.6) is very useful for deriving integral
representations of various special functions, e.g., Mittag-Leffler and Wright functions
(cf. (3.1) and (3.24) in Chapter 3). In practice, C can be deformed to facilitate the
1
numerical evaluation of the integral, which makes Γ(z) actually easier to compute
than Γ(z) itself [ST07].
For large argument, the following Stirling’s formula holds (for any > 0)
Γ(z) = (2π) 2 e−z z z− 2 (1 + O(z −1 )) as |z| → ∞, | arg(z)| < π − .

1 1
(2.7)
Further, for any two real numbers x and s, the following identity holds
lim (x s Γ(x))−1 Γ(x + s) = 1, (2.8)

x→∞
often known as Wendel’s limit [Wen48], due to his elegant proof using only Hölder’s
inequality and the recursion identity (2.2); see Exercise 2.4 for the detail.
A closely related function, the Beta function B(α, β) for α, β > 0 is defined by
∫ 1
B(α, β) = (1 − s)α−1 s β−1 ds. (2.9)
0
2.2 Riemann-Liouville Fractional Integral 21
It can be expressed by the Gamma function Γ(z) as

Γ(α)Γ(β)
B(α, β) = . (2.10)
Γ(α + β)
The (generalized) binomial coefficient C(α, β) is defined by
Γ(α + 1)
C(α, β) = ,
Γ(β + 1)Γ(α − β + 1)
provided β −1, −2, . . . and α − β −1, −2, . . .. For α, β ∈ N0 , it recovers the usual
binomial coefficient. The following identities hold for the binomial coefficients.
Lemma 2.1 Let j, k ∈ N0 , and α, β 0, −1, . . .. Then the following identities hold.
(−1) j C(α, j) = C( j − α − 1, j), (2.11)

C(α, k + j)C(k + j, k) = C(α, k)C(α − k, j), (2.12)

k
C(α, j)C(β, k − j) = C(α + β, k). (2.13)
j=0
Proof The identities (2.11) and (2.12) are direct from the definition. Indeed,
Γ(α + 1) Γ(k + j + 1)
C(α, k + j)C(k + j, k) =
Γ(k + j + 1)Γ(α − k − j + 1) Γ(k + 1)Γ( j + 1)
Γ(α + 1) Γ(α − k + 1)
=
Γ(k + 1)Γ(α − k + 1) Γ( j + 1)Γ(α − k − j + 1)
= C(α, k)C(α − k, j),
showing (2.12). The formula (2.13) can be found at [PBM86, p. 616, 4.2.15.13].
2.2 Riemann-Liouville Fractional Integral
The Riemann-Liouville fractional integral generalizes the well-known Cauchy’s

iterated integral formula. Let finite a, b ∈ R, with a < b, and we denote by D = (a, b).
Then for any integer k ∈ N, the k-fold integral a Ixk , based at the left end point x = a,
is defined recursively by a Ix0 f (x) = f (x) and
∫ x
k k−1
a I x f (x) = a Is f (s) ds, k = 1, 2, . . . .
a
We claim the following Cauchy’s iterated integral formula for k ≥ 1

∫ x
k 1
a I x f (x) = (x − s)k−1 f (s) ds (2.14)
(k − 1)! a
and verify it by mathematical induction. Clearly, it holds for k = 1. Now assume that
the formula (2.14) holds for some k ≥ 1, i.e.,
∫ x
k 1
a I x f (x) = (x − s)k−1 f (s) ds.
(k − 1)! a
Then by definition and changing integration order, we deduce
∫ x ∫ x∫ s
k+1 k 1
a Ix f (x) = a Is f (s) ds = (s − t)k−1 f (t) dtds
a (k − 1)! a a
∫ x∫ x ∫ x
1 1
= (s − t)k−1 ds f (t) dt = (x − s)k f (s) ds.
(k − 1)! a t k! a
This shows the formula (2.14). Likewise, there holds
∫ b
k 1
x Ib f (x) = (s − x)k−1 f (s) ds. (2.15)
(k − 1)! x
One way to generalize these formulas to any α > 0 is to use the Gamma function
Γ(z): the factorial (k −1)! in the formula can be replaced by Γ(k). Then for any α > 0,
by replacing k with α, and (k − 1)! with Γ(α), we arrive at the following definition of
Riemann-Liouville fractional integrals. It is named after Bernhard Riemann’s work
in 1847 [Rie76] and Joseph Liouville’s work in 1832 [Lio32], the latter of whom
was the first to consider the possibility of fractional calculus [Lio32].
Definition 2.1 For any f ∈ L 1 (D), the left-sided Riemann-Liouville fractional inte-
gral of order α > 0 based at x = a, denoted by a Ixα f , is defined by
∫ x
1
(a Ixα f )(x) = (x − s)α−1 f (s) ds, (2.16)
Γ(α) a
and the right-sided Riemann-Liouville fractional integral of order α > 0 based at
x = b, denoted by x Ibα f , is defined by
∫ b
1
(x Ibα f )(x) = (s − x)α−1 f (s) ds. (2.17)
Γ(α) x
When α = k ∈ N, they recover the k-fold integrals a Ixk f and x Ibk f . Further, we
define for x → a+ by
(a Ixα f )(a+ ) := lim+ (a Ixα f )(x),
x→a
if the right-hand side exists, and likewise (x Ibα f )(b− ). In case α = 0, we adopt the
convention a Ix0 f (x) = f (x), in view of the following result [DN68, (1.3)].
Lemma 2.2 Let f ∈ L 1 (D). Then for almost all x ∈ D, there hold
α
lim a Ix f (x) = lim+ x Ibα f (x) = f (x).
α→0+ α→0
∫h
Proof Let x ∈ D be a Lebesgue point of f , i.e., limh→0 h1 0 | f (x + s) − f (s)|ds = 0.
∫x
Then let Fx (t) = x−t f (s)ds, for any t ∈ [0, x − a]. Then it can be represented as
Fx (t) = t[ f (x) + ωx (t)], 0 ≤ t ≤ x − a,
where ωx (t) → 0 as t → 0. Thus for each > 0, there exists δ = δ() > 0 such that
|ωx (t)| < , ∀t ∈ (0, δ). (2.18)
Then by the definition of F and integration by parts, we deduce

∫ x−a
α 1
a I x f (x) := s α−1 f (x − s)ds
Γ(α) 0
∫
Fx (x − a) α − 1 x−a
= (x − a)α−1 − Fx (s)s α−2 ds
Γ(α) Γ(α) 0
Fx (x − a) α−1
= (x − a)α−1 − (x − a)α f (x)
Γ(α) Γ(α + 1)
∫ ∫ x−a
α−1 δ
− ωx (s)s α−1 ds + ωx (s)s α−1 ds .
Γ(α) 0 δ
Then it follows from (2.18) and the fact 0 is a pole of Γ(z) that limα→0+ | a Ixα f (x) −
f (x)| ≤ . Since > 0 is arbitrary, it proves the first identity for x. This completes
the proof since the set of Lebesgue points of f ∈ L 1 (D) is dense in D.
The integral operators a Ixα and x Ibα inherit certain properties from integer-order
integrals. Consider power functions (x − a)γ with γ > −1 (so that (x − a)γ ∈ L 1 (D)).
It can be verified directly that for α > 0 and x > a,
α Γ(γ + 1)
a I x (x − a)γ = (x − a)γ+α, (2.19)
Γ(γ + α + 1)
and similarly for x < b
α Γ(γ + 1)
x Ib (b − x)γ = (b − x)γ+α .
Γ(γ + α + 1)
Indeed, (2.19) follows from changing variables s = a + t(x − a) and (2.10):
∫ x ∫
α γ 1 α−1 γ (x − a)γ+α 1
I
a x (x − a) = (x − s) (s − a) ds = (1 − t)α−1 tγ dt
Γ(α) a Γ(α) 0
B(α, γ + 1) γ+α Γ(γ + 1) γ+α
= (x − a) = (x − a) ,
Γ(α) Γ(γ + α + 1)
Clearly, for any integer α ∈ N, these formulas recover the integral counterparts.
Example 2.1 We compute the Riemann-Liouville fractional integral 0 Ixα f (x), α > 0,
of the exponential function f (x) = eλx , λ ∈ R. By (2.19), we have
∞
(λx)k
∞ α k
0 Ix x
α
0 Ix f (x) = 0 Ixα = λk
k=0
Γ(k + 1) k=0 Γ(k + 1)
(2.20)

∞
x k+α ∞
(λx)k
k
= λ = xα = x α E1,α+1 (λx),
k=0
Γ(k + α + 1) k=0
Γ(k + α + 1)
zk
where the Mittag-Leffler function Eα,β (z) is defined by Eα,β (z) = ∞ k=0 Γ(kα+β) , for
z ∈ C, cf. (3.1) in Chapter 3. The interchange of summation and integral is legitimate
since Eα,β (z) is entire, cf. Proposition 3.1.
The following semigroup property holds for k-fold integrals: a Ixk a Ix f = a Ixk+ f ,
k, ∈ N0 , which holds also for a Ixα and x Ibα .
Theorem 2.1 For f ∈ L 1 (D), α, β ≥ 0, there holds
α β β α+β α β β α+β
a Ix a Ix f = a Ix a Ixα f = a Ix f, x Ib x Ib f = x Ib x Ibα f = x Ib f.
Proof The case α = 0 or β = 0 is trivial, since a Ix0 is the identity operator. It suffices
to consider α, β > 0. By changing integration order, we have
∫ x ∫ s
α β 1 α−1
(a Ix a Ix f )(x) = (x − s) (s − t)β−1 f (t) dtds
Γ(α)Γ(β) a a
∫ x ∫ x
1
= f (t) (x − s)α−1 (s − t)β−1 ds dt.
Γ(α)Γ(β) a t
Now by changing variables, the integral in the bracket can be simplified to

∫ x ∫ 1
(x − s)α−1 (s − t)β−1 ds = (x − t)α+β−1 (1 − s)α−1 s β−1 ds
t 0
α+β−1
= B(α, β)(x − t) .
By (2.10), we have
∫ x
β B(α, β) α+β
(a Ixα a Ix f )(x) = (x − s)α+β−1 f (s) ds = (a Ix f )(x).
Γ(α)Γ(β) a
The other identity follows similarly.
The next result collects several mapping properties for a Ixα , and similar results
hold for x Ibα . (i) indicates that a Ixα and x Ibα are bounded in L p (D) spaces, 1 ≤ p ≤ ∞.
Thus for any f ∈ L 1 (D), a Ixα f and x Ibα f of order α > 0 exists almost everywhere
(a.e.). The second part of (ii) is commonly known as the Hardy-Littlewood inequality.
The spaces AC(D) and C k,γ (D) are defined in Section A.1.1 in the appendix.
Theorem 2.2 The following mapping properties hold for a Ixα , and α > 0.
(i) a Ixα is bounded on L p (D) for any 1 ≤ p ≤ ∞.
(ii) For 1 ≤ p < α−1 , a Ixα is a bounded operator from L p (D) into L r (D) for
1 ≤ r < (1 − αp)−1 . If 1 < p < α−1 , then a Ixα is bounded from L p (D) to
p
L 1−α p (D).
(iii) For p > 1 and p−1 < α < 1 + p−1 or p = 1 and 1 ≤ α < 2, a Ixα is bounded from
1
L p (D) into C 0,α− p (D), and moreover, a Ixα f (0) = 0, for f ∈ L p (D).
(iv) a Ixα maps AC(D) into AC(D).
(v) If f ∈ L 1 (D) is nonnegative and nondecreasing then a Ixα f is nondecreasing.
Proof By Young’s inequality in (A.1), we deduce
α 1 (b − a)α
a Ix f L p (D) ≤ (x − a)α−1 L 1 (D) f L p (D) = f L p (D),
Γ(α) Γ(α + 1)
showing the assertion for a Ixα . The first part of (ii) follows similarly, and the proof
of the second part is technical and lengthy, and can be found in [HL28, Theorem 4].
To see (iii), consider first the case p > 1. Suppose first a + h ≤ x < x + h ≤ b. Then
Γ(α)(a Ixα f (x + h) − a Ixα f (x))

∫ x+h ∫ x
= (x + h − s)α−1 f (s)ds − (x − s)α−1 f (s)ds
a a
∫ x+h ∫ x
= (x + h − s)α−1 f (s)ds + ((x + h − s)α−1 − (x − s)α−1 ) f (s)ds
x x−h
∫ x−h
3
+ ((x + h − s)α−1 − (x − s)α−1 ) f (s)ds := Ii .
a i=1
By Hölder’s inequality, since (α−1)p > −(p−1) and ((p−1)−1 (α−1)p+1)p−1 (p−1) =
α − p−1 ,
∫ x+h 1 ∫ x+h p−1
p(α−1)
≤ c f L p (D) hα− p .
p p 1
|I1 | ≤ | f (s)| p ds (x + h − s) p−1 ds
x x
α− p1
This argument also yields |I2 | ≤ c f L p (D) h . Now for the term I3 , since
|(x + h − s)α−1 − (x − s)α−1 | ≤ h|1 − α|(x − s)α−2,
we have
∫ x−h
|I3 | ≤ h|1 − α| | f (s)|(x − s)α−2 ds
a
∫ x−h 1 ∫ x−h p−1
p p(α−2) p
p
≤ h|1 − α| | f (s)| ds (x − s) p−1 ds
a a
p−1
p−1
hα−
p 1
≤ |1 − α| p f L p (D) .
p + 1 − αp
Thus, the assertion follows. When x ≤ a + h, the estimate is direct. If p = 1, we have

α − 1 ≥ 0, and so ∫ x+h
|I1 | ≤ hα−1 | f (s)|ds,
x
and similarly
∫ x−h
|I3 | ≤ h(α − 1)hα−2 | f (s)|ds.
a
Thus the estimate also holds. This shows (iii). Next, if f ∈ AC(D), then f ∈ L 1 (D)
exists a.e. and f (x) − f (a) = (a Ix1 f )(x) for all x ∈ D. Then
α (x − a)α
a Ix f (x) = a Ixα a Ix1 f (x) + (a Ixα f (a))(x) = a Ix1 (a Ixα f )(x) + f (a) .
Γ(α + 1)
Both terms on the right-hand side belong to AC(D), showing (iv). Last,
∫
α (x − a)α 1
a I x f (x) = (1 − s)α−1 f ((x − a)s + a)ds,
Γ(α) 0
which is nondecreasing in x when f is nonnegative and nondecreasing.

Remark 2.1 The mapping properties of the Riemann-Liouville fractional integral
were systematically studied by Hardy and Littlewood [HL28, HL32]. (ii) can be
found in [HL28, Theorem 4]. They also showed that the result does not hold if p = 1
for any α ∈ (0, 1), and that if 0 < α < 1 and p = α−1 , a Ixα f is not necessarily
bounded. (iii) was proved in [HL28, Theorem 12], where it was also pointed out that
the result is false in the cases p > 1, α = p−1 and α = 1 + p−1 .
The next result asserts that a Ixα and x Ibα are adjoint to each other in L 2 (D). It
is commonly referred to as the “fractional integration by parts” formula [SKM93,
Corollary, p. 67].
Lemma 2.3 For any f ∈ L p (D), g ∈ L q (D), p, q ≥ 1 with p−1 + q−1 ≤ 1 + α
(p, q 1, when p−1 + q−1 = 1 + α), the following identity holds
∫ b ∫ b
α
g(x)(a Ix f )(x) dx = f (x)(x Ibα g)(x) dx.
a a
Proof For f , g ∈ L 2 (D), the proof is direct: by Theorem 2.2(i), for f ∈ L 2 (D),
∫b
α α
a I x f ∈ L (D), and hence a |g(x)||( a I x f )(x)| dx ≤ c f L 2 (D) g L 2 (D) . Now the
2
desired formula follows from Fubini’s Theorem. The general case follows from
Theorem 2.2(ii) instead.
The next result discusses the mapping to AC(D) [Web19a, Proposition 3.6]. It is
useful for characterizing fractional derivatives and solution concepts for fdes.
Lemma 2.4 Let f ∈ L 1 (D) and α ∈ (0, 1). Then a Ix1−α f ∈ AC(D) and a Ixα f (a) = 0
if and only if there exists g ∈ L 1 (D) such that f = a Ixα g.
2.3 Fractional Derivatives 27
Proof If there exists g ∈ L 1 (D) such that f = a Ixα g, then by Theorem 2.1,
1−α
a Ix f = a Ix1−α a Ixα g = a Ix1 g ∈ AC(D),
and since∫g ∈ L 1 (D), by the absolute continuity of Lebesgue integral, a Ix1−α f (a) =
x
limx→a+ a g(s)ds = 0. Conversely, suppose a Ix1−α f ∈ AC(D) and a Ix1−α f (a) = 0.
Let G(x) = a Ix1−α f so that G ∈ AC(D) and G(a) = 0. Then g := G exists for a.e.
x ∈ D with g ∈ L 1 (D), and G(x) = a Ix1 g. Furthermore, a Ixα G = a Ix1 f = a Ixα a Ix1 g =
1 α α 1 α
a I x a I x g. Since f , a I x g ∈ L (D), a I x f , a I x a I x g ∈ AC(D) and their derivatives
1 1
exist a.e. as L 1 (D) functions and are equal, that is, f = a Ixα g.
Last, we generalize the mean value theorem in calculus: if f ∈ C(D), g is

Lebesgue integrable on D and g does not change sign in D, then
∫ b ∫ b
f (x)g(x)dx = f (ξ) g(x)dx, (2.21)
a a
for some ξ ∈ (a, b). One fractional version is as follows [Die12, Theorem 2.1] (see
also [Die17]). The classical case is recovered by setting α = 1 and x = b.
Theorem 2.3 Let α > 0 and f ∈ C(D), and let g be Lebesgue integrable on D
and do not change its sign in D. Then for almost every x ∈ (a, b], there exists some
ξ ∈ (a, x) ⊂ D such that a Ixα ( f g)(x) = f (ξ)a Ixα g(x). If additionally, α ≥ 1 or
g ∈ C(D), then this result holds for every x ∈ (a, b].
∫x
Proof Note that a Ixα ( f g)(x) = Γ(α)
1
a
(x − s)α−1 f (s)g(s)ds. If α ≥ 1 and x ∈ (a, b],
α−1
then (x − s)α−1 is continuous in s. Hence, g̃(s) = (x−s) Γ(α) g(s) is integrable on [a, x]
and does not change sign on [a, x]. Thus by (2.21),
∫ x ∫ x
α
a I x ( f g)(x) = f (s)g̃(s)ds = f (ξ) g̃(s)ds = f (ξ)a Ixα g(x).
a a
If 0 < α < 1 and g is continuous, the same line of proof works. Finally, if 0 < α < 1
and g is only integrable, then one can still argue in a similar way, but the integrability
of g̃ holds only for almost all x ∈ D, cf. Theorem 2.2(i).
2.3 Fractional Derivatives
Now we discuss fractional derivatives, for which there are several different defini-
tions. We only discuss Riemann-Liouville and Djrbashian-Caputo fractional deriva-
tives, which represent two most popular choices in practice, and briefly mention
Grünwald-Letnikov fractional derivative. There are several popular textbooks with
extensive treatment on fractional derivatives [OS74, MR93, Pod99, KST06, Die10]
and also monographs [Dzh66, SKM93]. An encyclopedic treatment of Riemann-
Liouville fractional integral and derivatives is given in the monograph [SKM93],
and the book [Die10] contains detailed discussions on the Djrbashian-Caputo frac-
tional derivative. Like before, for any fixed a, b ∈ R, a < b, which are assumed to be
finite unless otherwise stated, we denote by D = (a, b), and by D = [a, b]. Further,
for k ∈ N, f (k) denotes the kth order derivative of f .
2.3.1 Riemann-Liouville fractional derivative
We begin with the definition of the Riemann-Liouville fractional derivative.

Definition 2.2 For f ∈ L 1 (D) and n − 1 < α ≤ n, n ∈ N, its left-sided Riemann-
Liouville fractional derivative of order α (based at x = a), denoted by RaDxα f , is
defined by
∫ x
R α dn n−α 1 dn
D
a x f (x) := ( a I f )(x) = (x − s)n−α−1 f (s) ds,
dx n x Γ(n − α) dx n a
if the integral on the right-hand side exists. Its right-sided Riemann-Liouville frac-
tional derivative of order α (based at x = b), denoted by RxDbα f , is defined by
∫ b
R α ndn n−α (−1)n dn
xDb f (x) := (−1) ( x I f )(x) = (s − x)n−α−1 f (s) ds,
dx n b Γ(n − α) dx n x
if the integral on the right-hand side exists.

Definition 2.2 does not give the condition for the existence of a Riemann-Liouville
fractional derivative. By Lemma 2.4, the condition a Ixn−α f ∈ AC n (D) is needed to
ensure RaDxα f ∈ L 1 (D), and similarly x Ibn−α f ∈ AC n (D) for RxDbα f ∈ L 1 (D). We
denote by
R α + R α
aDx f (a ) = lim+ aDx f (x),
x→a
if the limit on the right-hand side exists. The quantity RxDbα f (b− ) is defined similarly.
Due to the presence of a Ixn−α , RaDxα f is inherently nonlocal: the value of RaDxα f
at x > a depends on the values of f from a to x. The nonlocality dramatically
influences its analytical properties. We compute the derivatives RaDxα (x − a)γ , α > 0,
of the function (x − a)γ with γ > −1. The identities (2.19) and (2.2) yield
R α dn Γ(γ + 1)(x − a)γ+n−α Γ(γ + 1)(x − a)γ−α

aDx (x − a)γ = n
= . (2.22)
dx Γ(γ + 1 + n − α) Γ(γ + 1 − α)
One can draw a number of interesting observations. First, for α N, RaDxα f of the
constant function f (x) ≡ 1 (i.e., γ = 0) is not identically zero, since for x > a:
R α (x − a)−α
aDx 1 = . (2.23)
Γ(1 − α)
This fact is inconvenient for practical applications involving initial/boundary con-

ditions where the physical meaning of (2.23) then becomes unclear. Nonetheless,
when α ∈ N, since 0, −1, −2, . . . are the poles of Γ(z), the right-hand side of (2.23)
dn
vanishes, and thus it recovers the familiar identity dx n 1 = 0. Second, given α > 0,
for any γ = α − 1, α − 2, . . . , α − n, γ + 1 − α is a negative integer or zero, we have
R α
aDx (x − a)α−j ≡ 0, j = 1, . . . , n. (2.24)
In particular, for α ∈ (0, 1), (x − a)α−1 belongs to the kernel of the operator RaDxα and
plays the same role as a constant function for the first-order derivative. Generally, it
implies that if f , g ∈ L 1 (D) with a Ixn−α f , a Ixn−α g ∈ AC(D) and n − 1 < α ≤ n, then

n
R α
aDx f (x) = RaDxα g(x) ⇔ f (x) = g(x) + c j (x − a)α−j ,
j=1
where c j , j = 1, 2, . . . , n, are arbitrary constants. Hence, we have the following

useful result, where the notation (RaDxα−n f )(a+ ) is identified with (a Ixn−α f )(a+ ).
When α ∈ N, the result recovers the familiar Taylor polynomial.
Theorem 2.4 If n − 1 < α ≤ n, n ∈ N, f ∈ L 1 (D), and a Ixn−α f ∈ AC n (D), and
satisfies RaDxα f (x) = 0 in D, then it can be represented by

n R α−j
(D f )(a+ )
a x
f (x) = (x − a)α−j .
j=1
Γ(α − j + 1)
n
Proof By the preceding discussion, we have f (x) = j=1 c j (x − a)α−j . Applying
RD α−k , k = 1, . . . , n, to both sides gives
a x
(RaDxα−k f )(a+ ) = lim+ (RaDxα−k f )(x)

x→a

n
= c j lim+ RaDxα−k (x − a)α−j = ck Γ(α − k + 1),
x→a
j=1
in view of the identities (2.22) and (2.24), which gives directly the assertion.
Example 2.2 The integer-order derivative of eλx , λ ∈ R, is still a multiple of eλx .

Now we compute R0Dxα f of f (x) = eλx . From (2.22), we deduce
∞
(λx)k ∞ RD α x k
x
R α λx
0Dx e = R0Dxα = λk 0 (2.25)
k=0
Γ(k + 1) k=0
Γ(k + 1)

∞
λ k x k−α ∞
(λx)k
= = x −α = x −α E1,1−α (λx).
k=0
Γ(k + 1 − α) k=0
Γ(k + 1 − α)
Generally the right-hand side is no longer an exponential function if α N.

Example 2.3 For α > 0 and λ ∈ C, let f (x) = x α−1 Eα,α (λx α ). Then the Riemann-
Liouville fractional derivative R0Dxα f is given by
∞
(λx α )k ∞ RD α x (k+1)α−1
x
R α
0Dx f (x) = R0Dxα x α−1 = λk 0
k=0
Γ(kα + α) k=0
Γ(kα + α)

∞
λ k x kα−1
= = λx α−1 Eα,α (λx α ).
k=0
Γ(kα)
The interchange of the fractional derivative and summation is legitimate since Eα,β (z)
is an entire function, cf. Proposition 3.1 in Chapter 3. That is, x α−1 Eα,α (λx α ) is
invariant under R0Dxα , and is a candidate for an eigenfunction of the operator R0Dxα
(when equipped with suitable boundary conditions).
Remark 2.2 Note that in Examples 2.2 and 2.3, the obtained results depend on the
starting value a. For example, if we take a = −∞, then direct computation shows
R α λx
−∞ Dx e = λα eλx, λ > 0.
Next we examine analytical properties of RaDxα more closely. In calculus, both

integral and differential operators satisfy the commutativity and semigroup property.
This is also true for a Ixα . However, the operator RaDxα satisfies neither property. The
composition of two Riemann-Liouville fractional derivatives may lead to totally
unexpected results. Further, for f ∈ L 1 (D) and α, β > 0 with α + β = 1, the
following statements hold:
β
(i) If RaDxα (RaDx u) = f , then u(x) = a Ix1 f (x) + c0 + c1 (x − a)β−1 ;
β
(ii) If RaDx (RaDxα v) = f , then v(x) = a Ix1 f (x) + c2 + c3 (x − a)α−1 ;
R
(iii) If aDx w = f , then w(x) = a Ix1 f (x) + c4 ,
1
where ci, i = 0, . . . , 4 are arbitrary constants. Indeed, for (i), one first applies a Ixα to
β
both sides, and then a Ix . Thus two extra terms have to be included: (i) and (ii) include
two constants, whereas (iii) includes only one. Note that if we restrict f ∈ C(D),
then the singular term (x − a)α−1 (or (x − a)β−1 ) must disappear, indicating that a
composition rule might still be possible on suitable subspaces. Indeed, it does hold
γ
on the space a Ix (L p (D)), γ > 0 and p ∈ [1, ∞], defined by
γ p γ
a I x (L (D)) = { f ∈ L p (D) : f = a Ix ϕ for some ϕ ∈ L p (D)},
γ γ
and similarly the space x Ib (L p (D)). A function in a Ix (L p (D)) has the property that
its function value and a sufficient number of derivatives vanishing at x = a.
Lemma 2.5 For any α, β ≥ 0, there holds
R αR β α+β α+β
aDx aDx f = RaDx f, ∀ f ∈ a Ix (L 1 (D)).
α+β α+β
Proof Since f ∈ a Ix (L 1 (D)), f = a Ix ϕ for some ϕ ∈ L 1 (D). Now the desired
assertion follows directly from Theorem 2.6(i) below and Theorem 2.1.
There is no simple formula for the Riemann-Liouville fractional derivative of the

product of two functions. In the integer-order case, the Leibniz rule asserts

k
( f g)(k) = C(k, i) f (i) g (k−i) .
i=0
This formula is fundamental to many powerful tools in pde analysis. In contrast,

generally, for any α ∈ (0, 1), RaDxα ( f g) (RaDxα f )g + f RaDxα g, although there is a
substantially modified version. Hence, many useful tools derived from this, e.g.,
integration by parts formula, are either invalid or require substantial modification.
This result directly implies
R α
aDx ((x − a) f ) = (x − a)(RaDxα f ) + C(α, 1)(RaDxα−1 f ).
Theorem 2.5 Let f and g be analytic on D. Then for any α > 0 and β ∈ R,

∞
R α
aDx ( f g) = C(α, k)(RaDxα−k f )g (k),
k=0

∞
R α α−β−k β+k
aDx ( f g) = C(α, k + β)(RaDx f )(RaDx g).
k=−∞
Proof Since f , g are analytic, the product h = f g is also analytic. Hence, there
exists a power series representation of h based at x,

∞
h(k) (x)
h(s) = (−1)k (x − s)k ,
k=0
k!
and then applying the operator RaDxα termwise (which is justified by the uniform
convergence of the corresponding series), we obtain for n − 1 < α ≤ n, n ∈ N,
∫ x
∞
R α 1 dn (−1)k (x − s)k
aDx h(x) = (x − s)n−α−1 h(k) (x)ds
Γ(n − α) dx n a k=0
Γ(k + 1)
∞ ∫ x
dn (−1)k h(k) (x)
= n Γ(n − α)Γ(k + 1)
(x − s)k+n−α−1 ds
k=0
dx a
∞
dn (−1)k (x − a)k+n−α h(k) (x)
= .
k=0
dx n Γ(n − α)Γ(k + 1)(k + n − α)
Next we simplify the constant on the right-hand side. By the identity (2.2),
1 C(k + n − α − 1, k)
= .
Γ(n − α)Γ(k + 1)(k + n − α) Γ(k + n − α + 1)
Now by the combinatorial identity (2.11), we further deduce
(−1)k C(α − n, k)
= .
Γ(n − α)Γ(k + 1)(k + n − α) Γ(k + n − α + 1)
Consequently,

∞
dn (x − a)k+n−α h(k) (x)
R α
aDx h(x) = C(α − n, k) .
k=0
dx n Γ(n + k − α + 1)
Then by the standard Leibniz rule, we obtain

∞
n
(x − a) j+k−α h(k+j) (x)
R α
aDx h(x) = C(α − n, k) C(n, j)
k=0 j=0
Γ( j + k − α + 1)

∞
∞
(x − a) j+k−α h(k+j) (x)
= C(α − n, k)C(n, j) .
k=0 j=0
Γ( j + k − α + 1)
Now introducing a new variable = j + k of summation and changing the order of

summation lead to

∞

(x − a)−α h() (x)
R α
aDx h(x) = C(α − n, k)C(n, − k)
=0 k=0
Γ( − α + 1)

∞
(x − a)−α h() (x)
= C(α, ) , (2.26)
=0
Γ( − α + 1)
where we have used the identity (2.13). Applying the classical Leibniz rule, inter-
changing of summation order and the identity (2.12) give

∞
∞
(x − a)k+j−α
R α
aDx ( f g) = g (k) (x) C(α, k + j)C(k + j, k) f (j) (x)
k=0 j=0
Γ(k + j + 1 − α)

∞
∞
(x − a)k+j−α
= C(α, k)g (k) (x) C(α − k, j) f (j) (x) .
k=0 j=0
Γ(k + j + 1 − α)
Now using the identity (2.26) completes the proof of the first assertion. The second
assertion can be proved in a similar but more tedious manner, and thus we refer to
[SKM93, Section 15] or [Osl70] for a complete proof.
The next result gives the fundamental theorem of calculus for RaDxα : it is the left
inverse of a Ixα on L 1 (D), but generally not a right inverse.
Theorem 2.6 Let α > 0, n − 1 < α ≤ n, n ∈ N. Then
(i) For any f ∈ L 1 (D), RaDxα (a Ixα f ) = f .
(ii) If a Ixn−α f ∈ AC n (D), then

n−1
(x − a)α−j−1
αR α R α−j−1
a I x aDx f = f − aDx f (a+ ) .
j=0
Γ(α − j)
Proof (i) If f ∈ L 1 (D), by Theorem 2.2(i), a Ixα f ∈ L 1 (D), and by Theorem 2.1,
n−α α
a Ix a Ix f = a Ixn f ∈ AC n (D).
Upon differentiating, we obtain the desired assertion.

(ii) If a Ixn−α f ∈ AC n (D), by the characterization of the space AC n (D) (cf. Theorem
A.1 in the appendix), we deduce

n−1
(x − a) j
n−α
a Ix f = cj + a Ixn ϕ f , (2.27)
j=0
Γ( j + 1)
j−n+α
for some ϕ f ∈ L 1 (D), with c j = RaDx f (a+ ). Consequently,
αR α
a I x aDx f = a Ixα ϕ f . (2.28)
Now applying a Ixα to both sides of (2.27) and using Theorem 2.1, we obtain

n−1
(x − a)α+j
n
a Ix f = cj + a Ixα+n ϕ f .
j=0
Γ(α + j + 1)
Differentiating both sides n times gives

n−1
(x − a)α+j−n
f = cj + a Ixα ϕ f .
j=0
Γ(α + j − n + 1)
This together with (2.28) yields assertion (ii). Alternatively, by definition, we have
∫ x
αR α 1 dn
a I x aDx f (x) = (x − s)α−1 n (a Isn−α f )(s)ds
Γ(α) a ds
∫ x
d 1 dn
= (x − s)α n (a Isn−α f )(s)ds .
dx Γ(α + 1) a ds
Then applying integration by parts n times to the term in bracket gives
∫ x n
1 dn−j n−α (x − a)α−j
(x − s)α−n−1 a Isn−α f (s)ds − a I f (s)| s=a .
Γ(α − n) 0 j=1
ds n−j s Γ(α + 1 − j)
Note that the expression makes sense due to the conditions on f (x).
The next result gives a generalization of Theorem 2.6 [DN68, (1.13)–(1.15)].

Proposition 2.1 Let f ∈ L 1 (D).
Random documents with unrelated
content Scribd suggests to you:
1 1 1
360, macht 180, dividiere mit 250 in 180 macht 2
+ 5
+ 50
von einer
Elle. Eine Elle hat 7 Spannen, multipliziere mit 7: ihr Skd ist 5 15
Spannen.
Nr. 57. Eine Pyramide, 140 ist die uchatebet, 5 14 Spannen ihr
Skd,? die Pirems. Antwort: 93 13 .
Nr. 58. Pirems 93 13 , uchatebet 140,? Sqd. — Antwort: 5 14
wiederum. Die Rechnung enthält wieder einen groben Fehler des
Schreibers.
Nr. 59 a. Pirems ist 12, uchatebet 8? Sqd. Antwort wieder 5 14 .
und Nr. 59 b. Kontrollaufgabe, uchatebet gefragt, ganz
verworren abgeschrieben. Resultat 4 statt 8, wenn die eine Linie 12
und der Sqd 5 14 .
Nr. 60. Ein In von 15 Ellen an seinem Snti und 30 an seinem k3y
—n h r w. Lass mich seine Skd. wissen. Multipliziere 15; 12 davon ist
7 12 , multipliziere 7 12 mit 4 um 30 zu erhalten. Sein Rhi ist 4.... Das ist
sein Skd.
E i s e n l o h r bemerkt Rhi ist unklar, was jedoch ohne Einfluss
auf die Rechnung ist.
E i s e n l o h r und C a n t o r erklären nun die Pir—m—us als die
Seitenkante und die ucha tebet als die Diagonale des
Grundquadrates, während sie durch das Koptische gezwungen sind
die Kaienharu als die Höhe und die snti als die Grundlinie
aufzufassen; sie erklären also den Sekt in den fünf ersten Aufgaben
als den Cosinus des Neigungswinkels der Kante und Grundfläche und
in der letzten als die Cotangente des Böschungswinkels!
Dagegen wenden sich nun ganz unabhängig voneinander
R e v i l l o u t und B o r c h a r d t und schon W e y r trat ihnen bei,
beide zunächst vom Standpunkt des Steinhauers und Architekten;
beide bemerken, dass der Neigungswinkel für den Steinhauer ganz
wertlos.
Die Pyramide wurde in geraden Absätzen gebaut, die dann mit
mächtigen Steinplatten ausgefüllt wurden. Kannte der Arbeiter den
Böschungswinkel, der an allen Quadern derselbe, dann konnte er
jedes Stück richtig behauen; der Neigungswinkel ergab sich dann
ganz von selbst. (Figur.)
Eisenlohr erklärt Piremus als die aus der Säge heraustretende
und seqet leitet er von qd — ähnlich machen — ab und übersetzt es
mit Ähnlichkeit, besser wäre wohl Verhältnis. Revillout sagt, piremus
bedeutet hinausgehen in die Breite oder aus der Breite und beides
passt für die Höhe der Pyramide, die Linie, welche die Spitze mit der
Mitte der Grundlinie verbindet; uchatebet ist die Basis, und beide
Worte sind Synonyma für Kainharu und senti. C a n t o r noch in dem
Brief an Weyr und in der 2. Auflage mutet den Ägyptern die
Ungeheuerlichkeit zu zwei verschiedene Funktionen mit demselben
Namen bezeichnet zu haben. R e v i l l o u t und B o r c h a r d t sagen,
es sei stets die Cotangente des Böschungswinkels und nun erinnern
Sie sich daran, dass Ägypten aus zwei verschiedenen Ländern mit
verschiedener Sprache zusammengewachsen ist. Synonyma sind
häufig, wie wir aus analogen Gründen die ähnliche Erscheinung im
Englischen haben. Die Pyramide heisst smr und in, der Kreis Deben
und kd, der Vater heisst i f und atef, der König bjty und hk3 usw.
Die Höhe ist als Summe der Koordinaten.
Schichtenhöhen, ev. auch aus der
Schattenlänge und die Seite des Grundquadrats ganz direkt messbar.
Die Diagonale muss aber berechnet werden, und zwar mit dem
Pythagoras.
Eisenlohr hat die Winkel der Pyramiden γ und β (siehe Figur
S. 50) berechnet aus
1 1
cos β entweder 5 4 Sp oder 5 2 5 und damit
3 126 18
cos β = 4
oder = 175
= 25
und sie stimmen so ziemlich mit den Winkeln der erhaltenen

Pyramiden, was für den, der weiss, wie oft grobe Fehler der Schüler
geringe Fehler im Resultat geben, nicht wunderbar ist.
Aber Borchardt hat die Neigungswinkel aus der Cotangente
berechnet. Es sind dies Winkel von 54° 14′ 16″; 53° 7′ 48″ (kommt
4 mal vor) und in No. 60, 75° 57′ 50″.
Von diesen Neigungswinkeln gleicht der erste g e n a u bis auf
die Sekunde dem Winkel an der südlichen Steinpyramide von
Daschur (untere Hälfte), der zweite stimmt, ebenso haarscharf mit
dem von Petrie an Ort und Stelle gemessenen Winkel der zweiten
Pyramide von Giseh überein und der letzte ist ebenso g e n a u der
von Petrie 1892 nachgewiesene Winkel aus der Mastaba von
Meidum. Und dass die Werkleute wie die unsrigen nach solchen
Vorschriften, sogenannten Leeren, arbeiteten, das zeigen die von
Petrie aufgedeckten Winkelmauern an den Ecken der Mastaba No.
17 zu Meidum.
Diese Eckwinkelmauern zeigen wie die Erbauer der Mastaba sich
die anzulegende Neigung der Winkel g e n a u nach der in No. 60
gegebenen Vorschrift aufgezeichnet haben. Damit ist die Seqtfrage
entschieden. Sekt ist die Cotangente, rhi vermutlich nichts anderes
als die Ta n g e n t e, die also den Ägyptern auch schon bekannt war.
Das Wort seqd aber leitet Revillout von kd — bewegen ab und
aus dem hapt — Richtscheit, das ein unentbehrliches Werkzeug war;
seine aufrechtstehende Kathete ist 1 Elle, die untere ist in 7 Spannen
und 4 Finger geteilt, und eine Schnur wurde nach dem unteren
beweglichen Punkte geknüpft und gab dem Arbeiter damit durch den
sekd unmittelbar den Winkel, nach dem er seinen Stein zurichtete.
Ein ganz entscheidendes Moment liegt aber in der Natur der
Sache; der königliche Bauherr der Pyramide der sagt einfach, meine
Pyramide soll so und so viel im Geviert haben und so und so hoch
soll sie sein, die Ausführung überlässt er seinem Architekten.
Dass die Ägypter aber mit der Proportionalitäts- und
Ähnlichkeitslehre ganz vertraut waren, das zeigen ihre Abbildungen.
Sie teilten die Wand durch Linien in ein Netz von Quadraten, ganz
wie unsere Ingenieure ihr Zeichenpapier, und trugen in die einzelnen
Quadrate die Figuren in entsprechendem Massstab ein. In dem
sogenannten Grabe Belzoni zu Biban el Moluk ist ein solches
Coordinatensystem erhalten in einer unfertig gebliebenen
Grabkammer König Sety I., des gewaltigsten der Bamessiden nächst
seinem Sohne Ramses II.
Nun zeigen die Bilder und Inschriften, Darstellende Geometrie bei
ähnlich wie die japanischen, keine den Ägyptern.
Perspektive, und man nahm an, dass den
Ägyptern die Perspektive unbekannt gewesen sei. Aber vor etwa 10
Jahren wurden in Fayum, vom trockenen Wüstensand geschützt,
eine grosse Anzahl Totenmasken, Porträts der Verstorbenen,
gefunden, allerdings aus hellenistischer Zeit, die meisten
Handwerksarbeiten, aber auch eine ganze Anzahl Kunstwerke ersten
Ranges, die unsern besten Porträts nicht nachstehen. Und dazu noch
eins, was ich Ihnen nicht vorenthalten darf, oben auf dem Pylon des
Isis-Tempels von Phile, den Ptolemeus IX. Euergetes II. 150 v. Chr.
erneuert hat an der Stelle, von wo der Werkmeister seinen Bau am
besten übersehen konnte, sind in Stein geritzt zwei Zeichnungen
erhalten.
M. H. Nicht Menge, nicht Lambert, nicht Dürer sind die Urheber
der darstellenden Geometrie, hier sehen sie eine Werkzeichnung in
dem Sandstein der Plattform des Pylon, welche Borchardt 1878
aufgenommen hat, mit beigeschriebenen Massen, G r u n d r i s s und
A u f r i s s, und noch steht die Säule, welche genau danach
gearbeitet ist.
Resumieren wir die ägyptische Résumé.
Mathematik so weit wir jetzt davon wissen.
In der Arithmetik hatten sie eine entwickelte Fingerrechnung
und Rechenbretter, auf denen sie mit Steinen rechneten, kannten
alle vier Spezies mit ganzen und gebrochenen Zahlen, wussten mit
Gleichungen 1. und 2. Grades, arithmetischen und geometrischen
Reihen Bescheid und hatten Näherungsmethoden für die Ausziehung
der Quadratwurzeln.
In der Geometrie war ihre Konstruktions- oder Reisskunst hoch
entwickelt. 420 v. Chr. rühmte sich der grosse Demokrit, dass ihn in
der Reisskunst nicht einmal die ägyptischen Harpedonapten
überträfen; sie hatten eine sehr achtenswerte Quadratur des Kreises,
kannten Symmetrie und Proportion, waren mit der Kreisteilung
vertraut, hatten Ähnlichkeitslehre und Anfänge der Trigonometrie
und Elemente der darstellenden Geometrie.
II. Kapitel.
Babylonien — Assyrien.
Ich wende mich nun dem Zuge der semitischen

Völkerwanderung folgend nach dem uralten Kulturland, zwischen
den grossen Strömen Euphrat und Tigris, zum Zweistromland, dem
mâtu Pur Pur, nach Mesopotamien, Babylonien, Assyrien. Hier kam
zu den schon für Ägypten fliessenden Quellen noch B e r o s s o s
hinzu, und die Bibel in weit reichlicherem Masse. Berosus, ein
Babylonischer Priester des Bel der im 3. Jahre v. Chr. in griechischer
Sprache schrieb, hat sich als mit den Traditionen seines Volkes in
Mythos und Geschichte sehr vertraut erwiesen, und es ist zu
bedauern, dass von seinem grossen Werke nur Fragmente durch
Alexander Polyhistor und danach von Josephus und Eusebios
erhalten sind. Verdanken wir doch Berossos die Kunde von dem
Babylonischen Weltschöpfungsmythus, die Sintflut eingeschlossen,
der Quelle des mosaischen, eine Kunde, welche durch die Funde von
Kujundschik-Ninive so glänzend bestätigt und erweitert wurde. Aber
auch die Bibel hat sich als eine nicht zu unterschätzende
Geschichtsquelle erwiesen, und unter dem Einfluss der
Ausgrabungen in den letzten 30 Jahren ist »Babel und Bibel«
(P. D e l i t z s c h) zu einem Schlagwort geworden. Aber erst im
letzten Drittel des 19. Jahrhunderts gelang es durch Entzifferung der
rätselhaften Keilschrift, die Geschichte Vorderasiens auf urkundliche
Grundlage zu stellen. So bedeutend aber die Leistungen der Schüler
E b e r h a r d S c h r a d e r s im letzten Dezennium gewesen sind, so
sagt doch einer der berufensten unter ihnen P. J e n s e n: »Ein
jedes Werk von Assyriologen auch der besten ist und wird auf lange
noch vergleichbar bleiben einem Feld mit Hopfenstangen, von denen
sehr viele zwar annähernd oder durchaus korrekt und gerade, viele
aber, nach allen Richtungen hin, schief stehen.«
Im Gegensatz zu Ägypten, wo wir ein und dasselbe Volk bis zum
heutigen Tage vor uns haben, sind im Zweistromland zwei der Rasse
nach verschiedene Völker zu unterscheiden, die beide langsam
kulturell zusammengeschmolzen sind. Vom Nordosten her,
möglicherweise vom Altai und dem Pamirplateau kamen als
Nomaden in einzelnen Schwärmen die S u m e r e r, ein Volk, das bis
dato für sich steht, die sich vorzugsweise in Südbabylonien in Sumer
ansiedelten, vielleicht vom Meere aus in die Mündungen des Euphrat
und Tigris eindringend. Vom Nordwesten her in gleicher Weise die
S e m i t e n, die sich, zugleich oder früher, vorzugsweise in
Nordbabylonien, Accad (Agade) festsetzten. Naturgemäss mussten
beide Völker zusammenstossen, und in hin und her schwankenden
Kämpfen drangen Sumerer in Accad und Accader in Sumer ein, bis
seit C h a m m u r a b i die Sumerer endgültig den Semiten unterlagen,
die an den Beduinen Arabiens immer frischen Nachschub hatten.
Gehörte doch nach E d . M e y e r, welcher sich dabei stützt auf
R a n k e, Early Babyl. personal names (p. 33, Band VI, 1 des grossen
Hilprecht'schen Sammelwerkes über die Amerik. Ausgrabungen in
Nippur) und noch mehr auf die Monumente, Chammurabi selbst
einem solchen frischen Schwarm Amoriter Beduinen an.
Die sogen. Sumerische Frage gehörte zu Sumerische Frage.
den dunkelsten; während anfangs der
Siebziger die Sumerer als die Kulturträger, die Semiten als rohe
Nomadenhorden hingestellt wurden, hat später ein so bedeutender
Semitologe, wie Halévy, die ganze Existenz der Sumerer geleugnet
und ihre Schrift und Sprache für eine Art Stenographie der
Semitisch-Babylonischen erklärt. Gestützt auf die genaue
Untersuchung der ihm zugänglichen plastischen Denkmäler, hat
E d u a r d M e y e r in seiner Abhandlung »Sumerier und Semiten in
Babylonien« [Abh. d. Kön. Preuss. Akad. d. W. 1906 phil-hist.] die
Frage aufgehellt. An der Existenz der Sumerischen Sprache konnte,
wie Meyer mit Fug bemerkt, nach der Auffindung der griechischen
Übersetzungen bilinguer Syllabare, das sind Listen von Schriftzeichen
mit Angabe ihrer Sumerischen und Assyrischen Silben- und
Wortwerte, nicht mehr gezweifelt werden. Man vgl. die Abhandlung
von T. G . P i n c h e s in den Proc. Bib. Arch. 24, p. 108 und
A . H . S a y c e ibid. p. 120, in denen die Aspiration des p, k und t
durch die Griechische Übertragung konstatiert ist.
Die Rassenfrage wurde durch die Sumerer und Semiten.
bildlichen Darstellungen im wesentlichen auf
Grund der Ausgrabungen d e S a r z e c s, die von H e u z e y
vortrefflich ediert sind, und denen von Nippur, die seit 20 Jahren
ununterbrochen fortgesetzt sind, unzweifelhaft zugunsten eines
selbständigen Volks der Sumerer entschieden, wie es B e z o l d,
W i n k l e r, H i l p r e c h t etc. angenommen hatten. Abgesehen von
der Kleidung, dem sumerischen Mantel und dem semitischen bunten
Plaid, sind scharfe und stereotype Unterschiede vorhanden.
Zunächst zeichnen sich die Semiten wie noch heute durch üppig
wucherndes Bart- und Haupthaar aus, während die Sumerischen
Köpfe bis auf die Augenbrauen völlig ohne Haar sind. Die Nase ist
von der semitischen scharf verschieden, ebenso Mund, Backe und
Stirn. Auch die Frauenköpfe aus Tello sehen durchaus nicht semitisch
aus. »So lehren die Denkmäler mit unwiderleglicher Evidenz, dass es
zwei verschiedene Rassen in Babylonien gegeben hat, eine
semitische [vorzugsweise] im Norden, und eine nicht semitische
[vorzugsweise] im Süden, [die Sumerer]. Zu diesen beiden Rassen
kamen dann als drittes Element die Beduinischen Westsemiten
Chammurabis, die das Haupthaar kurz schneiden und die Lippen
rasieren.«
Die dritte Frage, die von M e y e r Anteil der Sumerer und der
naturgemäss nicht so entscheidend, wie die Semiten an der Kultur.
beiden ersten beantwortet wird, ist die Frage
nach dem Anteil der beiden Rassen an der Kultur. Da hat nun Meyer
nachgewiesen, dass die S u m e r e r d e r Z e i t G u d e a s (etwa
um 2600), i h r e G ö t t e r n i c h t m i t i h r e m e i g n e n
s u m e r i s c h e n Ty p u s , s o n d e r n i n G e s i c h t s b i l d u n g ,
Bart, Haar und Gewandung als Semiten gebildet
h a b e n. Danach haben auf religiösem Gebiete die Semiten
entschieden die Führung gehabt, wenn naturgemäss auch ihre
Religion durch die der Sumerer beeinflusst ist, bis sich eine
einheitliche Religion heranbildete. Meyer glaubt die Sagen von
Gilgamesch, dem Herkules der Babylonier, der Sintflut etc. den
Semiten zuweisen zu können, während besonders die Verbindung
der Götter mit den Sternen, insbesondere die Astrologie, der Hexen-
und Dämonenglauben sumerisch seien, der sich ja von Babylon aus
insbesondere durch das spätere Judentum und das Christentum über
die ganze Welt verbreitet hat.
Die Semiten scheinen auch auf dem Gebiet der Kunst die
Führenden gewesen zu sein, und sehr früh haben sie eine hohe
Stufe der Kunst erreicht, wie die unübertroffene Siegesstele des
Naramsin (s. u.) beweist (vgl. Abbildung).
Siegesstele des Naramsin.
Über einen Punkt aber herrscht unter den Assyriologen volle
Übereinstimmung, d i e E r f i n d u n g d e r B a b y l o n i s c h e n
Schrift, der Keilschrift, ist Eigentum der
S u m e r e r. Zwar ist die von H i l p r e c h t als sumerisch
angesprochene vorsargonische Periode Nippurs schriftlos, und wir
haben aus der Zeit wo in dieser Stadt, dem uralten
Stammesheiligtum der Babylonier, der Sumerische Sturmgott En-lil,
dessen Idiogramm später als Bel gelesen wird, seinen Kult hatte,
keine Tafeln mit Schriftzeichen gefunden, aber der Beweis liegt
darin, dass die semitischen Silbenzeichen ursprünglich sumerische
Worte bedeuten. Meyer weist mit Recht darauf hin, dass die Semiten
als Erfinder der Schrift, alle Konsonanten ihrer Sprache bezeichnet
hätten, und weist auf den entscheidenden Einfluss hin, den die
sumerische Schrift und Sprache auf das Semitische der Babylonier
für Phonetik und Satzbau geübt hat.
Durch die Ausgrabungen de Sarzecs Gudea und die
wissen wir, dass nach dem Tode der grossen Fürstpriester von Telloh.
Semitischen Fürsten Sargon und Naramsin
die Sumerer auch in Accad vorübergehend zur Macht gelangten in
dem Königreich von Sumer und Accad der Fürsten von Ur; wir
kennen durch die so erfolgreichen Ausgrabungen E . d e S a r z e c s
aus wunderbaren Statuen, denen leider der Kopf fehlte (vgl.
Abbildung) und einer Reihe von Schriften, genauer Vertonungen
ihren König oder richtiger Fürstpriester, pateïssi, denn nie nennt er
sich König, G u d e a; nach W i n k l e r war er Vasall des U r e n g u r
von Ur, König von Sumer und Accad, und Gudeas Vorgänger Urnina,
Entemena etc. Ihre Residenz war Schirpurla auch Lagasch, heute
Telloh geheissen; und die Urkunden aus jenen ältesten Zeiten sind
für die Entwicklung der Schrift ganz besonders wichtig. Der Plan und
der Massstab Gudeas (vgl. Abb. S. 62) ist für die Metrologie beinahe
unschätzbar; wie die p. 105 besprochene Arbeit B o r c h a r d t s
beweist, ist er zirka 3000 Jahr in Gültigkeit geblieben, und stimmt
nach der B o r c h a r d t'schen Messung mit L e h m a n n s Hypothesen
(p. 106) vortrefflich.
Gudea mit Plan und Massstab.
Plan der Gudeastatue, 1/2 der nat. Grösse.

Massstab der Gudeastatue, 1/2 der nat. Grösse.
Durch einen merkwürdigen Zufall ist uns Statuen des Gudea.
jetzt auch der K o p f G u d e a s bekannt
geworden. Der Nachfolger de Sarzecs in den Ausgrabungen von Tello
(Sirpurla), der Kapitän G . C r o s, fand unweit der Stelle, wo jener
einen prächtig gearbeiteten Kopf aus Diorit ausgegraben hatte, eine
kleine ganz disproportionierte Statue ohne Kopf, die laut Inschrift als
die der Gudea bezeichnet wurde, von ihm seinem speziellen
Schutzgott, dem er auch den neuen Tempel in Tello gebaut hatte,
dem Ningiszida, dem Sohn des Nin-a-zu (nach Meyer ein anderer
Name für den Götterkönig Anu, den Himmelsgott) gewidmet. L é o n
H e u z e y, der ausgezeichnete Leiter der Assyrischen Abteilung des
Louvre, bemerkte, dass die Brüche des Kopfes und des Torso zu
einander passten, er setzte den Kopf auf den Torso und ohne jeden
Kitt sass er fest (vgl. Rev. d'Assyr. Bd. VI, 1907 p. 19). Dadurch
besitzen wir jetzt 4 Köpfe des Gudea, darunter der von Hilprecht in
seinem Vortrag über die Ausgrabungen im Bêl-Tempel zu Nippur
S. 52 wiedergegebene »Marmorkopf von feinster Arbeit«. Die Köpfe
tragen sämtlich die sogenannte Kappe der Sumerischen Fürsten, die
wir bei Chammurabi (s. u.) wiederfinden, und drei davon den
Turban, der also uralt sumerischen Ursprungs ist. Die scheinbare
Plumpheit und Disproportioniertheit der Körper der Statuen aus Tello
hat Heuzey m. E. sehr zutreffend erklärt. Der Körper diente nur als
Sockel für den Kopf, falls der schwer zu bearbeitende Dioritblock für
eine ganze Statue zu klein war, und H e u z e y bemerkt sehr richtig,
dass unsere Büsten mit ihrer abgespalteten Brust den Sumerern, so
sonderbar vorgekommen waren, wie uns die ihren.
Kopf des Gudea, Federzeichnung nach dem Funde des
Cap. Cros.
Und von der entgegengesetzten Seite Semitische Einwanderung
her, wie heute ziemlich feststeht, von in Vorderasien.
Nordafrika her, drangen nomadische
Semitenschwärme, in verschiedene Volksstämme, richtiger Clane
gespalten in das reiche Zweistromland, und siedelten sich in der 13
Meridian breiten, paradiesisch fruchtbaren Ebene an. D e l i t z s c h
versetzt geradezu das Paradies in die Gegend von Babylon, den
Euphrat und Tigris nennt die Bibel selbst und die beiden andern
Ströme erklärt er für Kanäle, was nicht unmöglich, da die Babylonier
für Kanal und Fluss dasselbe Wort nâru haben. An der jetzigen
grauenhaften Verödung dieses Paradieses erklärt Delitzsch die
Türken für unschuldig, und sicher haben Beduinen und Islam vor
den Türken die Versandung der Kanäle und damit die Verödung des
Landes auf dem Gewissen. Wir hegen die begründete Hoffnung,
dass die deutsche Bagdadbahn und das deutsche Kapital in wenig
mehr als einem Menschenalter die jetzige Wüste wieder zu einem
grossen Garten umgeschaffen haben wird.
Die Unterwerfung der Sumerer gelang Sargon und Naramsin.
um so leichter, als sie keinen Grossstaat
hatten, sondern nur einzelne grosse Städte, in denen sich nach und
nach die Semiten ansiedeln. Die Städte standen unter sogenannten
Fürstpriestern, Pateissi, die sich gegenseitig unter einander
befehdeten, wie wir aus den Inschriften G u d e a s erfahren, und aus
dem von C r o s vor kurzem ausgegrabenen Bericht über die
Verwüstung Tellos durch Lugalzaggissi, den Pateissi der
Nachbarstadt Gishu, bis sie unter die Oberherrschaft Semitischer
»Grosskönige« gerieten, wie Tello unter die des grossen
Semitenfürsten S a r g o n I ., Besitzer von Argade (Accad), der von
Nordbabylonien, dem Lande Accad aus, auch Südbabylonien (Sumer)
unterwarf. Sargons und seines ebenfalls bedeutenden Sohnes
N a r a m s i n Existenz war lange sagenhaft, — die Moses-Mythe wird
auch von Sargon erzählt — bis Nabonahid und die Funde der
Amerikaner in N i p p u r, dem Sitz eines uralten Tempels des Bêl, ihre
historische Existenz bewiesen. Dort ist sogar der Stempel des Sargon
(vgl. Abb.) mit seinen altertümlichen Schriftzeichen gefunden
worden.
N a b o n a h i d, der letzte König von Babylon, war das, was wir
heute einen Romantiker nennen würden, seine Interessen wurzelten
in der Vorzeit, er wollte den uralten Dienst des Schamasch, der
Sonne, und des Sins, des Mondes, wiederherstellen und geriet so in
Konflikt mit der mächtigen Priesterschaft des Marduk-Bel in
Babylonien, deren Unterstützung Cyrus mehr für seinen Erfolg
verdankte als der Macht seiner Waffen. Im Grundstein des Tempels
von Sippar, den Nabonid erneuern wollte, fand er die Urkunde
Naramsins, des Sohnes des Sar-u-ukin.
Die Gelehrten des Königs berechneten
nach den Königslisten die
Regierungszeit des Naramsin auf 3200
Jahre früher, wodurch Sargon auf 3800
v. Chr. gerückt wurde, und mit ihm
Gudea. Trotz mancher Bedenken,
welche gegen dieses hohe Alter
geltend gemacht wurden, insbesondere
von H. Winkler und
C . F. L e h m a n n, nahm doch noch
B e z o l d 1903 diese Daten als richtig
an. Aber der Fund der neuen
Königsliste von Nippur, aus dem Ende
des 3. Jahrtausend der Schrift nach,
durch H i l p r e c h t 1906 im XX. Bd. der
Berichte publiziert und interpretiert, bewies, dass Lehmann mit
seiner Vermutung, dass die Gelehrten des Nabonid sich um etwa 800
Jahre geirrt hatten, im Recht war und die neue Chronologie von
L . W. K i n g (Chronicles conc. early Babyl. kings 2 vol 1907) setzt
Sargon von Akkad auf 2500 v. Chr. auf Grund der Arbeiten
H . R a n k e s.
Über die Chronologie sei gleich hier Babylonisch-Assyrische
bemerkt, dass der Hang der Babylonier zum Chronologie.
genauen Datieren, insbesondere auch die
zahllosen Geschäftsurkunden, die wir von Gudea bis Nabonid
besitzen, uns über die Chronologie der Assyrer weit besser als über
die der Ägypter unterrichtet haben. In Kürze werden uns die
Ausgrabungen, besonders die der Pennsylvania Universität in Nippur
bis ins 4. Jahrtausend hinein eine völlig gesicherte Zeitfolge der
Geschichte gewähren, von Chammurabi bis Kyros, von 2000 bis 539
steht sie schon jetzt auf sicherem Boden. Vom 15. Jahrhundert bis
zum Jahr 1000 können wir uns auf die sogen.
s y n c h r o n i s t i s c h e Geschichte stützen. Nach H . W i n k l e r (die
Keilinschr. u. das alte Test. 3. Aufl. 1903 p. 47) ist es ein Dokument,
in welchem A d a d - n i r a r i III. von Assyrien (812–783) die
Vereinigung Assyriens und Babyloniens als im Interesse beider
Völker hinstellt, nach B e z o l d ein Staatsvertrag beider Länder.
Jedenfalls wird darin in Kürze die Geschichte beider Länder
chronologisch erzählt. Die synchron. Geschichte ist immerhin nicht
ganz einwandfrei, sie enthält gewissermassen den persönlichen
Fehler Adad-niraris. Von diesen sind für Assyrien die
E p o n y m e n k a n o n e s, für Babylonien die K ö n i g s l i s t e n frei.
Das Jahr wurde von Adad-nirari II., etwa um 900 an, zunächst nach
dem die Regierung antretenden Herrscher und dann der Reihe
gemäss, nach den höchsten Beamten benannt, wie in Athen nach
den Archonten. Beide Listen sind Chroniken zum Zweck genauer
Datierung von Rechtshandlungen. Die Vergleichbarkeit des Kanons
mit unserer Zeitrechnung wurde möglich durch Erwähnung der
Sonnenfinsternis im Monat Sivan bei Gelegenheit eines Aufstands
gegen Assur-daja. Die Astronomische Berechnung ergab den 15.
Juni 763. Eine weitere Kontrolle ergab dann der völlig zuverlässige
Kanon des grossen Astronomen Ptolemaios (vgl. Hellas), der uns
hilft bis zur S e l e u c i d e n-Ära (Berossos), deren Beginn zwischen
312 und 311 schwankt und die Arsaciden-Ära von 248, welche
neben der Seleucidenära hergeht.
Die Semiten überschwemmten ganz Westasien, längs der Küste
des Mittelmeeres zogen die Phönizier, besser Kanaanäer, zu denen
die Chabiri, die wir jetzt als Hebräer bezeichnen, gehören, die, wie
es scheint, noch im Anfange der historischen Zeit nicht sesshaft
waren, und erst zur Zeit Chinatôns ihre Stammesgenossen angriffen.
Arvat, Byblos und vor allem Sydon und Tyrus sind Städte der
Phönizier. Die zweite Sammelgruppe der Beduinenschwärme bilden
die Aramäer, mit dem Hauptzweig der Syrer, die südlich von den
Kanaanäern hielten und sich weit nach Norden und Osten
vorschoben. Hier kam es nur in Damaskus, der alt berühmten noch
heute blühenden Handelsstadt zu einer Staatenbildung. Am
ausgedehntesten war die Wanderung des an Zahl stärksten dritten
Zweiges, der Babylonier und Assyrer, die sprachlich und genealogisch
nahe verwandt sind. Doch sind nach den Abbildungen die Babylonier
weit stärker mit den Sumerern blutgemischt als die Assyrer.
Die Assyrer sind sprachlich und auch Geschichte der Babylonier
dem Rassentypus nach mit den Babyloniern und Assyrer.
so nahe verwandt, dass die Annahme ihrer
Abzweigung von diesen, etwa um 1150, nach einem siegreichen
Einfall der Elamiten, sehr wahrscheinlich ist. Sie waren ein Krieger-
und Herrenvolk, das den Priestern einen weit geringeren Einfluss
einräumte als die Babylonier. Ihre Kämpfe, wie die der Babylonier,
gelten, wie leicht begreiflich ist, dem Bestreben, sich die grossen
Handelsstrassen nach Indien und nach dem Kulturzentrum, dem
Mittelmeerbecken offen zu halten. Wird ihnen, durch das
Aufkommen einer nicht semitischen Grossmacht ein Handelsweg im
Westen verlegt, so erkämpfen sie sich einen neuen im Osten. Sehr
bald gingen sie gegen Babylonien aggressiv vor, und der grausame
aber tüchtige A s s u r n a s s i r p a l bringt Babylon völlig unter seinen
Einfluss. Der eigentliche Begründer der Assyrischen Weltmacht
T i g l a t P i l e s e r III. besteigt dann 744 unter dem Namen Pulu
(Phul der Bibel) den Thron Babels und nennt sich König von Sumer
und Accad. Diese Glanzzeit Assyriens hält unter Sargon II. und
seinem Sohn S a n h e r i b an, aber kurz nachdem Sanherib Babylon
zerstört hatte (689) und nach der erfolgreichen Regierung
A s s u r b a n i p a l s (Sardanapal) wird auch Ninive, die Residenz seit
Sanherib von den Medern unter Kyaxares zerstört und zwar weit
gründlicher als Babel.
Bis an die Hochebene Mediens in Nordosten, Elams oder Susa in
Südosten, im Süden bis an die Sümpfe der Mündung des Euphrat
und Tigris in den persischen Busen drangen die Semiten, auch hier
zunächst kein Grossstaat, sondern Städte, die das Stammesheiligtum
bargen als Zentren des Kultus, des Marktverkehrs und Sitz der
Fürsten. Nach Agade und Sirpurla nenne ich Kis, Ur (deren Fürsten
sich seit U r e n g u r Könige der vier Weltgegenden nannten und
Nordbabylonien in Abhängigkeit brachten), Nippur, Larsam und
Babel, die mehr oder minder zentrale Bedeutung gewannen bis
Chammurabi (vielleicht der Amraphel der Bibel) Babel zur Hauptstadt
des Grossstaats Babylon machte, der nun Nordbabylonien (Accad)
und Südbabylonien (Sumer) durch Eroberung von Larsam im Süden
und Absetzung des dortigen Königs einte.
Babel war eigentlich eine Doppelstadt, an einem Ufer Babel —
das Tor Gottes, am andern Borsippa (Birs) — die Stadt des
Mondgottes Sin, dessen Kult in Sumer, insbesondere in Ur blühte,
während in Nordbabylonien der Dienst der Sonne (Schamasch und
Marduk) in den Vordergrund trat.
Ḫammurabi empfängt von Schamasch seine Gesetze.

Wir kennen C h a m m u r a b i wie Chammurabi.
wenige Fürsten des Altertums, und wenige
Regenten dürften ihn in alter und neuer Zeit an Kraft und Weisheit,
und wenn wir seinen Gesichtszügen (s. Abb.) und den zahlreichen
Rechtsschriften Glauben schenken, auch an Gerechtigkeit und Milde
übertroffen haben. Was er für die Stadt Babel getan, berichtet er uns
selbst sumerisch und babylonisch: »Chammurabi, der mächtige
König, der König von Babylon, der König der vier Weltgegenden, der
Begründer des Landes, der König, dessen Taten dem Fleische des
Gottes Schamasch und des Gottes Marduk wohltun, bin ich. Die
Spitze der Mauer von Sippar habe ich mit Erdreich wie einen Berg
erhöht, mit Rohrgeflecht habe ich sie umgeben. Den Euphrat grub
ich ab gen Sippar zu und liess einen Damm dafür aufwerfen.
Chammurabi, der Begründer des Landes, dessen Taten etc. wohltun,
bin ich. Sippar und Babel habe ich auf immerdar zu behaglichen
Wohnstätten gemacht. Chammurabi, der Günstling des Gottes
Schamasch, der Liebling des Gottes Marduk bin ich. Was seit uralten
Tagen kein König dem Herrn der Stadt (dem Schutzgott) gebaut hat,
das habe ich für Schamasch, meinen Herrn, grossartig ausgeführt.«
Chammurabi.
Hatte C . B e z o l d in Ninive und Codex des Ḫammurabi.
Babylon schon C h a m m u r a b i in der eben
zitierten Weise gewürdigt, so wurde die Gestalt dieses grossen
Fürsten in noch weit helleres Licht gerückt durch die Erfolge der
französischen Ausgrabung unter G . d e M o r g a n in Susa, der
Hauptstadt von Elam. In drei Stücken wurde dort im Dezember 1901
und Januar 1902 die Standsäule mit der Gesetzsammlung
Ḫammurabis gefunden, welche 1903 von V. S c h e i l zum ersten
Male ediert und in französischer Sprache erklärt wurde und 1904 von
H . W i n k l e r deutsch und von R . H a r p e r englisch ebenfalls
1904, und vom juristischen Standpunkt von J . K ö h l e r und
E . P e i s e r 1904. Der Codex Hammurabis steht auf einer ethischen
Höhe, welche dem mosaischen vom Sinai nichts nachgibt, und ist
das erste uns erhaltene Corpus juris. Sie genoss, Winkler zufolge,
viele Jahrhunderte das höchste Ansehen — wie die Gesetze des
Moses sind sie von Gott gegeben, das Bild der Säule zeigt, wie der
König die Gesetze von Schamasch empfängt, leider ist das Antlitz
des Königs, der Kappe und Stab trägt, verstümmelt, der Sonnengott
ist mit T u r b a n und Faltenrock bekleidet — sie hat das griechische
Recht, dieses das römische und dieses das unsrige in hohem Grade
beeinflusst. Die Strafe ist natürlich wie bei den Hebräern und
Römern Vergeltung, bei Sittlichkeitsvergehen Abschreckung. Im
Zivilprozess spielt der Eid, grade wie bedauerlicherweise noch heute,
eine hervorragende Rolle. Die Sammlung weist der Frau eine
rechtliche Stellung an, welche sie noch heute in der Türkei nicht
errungen hat, sie schränkt die väterliche Gewalt, ich nenne nur §
168, die Ausweisung des Sohnes betreffend, erheblich ein, und das
Erbrecht ist in sehr zu billigender Weise geregelt, denn auch hier ist
die Frau und die Tochter geschützt. Das Handelsrecht hat er wohl
kaum modifizieren können, denn das war ja zugleich international,
aber das sogenannte Sumerische Familienrecht zeigt, dass dieser
Schutz der weiblichen Familienglieder so recht dem eigenen Sinn des
grossen Königs entsprungen ist. Und so können wir den Worten, mit
denen er auf der Säule sich seiner Taten nach orientalischer Sitte
rühmt — Einleitung und Schluss — wohl Glauben schenken. Die
Stele kam nach Susa als Trophäe zugleich mit anderen wichtigen
steinernen Urkunden im 12/11 Jahr v. Chr., als die Elamiten unter
Sutruk-Nahunte Sippar und Babylonien erobert hatten. Es sei hier
auch erwähnt, dass von dem Kampfe Abrahams zur Befreiung Lots
auch eine Urkunde Chammurabis berichten soll. Die Stele mit der
Gesetzsammlung zeigt am Anfang das Relief, welches die Übergabe
des Codex an den König durch Schamasch schildert, das Relief ist
verstümmelt; (Abbild. S. 69) die Legende ist um so deutlicher.
Die Geschichte Babyloniens und Babylonisch-Assyrische
Assyriens kann ich hier nicht erzählen, sie ist Kultur.
z. T. in der Bibel und bei Herodot und später
bei Arrian, Diodor, und vor allem bei Berossos etc. wenigstens von
2000 ab erzählt; sie ist jetzt bis 4000 v. Chr. so ziemlich aufgehellt;
sie wurde in grossen Zügen durch die verschiedenen Schichten der
einwandernden nomadischen Semitenschwärme und durch die
geographische Lage im einzelnen bedingt. Nach Westen und
Südosten Kämpfe mit den Aramäern und weiter nördlich mit den
Kanaanäern, Phöniziern und Hebräern, die an dem nahen Ägypten
Rückendeckung hatten. Im nördlichen Syrien auch Kämpfe mit dem
uralten vermutlich von Kappadocien her eingedrungenen vielleicht
indogermanischen Stamm der Cheti oder Hetiter, die sich später mit
den Hebräern vermischt haben und mit den Mitani, die noch ziemlich
rätselhaft sind. Im Norden, Osten oder Südosten ist es die
indogermanische Wanderung, die unausgesetzt das babylonisch-
assyrische Reich bedroht; im Norden zusammengefasst als Skythen,
im Osten die Meder, in Südosten die Elamiter mit der Hauptstadt
Susa. Im Süden wieder hemmten die Chaldäer, die im sogenannten
neubabylonischen Reiche nach jahrhundertelangen Kämpfen
schliesslich die Herrschaft an sich rissen. Und hinter den Medern und
Elamitern wieder Indogermanen, deren bedeutsamster Stamm, die
Perser, das ganze babylonische Reich zerstörten.
Die Erschliessung d e r Grotefend und die
b a b y l o n i s c h - a s s y r i s c h e n K u l t u r Entzifferung der Keilschrift.
verdanken wir in erster Linie dem Lehrer am
Gymnasium zu Göttingen: G e o r g F r i e d r i c h G r o t e f e n d. Aus
den Ruinen von Persepolis, der von Alexander dem Grossen in der
Trunkenheit in Brand gesteckten Hauptstadt Persiens, waren im
Laufe der Zeit einige Inschriften in eigentümlichen keilförmigen
Zeichen bekannt geworden, und C a r s t e n N i e b u h r, der Vater
des berühmten Historikers hatte 1770 äusserst sorgfältige und
ausführliche Kopien mitgebracht, welche die allgemeine
Aufmerksamkeit auf die Keilschrift lenkten; er hatte auch schon
bemerkt, dass die Inschriften drei verschiedenen Schriftsystemen
angehörten und von links nach rechts zu lesen waren. Zufällig wurde
Grotefend auf einem Spaziergang im Juli 1802 veranlasst, sich mit
der Entzifferung zu beschäftigen und schon am 4. September 1802
legte er die Resultate seiner Forschung der Göttinger gelehrten
Gesellschaft vor. Er ging davon aus, dass die in drei verschiedenen
Keilschriften und also auch wohl in drei verschiedenen Sprachen
verfassten Inschriften von den Erbauern der Paläste, den persischen
Achämeniden Darius, Xerxes, Artaxerxes etc. herrührten; dass also
vermutlich die erste der drei Sprachen die persische, dass die Texte
wahrscheinlich auch die Namen der Könige enthielten, dass endlich
die Schrift des ersten Systems wegen der geringen Anzahl der
Zeichen eine Buchstabenschrift sein musste; danach verglich
Grotefend die ihm aus der Bibel und den Klassikern und aus der
Zendsprache in den heiligen Büchern Zarathustras bekannten Namen
dieser Könige auf ihre Länge und die Wiederkehr gewisser Zeichen
und kam zu folgendem Schluss: Eine häufig wiederkehrende Gruppe
von Zeichen musste König oder verdoppelt König der Könige
bedeuten, und in den dieser Gruppe vorangehenden Zeichen war der
Name des Königs enthalten; so fand er Darius oder vielmehr die
altpersische Form Dārheūsch, und ein zweiter Name liess sich als
Xerxes-Khschêrsche, ein dritter als Hystaspes-Gôschtaspähe deuten
und ebenso bekam er das Wort Sohn heraus. Die Göttinger gelehrte
Gesellschaft verfuhr mit der Abhandlung Gr. ähnlich wie die dänische
mit der Kaspar Wessels über die geometrische Darstellung der
Complexen Zahlen und die Pariser Akademie mit A b e l s grösster
Arbeit: sie lehnte es ab, die Abhandlung zu veröffentlichen. »Erst
neunzig Jahre später (1893) ist seine Originalabhandlung von Prof.
Wilhelm Meyer in Göttingen wieder aufgefunden und in den
»Gelehrten Nachrichten« der Akademie veröffentlicht worden.«
(H . V. H i l p r e c h t, die Ausgrabungen in Assyrien und Babylonien
1904).
Aber die Entdeckungen Grotefends wurden vor dem Schicksal
der Wessel'schen und Abel'schen bewahrt, dadurch dass sie
Aufnahme fanden in das s. Z. epochemachende Werk von
A . H e e r e n, Ideen über Politik, den Verkehr und den Handel der
alten Welt 4. Aufl. I, 2 S. 345. So war die Grundlage geschaffen, auf
der dann die anderen, ich nenne B e n f e y, H i n k s, O p p e r t,
S p i e g e l weitergebaut haben, so dass jetzt die bisher bekannten
derartigen Texte, mit voller Sicherheit gelesen werden.
In der zweiten Schrift entdeckten N o r r i s und O p p e r t eine
aus Silbenzeichen und einigen Wortzeichen konstruierte Schrift, in
der, wie heute feststeht, die susische oder elamitische Sprache
ausgedrückt wurde; sie enthält gegen 100 Zeichen.
Weit grössere Schwierigkeit bot das dritte System, das über 300
verschiedene Keilschriftzeichen enthielt. Die Entzifferung war schwer
möglich und sie gelang Grotefend nicht. Da entdeckte J a m e s
R i c h, ein geborener Franzose, aber Resident der ostindischen
Kompagnie in Bagdad im Jahre 1820–21 gegenüber der blühenden
Handelsstadt Mossul (Musselin) auf dem linken Tigrisufer die Ruinen
von Ninive und fand zahlreiche Inschriften des dritten Systems.
Bemerkenswert ist es, dass schon im 12. Jahrhundert der spanische
Rabbi B e n j a m i n v o n T u d e l a den Ort von Ninive bestimmt
bezeichnete.
Fast gleichzeitig wurde die sogenannte grosse Dariusinschrift,
eine sehr lange dreisprachige Inschrift am Felsen von Behistun, einer
100 Meter steilen Felswand, an der Grenze des alten Mediens
gefunden und 1835 von H e n r y R a w l i n s o n vermittelst hoher
Leitern auf ungeheueren Papierabklatschen aufgenommen unter
grosser Lebensgefahr —, man nennt die Dariusschrift den
Babylonischen Stein von Rosette —. Von nun ab wuchs die Menge
der ausgegrabenen Inschriften rapide, besonders durch die Arbeiten
von Sir H e n r y L a y a r d und R a s s a m, im Auftrage des British
Museum, in Nimrud, 25 Kilometer von Mossul, die alte Residenz
K e l a c h.
Im Jahre 1881 entdeckte H o r m u z Die wichtigsten
R a s s a m die Ruinen von Sippar. R. hatte Ausgrabungen.
schon 1878 in Balawat, die für die assyrische
Kunst- und Kulturgeschichte gleich wichtigen Bronzetüren
Salmanassars II. gefunden. Von grösster Bedeutung sind die
Ausgrabungen der Franzosen in Tello gewesen, schon dadurch dass
die wunderbaren Funde E . d e S a r z e c s Franzosen, Engländer,
Amerikaner, Deutsche, ja selbst die hohe Pforte zu weiteren Arbeiten
anspornte. Vor de Sarzec hatten schon im Auftrage der französischen
Regierung B o t t a u n d P l a c e in Korsabad den Palast Sargons II.
gefunden und mit Glück gearbeitet, und den Grund zu der grossen
Sammlung im Louvre gelegt.
D e S a r z e c s »Découvertes en Chaldée« von L é o n
H e u z e y 1868 auf Kosten der Regierung herausgegeben, wie schon
die Prachtwerke, welche über Bottas und Places Arbeiten
berichteten: Monument de Ninive découvert et décrit par
E . B o t t a, mesuré et dessiné, par E . F l a n d i n, Paris 1846–50
und V. P l a c e, Ninive et l'Assyrie 1866–69, haben der modernen
Assyriologie den stärksten Impuls gegeben. Die Franzosen setzen die
Ausgrabungen von Tello bis heute fort, daneben hat die Expedition
nach Elam (Susa) unter D e M o r g a n, deren Resultate der
hochverdiente V. S c h e i l mitgeteilt hat, u. a. den Kodex des
Chammurabi aufgefunden. Die Engländer ihrerseits haben fleissig
unter Budge und King in Kujundschik, das Layard seinerzeit den
Franzosen weggenommen, gearbeitet. Die Deutsche
O r i e n t g e s e l l s c h a f t arbeitet seit 1899 unter R . K o l d w e y
und L . B o r c h a r d t mit grossem Erfolg in Babylon und besonders
in Assur. Aber mit den Riesensummen, welche der Staat
Pennsylvanien und seine Universität Philadelphia auf die
Ausgrabungen in Nippur verwandt hat, ist keine Konkurrenz möglich.
Von den Leitern J . P. P e t e r s, H . V. H i l p r e c h t,
J . H . H a y n e s ist besonders der Deutsche Hilprecht der
eigentliche Assyriologe, unter dessen Leitung die Excavations in
Assyria and Babylonia die Resultate der seit 1879 bis jetzt
fortgesetzten Ausgrabungen der Mit- und Nachwelt zugänglich
machen.
Es gelang vier grossen Forschern Die Keilschrift.
R a w l i n s o n, O p p e r t, D e S a u l c y und
dem scharfsinnigen Irländer H i n k s die dritte Schrift und die
Sprache zu entziffern. Die Schrift war eine Verbindung von Wort und
Silbenzeichen, die Sprache eine der arabischen und hebräischen
nahe verwandte, es war die babylonisch-assyrische Sprache. Die
Schrift war ursprünglich eine ziemlich rohe Bilderschrift, zeigt aber
schon in ihren ältesten Formen das Bestreben, Bogen durch Striche
zu ersetzen, aus denen sich dann die Keilschrift entwickelte. So sind
z. B. die ältesten Formen für »Stern«, »Sonne«, »Rohrpflanze«:
für , später
und weiterhin vereinfacht:
und analog haben sich aus den Bildern für Fuss und Weib
die betreffenden Keilschriftzeichen entwickelt.
Diese Keilschriftzeichen lassen sich im wesentlichen auf drei
Grundelemente: den horizontalen Keil , den vertikalen Keil und
den schrägen Keil zurückführen, selten sind die umgekehrten
Keile, der Winkelhaken ist wohl aus Vereinigung zweier Keile
hervorgegangen. Die Keile konnten durch Wiederholung, Neben- und
Übereinanderstellung und Kreuzung zu den mannigfachsten, oft
äusserst komplizierten Gruppen vereinigt, sowohl Worte als Silben im
Assyrischen bezeichnen. Dabei zeigte sich aber eine anfangs
äusserst rätselhafte Erscheinung, die sogenannte Polyphonie.
Dasselbe Zeichen bedeutet sehr oft ein oder mehrere Worte und
daneben noch ein oder mehrere Silben. So bedeutet das Zeichen
nicht nur »Stern«, assyrisch Kakkabu, sondern auch Himmel
schami und Gott ilu und hatte die Silbenwerte an und il. Das Zeichen
hatte nicht nur die Wortbedeutungen »Land« (matu) »Berg«
(schadu), erreichen, erobern Kaschādu; aufgehen (von der Sonne,
napāchu), sondern konnte auch ausserdem als Silbenzeichen in
seinen verschiedenen Zusammenstellungen mit andern Zeichen noch
kur, mad, mat, schad, schat, lat, nad, nat, kin oder gin gelesen
werden.
Das Rätsel löste sich mit einem Schlage als R a w l i n s o n aus
einer Anzahl sehr alter Keilschrifttexte eine neue Sprache in genau
derselben Schrift entdeckte, die Sprache der Sumerer.
Die Beduinenhorden der Babylonier hatten sich mit dem Lande
zugleich der S c h r i f t der Sumerer bemächtigt, der Himmel
hiess sumerisch an, hoch und wurde im Babylonischen Zeichen für
den Begriff Himmel und für die Lautsilbe an, Wortzeichen und
Determinativ für Gott und ebenso wurde Land; Berg, sumerisch
kur als Wortzeichen und Determinativ für Land und Berg und
Silbenzeichen gebraucht.
Diese Erklärung wurde später durch die Auffindung einer
grossen Menge zweisprachiger Texte, babylonisch und sumerisch, in
derselben Schrift bestätigt. (E . B e z o l d: Ninive und Babylon,
Monographien zur Weltgeschichte XVIII 1903.)
Über die Entwicklung der Schrift oder Entwicklung der Keilschrift
den Ursprung der Keilinschriften hat F r. nach Delitzsch.
D e l i t z s c h, dem wir Wörterbuch und
Grammatik des Assyrischen verdanken, 1897 ein Werk veröffentlicht,
das, mögen auch Einzelheiten verbesserungsfähig sein, die
Prinzipien völlig einleuchtend festlegt, nach denen die Sumerischen
Priesterfürsten die Schrift als Verbindung von Wortzeichen —
Idiogrammen — und Silbenzeichen geschaffen haben. Und wenn die
Schrift planmässig mittelst weniger aber wirksamer
Grundgedanken aus der Bilderschrift entstanden ist, so wird damit
auch meine Ansicht, dass das Z a h l s y s t e m eine planmässige und
mit Überlegung ausgeführte Schöpfung derselben Gelehrten ist, im
höchsten Grade wahrscheinlich. Gestützt auf die Formen der Schrift
aus Telloh und die noch älteren aus Nippur, die Geierstele, die Vase
Entemenàs, die Vase Lugat-šug-engur, welche sicher bis gegen 4000
(3700) heraufreicht, und, anknüpfend an des grossen 1905
verstorbenen Jules Oppert Expédition en Mésopotamie 1859 Kap. I,
schied D. zunächst 37 Urzeichen aus, welche sich aus 21 Urbildern
und 16 Urmotiven zusammensetzen. Ich gebe hier die wichtigsten
an: Stern etc., Sonne, aufgehend, Tag, Licht, hell sein,
untergehende Sonne, schwach werden, niedergehen.
Zunehmender Mond (Horn), zunehmen, voll werden, schwinden,
zurückkehren (abnehmender Mond), penis = Mann, männlich,
Mann, Diener, (volva) = Weib, Auge aus ; Hand,
(Fuss) gehen, stehen. Herz, Ochse, Werkzeug zum
Öffnen, daher öffnen, auflösen, Tod, Netz, Geflecht, Gefüge,
Umschliessung, Raum, Kreis (aus ), das Richtungsmotiv,
dessen Ecken die 4 Kardinalpunkte und dessen Axe die Nord-Südlinie
verbildlicht; oder Spitze, daher Gebirge, Kopf,
Bogen, Kurve etc.
Aus diesen Grundelementen werden dann durch
Zusammensetzung gleicher oder verschiedener Zeichen beliebig viele
neue Wortzeichen abgeleitet, welche sich häufig als Definitionen der
dargestellten Begriffe erweisen und auf die Psyche und die Kultur
des Volkes der Sumerer ein so helles Schlaglicht werfen, dass D.
daraufhin den Versuch wagen konnte, ihren Kulturzustand zur Zeit
der Schrifterfindung zu rekonstruieren.
Die Verdoppelung, im Altbabylonischen auch als Kreuzung
sichtbar gemacht, dient zunächst als Pluralzeichen und Iterativum
wie das hebräische Piël, dann aber auch zu Neubildungen. Aus
geben wird durch hinzugeben, addieren tab, dap; aus
gross (nun-rabû) wird Herr d. i. Grösster (Grossmann der
Hottentotten), mit doppelten Zeichen des Umschliessens wird die
Summe bezeichnet: entwickelt zu . Für die
Zusammensetzung ungleicher Zeichen greife ich aus den Beispielen
von D. die folgenden heraus: berufen, erwählen = Auge + werfen,
König = gross + Mensch, Hirt, bei Gudea = Stab + Träger.
Fügte man in das Zeichen für Mund das Zeichen für Brot ein, so
erhielt man: essen, und das eingefügte (Wasser) ergab trinken
und tränken. Die »Schlacht« wird dargestellt als »Handwerk des
Kriegers«, der Regen als gleich Wasser des Himmels, die
Tränen als Wasser des Auges ; Vater als Schützer des
Hauses zu erklären unter Hinweis auf das entsprechende lateinische
pater familias scheint allerdings zweifach fehlerhaft, insofern das
Zeichen im Haus den Feind bedeutet und das sanscrit paṭar schützen
mit piter Vater gar nichts zu tun hat. Die Verkürzung des a zu i in
Jupiter und der Komposition (z. B. suscipio) ist eine ganz spez.
lateinische Eigentümlichkeit. Eins der schlagendsten Beispiele ist
Mond oder Monat, das durch Tag und 30 bezeichnet wird; und
also .
Ein ebenso einfaches wie weittragendes Die Gunierung.
Mittel der Weiterbildung ist die von den
Babylonisch-Assyrischen Grammatikern gunû, d. i. Beschwerung,
genannte Steigerung. Sie besteht in der Hinzufügung von 4 Strichen
oder Keilen, d. h. also Paare von Paaren, die aus Rücksicht auf den
Raum mitunter auf drei reduziert werden. So wird aus Wohnung,
Wohnraum durch Gunierung Palast, Residenz, Grossstadt, und
damit das Determinativ für die Sitze der Pateissi. Aus dem Bilde
des Unterschenkels mit Fuss, das zugleich gehen, stehen, stellen etc.
bedeutet, wird durch Gunierung »Fundament«. Zu den von den
Babylonischen Grammatikern, insbesondere von dem so äusserst
wichtigen Syllabar b der Bibliothek Sardanapals (s. u.) gegebenen
hat D. eine ganze Reihe neuer Gunû Idiogramme abgeleitet, von
denen ich erwähne das Schwert als grosser Dolch; der Vollmond ist
der gunierte Mond, d. h. der grosse, volle, Mond, die Monatsmitte,
die vom Neulicht (s. u.) gezählt wurde und dann Mitte schlechtweg,
archaisch , und das Neulicht selbst wird als der g r o s s e Eingang
des Tages oder als Anfang einer Tagesreihe guniert geschrieben. Es
ist D. gelungen, für einen sehr grossen Teil der Idiogramme meist
recht einleuchtende Ableitungen zu geben, auf Grund derer er es
eben wagen konnte ein Bild des Kulturzustandes der Sumerer nach
Erfindung der Schrift zu geben. Und selbst Erklärung wie die des
Zeichen für Mensch als des auf das Antlitz geworfenen
Knechts oder »H u n d e s« der Götter sind in Anbetracht, dass es
Priester waren, welche die Schrift erfanden, nicht unglaubwürdig,
und recht einleuchtend ist die Erklärung für Ehemann oder Frau als
Verbindung von und durch das Vereinigungszeichen p.
161 (vgl. Abb.).

9611

Uploaded by

Copyright:

Available Formats

9611

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

9611

Uploaded by

Copyright:

Available Formats

Full download ebooks at ebookmeta.

Fractional Differential Equations - An Approach

Download more ebook from https://ebookmeta.com

Fractional Stochastic Differential Equations -

Numerical Methods for Fractal-Fractional Differential

Methods of Mathematical Modelling Fractional

Natures Patterns and the Fractional Calculus Fractional

Differential Equations: A Linear Algebra Approach 1st

Stochastic Models for Fractional Calculus Mark M.

Fractional Inequalities In Banach Algebras 1st Edition

An Elementary Course on Partial Differential Equations

More information about this series at http://www.springer.com/series/34

ISSN 0066-5452 ISSN 2196-968X (electronic)

Mathematics Subject Classiﬁcation: 26A33, 33E12, 35R11, 34A08, 35R30

Fractional differential equations (FDES), i.e., differential equations involving

mathematical theory, especially distinct features when compared with the

London, UK Bangti Jin

3 Mittag-Leffler and Wright Functions . . . . . . . . . . . . . . . . . . . . . . . . 59

Part II Fractional Ordinary Differential Equations

4 Cauchy Problem for Fractional ODEs . . . . . . . . . . . . . . . . . . . . . . . 97

5 Boundary Value Problem for Fractional ODEs . . . . . . . . . . . . . . . . 137

Part III Time-Fractional Diffusion

6 Subdiffusion: Hilbert Space Theory . . . . . . . . . . . . . . . . . . . . . . . . . 175

7 Subdiffusion: Hölder Space Theory . . . . . . . . . . . . . . . . . . . . . . . . . 277

A Mathematical Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 331

CðzÞ Gamma function

Rþ the set ð0; 1Þ

The classical diﬀusion equation

∂t u(x, t) − κΔu(x, t) = 0 (1.1)

is widely used to describe transport phenomena observed in nature, where u denotes

etc. denote Djrbashian-Caputo / Riemann-Liouville fractional derivatives, etc. in

1.1 Random Walk on a Lattice

In the one-dimensional case, such a process can be modeled by

p j (t + Δt) = 12 p j−1 (t) + 12 p j+1 (t),

p(x, t + Δt) = p(x, t) + Δt∂t p(x, t) + O((Δt)2 ),

κ= lim (2Δt)−1 (Δx)2 (1.3)

xn = Δx1 + Δx2 + . . . + Δxn

pn (n− 2 ξ) = (1 − n−1 κtξ 2 + O(n−2 ))n

(1.2), and is given by the Gaussian

1.2 Continuous Time Random Walk

E[Δtn ] = τ, E[Δtn2 ] = τ 2, E[Δxn ] = 0, E[Δxn2 ] = σ 2,

Then the characteristic

The characteristic function for the survival probability Ψ is

Ψ(z) = z−1 − z −1 ψ(z).

Then the probability of taking exactly n steps up to time t is given by

The convolution rule for Laplace transform yields

χn (z) = ψn (z)Ψ(z) = ψ(z)n z −1 (1 − ψ(z)).

We denote the Fourier-Laplace transform of p(x, t) by

Then p(ξ, z) is given by

Since λ and ψ are pdfs,

 ψ(z)| < 1 so the series in (1.4) is absolutely convergent:

Diﬀerent types of ctrw processes can be classiﬁed according to the characteristic

using the normalization condition (1.7), we have

ψτ (t) = τ −1 ψ(τ −1 t) and λσ (x) = σ −1 λ(σ −1 x). (1.8)

Now consider the limit of p(ξ, z; σ, τ) as τ, σ → 0+ . Then simple computation shows

ψτ (z) = ψ(τz) and σ (ξ) = λ(σξ).

ψ(z) = ψ(0) + ψ (0)z + . . . = 1 − z + O(z 2 ) as z → 0.

Using the last two relations, simple algebraic manipulation gives

Now in the formula (1.9), we let σ, τ → 0+ while keeping the relation

with Bα = aα−1 Γ(1 − α).

pn (n− 2 ξ) = (1 − n−1 κtξ 2 + O(n−2 ))n

Then p(ξ, z) is given by

ψ(z)| < 1 so the series in (1.4) is absolutely convergent:

Now consider the limit of p(ξ, z; σ, τ) as τ, σ → 0+ . Then simple computation shows

ψτ (z) = ψ(τz) and σ (ξ) = λ(σξ).

Once again, in (1.13), we ﬁnd the Fourier-Laplace transform p(ξ, z) by letting

p(ξ, t) = Eα,1 (−κα t α ξ 2 ),