Example 1.  An arrow is shot upward from the origin with an initial velocity of  300 ft/sec.  Assume that there is no air resistance and use the model  
        [Graphics:Images/ProjectileMotionMod_gr_38.gif].  
Find the velocity and position as a function of time.  Find the ascent time, the descent time, maximum height, and the impact velocity.

Solution 1.

First, compute the solution using the Runge-Kutta method for second order D.E.'s.

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



[Graphics:../Images/ProjectileMotionMod_gr_40.gif]
[Graphics:../Images/ProjectileMotionMod_gr_41.gif]
[Graphics:../Images/ProjectileMotionMod_gr_42.gif]
[Graphics:../Images/ProjectileMotionMod_gr_43.gif]
[Graphics:../Images/ProjectileMotionMod_gr_44.gif]
[Graphics:../Images/ProjectileMotionMod_gr_45.gif]
[Graphics:../Images/ProjectileMotionMod_gr_46.gif]
[Graphics:../Images/ProjectileMotionMod_gr_47.gif]
[Graphics:../Images/ProjectileMotionMod_gr_48.gif]
[Graphics:../Images/ProjectileMotionMod_gr_49.gif]
[Graphics:../Images/ProjectileMotionMod_gr_50.gif]

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

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

Now we can plot the solution.

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


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

[Graphics:../Images/ProjectileMotionMod_gr_54.gif]
[Graphics:../Images/ProjectileMotionMod_gr_55.gif]
[Graphics:../Images/ProjectileMotionMod_gr_56.gif]
[Graphics:../Images/ProjectileMotionMod_gr_57.gif]
[Graphics:../Images/ProjectileMotionMod_gr_58.gif]

Second, compute the analytic solution to the D.E.'s.

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


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

[Graphics:../Images/ProjectileMotionMod_gr_61.gif]
[Graphics:../Images/ProjectileMotionMod_gr_62.gif]
[Graphics:../Images/ProjectileMotionMod_gr_63.gif]
[Graphics:../Images/ProjectileMotionMod_gr_64.gif]
[Graphics:../Images/ProjectileMotionMod_gr_65.gif]
[Graphics:../Images/ProjectileMotionMod_gr_66.gif]
[Graphics:../Images/ProjectileMotionMod_gr_67.gif]
[Graphics:../Images/ProjectileMotionMod_gr_68.gif]
[Graphics:../Images/ProjectileMotionMod_gr_69.gif]
[Graphics:../Images/ProjectileMotionMod_gr_70.gif]
[Graphics:../Images/ProjectileMotionMod_gr_71.gif]
[Graphics:../Images/ProjectileMotionMod_gr_72.gif]

Compare the Runge-Kutta solution with the analytic solution.

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


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

[Graphics:../Images/ProjectileMotionMod_gr_75.gif]
[Graphics:../Images/ProjectileMotionMod_gr_76.gif]

Notice that the maximum altitude will occur when the time is near t = 9,
and the arrow will hit the ground when the time is near t = 19.

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


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

[Graphics:../Images/ProjectileMotionMod_gr_79.gif]
[Graphics:../Images/ProjectileMotionMod_gr_80.gif]
[Graphics:../Images/ProjectileMotionMod_gr_81.gif]
[Graphics:../Images/ProjectileMotionMod_gr_82.gif]
[Graphics:../Images/ProjectileMotionMod_gr_83.gif]
[Graphics:../Images/ProjectileMotionMod_gr_84.gif]
[Graphics:../Images/ProjectileMotionMod_gr_85.gif]
[Graphics:../Images/ProjectileMotionMod_gr_86.gif]
[Graphics:../Images/ProjectileMotionMod_gr_87.gif]

In this model, the ascent time is equal to the descent time.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004