Lyapunov exponents play an important role among the invariant characteristics of dynamical systems. Analysis of the spectrum of Lyapunov exponents is widely used to study the complex dynamics in systems of ordinary differential equations and models that can be reduced to maps. In this paper we consider the question of numerical evaluation of Lyapunov exponents for delay differential equations. We used a new algorithm with FFT, tested it on Hutchinson equation with known Lyapunov exponents spectrum and compared with results of the old algorithm. The so-called quasi-stable behavior is typical for a number of gene networks models and neuronal associations that have been studied recently. The phenomenon of cycle (k-dimensional torus) quasistability in dynamical system is characterized by the fact that some of its multipliers are asymptotically close to the unit circle, and the remaining multipliers are modulo less than one (with the exception of a simple root 1 (k such roots)). In some cases, it is possible to prove the existence and give an asymptotic estimate of the multipliers of considered system by using the large parameter methods. However, if these methods are not applicable, it is necessary to obtain a tool for the numerical estimation of multipliers. Such a tool is provided by algorithms for Lyapunov exponents estimation. Given that the equations with delay are often applied in the models of neural and gene networks, the algorithm for Lyapunov exponents estimation for such systems will be in demand.