RK4を使って一階連立微分方程式を解くためのモジュール。 サブルーチンset_funcでのを設定し、 関数nextで時刻tの値から時刻t+hの値を算出する。 FORTRAN 2003以降でないとコンパイルできないと思う。 module runge_kutta_4 implicit none interface functio…