回答編集履歴
1
コード修正
answer
CHANGED
|
@@ -3,16 +3,15 @@
|
|
|
3
3
|
import numpy as np
|
|
4
4
|
from scipy.integrate import odeint
|
|
5
5
|
|
|
6
|
-
def func(y, t):
|
|
6
|
+
def func(y, t, m, g, k):
|
|
7
|
-
m = 10
|
|
8
|
-
g = 9.8
|
|
9
|
-
k = 0.1
|
|
10
|
-
dydt = g - (k * y) / m
|
|
7
|
+
dydt = (g - (k * y) / m)
|
|
11
8
|
return dydt
|
|
12
|
-
|
|
9
|
+
m = 10
|
|
10
|
+
g = 9.8
|
|
11
|
+
k = 0.1
|
|
13
12
|
t_list = np.arange(0.0, 1000, 1)
|
|
14
13
|
y_init = 0.0 #初期値
|
|
15
|
-
y_list = odeint(func, y_init, t_list)
|
|
14
|
+
y_list = odeint(func, y_init, t_list, args=(m,g,k))
|
|
16
15
|
print(f"v(t=100) is {y_list[100]}, v(t=200) is {y_list[200]}")
|
|
17
16
|
# v(t=100) is [619.47815088], v(t=200) is [847.37142731]
|
|
18
17
|
```
|