回答編集履歴
1
コード修正
test
CHANGED
|
@@ -8,25 +8,23 @@
|
|
|
8
8
|
|
|
9
9
|
|
|
10
10
|
|
|
11
|
-
def func(y, t):
|
|
11
|
+
def func(y, t, m, g, k):
|
|
12
12
|
|
|
13
|
-
m = 10
|
|
14
|
-
|
|
15
|
-
g = 9.8
|
|
16
|
-
|
|
17
|
-
k = 0.1
|
|
18
|
-
|
|
19
|
-
dydt = g - (k * y) / m
|
|
13
|
+
dydt = (g - (k * y) / m)
|
|
20
14
|
|
|
21
15
|
return dydt
|
|
22
16
|
|
|
17
|
+
m = 10
|
|
23
18
|
|
|
19
|
+
g = 9.8
|
|
20
|
+
|
|
21
|
+
k = 0.1
|
|
24
22
|
|
|
25
23
|
t_list = np.arange(0.0, 1000, 1)
|
|
26
24
|
|
|
27
25
|
y_init = 0.0 #初期値
|
|
28
26
|
|
|
29
|
-
y_list = odeint(func, y_init, t_list)
|
|
27
|
+
y_list = odeint(func, y_init, t_list, args=(m,g,k))
|
|
30
28
|
|
|
31
29
|
print(f"v(t=100) is {y_list[100]}, v(t=200) is {y_list[200]}")
|
|
32
30
|
|