Example 1.  Solve the second order  I.V.P.  

    [Graphics:Images/PendulumMod_gr_6.gif]  with  [Graphics:Images/PendulumMod_gr_7.gif],   [Graphics:Images/PendulumMod_gr_8.gif].    

Use the Runge-Kutta method to compute the solution over the interval
[Graphics:Images/PendulumMod_gr_9.gif].  

Solution 1.

The D. E. can be written as   [Graphics:../Images/PendulumMod_gr_10.gif]  where  [Graphics:../Images/PendulumMod_gr_11.gif].  
We use the substitution [Graphics:../Images/PendulumMod_gr_12.gif]  and make the second order I.V.P. into a system of first order I.V.P.'s  

    [Graphics:../Images/PendulumMod_gr_13.gif],  
and  
    [Graphics:../Images/PendulumMod_gr_14.gif],  
where  
    [Graphics:../Images/PendulumMod_gr_15.gif],  
and  
     [Graphics:../Images/PendulumMod_gr_16.gif],   [Graphics:../Images/PendulumMod_gr_17.gif].  

 

[Graphics:../Images/PendulumMod_gr_18.gif]
[Graphics:../Images/PendulumMod_gr_19.gif] [Graphics:../Images/PendulumMod_gr_20.gif] [Graphics:../Images/PendulumMod_gr_21.gif] [Graphics:../Images/PendulumMod_gr_22.gif]

Compute the Runge-Kutta solution.

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

The solution we seek is the first coordinate in the 2D system.

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

Now we can plot the solution.

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

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

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2003