Open Access
2005 A Software Package for the Numerical Integration of ODEs by Means of High-Order Taylor Methods
Àngel Jorba, Maorong Zou
Experiment. Math. 14(1): 99-117 (2005).

Abstract

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 paper discusses 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.

Citation

Download Citation

Àngel Jorba. Maorong Zou. "A Software Package for the Numerical Integration of ODEs by Means of High-Order Taylor Methods." Experiment. Math. 14 (1) 99 - 117, 2005.

Information

Published: 2005
First available in Project Euclid: 30 June 2005

zbMATH: 1108.65072
MathSciNet: MR2146523

Rights: Copyright © 2005 A K Peters, Ltd.

Vol.14 • No. 1 • 2005
Back to Top