質問編集履歴
1
コードを追加しました。
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -12,10 +12,30 @@
|
|
|
12
12
|
|
|
13
13
|
### 該当のソースコード
|
|
14
14
|
|
|
15
|
-
```
|
|
15
|
+
```python
|
|
16
16
|
ソースコード
|
|
17
|
-
```
|
|
17
|
+
```import numpy as np
|
|
18
18
|
|
|
19
|
+
def walk(n_walker=10, nt=10):
|
|
20
|
+
pos = np.zeros((nt+1,n_walker),dtype=np.int32)
|
|
21
|
+
irand = 2*(np.random.randint(0,2,(nt,n_walker)))-1
|
|
22
|
+
if irand==-1:
|
|
23
|
+
pos=pos-1
|
|
24
|
+
else:
|
|
25
|
+
pos=pos+1
|
|
26
|
+
for walker in range(n_walker):
|
|
27
|
+
for step in range(nt):
|
|
28
|
+
|
|
29
|
+
|
|
30
|
+
return pos(sum(step)for step in n_waker)
|
|
31
|
+
|
|
32
|
+
nt = 500
|
|
33
|
+
n_walker = 1000
|
|
34
|
+
pos = walk(n_walker=n_walker,nt=nt)
|
|
35
|
+
print('Number of walkers',n_walker)
|
|
36
|
+
print('Average position:',pos.mean(n_walker))
|
|
37
|
+
print('Variance:',pos.var(n_walker))
|
|
38
|
+
|
|
19
39
|
### 試したこと
|
|
20
40
|
|
|
21
41
|
forやifを使ってみたのですが、エラーが出てしまいできませんでした。
|