Базовый результат уравнений с запаздыванием такой (пусть

и все функции хорошие): если зафиксировать поведение искомой функции на отрезке
![$[-c,0]$ $[-c,0]$](https://dxdy-04.korotkov.co.uk/f/b/e/f/bef58064d1766256313f62e90129d23782.png)
, то существует единственное решение с данным фиксированным началом и удовлетворяющее уравнению на отрезке
![$[0,T]$ $[0,T]$](https://dxdy-03.korotkov.co.uk/f/a/a/b/aab0f08201b211261f795050337fa8df82.png)
, причем его можно неограниченно продолжать вправо. А вот влево продолжать не всегда удается.
Если же Вы рассматриваете уравнение на всей оси, то решения в точности соответствуют начальным данным (=функциям на отрезке
![$[-c,0]$ $[-c,0]$](https://dxdy-04.korotkov.co.uk/f/b/e/f/bef58064d1766256313f62e90129d23782.png)
), для которых решение можно продолжить не только вправо, но и влево. Наверное, такие начальные данные стоит искать на аттракторе системы, по аналогии с задачами математической физики, где решения, вообще говоря, не определяются в отрицательные моменты времени, зато таким свойством обладают решения составляющие аттрактор соответствующей динамической системы.
Именно для этого уравнения при положительном c продолжать можно влево, а вправо - не всегда.
-- 10 окт 2017, 11:52 --Вправо такое просто

. То есть производная зависит от прошлых, а не будущих значений. В данном случае проще всего шагами длиной c шагать. Задать значения на [-1, 0] и рассматривать уравнение, как линейное неоднородное

, где g - значения функции f на [-1, 0]. Получится

, C выбирается так, чтобы не было разрыва функции, а разрывы производной вообще-то не исключаются. Потом опять шаг на c, используя вычисленное значение и так далее
(Оффтоп)
до полного удовлетворения
Ну и к более сложным уравнениям техника аналогичная.
Для линейных есть теория, сходная с таковой для обыкновенных линейных ДУ, только вместо полинома появляется квазиполином с экспонентами и с, вообще говоря, бесконечным количеством корней.