Я анализирую устойчивость системы, моделируемой дифференциальными уравнениями с запаздыванием нейтрального типа.Единственная задержка появляется во всех нулевых, первых и вторых производных функции смещения:
d^2(x(t))/dt^2 = f(x(t-tau), d(x(t-tau))/dt, d^2(x(t-tau))/dt^2)
Я использую dde23 для численного моделирования, но есть проблема. Для ввода функций с задержкой я использую: xlag=Z(:,1)
xlag(1)=x(t-tau)
xlag(2)=d(x(t-tau))/dt
а как же d^2(x(t-tau))/dt^2?Я думаю, что это не xlag(3).