Pythonで円周率を計算!様々な方法を徹底解説
Pythonを使って円周率を推定するためにモンテカルロ法を利用することができます。具体的な手順は次のとおりです。
- 正方形の領域内にランダムに多数の点を生成する。
- 単位円内にある点の数を数える
- 単位円の面積(つまりπ/4)と単位正方形の面積の比を計算し、πの近似値を得る。
以下是一倅简单的Python代码示例:
import random
def estimate_pi(num_points):
points_inside_circle = 0
for _ in range(num_points):
x = random.uniform(0, 1)
y = random.uniform(0, 1)
if x**2 + y**2 <= 1:
points_inside_circle += 1
pi_estimate = 4 * points_inside_circle / num_points
return pi_estimate
num_points = 1000000
pi_approx = estimate_pi(num_points)
print("Approximated value of pi:", pi_approx)
上記の例では、100万個の点を生成し、単位円内に落ちる点の数を数えることで円周率の値を推定しました。より正確な推定結果を得るためには、num_pointsの値を調整してみてください。