The following is a code I'm writing to solve the Fermi-Pasta-Ulam-Tsingou problem to Model the behavior of a chain of particles.
the thing that the equation depends on the initial condition of the nearest neighbors in addition the main particles:
I have defined the variables before and I have issue with working Odeint and vectors:
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] ) ))]