Численное решение ODE с фиксированным условием в последний раз с Python

Question

Я пытаюсь решить ODE вида

x'=f(x), x(T)=x_{T}

в python, используя scipy.integrate.odeint.

Однако это имеет вид: scipy.integrate.odeint(func, y0, t, ...), где y0 - вектор начальных условий.У меня нет начального вектора, вместо этого у меня есть вектор для последнего времени T. Как мне изменить подпрограмму scipy.integrate.odeint?

Заранее спасибо

Answer 1

odeint может интегрироваться «назад».Задайте известное значение x(T) в качестве «начального» условия, а затем задайте значения t в порядке убывания от T до того места, где вы хотите остановить решение.

См. thisответ для примера.