6533b85bfe1ef96bd12ba151

RESEARCH PRODUCT

Reprint of: Approximate Taylor methods for ODEs

Giovanni RussoPep MuletAntonio BaezaDavid ZoríoSebastiano Boscarino

subject

ODE integratorsGeneral Computer ScienceTaylor methodsMathematicsofComputing_NUMERICALANALYSISGeneral EngineeringOdeFunction (mathematics)Present procedure01 natural sciences010101 applied mathematicsFaà di Bruno's formulasymbols.namesakeTaylor seriessymbolsApplied mathematicsOrder (group theory)0101 mathematicsMathematics

description

Abstract A new method for the numerical solution of ODEs is presented. This approach is based on an approximate formulation of the Taylor methods that has a much easier implementation than the original Taylor methods, since only the functions in the ODEs, and not their derivatives, are needed, just as in classical Runge–Kutta schemes. Compared to Runge–Kutta methods, the number of function evaluations to achieve a given order is higher, however with the present procedure it is much easier to produce arbitrary high-order schemes, which may be important in some applications. In many cases the new approach leads to an asymptotically lower computational cost when compared to the Taylor expansion based on exact derivatives. The numerical results that are obtained with our proposal are satisfactory and show that this approximate approach can attain results as good as the exact Taylor procedure with less implementation and computational effort.

https://doi.org/10.1016/j.compfluid.2018.03.058