A software package for the numerical integration of ODEs by means of high-order Taylor methods. (English) Zbl 1108.65072

Jorba, Àngel; Zou, Maorong

Exp. Math. 14, No. 1, 99-117 (2005).

Summary: This paper revisits the Taylor method for the numerical integration of initial value problems of ordinary differential equations (ODEs). The main goal is to present a computer program that outputs a specific numerical integrator for a given set of ODEs. The generated code includes a function to compute the jet of derivatives of the solution up to a given order plus adaptive selection of order and step size at run time. The package provides support for several extended precision arithmetics, including user-defined types.
The authors discuss the performance of the resulting integrator in some examples, showing that it is very competitive in many situations. This is especially true for integrations that require extended precision arithmetic. The main drawback is that the Taylor method is an explicit method, so it has all the limitations of these kind of schemes. For instance, it is not suitable for stiff systems.

MSC:

65L05	Numerical methods for initial value problems involving ordinary differential equations
34A34	Nonlinear ordinary differential equations and systems
65Y15	Packaged methods for numerical algorithms

Keywords:

numerical examples; automatic differentiation; software; Taylor method; initial value problems; computer program; package; extended precision arithmetics

Software:

Taylor; ATOMFT; ADIFOR