Example
7. Solve
, by
using the Runge-Kutta method to solve the system of D.
E.'s
,
and
.
Solution 7.
Enter the functions
and
and
form the vector function
.
Compute the Runge-Kutta solutions.
First, graph four spirals around the origin.
![[Graphics:../Images/PendulumMod_gr_198.gif]](../Images/PendulumMod_gr_198.gif)
![[Graphics:../Images/PendulumMod_gr_200.gif]](../Images/PendulumMod_gr_200.gif)
![[Graphics:../Images/PendulumMod_gr_202.gif]](../Images/PendulumMod_gr_202.gif)
![[Graphics:../Images/PendulumMod_gr_204.gif]](../Images/PendulumMod_gr_204.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_206.gif]](../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_208.gif]](../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_210.gif]](../Images/PendulumMod_gr_210.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_212.gif]](../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
instead of
.
![[Graphics:../Images/PendulumMod_gr_216.gif]](../Images/PendulumMod_gr_216.gif)
![[Graphics:../Images/PendulumMod_gr_218.gif]](../Images/PendulumMod_gr_218.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_220.gif]](../Images/PendulumMod_gr_220.gif)
Next, graph four spirals around the point
.
![[Graphics:../Images/PendulumMod_gr_223.gif]](../Images/PendulumMod_gr_223.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_225.gif]](../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
instead of
,
and the time interval is larger than before.
![[Graphics:../Images/PendulumMod_gr_229.gif]](../Images/PendulumMod_gr_229.gif)
![[Graphics:../Images/PendulumMod_gr_231.gif]](../Images/PendulumMod_gr_231.gif)
Combine the above graphs.
![[Graphics:../Images/PendulumMod_gr_233.gif]](../Images/PendulumMod_gr_233.gif)
(c) John H. Mathews 2004