Author 
Comment/Response 
James

07/20/13 3:25pm
I'm trying to figure out how to successfully write this code, but I can't seem to figure it out.
The x's are points, and Z is a coupling that is like some liquid, so it isn't a point. What could be the problem with my code?
N0 = 4; \[Mu] = 2; \[Rho] = 1; d = 1; n = 1;
eqone = Flatten[
Join[
Table[x[i]''[t] == \[Epsilon]*(1  x[i][t]^2)*x[i]'[t]  x[i][t] +
d*Z[t]  x[i][t], {i, 0,
N0}], {Z'[t] == \[Rho]*d/n (Sum[x[i][t]  Z[t], {i, 0, N0}]),
x[0][0] == 1, x[0]'[0] == .5, x[N0/2][0] == 0, x[N0/2]'[0] == 0},
Table[x[i]'[0] == 0, {i, 1, N0/2  1}],
Table[x[i][0] == 0, {i, 1, N0/2  1}],
Table[x[i]'[0] == 0, {i, N0/2 + 1, N0}],
Table[x[i][0] == 0, {i, N0/2 + 1, N0}]]]
eqtwo = NDSolve[eqone, {Table[x[i], {i, 0, N0}], Z}, {t, 100}]
URL: , 
