Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: Could someone else verify that an example from the Numerical

• To: mathgroup at smc.vnet.net
• Subject: [mg127321] Re: Could someone else verify that an example from the Numerical
• From: "Nasser M. Abbasi" <nma at 12000.org>
• Date: Mon, 16 Jul 2012 05:00:08 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: mathgroup-newout@smc.vnet.net
• Delivered-to: mathgroup-newsend@smc.vnet.net
• References: <jttv31$gut$1@smc.vnet.net>
• Reply-to: nma at 12000.org

On 7/15/2012 3:30 AM, W Craig Carter wrote:

<snip>

Ah! But you changed the code from the pdf file !

You wrote

> GraphicsGrid[Partition[Table[
>     NDSolve`Iterate[state, tau];
>     Plot3D[
>      Evaluate[
>       u[tau, x, y] /.
>        NDSolve`ProcessSolutions[state, "Forward"]], {x, -L, L}, {y, -L,
>        L}, PlotRange -> {-1/4, 1/4}],
>     {tau, 0, 20, 5}], 2]]
>

Notice the limit you have on the tau. You wrote

       {tau, 0, 20, 5}

while it should be

       {tau, 0, 20., 5.}

This is very important here, due to the replacement rule above
that wants to pattern match a solution obtained from NDSolve`ProcessSolutions
which will be numeric! Hence, what happened is that no match was
found and that is why the plots were coming out empty.

If you just change this line and put a "." there, it will work :)

regards,
--Nasser

• Follow-Ups:
  • Re: Could someone else verify that an example from the Numerical
    • From: W Craig Carter <ccarter@MIT.EDU>