Python，Scipy：Odeint对标量变量的无效索引

Question

以下是我正在编写的代码，用于解决Fermi-Pasta-Ulam-Tsingou问题以模拟粒子链的行为。 除了主要粒子之外，等式还取决于最近邻居的初始条件：

之前我已经定义了变量，但我遇到了使用Odeint和vector的问题：

def FPUT(y, t):
    for j in N:
        u1, u2 = y
        dydt = [u2[j], -( k( u1[j+2] + u1[j] - 2*u1[j+1]) * ( 1 + alpha( u1[j+2] - u1[j] ) ))]
return u1,dydt

y0 = [0,0] 
t = np.linspace(0, 10, 101)

。

sol = odeint(FPUT, y0, t)

运行代码时出现以下错误：

IndexError：标量变量的索引无效。

error description screenshot

请帮助，我对如何解决这个问题一无所知。

如果您需要更多信息：https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem