回答編集履歴
2
edit
test
CHANGED
@@ -147,3 +147,21 @@
|
|
147
147
|
plt.show()
|
148
148
|
|
149
149
|
```
|
150
|
+
|
151
|
+
|
152
|
+
|
153
|
+
---
|
154
|
+
|
155
|
+
|
156
|
+
|
157
|
+
Runge-Kuttaの式に関しては、
|
158
|
+
|
159
|
+
```python
|
160
|
+
|
161
|
+
k1=y0*dt
|
162
|
+
|
163
|
+
l1=u0*dt
|
164
|
+
|
165
|
+
```
|
166
|
+
|
167
|
+
が多分正しくない気がします。
|
1
edit
test
CHANGED
@@ -1,3 +1,9 @@
|
|
1
|
+
デバッグに伴って、逆にわかりにくくなっているので、落ち着いて書き直したほうが早いかもしれません。
|
2
|
+
|
3
|
+
実装自体が目的と想定して、正しい値を出すであろうscipyの結果との比較を入れました。
|
4
|
+
|
5
|
+
|
6
|
+
|
1
7
|
```python
|
2
8
|
|
3
9
|
import numpy as np
|