Я пишу сценарий Python для решения определенного дифференциального уравнения с помощью SciPy. Я только что скопировал пример на странице do c:
def pend(y,t,b,c):
theta, omega = y
dydt = [omega, -b*omega -c*np.sin(theta)]
return dydt
b = 0.25
c = 5.0
y0 = [np.pi-0.1,0.0]
t = np.linspace(0,10,101)
sol = odeint(pend, y0, t, args = (b,c))
plt.plot(t,sol[:,1])
plt.plot(t,sol[:,0])
Это прекрасно работает и все, но когда я пытаюсь использовать систему Lotka-Volterra, код не работает:
def F(t,n,a,b,c,d):
x, y = n
deriv = [a*x-b*x*y,c*x*y-d*y]
return deriv
t = np.linspace(0,100,100)
a = 1.1
b= 0.4
c = 0.1
d = 0.4
n0 = [10,10]
sol = odeint(F,n0,t,args = (a,b,c,d))
Это возвращает TypeError
<ipython-input-14-ea2a41feaef2> in F(t, n, a, b, c, d)
1 def F(t,n,a,b,c,d):
----> 2 x, y = n
3 deriv = [a*x-b*x*y,c*x*y-d*y]
4 return deriv
5
TypeError: cannot unpack non-iterable float object
Может ли кто-нибудь помочь мне увидеть, что мне не хватает? В частности, почему пример кода работает, если второй код написан с той же структурой. Спасибо!