使用Python在最后一次固定条件下的ODE的数值解

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 answer。