Runge-Kutta-Gill法の各ステップの係数のチートシートを探していて次の論文に行き当たったのでメモ。 CiNii 論文 - Runge-Kutta-Gill法について 伊理 正夫, 松谷 泰行, 情報処理, 8(2), pp.103-107 (1967). ソースコードを見ると k=h*f(i) qi=q(i) r=0.5*K-q…