落体の運動方程式をEuler法により数値的に解くFORTRANプログラム。 program ode implicit none double precision, parameter :: g = 9.80665 ! gravitational acceleration double precision, parameter :: t0 = 0 ! initial t double precision, parameter …