Есть ли у scipysolve_ivp метод, аналогичный методу ode set_f_params? ⇐ Python

[Гость]

Для метода ode scipy существует функция set_f_params для обновления любых предоставленных аргументов: используйтеsolve_ivp с различными параметрами. Это то, что я ищу с помощью методаsolve_ivp, но не смог найти воспроизведение. Кто-нибудь знает, как это сделать?
По сути, я распространяю вектор состояния y, и я хочу вычислить изменение дополнительных аргументов, которые я передаю вsolve_ivp, затем я хочу снова передать эти измененные аргументы вsolve_ivp. Для тех, кому интересно, это относится к распространению двух естественных спутников, и я хочу обновить коэффициенты гравитационного поля во время распространения, но мне это очень трудно сделать.
Чтобы быть более конкретно, точный набор уравнений, который я пытаюсь реализовать, выглядит следующим образом:
Набор связанных дифференциальных уравнений
Я могу реализовать уравнение. 5.20, 5.21 и второе уравнение 5.22 штрафа в состоянии y. Затем я могу вызватьsolve_ivp, и он просто вернет мне решение. Однако я не понимаю, как на основе результирующего $\Delta Z^\nu {l,m}$ для каждого шага интегрирования я могу обновить $\Delta Z{l,m}$ соответственно, используя левое уравнение в 5.22. Эта новая дельта $Z_{l,m}$ будет использоваться в правой части уравнений 5.20 и 5.21, поэтому важно, чтобы я делал все это одновременно.
Я бы сделал ценю любую помощь

Подробнее здесь: https://stackoverflow.com/questions/781 ... t-f-params