Numerical detection of instability regions for delay models with delay-dependent parameters

In this paper we are interested in gaining local stability insights about the interior equilibria of delay models arising in biomathematics. The models share the property that the corresponding characteristic equations involve delay-dependent coefficients. The presence of such dependence requires the use of suitable criteria which usually makes the analytical work harder so that numerical techniques must be used. Most existing methods for studying stability switching of equilibria fail when applied to such a class of delay models. To this aim, an efficient criterion for stability switches was recently introduced in [E. Beretta, Y. Kuang, Geometric stability switch criteria in delay differential systems with delay dependent parameters, SIAM J. Math. Anal. 33 (2002) 1144–1165] and extended [E. Beretta, Y. Tang, Extension of a geometric stability switch criterion, Funkcial Ekvac 46(3) (2003) 337–361]. We describe how to numerically detect the instability regions of positive equilibria by using such a criterion, considering both discrete and distributed delay models.