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.