“…For this reason, in the last years the interest in the study of the dynamics of delay models has been increasing and important challenges, in particular numerical, have been identified. Indeed, delay equations describe infinite-dimensional dynamical systems, and theoretical results should be complemented with efficient numerical methods to approximate solutions of initial value problems [3,4,5,7,17,15,16,19,44,43], boundary value problems [48,49,50], and to investigate the stability of equilibria and periodic solutions [10,11,12,13,14,47,52,53,66]. In applications the attention is focused not only on the approximation of the dynamical properties for some given parameter values, but also on how such properties change when varying some parameters.…”