problem with NDSolve::dvnoarg:
- To: mathgroup at smc.vnet.net
- Subject: [mg116529] problem with NDSolve::dvnoarg:
- From: tarun dutta <tarunduttaz at gmail.com>
- Date: Fri, 18 Feb 2011 04:36:07 -0500 (EST)
n = 5;
p = 5/1000;
q = 1/10;
c[-1][t] = 0; d[-1][t] = 0;
c[n + 1][t] = 0; d[n + 1][t] = 0;
eqn1 = Table[
c[i]'[t] == ((1/2)*i (i - 1) - q*i)*c[i][t] -
p*(Sqrt[i]*
c[i - 1][t]*(Sum[Sqrt[i]*d[i - 1][t]*d[i][t], {i, 0, n}]) +
Sqrt[i + 1]*
c[i + 1][
t]*(Sum[Sqrt[i + 1]*d[i - 1][t]*d[i][t], {i, 0, n}])), {i,
0, n}];
eqn2 = Table[
d[i]'[t] == ((1/2)*i (i - 1) - q*i)*d[i][t] -
p*(Sqrt[i]*
d[i - 1][t]*(Sum[Sqrt[i]*c[i - 1][t]*c[i][t], {i, 0, n}]) +
Sqrt[i + 1]*
d[i + 1][
t]*(Sum[Sqrt[i + 1]*c[i - 1][t]*c[i][t], {i, 0, n}])), {i,
0, n}];
eqn3 = Table[Sum[(c[i]^2)[t], {i, 0, n}] == 1, {i, 1}];
eqn4 = Table[Sum[(d[i]^2)[t], {i, 0, n}] == 1, {i, 1}];
eqns = Flatten[Join[eqn1, eqn2, eqn3, eqn4]];
bcs = {{c[0][0] == d[0][0] == 0.004567},
Table[c[i][0] == 0.0000000034, {i, 1, 5}],
Table[d[i][0] == 0.0000000034, {i, 1, 5}]};
var = Join[Table[c[i], {i, 0, 5}], Table[d[i], {i, 0, n}]];
sol = NDSolve[{eqns, bcs}, var, {t, 1/10}]
NDSolve::dvnoarg:
I can not understand the meaning of this kind of error,,
can anyone explain and give some valuable insight.
regards,
tarun