回答編集履歴
1
追記
answer
CHANGED
|
@@ -24,4 +24,27 @@
|
|
|
24
24
|
```python3
|
|
25
25
|
# 描画、これが重い
|
|
26
26
|
plot(U,(x,-0.005,0.005))
|
|
27
|
-
```
|
|
27
|
+
```
|
|
28
|
+
|
|
29
|
+
|
|
30
|
+
追記
|
|
31
|
+
係数を簡略化して置換積分すると、どうにか計算可能。
|
|
32
|
+
|
|
33
|
+
```python3
|
|
34
|
+
from sympy import Symbol, integrate, cos, sin, sqrt
|
|
35
|
+
x = Symbol('x')
|
|
36
|
+
s = Symbol('s')
|
|
37
|
+
|
|
38
|
+
t = Symbol('t')
|
|
39
|
+
a = 1 # 9.929180321080256*10**6
|
|
40
|
+
b = 0 # 0.0625
|
|
41
|
+
# t = b+(x-s)**2.0
|
|
42
|
+
|
|
43
|
+
dsdt = -1/2/sqrt(t) # 置換積分用 ds/dt
|
|
44
|
+
|
|
45
|
+
expr=cos(a*sqrt(t))/t*dsdt
|
|
46
|
+
S=integrate(expr, t)
|
|
47
|
+
S
|
|
48
|
+
```
|
|
49
|
+
|
|
50
|
+
|