MathGroup Archive 2012

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

Search the Archive

Could someone else verify that an example from the Numerical Solutions to Differential Equations isn't working.

I am looking at the example on page 348 of the tutorial for a numerical solution to the Sine-Gordon equation.
It doesn't work for me for M8.0.4, MacBook, MacOs 10.7.4. Can someone else verify that this doesn't work?

Here is the example:

L = 10; (*the tutorial has L=-10, but I think that is a typo. In any case, the example doesn't work with + or - *)

state = First[NDSolve`ProcessEquations[{D[u[t, x, y], t, t] ==
     D[u[t, x, y], x, x] + D[u[t, x, y], y, y] - Sin[u[t, x, y]],
    u[0, x, y] == Exp[-(x^2 + y^2)],
    Derivative[1, 0, 0][u][0, x, y] == 0,
    u[t, -L, y] == u[t, L, y], u[t, x, -L] == u[t, x, L]}, u,
   t, {x, -L, L},
   {y, -L, L},
   Method -> {"MethodOfLines",
     "SpatialDiscretization" -> {"TensorProductGrid",
       "DifferenceOrder" -> "Pseudospectral"}}]]

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

Thanks, Craig

W Craig Carter
Professor of Materials Science, MIT

  • Prev by Date: Re: Quit and Restart kernel quickly?
  • Next by Date: Re: Algorithm Analysis Course: Should I use Mathematica for projects?
  • Previous by thread: Re: Using Mathematica to Plot Business Business Data By
  • Next by thread: Re: Could someone else verify that an example from the