質問編集履歴
11
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -62,9 +62,7 @@
|
|
62
62
|
|
63
63
|
U = []
|
64
64
|
|
65
|
-
for i in range(1000
|
65
|
+
for i in range(1,10001):
|
66
|
-
|
67
|
-
|
68
66
|
|
69
67
|
x = np.random.rand(i)
|
70
68
|
|
@@ -94,11 +92,7 @@
|
|
94
92
|
|
95
93
|
```
|
96
94
|
|
97
|
-
ipython-input-5-b0ad3ca2b10a>:11: RuntimeWarning: invalid value encountered in long_scalars
|
98
|
-
|
99
|
-
p = 4*m/i
|
100
|
-
|
101
|
-
![イメージ説明](
|
95
|
+
![イメージ説明](393b6a17edf302ab22a4c7bce19c547e.png)
|
102
96
|
|
103
97
|
|
104
98
|
|
10
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -62,7 +62,9 @@
|
|
62
62
|
|
63
63
|
U = []
|
64
64
|
|
65
|
-
for i in range(
|
65
|
+
for i in range(10000):
|
66
|
+
|
67
|
+
|
66
68
|
|
67
69
|
x = np.random.rand(i)
|
68
70
|
|
@@ -72,7 +74,7 @@
|
|
72
74
|
|
73
75
|
m = np.sum(kyori <= 1)
|
74
76
|
|
75
|
-
p = 4*m/
|
77
|
+
p = 4*m/i
|
76
78
|
|
77
79
|
U.append(p - np.pi)
|
78
80
|
|
@@ -92,11 +94,11 @@
|
|
92
94
|
|
93
95
|
```
|
94
96
|
|
97
|
+
ipython-input-5-b0ad3ca2b10a>:11: RuntimeWarning: invalid value encountered in long_scalars
|
95
98
|
|
99
|
+
p = 4*m/i
|
96
100
|
|
97
|
-
![イメージ説明](
|
101
|
+
![イメージ説明](a7d4e7433b76550914376d40af356654.png)
|
98
|
-
|
99
|
-
直線になってしまいました。
|
100
102
|
|
101
103
|
|
102
104
|
|
9
説明の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -68,7 +68,7 @@
|
|
68
68
|
|
69
69
|
y = np.random.rand(i)
|
70
70
|
|
71
|
-
kyori = np.sqrt(x*
|
71
|
+
kyori = np.sqrt(x*x+y*y)
|
72
72
|
|
73
73
|
m = np.sum(kyori <= 1)
|
74
74
|
|
8
説明の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -14,7 +14,7 @@
|
|
14
14
|
|
15
15
|
|
16
16
|
|
17
|
-
kyori = np.sqrt(x*
|
17
|
+
kyori = np.sqrt(x*x+y*y)
|
18
18
|
|
19
19
|
m = np.sum(kyori <= 1)
|
20
20
|
|
@@ -97,3 +97,9 @@
|
|
97
97
|
![イメージ説明](4a840cb42dd6fe3481f8276bd230eb8e.png)
|
98
98
|
|
99
99
|
直線になってしまいました。
|
100
|
+
|
101
|
+
|
102
|
+
|
103
|
+
(追記)
|
104
|
+
|
105
|
+
計算式を間違えていましたx*y+y*yをx*x+y*y
|
7
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -58,19 +58,15 @@
|
|
58
58
|
|
59
59
|
import matplotlib.pyplot as plt
|
60
60
|
|
61
|
-
|
61
|
+
n = 10000
|
62
|
-
|
63
|
-
|
64
62
|
|
65
63
|
U = []
|
66
64
|
|
67
65
|
for i in range(n):
|
68
66
|
|
69
|
-
|
67
|
+
x = np.random.rand(i)
|
70
68
|
|
71
|
-
x = np.random.rand(n)
|
72
|
-
|
73
|
-
y = np.random.rand(
|
69
|
+
y = np.random.rand(i)
|
74
70
|
|
75
71
|
kyori = np.sqrt(x*y+y*y)
|
76
72
|
|
@@ -96,10 +92,8 @@
|
|
96
92
|
|
97
93
|
```
|
98
94
|
|
99
|
-
<function matplotlib.pyplot.show(*args, **kw)>
|
100
95
|
|
101
96
|
|
97
|
+
![イメージ説明](4a840cb42dd6fe3481f8276bd230eb8e.png)
|
102
98
|
|
103
|
-
|
99
|
+
直線になってしまいました。
|
104
|
-
|
105
|
-
![イメージ説明](cd004496425c8e035ea1e6d4f0291df6.png)
|
6
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -101,3 +101,5 @@
|
|
101
101
|
|
102
102
|
|
103
103
|
書き直しました。がグラフがすごい荒れてます
|
104
|
+
|
105
|
+
![イメージ説明](cd004496425c8e035ea1e6d4f0291df6.png)
|
5
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -96,4 +96,8 @@
|
|
96
96
|
|
97
97
|
```
|
98
98
|
|
99
|
+
<function matplotlib.pyplot.show(*args, **kw)>
|
100
|
+
|
101
|
+
|
102
|
+
|
99
103
|
書き直しました。がグラフがすごい荒れてます
|
4
説明の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -96,4 +96,4 @@
|
|
96
96
|
|
97
97
|
```
|
98
98
|
|
99
|
-
|
99
|
+
書き直しました。がグラフがすごい荒れてます
|
3
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -60,29 +60,27 @@
|
|
60
60
|
|
61
61
|
|
62
62
|
|
63
|
-
n = 10000
|
64
63
|
|
65
|
-
x = np.random.rand(n)
|
66
64
|
|
67
|
-
y = np.random.rand(n)
|
68
|
-
|
69
|
-
|
65
|
+
U = []
|
70
66
|
|
71
67
|
for i in range(n):
|
68
|
+
|
69
|
+
n = 10000
|
70
|
+
|
71
|
+
x = np.random.rand(n)
|
72
|
+
|
73
|
+
y = np.random.rand(n)
|
72
74
|
|
73
75
|
kyori = np.sqrt(x*y+y*y)
|
74
76
|
|
75
77
|
m = np.sum(kyori <= 1)
|
76
78
|
|
77
|
-
print('原点からの距離が1以内の点の数:',m)
|
78
|
-
|
79
79
|
p = 4*m/n
|
80
80
|
|
81
|
-
|
81
|
+
U.append(p - np.pi)
|
82
82
|
|
83
|
-
|
84
|
-
|
85
|
-
plt.plot(
|
83
|
+
plt.plot(U)
|
86
84
|
|
87
85
|
plt.title('Error/Number of points')
|
88
86
|
|
@@ -98,76 +96,4 @@
|
|
98
96
|
|
99
97
|
```
|
100
98
|
|
101
|
-
ValueError Traceback (most recent call last)
|
102
|
-
|
103
|
-
<ipython-input-8-7a8c12fc9b9e> in <module>
|
104
|
-
|
105
|
-
7 x = []
|
106
|
-
|
107
|
-
8 for i in range(n):
|
108
|
-
|
109
|
-
----> 9 kyori = np.sqrt(x*y+y*y)
|
110
|
-
|
111
|
-
10 m = np.sum(kyori <= 1)
|
112
|
-
|
113
|
-
11 print('原点からの距離が1以内の点の数:',m)
|
114
|
-
|
115
|
-
|
116
|
-
|
117
|
-
ValueError: operands could not be broadcast together with shapes (0,) (10000,)
|
118
|
-
|
119
|
-
|
120
|
-
|
121
|
-
|
122
|
-
|
123
|
-
ここまで書き直しました
|
124
|
-
|
125
|
-
```
|
126
|
-
|
127
|
-
import numpy as np
|
128
|
-
|
129
|
-
import matplotlib.pyplot as plt
|
130
|
-
|
131
|
-
|
132
|
-
|
133
|
-
n = 10000
|
134
|
-
|
135
|
-
x = np.random.rand(n)
|
136
|
-
|
137
|
-
y = np.random.rand(n)
|
138
|
-
|
139
|
-
|
99
|
+
もう少しです
|
140
|
-
|
141
|
-
for i in range(10000):
|
142
|
-
|
143
|
-
kyori = np.sqrt(x*y+y*y)
|
144
|
-
|
145
|
-
m = np.sum(kyori <= 1)
|
146
|
-
|
147
|
-
p = 4*m/n
|
148
|
-
|
149
|
-
U.append(int(p - np.pi))
|
150
|
-
|
151
|
-
|
152
|
-
|
153
|
-
plt.plot(U,)
|
154
|
-
|
155
|
-
plt.title('Error/Number of points')
|
156
|
-
|
157
|
-
plt.xlabel('Number of points')
|
158
|
-
|
159
|
-
plt.ylabel('Error')
|
160
|
-
|
161
|
-
plt.grid(True)
|
162
|
-
|
163
|
-
plt.show
|
164
|
-
|
165
|
-
コード
|
166
|
-
|
167
|
-
```
|
168
|
-
|
169
|
-
<function matplotlib.pyplot.show(*args, **kw)>
|
170
|
-
|
171
|
-
|
172
|
-
|
173
|
-
ここから縦軸が何なのか、リスト化して折れ線グラフにする方法がわかりません。上の条件どおりリスト化できてるかも怪しいです。
|
2
実行結果の追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -166,4 +166,8 @@
|
|
166
166
|
|
167
167
|
```
|
168
168
|
|
169
|
+
<function matplotlib.pyplot.show(*args, **kw)>
|
170
|
+
|
171
|
+
|
172
|
+
|
169
173
|
ここから縦軸が何なのか、リスト化して折れ線グラフにする方法がわかりません。上の条件どおりリスト化できてるかも怪しいです。
|
1
説明の追加
test
CHANGED
@@ -1 +1 @@
|
|
1
|
-
リスト化、
|
1
|
+
リスト化、折れ線グラフの表し方
|
test
CHANGED
@@ -115,3 +115,55 @@
|
|
115
115
|
|
116
116
|
|
117
117
|
ValueError: operands could not be broadcast together with shapes (0,) (10000,)
|
118
|
+
|
119
|
+
|
120
|
+
|
121
|
+
|
122
|
+
|
123
|
+
ここまで書き直しました
|
124
|
+
|
125
|
+
```
|
126
|
+
|
127
|
+
import numpy as np
|
128
|
+
|
129
|
+
import matplotlib.pyplot as plt
|
130
|
+
|
131
|
+
|
132
|
+
|
133
|
+
n = 10000
|
134
|
+
|
135
|
+
x = np.random.rand(n)
|
136
|
+
|
137
|
+
y = np.random.rand(n)
|
138
|
+
|
139
|
+
U = []
|
140
|
+
|
141
|
+
for i in range(10000):
|
142
|
+
|
143
|
+
kyori = np.sqrt(x*y+y*y)
|
144
|
+
|
145
|
+
m = np.sum(kyori <= 1)
|
146
|
+
|
147
|
+
p = 4*m/n
|
148
|
+
|
149
|
+
U.append(int(p - np.pi))
|
150
|
+
|
151
|
+
|
152
|
+
|
153
|
+
plt.plot(U,)
|
154
|
+
|
155
|
+
plt.title('Error/Number of points')
|
156
|
+
|
157
|
+
plt.xlabel('Number of points')
|
158
|
+
|
159
|
+
plt.ylabel('Error')
|
160
|
+
|
161
|
+
plt.grid(True)
|
162
|
+
|
163
|
+
plt.show
|
164
|
+
|
165
|
+
コード
|
166
|
+
|
167
|
+
```
|
168
|
+
|
169
|
+
ここから縦軸が何なのか、リスト化して折れ線グラフにする方法がわかりません。上の条件どおりリスト化できてるかも怪しいです。
|