Системы уравнений с запаздываниями широко применяются в различных областях современного математического моделирования. В ходе разработки структуры математической модели и идентификации ее параметров приходится многократно решать задачу Коши для подобных систем. В случае высокой размерности системы, а также при условии жесткости задачи численное решение уравнений с запаздываниями может требовать значительных вычислительных и временных затрат. Таким образом, разработка и реализация эффективных алгоритмов численного решения различных классов уравнений с запаздывающими аргументами является актуальной задачей. В настоящей статье представлена модифицированная версия программного комплекса DIFSUBDEL, в которой реализованы методы численного решения дифференциальных уравнений с запаздываниями на основе линейных многошаговых методов. Переработанная версия разработана с применением принципов структурного программирования и является значительно более удобной в эксплуатации, чем исходная, а также обладает свойством потокобезопасности, что позволяет использовать комплекс в качестве блока в системах, основанных на технологиях параллельного программирования с общей памятью. Был проведен сравнительный анализ производительности переработанной системы DIFSUBDEL c другими существующими программными реализациями численных методов решения дифференциальных уравнений с запаздыванием и показана ее эффективность.