Asymptotic and numerical stability of systems of neutral differential equations with many delays

We are concerned with the asymptotic stability of a system of linear neutral differential equations with many delays in the form y′(t)=Ly(t)+∑i=1dMiy(t−τi)+∑i=1dNiy′(t−τi), where L,Mi,Ni∈CN×N(i=1,2,…,d) are constant complex matrices, τi>0(i=1,2,…,d) are constant delays and y(t)=(y1(t),y2(t)…yN(t))T is an unknown vector-valued function for t>0t>0. We first establish a new result for the distribution of the roots of its characteristic function, next we obtain a sufficient condition for its asymptotic stability and then we investigate the corresponding numerical stability of linear multistep methods applied to such systems. One numerical example is given to testify our numerical analysis.