MathGroup Archive 2005

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

Search the Archive

Why am I getting this error?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg60341] Why am I getting this error?
  • From: ksmath <kdsimon at uiuc.edu>
  • Date: Tue, 13 Sep 2005 06:07:51 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

I am attempting to solve a system with 3 pde's and I keep getting this error even though I think all of my equations have independent variables:

NDSolve::dvnoarg: The function v appears with no arguments.

If anyone could find the mistake in my code I would really appreciate it.  C,V, and R are the dependent variables, and t (time) and x(distance-does not depend on time) are the independent variables.  All other variables are constants.

Thank you for any assistance you can provide me.

\!\(\(rRHS\  = \ \(d\_r\) D[r[t, 
      x], {x, 2}] + β \((1 - \((r[t, x] - r\_base)\)/\((k\_r - 
        r\_base)\))\) \((r[t, x] - r\_base)\);\)\[IndentingNewLine]
  \(vRHS\  = \ \(d\_v\) D[v[t, x], {x, 2}] - Ï?\ v[t, x]\  - 
            η\ v[t, x]\ c[t, x]/\((γ\_1 + v[t, x])\) + If[
    q\_large < q\_min, λ \((q\_large - q\_min)\), 0];\)\[IndentingNewLine]
  \(\(cRHS\  = \ D[c[t, x]\ v[t, x]\ D[v[t, x], x], x] + If[v[
        t, x] â?¥ v\_g, Ï?\ v[t, x]\ c[t, 
        x]\ /\((γ\_2 + 
            v[t, x])\), Ï? \((v[t, x] - 
              v\_d)\) c[t, x]];\)\(\[IndentingNewLine]\)
  \)\[IndentingNewLine]
  \(reqn\  = D[r[t, x], t]\  == \ rRHS;\)\[IndentingNewLine]
  \(veqn\  = \ D[v[t, x], t]\  == \ vRHS;\)\[IndentingNewLine]
  \(ceqn\  = \ D[c[t, x], t]\  == \ cRHS;\)\)

BC1 = Derivative[0, 1][r][t, 0] == 0;
BC2 = Derivative[0, 1][r][t, 1] == 0;
BC3 = Derivative[0, 1][v][t, 0] == 0;
BC4 = Derivative[0, 1][v][t, 1] == 0;
BC5 = c[t, 0] == 1;
BC6 = c[t, 1] v[t, 1] Derivative[0, 1][v][t, 1] == 0;

\!\(\(c\_init = \(c\_o\) Exp \((\(-150\) x\^2)\) == c[0, 
          x];\)\[IndentingNewLine]
  \(r\_init = \((K\_r - 
          r\_base)\) Exp \((\(-100\) x\^2)\) + r\_base == r[0, x];\)\
\[IndentingNewLine]
  \(v\_init = \((v\_g - v\_d)\)/\((2  v\_o)\) == v[0, x];\)\)

\!\(\(solution\  = \ NDSolve[{reqn, \ veqn, \ ceqn, q1\_i,
         q2\_i, \ c\_init, \ v\_init, \ 
          r\_init, \ BC1, \ BC2, \ BC3, \ BC4, \ BC5, \ BC6}, {r[t, x], v[t, 
          x], c[t, x], q1[t, x], q2[t, x]}, {t, 0, 1}, {x, 0, 1}];\)\)

Any help is much appreciated.

Katie


  • Prev by Date: Re: smooth eigenvalues and eigenvectors as a function of frequency
  • Next by Date: Re: Why this function does not return a single value
  • Previous by thread: Re: Draw3D ColorFunction->Hue
  • Next by thread: Re: Why am I getting this error?