回答編集履歴
1
グラフを表示できるようにコードを改善
answer
CHANGED
|
@@ -26,7 +26,10 @@
|
|
|
26
26
|
|
|
27
27
|
これをpythonでコード化します。
|
|
28
28
|
|
|
29
|
+
依頼を受けてグラフを表示できるように修正しました。
|
|
30
|
+
|
|
29
31
|
```python
|
|
32
|
+
import matplotlib.pyplot as plt
|
|
30
33
|
|
|
31
34
|
r = 0.01
|
|
32
35
|
k = 30.0
|
|
@@ -34,12 +37,17 @@
|
|
|
34
37
|
dt = 1.0
|
|
35
38
|
|
|
36
39
|
t = 0
|
|
40
|
+
x = [t]
|
|
41
|
+
y = [n]
|
|
37
|
-
while t <=
|
|
42
|
+
while t <= 1000:
|
|
38
43
|
next_n = n + dt * r * ( 1.0 - n / k ) * n
|
|
44
|
+
t = t + 1
|
|
45
|
+
x.append(t)
|
|
39
|
-
|
|
46
|
+
y.append(next_n)
|
|
40
47
|
n = next_n
|
|
41
|
-
t = t + 1
|
|
42
48
|
|
|
49
|
+
plt.plot(x,y)
|
|
50
|
+
plt.show()
|
|
43
51
|
```
|
|
44
52
|
|
|
45
53
|
このプログラムで計算すると、900秒ぐらいでnが30に収束することがわかります。
|