Numerical Solution of Differential Equations

This book is a tutorial written by researchers and developers behind the FEniCS Project and explores an advanced, expressive approach to the development of mathematical software. The presentation spans mathematical background, software design and the use of FEniCS in applications. Theoretical aspects are complemented with computer code which is available as free/open source software. The book begins with a special introductory tutorial for beginners. Following are chapters in Part I addressing fundamental aspects of the approach to automating the creation of finite element solvers. Chapters in Part II address the design and implementation of the FEnicS software. Chapters in Part III present the application of FEniCS to a wide range of applications, including fluid flow, solid mechanics, electromagnetics and geophysics.

The numerical analysis of stochastic differential equations differs significantly from that of ordinary differential equations due to peculiarities of stochastic calculus. This book provides an introduction to stochastic calculus and stochastic differential equations, in both theory and applications, emphasising the numerical methods needed to solve such equations. It assumes of the reader an undergraduate background in mathematical methods typical of engineers and physicists, though many chapters begin with a descriptive summary.

This is a textbook for a one semester course on numerical analysis for senior undergraduate or beginning graduate students with no previous knowledge of the subject. The prerequisites are calculus, some knowledge of ordinary differential equations, and knowledge of computer programming using Fortran. Normally this should be half of a two semester course, the other semester covering numerical solution of linear systems, inversion of matrices and roots of polynomials. Neither semester should be a prerequisite for the other.

Domain decomposition methods are divide and conquer methods for the parallel and computational solution of partial differential equations of elliptic or parabolic type. They include iterative algorithms for solving the discretized equations, techniques for non-matching grid discretizations and techniques for heterogeneous approximations.