質問編集履歴
1
不明点,エラーコードを追記しました
test
CHANGED
File without changes
|
test
CHANGED
@@ -1,4 +1,4 @@
|
|
1
|
-
###
|
1
|
+
### 質問内容
|
2
2
|
|
3
3
|
|
4
4
|
|
@@ -6,7 +6,11 @@
|
|
6
6
|
|
7
7
|
r,θを変数とし,各座標の値を代入して,その場の速度を算出してベクトル図を描ければと考えましたが上手くいきませんでした.
|
8
8
|
|
9
|
+
|
10
|
+
|
11
|
+
下記エラーコードが出たため,その対策も調べ確認してみましたが,上手くいかずご質問させていただきました.
|
12
|
+
|
9
|
-
|
13
|
+
アドバイスなど頂けるとありがたいです.
|
10
14
|
|
11
15
|
|
12
16
|
|
@@ -225,3 +229,69 @@
|
|
225
229
|
plt.show()
|
226
230
|
|
227
231
|
```
|
232
|
+
|
233
|
+
|
234
|
+
|
235
|
+
### エラーコード
|
236
|
+
|
237
|
+
|
238
|
+
|
239
|
+
```
|
240
|
+
|
241
|
+
ValueError Traceback (most recent call last)
|
242
|
+
|
243
|
+
<ipython-input-1-ee1d8f65e98c> in <module>
|
244
|
+
|
245
|
+
86
|
246
|
+
|
247
|
+
87 A_M = M(lambda r: r2 / (r * cmath.sin(beta2)), r2, r1)
|
248
|
+
|
249
|
+
---> 88 A_F_rc = F_rc(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (r / radius)**(n + 1), radius, r3)
|
250
|
+
|
251
|
+
89 A_G_rc = G_rc(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (radius / r)**(n - 1), r3, radius)
|
252
|
+
|
253
|
+
90 A_F_r2 = F_r2(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (r / r2)**(n + 1), r3, r2)
|
254
|
+
|
255
|
+
|
256
|
+
|
257
|
+
<ipython-input-1-ee1d8f65e98c> in F_rc(func, a, b, **kwargs)
|
258
|
+
|
259
|
+
22 def imag_func(r):
|
260
|
+
|
261
|
+
23 return scipy.imag(func(r))
|
262
|
+
|
263
|
+
---> 24 real_integral = quad(real_func, a, b, **kwargs)
|
264
|
+
|
265
|
+
25 imag_integral = quad(imag_func, a, b, **kwargs)
|
266
|
+
|
267
|
+
26 return (real_integral[0] + 1j*imag_integral[0], real_integral[1:], imag_integral[1:])
|
268
|
+
|
269
|
+
|
270
|
+
|
271
|
+
~/opt/anaconda3/lib/python3.7/site-packages/scipy/integrate/quadpack.py in quad(func, a, b, args, full_output, epsabs, epsrel, limit, points, weight, wvar, wopts, maxp1, limlst)
|
272
|
+
|
273
|
+
336
|
274
|
+
|
275
|
+
337 # check the limits of integration: \int_a^b, expect a < b
|
276
|
+
|
277
|
+
--> 338 flip, a, b = b < a, min(a, b), max(a, b)
|
278
|
+
|
279
|
+
339
|
280
|
+
|
281
|
+
340 if weight is None:
|
282
|
+
|
283
|
+
|
284
|
+
|
285
|
+
ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()
|
286
|
+
|
287
|
+
|
288
|
+
|
289
|
+
```
|
290
|
+
|
291
|
+
|
292
|
+
|
293
|
+
### 出力結果
|
294
|
+
|
295
|
+
|
296
|
+
|
297
|
+
![出力結果](1439a0726d7a092ba158827903ae8c6b.png)
|