Example 7.   Solve  [Graphics:Images/PendulumMod_gr_186.gif],  by using the Runge-Kutta method to solve the system of D. E.'s  
        [Graphics:Images/PendulumMod_gr_187.gif] ,  
and
        [Graphics:Images/PendulumMod_gr_188.gif].    

Solution 7.

Enter the functions  [Graphics:../Images/PendulumMod_gr_189.gif]  and  [Graphics:../Images/PendulumMod_gr_190.gif]  and form the vector function  [Graphics:../Images/PendulumMod_gr_191.gif].  

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


[Graphics:../Images/PendulumMod_gr_193.gif]
[Graphics:../Images/PendulumMod_gr_194.gif]
[Graphics:../Images/PendulumMod_gr_195.gif]
[Graphics:../Images/PendulumMod_gr_196.gif]

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

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

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

 

 

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

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

 

 

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

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

 

 

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

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

 

 

Combine the above graphs.

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

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

 

 

Second, graph a curves where the initial velocity is large enough to make the pendulum turn twice about its pivot.

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

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

 

 

Third, graph a curves where the initial velocity is large enough to make the pendulum turn three about its pivot.

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

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

 

 

Combine the above graphs.

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

[Graphics:../Images/PendulumMod_gr_212.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_213.gif] instead of [Graphics:../Images/PendulumMod_gr_214.gif].

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

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

 

 

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

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

 

 

Combine the above graphs.

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

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

 

 

Next, graph four spirals around the point [Graphics:../Images/PendulumMod_gr_221.gif].

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

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

 

 

Combine the above graphs.

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

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

 

 

Next, 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_226.gif] instead of [Graphics:../Images/PendulumMod_gr_227.gif], and the time interval is larger than before.

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

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

 

 

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

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

 

 

Combine the above graphs.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004