Estimations on numerically stable step-size for neutral delay differential systems with multiple delays (Q1300767)

The authors present an asymptotic stability analysis for explicit linear multistep and explicit Runge-Kutta methods for systems of linear neutral delay differential equations of the form \[ \dot u(t)= Lu(t)+ \sum^m_{j= 1} [M_ju(t- \tau_j)+ N_j\dot u(t- \tau_j)],\quad t\geq 0,\tag{1} \] where \(L\), \(M_j\) and \(N_j\) are constant complex \(d\times d\) matrices and the delays satisfy \(0< \tau_1<\cdots< \tau_m\). The derivation of bounds for numerically stable stepsizes for these methods is based on the logarithmic norm and the spectral radius of a certain matrix associated with the given matrices in (1). The results of this paper extend recent ones by \textit{K. J. in't Hout} [IMA J. Numer. Anal. 17, No. 1, 17-27 (1997; Zbl 0867.65046)], and by \textit{G. Hu}, \textit{G. Hu} and \textit{S. A. Meguid} [IMA J. Numer. Anal. 19, No. 3, 349-356 (1999; reviewed above)]. Two numerical examples, with \(m= 2\) and \(d= 2,3\), are used to verify the theoretical results.