Re: NDSolve

*To*: mathgroup at smc.vnet.net*Subject*: [mg66002] Re: NDSolve*From*: "D.C.Grady at gmail.com" <D.C.Grady at gmail.com>*Date*: Wed, 26 Apr 2006 04:38:16 -0400 (EDT)*References*: <e2kq0u$77u$1@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

Patrick, When you use NDSolve to generate a solution, Mathematica will give you a replacement rule rather than assigning the solution directly to the variable name you're using. You can, for example, plot your solution by doing something like this: sol = NDSolve[ equations, Q, {x, 0.03, 0.06}, {t, 0.001, 1}] Plot3D[Q[x,t] /. sol, {x, 0.03, 0.06}, {t, 0.001, 1}] The documentation for the NDSolve command has many examples. In the code you provide, you twice write D[S, x]. Since S is itself an expression in terms of x, should this be D[S, t]?