Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'ND Solve problem' topicStudent Support Forum > General > "ND Solve problem"

Next Comment >Help | Reply To Topic
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: ,

Subject (listing for 'ND Solve problem')
Author Date Posted
ND Solve problem James 07/20/13 3:25pm
Re: ND Solve problem Bill Simpson 07/20/13 11:24pm
Re: Re: ND Solve problem James 07/23/13 1:01pm
Re: ND Solve problem Peter Pein 07/21/13 06:58am
Next Comment >Help | Reply To Topic