Solve_ivp 设置具有序列的数组元素。请求的数组在 1 维之后具有不均匀的形状。检测到的形状为 (4,)+不均匀

Question

我将solve_ivp 用于4 个ODE 方程的系统，但它给了我这个错误。可能是因为我需要的solve_ivp中的y0是4项（初始值：x0，y0，w0，z0），我像数组一样输入它。如何在求解器中输入这些初始值？如果我尝试像列表一样放置会发生相同的错误。 我写了一个我的代码示例（写完太长了）

def funct(iv,t):
    a=iv[0]
    b=iv[1]
    c=iv[2]
    c=iv[3]
    # ODE
    dxdt=...
    dydt=...
    dwdt=...
    dzdt=...
    ODES=[dxdt,dydt,dwdt,dzdt]
    return ODES
x0=1
y0=2
w0=3 
z0=4 
initial_values=np.array([x0,y0,w0,z0]) #*******
t_s = np.linspace(0,100,1001)
#Integración del modelo
inputs = (D,E,F,G,H,J,K)
solut = odeint(funct, initial_values, t_s)


 

Answer 1

根据文档：https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.solve_ivp.html

solut = solve_ivp(funct, t_s, initial_values)

这将是正确的顺序。