Я использую solve_ivp из scipy.integrate для решения проблемы начальных значений y'(t) = fun(t,y) с y(t0)=y0. Справка: здесь , из которой видно, что наиболее простой c вызов функции имеет вид:
solve_ivp(fun,[t0,tf],y0)
Как указано здесь, y0=y(t0). Я уже понял, что можно интегрировать в обратном направлении во времени, введя tf, который меньше t0. Однако я не знаю, как я могу интегрировать и назад, и вперед во времени в одном go.
Т.е. я хочу интегрировать в интервале [ti,tf], где ti<t0<tf и где еще y0=y(t0) как и раньше. Это возможно? Или мне нужно решить дважды, а потом собрать воедино решения?
Спасибо за предложения!