A three-stage, VSVO, Hermite-Birkhoff-Taylor, ODE solver

The new variable-step, variable-order, ODE solver, HBT(p) of order p, presented in this paper, combines a three-stage Runge-Kutta method of order 3 with a Taylor series method of order p-2 to solve initial value problems , where y:R→R^d and f:R×R^d→R^d. The order conditions satisfied by HBT(p) are formulated and they lead to Vandermonde-type linear algebraic systems whose solutions are the coefficients in the formulae for HBT(p). A detailed formulation of variable-step HBT(p) in both fixed-order and variable-order modes is presented. The new method and the Taylor series method have similar regions of absolute stability. To obtain high-accuracy results at high order, this method has been implemented in multiple precision.     

Variable-step variable-order 3-stage Hermite-Birkhoff-Obrechkoff DDE solver of order 4 to 14

This article presents a solver for delay differential equations (DDEs) called HBO414DDE based on a hybrid variable-step variable-order 3-stage Hermite-Birkhoff-Obrechkoff ODE solver of order 4 to 14. The current version of our method solves DDEs with state dependent, non-vanishing, small, vanishing and asymptotically vanishing delays, except neutral type and initial value DDEs. Delayed values are computed using Hermite interpolation, small delays are dealt with by extrapolation, and discontinuities are located by a bisection method. HBO414DDE was tested on several problems and results were compared with those of known solvers like SYSDEL and the recent Matlab DDE solver ddesd and statistics show that it gives, most of the time, a smaller relative error than the other solvers for the same number of function evaluations.