RK45 ode求解器python3

请如何使用此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

第一个论点应该是什么?

qijiannian 回答:RK45 ode求解器python3

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

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.tsolution.y中找到解决方案成分,其中datei <- scan("https://isor.univie.ac.at/about-us/People.html",what = list(name = "",room ="",email = "",website = "")) datei 是每个状态成分的时间序列的元组。

本文链接:https://www.f2er.com/3049082.html

大家都在问