Example 1.  Solve the Lorenz I. V. P.  
        [Graphics:Images/LorenzAttractorMod_gr_17.gif]        and      [Graphics:Images/LorenzAttractorMod_gr_18.gif]  
        [Graphics:Images/LorenzAttractorMod_gr_19.gif]        and      [Graphics:Images/LorenzAttractorMod_gr_20.gif]  
    
    [Graphics:Images/LorenzAttractorMod_gr_21.gif]        and      [Graphics:Images/LorenzAttractorMod_gr_22.gif]  

Solution 1.

Enter the vector function  [Graphics:../Images/LorenzAttractorMod_gr_23.gif].  

[Graphics:../Images/LorenzAttractorMod_gr_24.gif]
[Graphics:../Images/LorenzAttractorMod_gr_25.gif]



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

[Graphics:../Images/LorenzAttractorMod_gr_27.gif]
[Graphics:../Images/LorenzAttractorMod_gr_28.gif]
[Graphics:../Images/LorenzAttractorMod_gr_29.gif]

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

Compute the Runge-Kutta solution. Since this curve is known to twist a lot in 3D space, it will require about 1000 subintervals to get good resolution.

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


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

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

[Graphics:../Images/LorenzAttractorMod_gr_34.gif]
[Graphics:../Images/LorenzAttractorMod_gr_35.gif]
[Graphics:../Images/LorenzAttractorMod_gr_36.gif]

[Graphics:../Images/LorenzAttractorMod_gr_37.gif]
[Graphics:../Images/LorenzAttractorMod_gr_38.gif]

We can plot the projections on the coordinate planes.

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

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

[Graphics:../Images/LorenzAttractorMod_gr_41.gif]
[Graphics:../Images/LorenzAttractorMod_gr_42.gif]


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

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

[Graphics:../Images/LorenzAttractorMod_gr_46.gif]
[Graphics:../Images/LorenzAttractorMod_gr_47.gif]


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

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

[Graphics:../Images/LorenzAttractorMod_gr_51.gif]
[Graphics:../Images/LorenzAttractorMod_gr_52.gif]


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

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

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


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

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

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


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

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

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

We can plot the coordinate functions.

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

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

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


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

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

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


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

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

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


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

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

[Graphics:../Images/LorenzAttractorMod_gr_75.gif]
[Graphics:../Images/LorenzAttractorMod_gr_76.gif]
[Graphics:../Images/LorenzAttractorMod_gr_77.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004