MathGroup Archive 2005

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

Search the Archive

Re: ndsolve error message

  • To: mathgroup at smc.vnet.net
  • Subject: [mg56653] Re: [mg56648] ndsolve error message
  • From: Chris Chiasson <chris.chiasson at gmail.com>
  • Date: Mon, 2 May 2005 01:32:36 -0400 (EDT)
  • References: <200505010713.DAA16130@smc.vnet.net>
  • Reply-to: Chris Chiasson <chris.chiasson at gmail.com>
  • Sender: owner-wri-mathgroup at wolfram.com

Ack!

Apparently the formatting for the script N and the chi character were
messed up as part of the email transmission process...

Here is the input with those characters replaced by NN and XX, respectively:

NDSolve[NNA[z,
        t]\[Equal](NNA[z,t]+NNB[z,t]+NNC[z,t])*
          XXA[z,t]-(0.0003610152163600286*
              Derivative[1,0][XXA][z,
                t])/((0.3040465336986682*XXC[z,t])/(XXA[z,t]+
                    XXB[z,t])+(0.39170805107300855*XXB[z,t])/(XXA[z,t]+
                    XXC[z,t]))&&
    NNB[z,t]\[Equal](NNA[z,t]+NNB[z,t]+NNC[z,t])*
          XXA[z,t]-(0.0003610152163600286*
              Derivative[1,0][XXB][z,
                t])/((0.2905456429313062*XXC[z,t])/(XXA[z,t]+
                    XXB[z,t])+(0.39170805107300855*XXA[z,t])/(XXB[z,t]+
                    XXC[z,t]))&&
    NNC[z,t]\[Equal](NNA[z,t]+NNB[z,t]+NNC[z,t])*
          XXA[z,t]-(0.0003610152163600286*
              Derivative[1,0][XXC][z,
                t])/((0.2905456429313062*XXB[z,t])/(XXA[z,t]+
                    XXC[z,t])+(0.3040465336986682*XXA[z,t])/(XXB[z,t]+
                    XXC[z,t]))&&
    0.0003610152163600286*(XXA[z,t]*Derivative[0,1][P][z,t]+
            P[z,t]*Derivative[0,1][XXA][z,t]+2769.9663468*
              Derivative[1,0][NNA][z,t])\[Equal]0&&
    0.0003610152163600286*(XXB[z,t]*Derivative[0,1][P][z,t]+
            P[z,t]*Derivative[0,1][XXB][z,t]+2769.9663468*
              Derivative[1,0][NNB][z,t])\[Equal]0&&
    0.0003610152163600286*(XXC[z,t]*Derivative[0,1][P][z,t]+
            P[z,t]*Derivative[0,1][XXC][z,t]+2769.9663468*
              Derivative[1,0][NNC][z,t])\[Equal]0&&
    XXA[z,t]+XXB[z,t]+XXC[z,t]\[Equal]1&&
    XXA[z,0]\[Equal]XXA[0,t]\[Equal]0.13238383827902114&&
    XXB[z,0]\[Equal]XXB[0,t]\[Equal]0.6681373176389893&&
    XXC[z,0]\[Equal]XXC[0,t]\[Equal]0.1994788440819895&&
    P[z,0]\[Equal]P[0,t]\[Equal]100000&&
    NNA[1/4000,t]\[Equal]0.011156006364769414&&NNB[1/4000,t]\[Equal]0&&
    NNC[1/4000,t]\[Equal]0,{P,NNA,NNB,NNC,XXA,XXB,XXC},{z,0,1/4000},{t,0,2}]

Sorry about the multiple emails,

