回答編集履歴
2
コード修正
answer
CHANGED
|
@@ -8,7 +8,7 @@
|
|
|
8
8
|
|
|
9
9
|
xs = np.linspace(0, 10, 11)
|
|
10
10
|
f = lambda z: kv(5/3,z)
|
|
11
|
-
fs = [quad(f,x,np.inf) for x in xs]
|
|
11
|
+
fs = [quad(f,x,np.inf)[0]*x for x in xs]
|
|
12
12
|
|
|
13
13
|
fig = plt.figure()
|
|
14
14
|
ax = fig.add_subplot(1,1,1)
|
|
@@ -18,4 +18,4 @@
|
|
|
18
18
|
ax.set_ylabel('F')
|
|
19
19
|
plt.show()
|
|
20
20
|
```
|
|
21
|
-
|
1
コード修正
answer
CHANGED
|
@@ -9,6 +9,13 @@
|
|
|
9
9
|
xs = np.linspace(0, 10, 11)
|
|
10
10
|
f = lambda z: kv(5/3,z)
|
|
11
11
|
fs = [quad(f,x,np.inf) for x in xs]
|
|
12
|
+
|
|
13
|
+
fig = plt.figure()
|
|
14
|
+
ax = fig.add_subplot(1,1,1)
|
|
12
|
-
|
|
15
|
+
ax.plot(xs,fs)
|
|
16
|
+
ax.set_title('test')
|
|
17
|
+
ax.set_xlabel('x')
|
|
18
|
+
ax.set_ylabel('F')
|
|
19
|
+
plt.show()
|
|
13
20
|
```
|
|
14
|
-
|