Ageneral and easy-to-code numerical method based on radial basis functions
(RBFs) collocation is proposed for the solution of delay differential equations
(DDEs). It relies on the interpolation properties of infinitely smooth RBFs, which allow
for a large accuracy over a scattered and relatively small discretization support.
Hardy's multiquadric is chosen as RBF and combined with the Residual Subsampling
Algorithm of Driscoll and Heryudono for support adaptivity. The performance
of the method is very satisfactory, as demonstrated over a cross-section of
benchmark DDEs, and by comparison with existing general-purpose and specialized
numerical schemes for DDEs.