On 5/1/05, Chris Chiasson <chris.chiasson at gmail.com> wrote:
> Hi MathGroup gurus,
> 
> I get the following error from my input to NDSolve, but I am not sure
> if I am actually violating the rule it says I am:
> 
> Boundary values may only be specified for one independent variable.
> Initial values may only be specified at one value of the other
> independent variable.
> 
> ***Do you think I have violated the rule? What should I change?***
> 
> the input was as follows: (copy it into Mathematica and hit
> ctrl+shift+n to get a more readable form)
> 
> NDSolve[\[ScriptCapitalN]A[z, t] == (\[ScriptCapitalN]A[z, t] +
> \[ScriptCapitalN]B[z, t] + \[ScriptCapitalN]C[z, t])*÷A[z, t] -
>     (0.0003610152163600286*Derivative[1, 0][÷A][z, t])/
>      ((0.3040465336986682*÷C[z, t])/(÷A[z, t] + ÷B[z, t]) +
>       (0.39170805107300855*÷B[z, t])/(÷A[z, t] + ÷C[z, t])) &&
>   \[ScriptCapitalN]B[z, t] == (\[ScriptCapitalN]A[z, t] +
> \[ScriptCapitalN]B[z, t] + \[ScriptCapitalN]C[z, t])*÷A[z, t] -
>     (0.0003610152163600286*Derivative[1, 0][÷B][z, t])/
>      ((0.2905456429313062*÷C[z, t])/(÷A[z, t] + ÷B[z, t]) +
>       (0.39170805107300855*÷A[z, t])/(÷B[z, t] + ÷C[z, t])) &&
>   \[ScriptCapitalN]C[z, t] == (\[ScriptCapitalN]A[z, t] +
> \[ScriptCapitalN]B[z, t] + \[ScriptCapitalN]C[z, t])*÷A[z, t] -
>     (0.0003610152163600286*Derivative[1, 0][÷C][z, t])/
>      ((0.2905456429313062*÷B[z, t])/(÷A[z, t] + ÷C[z, t]) +
>       (0.3040465336986682*÷A[z, t])/(÷B[z, t] + ÷C[z, t])) &&
>   0.0003610152163600286*(÷A[z, t]*Derivative[0, 1][P][z, t] +
>      P[z, t]*Derivative[0, 1][÷A][z, t] +
>      2769.9663468*Derivative[1, 0][\[ScriptCapitalN]A][z, t]) == 0 &&
>   0.0003610152163600286*(÷B[z, t]*Derivative[0, 1][P][z, t] +
>      P[z, t]*Derivative[0, 1][÷B][z, t] +
>      2769.9663468*Derivative[1, 0][\[ScriptCapitalN]B][z, t]) == 0 &&
>   0.0003610152163600286*(÷C[z, t]*Derivative[0, 1][P][z, t] +
>      P[z, t]*Derivative[0, 1][÷C][z, t] +
>      2769.9663468*Derivative[1, 0][\[ScriptCapitalN]C][z, t]) == 0 &&
>   ÷A[z, t] + ÷B[z, t] + ÷C[z, t] == 1 &&
>   ÷A[z, 0] == ÷A[0, t] == 0.13238383827902114 &&
>   ÷B[z, 0] == ÷B[0, t] == 0.6681373176389893 &&
>   ÷C[z, 0] == ÷C[0, t] == 0.1994788440819895 &&
>   P[z, 0] == P[0, t] == 100000 && \[ScriptCapitalN]A[1/4000, t] ==
>    0.011156006364769414 && \[ScriptCapitalN]B[1/4000, t] == 0 &&
> \[ScriptCapitalN]C[1/4000, t] == 0,
>   {P, \[ScriptCapitalN]A, \[ScriptCapitalN]B, \[ScriptCapitalN]C, ÷A,
> ÷B, ÷C}, {z, 0, 1/4000}, {t, 0, 2}]
> 
> Thank you for your valuable time,
> --
> Chris Chiasson
> http://chrischiasson.com
> Kettering University
> Mechanical Engineering
> Graduate Student
> 1 810 265 3161
> 
> 


-- 
Chris Chiasson
http://chrischiasson.com
Kettering University
Mechanical Engineering
Graduate Student
1 810 265 3161


  • Prev by Date: Re: books on writing packages
  • Next by Date: Re: Plot3D attempt.
  • Previous by thread: ndsolve error message
  • Next by thread: Re: Can't Get a numric answer for this