ルンゲクッタ法で2階連立微分方程式を解く - きしだのHatena今回は、空気抵抗があるときの物体の軌跡を描いてみました。下のプログラムでは、実際にはアニメーションさせてます。 ※追記 2023/4/30 動画を投稿しました 空気抵抗がある中でモノを投げるシミュレーションhttps://t.co/zd6SCX8qE1 pic.twitter.com/IM2HbKomzc— きしだൠ(K1S) (@kis) 2023年4月30日 2階以上の微分方程式は、 を、 と置くことで という1階連立微分方程式に分解して、ルンゲクッタ法を使います。 このとき物体の位置x,yの微分方程式は、空気抵抗の係数をT、重力加速度をgとすると となります。 ここで と置くことで のような１階微分方程式に変形して計算を行います。 プログラム中では、funcw・funcvメソッドが変形後の式にあたります。 こうやって2階微分方程式が解けるようになると、いろいろな物理現象がシミュレート
