Example
6. Solve
, by
using the Runge-Kutta method to solve the system of D.
E.'s
,
and
.
Solution 6.
Enter the functions
and
and
form the vector function
.
Compute the Runge-Kutta solutions.
First, graph four closed loops around the origin.
![[Graphics:../Images/PendulumMod_gr_153.gif]](../Images/PendulumMod_gr_153.gif)
![[Graphics:../Images/PendulumMod_gr_155.gif]](../Images/PendulumMod_gr_155.gif)
![[Graphics:../Images/PendulumMod_gr_157.gif]](../Images/PendulumMod_gr_157.gif)
![[Graphics:../Images/PendulumMod_gr_159.gif]](../Images/PendulumMod_gr_159.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_161.gif]](../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_163.gif]](../Images/PendulumMod_gr_163.gif)
![[Graphics:../Images/PendulumMod_gr_165.gif]](../Images/PendulumMod_gr_165.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_167.gif]](../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
instead of
.
![[Graphics:../Images/PendulumMod_gr_171.gif]](../Images/PendulumMod_gr_171.gif)
![[Graphics:../Images/PendulumMod_gr_173.gif]](../Images/PendulumMod_gr_173.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_175.gif]](../Images/PendulumMod_gr_175.gif)
Next, graph four closed loops around the point
.
![[Graphics:../Images/PendulumMod_gr_178.gif]](../Images/PendulumMod_gr_178.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_180.gif]](../Images/PendulumMod_gr_180.gif)
Next, graph four closed loops around the point
.
![[Graphics:../Images/PendulumMod_gr_183.gif]](../Images/PendulumMod_gr_183.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_185.gif]](../Images/PendulumMod_gr_185.gif)
(c) John H. Mathews 2004