古典的ルンゲ＝クッタ法のPythonプログラムをWikiPediaの記述と同じ変数名で書いてみた - あらきけいすけのメモ帳

教育用の覚書 Python だと汎用の Runge-Kutta 法のルーチンが６行で書けてしまう*1。以下のコードではルンゲクッタ法ルーチン rk4() を書き換えることなく、1階, 2階, 3階の常微分方程式(ordinary differential equation(s), ODE)を解いている*2*3。関数 rk4() で使っている変数名はRunge–Kutta methods - Wikipediaに揃えた。C言語*4で数値計算をコーディング/教育する気力が一気に失われてしまったｗｗｗｗｗｗ 次のコンプリートコードは次の微分方程式を積分して、あらかじめ求めておいた解の関数の値と比較している*5： class getExpt(): [1階] (比較する解：); class getSint(): [2階] (比較する解：); [1階連立ODE] class getExptSint()

[arakik10]

Runge-Kutta法をPythonで書いてみたらWikipediaの記述どおりの式を書いて6行(インターフェース1行,本体4行,戻り値1行)で済んでしまった。Python以外でコーディングする気力が失せてしまうｗｗｗｗｗｗ

[samu_i]

これでRungeKuttaをかく演習は駆逐されるね。やったね

[dorokei]

めっちゃ懐かしいルンゲクッタ。なんなのかは思い出せない