| Author |
Comment/Response |
Michael
|
 12/02/12 09:22am
Put the initialization of constants first:
imagesize = 350; n = 2; xi = 3; ti = 0; ts = 1000; m = .3; tp = 2;
l = -1; ph = 1.9; phd = 1.7; \[Kappa] = 2
and I would avoid SetDelayed (:=) and use an ordinary =.
Next change the typos ai to xi in the last cell.
The Method names are usually supposed to be strings, but they make work as symbols. I changed it to
Method -> {"StiffnessSwitching",
Method -> {"ExplicitRungeKutta", Automatic}}
True is not a Graphics directive, so change it to one or to {}.
Then it all works with only a couple of warnings.
NDSolve often gives warnings, because DEs can have all sorts of mathematical difficulties. You can turn them off with the function Quiet[], but you probably want them. You might need to check if they are indicating an intrinsic problem with the system or a problem with your setup.
URL: , |
|
