Example 4.  Solve the Lorenz I. V. P.
        [Graphics:Images/LorenzAttractorMod_gr_183.gif]  
        [Graphics:Images/LorenzAttractorMod_gr_184.gif]  
    
    [Graphics:Images/LorenzAttractorMod_gr_185.gif]  
Use Mathematica's NDSolve procedure.  

Solution 4.

Caution. The syntax must be followed carefully.  The symbol "[Graphics:../Images/LorenzAttractorMod_gr_186.gif] is the boolean equal, you can use two ordinary equal signs next to each other if you prefer [Graphics:../Images/LorenzAttractorMod_gr_187.gif]".  

Mathematica's numerical D.E. solver.
  It is more robust than the Runge-Kutta method, we can solve the D.E.'s over the larger interval  [Graphics:../Images/LorenzAttractorMod_gr_188.gif].  

Caution.  If you do not use the option  "[Graphics:../Images/LorenzAttractorMod_gr_189.gif]" then Mathematica will give you the error message:

NDSolve::"mxst": "Maximum number of 1000 steps reached at the point t == 11.146875180780313`."

After you solve the problem successfully, you should make copies of the commands and delete the MaxSteps option so you can see "numerical analysis in action."  (Or the lack of it, because no one wants to report a bad answer.)

[Graphics:../Images/LorenzAttractorMod_gr_190.gif]
[Graphics:../Images/LorenzAttractorMod_gr_191.gif]



[Graphics:../Images/LorenzAttractorMod_gr_192.gif]

[Graphics:../Images/LorenzAttractorMod_gr_193.gif]

[Graphics:../Images/LorenzAttractorMod_gr_194.gif]

[Graphics:../Images/LorenzAttractorMod_gr_195.gif]

[Graphics:../Images/LorenzAttractorMod_gr_196.gif]

[Graphics:../Images/LorenzAttractorMod_gr_197.gif]

[Graphics:../Images/LorenzAttractorMod_gr_198.gif]

[Graphics:../Images/LorenzAttractorMod_gr_199.gif]

 

 

[Graphics:../Images/LorenzAttractorMod_gr_200.gif]

[Graphics:../Images/LorenzAttractorMod_gr_201.gif]

[Graphics:../Images/LorenzAttractorMod_gr_202.gif]
[Graphics:../Images/LorenzAttractorMod_gr_203.gif]
[Graphics:../Images/LorenzAttractorMod_gr_204.gif]
[Graphics:../Images/LorenzAttractorMod_gr_205.gif]
[Graphics:../Images/LorenzAttractorMod_gr_206.gif]

[Graphics:../Images/LorenzAttractorMod_gr_207.gif]

[Graphics:../Images/LorenzAttractorMod_gr_208.gif]

[Graphics:../Images/LorenzAttractorMod_gr_209.gif]

 

 

[Graphics:../Images/LorenzAttractorMod_gr_210.gif]

[Graphics:../Images/LorenzAttractorMod_gr_211.gif]

[Graphics:../Images/LorenzAttractorMod_gr_212.gif]
[Graphics:../Images/LorenzAttractorMod_gr_213.gif]
[Graphics:../Images/LorenzAttractorMod_gr_214.gif]
[Graphics:../Images/LorenzAttractorMod_gr_215.gif]
[Graphics:../Images/LorenzAttractorMod_gr_216.gif]

[Graphics:../Images/LorenzAttractorMod_gr_217.gif]

[Graphics:../Images/LorenzAttractorMod_gr_218.gif]

[Graphics:../Images/LorenzAttractorMod_gr_219.gif]

 

 

[Graphics:../Images/LorenzAttractorMod_gr_220.gif]

[Graphics:../Images/LorenzAttractorMod_gr_221.gif]

[Graphics:../Images/LorenzAttractorMod_gr_222.gif]
[Graphics:../Images/LorenzAttractorMod_gr_223.gif]
[Graphics:../Images/LorenzAttractorMod_gr_224.gif]
[Graphics:../Images/LorenzAttractorMod_gr_225.gif]
[Graphics:../Images/LorenzAttractorMod_gr_226.gif]

[Graphics:../Images/LorenzAttractorMod_gr_227.gif]

[Graphics:../Images/LorenzAttractorMod_gr_228.gif]

[Graphics:../Images/LorenzAttractorMod_gr_229.gif]

 

 

[Graphics:../Images/LorenzAttractorMod_gr_230.gif]

[Graphics:../Images/LorenzAttractorMod_gr_231.gif]

[Graphics:../Images/LorenzAttractorMod_gr_232.gif]
[Graphics:../Images/LorenzAttractorMod_gr_233.gif]
[Graphics:../Images/LorenzAttractorMod_gr_234.gif]
[Graphics:../Images/LorenzAttractorMod_gr_235.gif]
[Graphics:../Images/LorenzAttractorMod_gr_236.gif]

[Graphics:../Images/LorenzAttractorMod_gr_237.gif]

[Graphics:../Images/LorenzAttractorMod_gr_238.gif]

[Graphics:../Images/LorenzAttractorMod_gr_239.gif]

 

 

 Return to Numerical Methods - Numerical Analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004