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

MathGroup Archive 2012

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

Search the Archive

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



  • Prev by Date: IsoWeek Function
  • Next by Date: Re: IsoWeek Function
  • Previous by thread: Re: Could someone else verify that an example from the Numerical
  • Next by thread: Re: Could someone else verify that an example from the Numerical