In this paper, a class of Runge-Kutta methods for solving neutral delay differential equations (NDDEs) is proposed, which was first introduced by Bassenne et al. (J. Comput. Phys. 424, 109847, 2021) for ODEs. In the study, the explicit Runge-Kutta method is multiplied by an operator, which is a Time-Accurate and highly-Stable Explicit operator (TASE-RK), resulting in higher stability than explicit RK. Recently, the multi-parameter TASE-W method was extended by González-Pinto et al. (Appl. Numer. Math. 188, 129–145, 2023). We generalized TASE-RK and TASE-W to NDDEs for the first time. Then, by applying the argument principle, sufficient conditions for delay-dependent stability of TASE-RK and TASE-W combined with Lagrange interpolation for NDDEs are investigated. Finally, numerical examples are carried out to verify the theoretical results.
