F.A. Harris

Jan. 12, 2000

1.) Problem 3.1, parts a, b, c, and e. See SUB Euler_Richardson (below) for part c. For part c, also compare your result with the analytical result. For part e, get a copy of plot_axis from your instructor. Use this as a beginning of a library (See Appendix 3A). NOTE: Your library file must be in the same directory as your True BASIC program.

SUB Euler_Richardson(y,v,a,g,t,dt)

LET a = -g

LET vmid = v + a*dt/2 ! velocity at midpoint

LET ymid = y + v*dt/2 ! position at midpoint

LET amid = -g ! acceleration at midpoint

LET v = v + amid*dt

LET y = y + vmid*dt

LET t = t + dt

END SUB

2.) Problem 3.3, part a only. For help on modifying your program, see PROGRAM Drag. Do not use PROGRAM Drag directly because it does something special with t0.

SUMMARY: For this week's summary, describe the numerical solution of a freely falling body near the surface of the earth, the three algorithms discussed in class, and their relevance to this problem.