Example 6.   Solve   [Graphics:Images/PendulumMod_gr_141.gif],  by using the Runge-Kutta method to solve the system of D. E.'s  
        [Graphics:Images/PendulumMod_gr_142.gif] ,  
and
        [Graphics:Images/PendulumMod_gr_143.gif].    

Solution 6.

Enter the functions  [Graphics:../Images/PendulumMod_gr_144.gif]  and  [Graphics:../Images/PendulumMod_gr_145.gif]  and form the vector function  [Graphics:../Images/PendulumMod_gr_146.gif].  

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


[Graphics:../Images/PendulumMod_gr_148.gif]
[Graphics:../Images/PendulumMod_gr_149.gif]
[Graphics:../Images/PendulumMod_gr_150.gif]
[Graphics:../Images/PendulumMod_gr_151.gif]

Compute the Runge-Kutta solutions.
First, graph four closed loops around the origin.

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

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

 

 

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

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

 

 

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

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

 

 

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

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

 

 

Combine the above graphs.

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

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

 

 

Second, graph two curves where the initial velocity is large enough to make the pendulum continue turning about its pivot.

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

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

 

 

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

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

 

 

Combine the above graphs.

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

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

 

 

Third, graph two curves where the initial velocity is large enough to make the pendulum continue turning about its pivot but in the opposite direction.
Remark.  In order to keep these curves from plotting in quadrant III, we have started them at [Graphics:../Images/PendulumMod_gr_168.gif] instead of [Graphics:../Images/PendulumMod_gr_169.gif].

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

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

 

 

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

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

 

 

Combine the above graphs.

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

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

 

 

Next, graph four closed loops around the point [Graphics:../Images/PendulumMod_gr_176.gif].

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

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

 

 

Combine the above graphs.

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

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

 

 

Next, graph four closed loops around the point [Graphics:../Images/PendulumMod_gr_181.gif].

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

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

 

 

Combine the above graphs.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004