Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

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


  • Prev by Date: Re: Another point about Mathematica 8.0
  • Next by Date: Re: get rid of I in the result of an integral
  • Previous by thread: Re: Emacs with Mathematica
  • Next by thread: Re: problem with NDSolve::dvnoarg: