RK45 ode求解器python3

Question

请如何使用此function？

我有一个函数，用于返回右侧微分方程right(t, x0)的列表。

import scipy.integrate as ode
t, r =ode.RK45(right(t, x0), t0, tmax, dt)

错误：

TypeError: 'list' object is not callable

第一个论点应该是什么？

Answer 1

因此，您将函数作为对象而不是函数值传递

solver = ode.RK45(right)
solver.set_initial(y0,t0)

然后在solver.integrate上循环计算所需的集成步骤，请参阅文档中的示例。

如果您有scipy的任何最新版本，请使用较新的solve_ivp界面，您可以在一个呼叫中提供所有内容

solution = solve_ivp(right, [t0, tmax], y0, method = "RK45", t_eval = np.arange(t0,tmax,dt))

并在solution.t和solution.y中找到解决方案成分，其中datei <- scan("https://isor.univie.ac.at/about-us/People.html", what = list(name = "", room ="", email = "", website = "")) datei 是每个状态成分的时间序列的元组。