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