A new efficient algorithm for the computation of the stability chart of linear time delay systems is proposed and tested on several examples. The stability chart is obtained by investigating the 2d-parameter space by a first coarse square grid which is then adaptively refined by triangulation to match the desired tolerance. This leads to a considerable reduction in computational cost with respect to investigate a uniform fine square grid. Stability of each point is determined by approximating the rightmost characteristic root real part via a numerical scheme recently developed by the authors and based on pseudospectral differencing methods. A Matlab code is included in appendix.
