Python 微分方程式

Question

前提・実現したいこと

mdv/dt = mg - k*v
上の微分方程式をt=100,200の時のvの値をpythonで解きたいのですが、どのようにすればよいか教えていただきたいです。

発生している問題・エラーメッセージ

下のようなエラー？が発生してしまいました。
<ipython-input-31-6b2f8cf73b86>:6: RuntimeWarning: invalid value encountered in true_divide
  return g - (k * y) / m
C:\Users\s317048\Anaconda3\envs\numpy\lib\site-packages\scipy\integrate\odepack.py:248: ODEintWarning: Excess accuracy requested (tolerances too small). Run with full_output = 1 to get quantitative information.
  warnings.warn(warning_msg, ODEintWarning)

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint

def func(y,m,g,k,t):
    return g - (k * y) / m

m = 10
g = 9.8
k = 0.1
y = 0
t = np.arange(0, 100, 0.1)

v = odeint(func,y,t,args=(m,g,k))

Comments

[jeanbiego]

>上の微分方程式をt=100,200の時のvの値
積分型に直して、初期速度0としてt=100,200のときのvの値が知りたいみたいなことですか？

[退会済みユーザー]

初期値はt=0のときv=0とします。またm=10, g=9.8 k=0.1とします。時間間隔（差分時間）は任意とします。t=100とt=200 の時のvの値を求めてください。
\begin{eqnarray}
m\frac {d v}{d t} = mg - kv
\end{eqnarray}

文章的にはこのようになります。

Answer 1

下記でいかがでしょう。

import numpy as np
from scipy.integrate import odeint

def func(y, t, m, g, k):
    dydt = (g - (k * y) / m)
    return dydt
m = 10
g = 9.8
k = 0.1
t_list = np.arange(0.0, 1000, 1)
y_init = 0.0  #初期値
y_list = odeint(func, y_init, t_list, args=(m,g,k))
print(f"v(t=100) is {y_list[100]}, v(t=200) is {y_list[200]}")
# v(t=100) is [619.47815088], v(t=200) is [847.37142731]

Comments

[退会済みユーザー]

ありがとうございます。