Sub-ODEs Simplify Taylor Series Algorithms for Ordinary Differential Equations
Abstract
A Taylor method for solving an ordinary differential equation initial-value problem , , computes the Taylor series (TS) of the solution at the current point, truncated to some order, and then advances to the next point by summing the TS with a suitable step size. A standard ODE method (e.g. Runge-Kutta) treats function as a black box, but a Taylor solver requires to be preprocessed into a code-list of elementary operations that it interprets as operations on (truncated) TS. The trade-off for this extra work includes arbitrary order, typically enabling much larger step sizes. For a standard function, such as , this means evaluating , where are TS. The sub-ODE method applies the ODE , obeyed by , to in-line this operation as . This gives economy of implementation: each function that satisfies a simple ODE goes into the "Taylor library" with a few lines of code--not needing a separate recurrence relation, which is the typical approach. Mathematically, however, the use of sub-ODEs generally transforms the original ODE into a differential-algebraic system, making it nontrivial to ensure a sound system of recurrences for Taylor coefficients. We prove that, regardless of how many sub-ODEs are incorporated into , this approach guarantees a sound system. We introduce our sub-ODE-based Matlab ODE solver and show that its performance compares favorably with solvers from the Matlab ODE suite.
Cite
@article{arxiv.2503.21078,
title = {Sub-ODEs Simplify Taylor Series Algorithms for Ordinary Differential Equations},
author = {Nedialko S. Nedialkov and John D. Pryce},
journal= {arXiv preprint arXiv:2503.21078},
year = {2025}
}
Comments
25 pages