RK4モジュールを使って落体の運動方程式。 せっかく計算本体を別モジュールにしたにも関わらず記述量が変わっていないのがなんとも。 program ode use runge_kutta_4 implicit none double precision, parameter :: g = 9.80665 ! gravitational acceleratio…