我有一个微分方程系统,需要计算雅可比矩阵。以下代码抛出AttributeError: 'list' object has no attribute 'ravel'
。我错过了什么?
import numpy as np
import numdifftools as ndt
def rhs(z, t=0):
x,y = z
xdot = (x/5 + y)*(-x**2+1)
ydot = -x*(-y**2+1)
return [xdot, ydot]
Jfun = ndt.Jacobian(rhs)
Jfun([1,1])
答案 0 :(得分:8)
只是做:
return np.array([xdot, ydot])
代替。这应该有用......