Я хотел бы решить пару вложенных дифференциальных уравнений в R. Предположим, что x '= x + a с заданным x (T) и y' = y + x + b с заданным y (0), где a и б постоянны. Первое уравнение не зависит от второго, и оно может быть решено с помощью deSolve путем обращения времени, поскольку оно имеет терминальное условие. Второе уравнение должно быть решено с использованием решения первого уравнения и со временем, продвигающимся вперед, поскольку оно имеет начальное условие. Есть ли способ сохранить значение решения для x таким образом, чтобы его можно было использовать для решения уравнения y?
Я попытался использовать X <- out [, "X"], где out - это выход для решения функции для уравнения x, и подключить X к решению функции для y, но я получаю следующее сообщение: </p>
Число производных, возвращаемых функцией func (), должно равняться длине вектора начальных условий.
Я понимаю природу проблемы, но не знаю, как от этого избавиться.