Pythonの図形の書き方を教えて頂きたいです

Question

モンテカルロ法による円周率の算出をPython により行いたいです。
そこで、一辺の長さが1の正方形に半径1の1/4の円が入っているとします。
この中に次々とダーツ(点)を投げます。
この時の様子を可視化したいのですが、
1/4の円を描画し、さらにその上から点を描画する方法を教えて頂きたいです。

Comments

[otn]

質問タイトルは、「円周率の算出方法」じゃなくて「図形の書き方」では？

Answer 1

matplotlibにて、円（扇形）、正方形を描画するにはmatplotlib.patchesAPIをお使いください

https://matplotlib.org/3.1.1/api/patches_api.html

import matplotlib.pyplot as plt
import matplotlib.patches as pat

ax = plt.subplot()

# 中心:(0,0),半径:1,0の扇形（0度～90度）を生成
# Face Color:blue, Alpha値:0.3
w = pat.Wedge((0,0), 1, 0, 90, fc='blue', alpha=0.3)
# 扇形をグラフに追加
ax.add_patch(w)

# 始点:(0,0),幅:1,高さ:1 の正方形を生成
# Face Color:なし(透明),Edge Color:Black, Line Width:2
r = pat.Rectangle((0,0), 1, 1, fc='none', ec='black', lw=2)
# 正方形をグラフに追加
ax.add_patch(r)

# (0.5, 0.5)の点を描画
ax.plot(0.5, 0.5, '.r')

# 描画範囲を設定
ax.set_xlim(-0.1, 1.1)
ax.set_ylim(-0.1, 1.1)

# 正方ピクセルに設定
ax.set_aspect('equal')

plt.show()