Scipyのsolve_ivpで条件を満たすまでシミュレーションを行う - Qiita

概要 scipyライブラリのsolve_ivpのevent機能について試すために自由落下のシミュレーションを行います。 今回は地面に到着した瞬間までシミュレーションを行います。 式 自由落下の式は次のようになります。 ここで、物体の高さを$y$、速度を$v$とします。 $m$と$g$は定数です。 $y=0$のとき地面と衝突したと考えます。 solve_ivpのeventはfloatを返す関数で定義され、eventの発生時に0を返します。 今回はそのまま物体の高さ$y$を返します。 プログラム 上記の式を用いてシミュレーションプログラムを作成します。 from scipy.integrate import solve_ivp import numpy from matplotlib import pyplot # 物理定数 m = 1 g = 1 def free_fall(t, y): '